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

Inhibition of the voltage-gated potassium channel Kv1.5 by hydrogen sulfide attenuates remodeling through S-nitrosylation-mediated signaling

The voltage-gated K+ channel plays a key role in atrial excitability, conducting the ultra-rapid rectifier K+ current (IKur) and contributing to the repolarization of the atrial action potential. In this study, we examine its regulation by hydrogen sulfide (H2S) in HL-1 cardiomyocytes and in HEK293 cells expressing human Kv1.5. Pacing induced remodeling resulted in shorting action potential duration, enhanced both Kv1.5 channel and H2S producing enzymes protein expression in HL-1 cardiomyocytes. H2S supplementation reduced these remodeling changes and restored action potential duration through inhibition of Kv1.5 channel. H2S also inhibited recombinant hKv1.5, lead to nitric oxide (NO) mediated S-nitrosylation and activated endothelial nitric oxide synthase (eNOS) by increased phosphorylation of Ser1177, prevention of NO formation precluded these effects. Regulation of Ikur by H2S has important cardiovascular implications and represents a novel and potential therapeutic target.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Optimal entrainment for removal of pinned spiral waves

Cardiac fibrillation is caused by self-sustaining spiral waves that occur in the myocardium, some of which can be pinned to anatomical obstacles, making them more difficult to eliminate. A small electrical stimulation is often sufficient to unpin these spirals but only if it is applied during the vulnerable unpinning window. Even if these unpinning windows can be inferred from data, when multiple pinned spirals exist, their unpinning windows will not generally overlap. Using phase-based reduction techniques, we formulate and solve an optimal control problem to yield a time-varying external voltage gradient that can synchronize a collection of spiral waves that are pinned to a collection of heterogeneous obstacles. Upon synchronization, the unpinning windows overlap so that they can be simultaneously unpinned by applying an external voltage gradient pulse at an appropriate moment. Numerical validation is presented in bidomain model simulations. Results represent a proof-of-concept illustration of the proposed unpinning strategy which explicitly incorporates heterogeneity in the problem formulation and requires no real-time feedback about the system state.

Panoramic Endocardial Optical Mapping Demonstrates Serial Rotors Acceleration and Increasing Complexity of Activity During Onset of Cholinergic Atrial Fibrillation

Background

Activation during onset of atrial fibrillation is poorly understood. We aimed at developing a panoramic optical mapping system for the atria and test the hypothesis that sequential rotors underlie acceleration of atrial fibrillation during onset.

Methods and Results

Five sheep hearts were Langendorff perfused in the presence of 0.25 µmol/L carbachol. Novel optical system recorded activations simultaneously from the entire left and right atrial endocardial surfaces. Twenty sustained (>40 s) atrial fibrillation episodes were induced by a train and premature stimuli protocol. Movies obtained immediately (Initiation stage) and 30 s (Early Stabilization stage) after premature stimulus were analyzed. Serial rotor formation was observed in all sustained inductions and none in nonsustained inductions. In sustained episodes maximal dominant frequency increased from (mean±SD) 11.5±1.74 Hz during Initiation to 14.79±1.30 Hz at Early Stabilization (P<0.0001) and stabilized thereafter. At rotor sites, mean cycle length (CL) during 10 prerotor activations increased every cycle by 0.53% (P=0.0303) during Initiation and 0.34% (P=0.0003) during Early Stabilization. In contrast, CLs at rotor sites showed abrupt decreases after the rotors appearances by a mean of 9.65% (P<0.0001) during both stages. At Initiation, atria‐wide accelerations and decelerations during rotors showed a net acceleration result whereby post‐rotors atria‐wide minimal CL (CLmin) were 95.5±6.8% of the prerotor CLmin (P=0.0042). In contrast, during Early Stabilization, there was no net acceleration in CLmin during accelerating rotors (prerotor=84.9±11.0% versus postrotor=85.8±10.8% of Initiation, P=0.4029). Levels of rotor drift distance and velocity correlated with atria‐wide acceleration. Nonrotor phase singularity points did not accelerate atria‐wide activation but multiplied during Initiation until Early Stabilization. Increasing number of singularity points, indicating increased complexity, correlated with atria‐wide CLmin reduction (P<0.0001).

Conclusions

Novel panoramic optical mapping of the atria demonstrates shortening CL at rotor sites during cholinergic atrial fibrillation onset. Atrial fibrillation acceleration toward Early Stabilization correlates with the net result of atria‐wide accelerations during drifting rotors activity.

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

The impact of wall thickness and curvature on wall stress in patient-specific electromechanical models of the left atrium

The left atrium (LA) has a complex anatomy with heterogeneous wall thickness and curvature. The anatomy plays an important role in determining local wall stress; however, the relative contribution of wall thickness and curvature in determining wall stress in the LA is unknown. We have developed electromechanical finite element (FE) models of the LA using patient-specific anatomical FE meshes with rule-based myofiber directions. The models of the LA were passively inflated to 10mmHg followed by simulation of the contraction phase of the atrial cardiac cycle. The FE models predicted maximum LA volumes of 156.5 mL, 99.3 mL and 83.4 mL and ejection fractions of 36.9%, 32.0% and 25.2%. The median wall thickness in the 3 cases was calculated as [Formula: see text] mm, [Formula: see text] mm, and [Formula: see text] mm. The median curvature was determined as [Formula: see text] [Formula: see text], [Formula: see text], and [Formula: see text]. Following passive inflation, the correlation of wall stress with the inverse of wall thickness and curvature was 0.55-0.62 and 0.20-0.25, respectively. At peak contraction, the correlation of wall stress with the inverse of wall thickness and curvature was 0.38-0.44 and 0.16-0.34, respectively. In the LA, the 1st principal Cauchy stress is more dependent on wall thickness than curvature during passive inflation and both correlations decrease during active contraction. This emphasizes the importance of including the heterogeneous wall thickness in electromechanical FE simulations of the LA. Overall, simulation results and sensitivity analyses show that in complex atrial anatomy it is unlikely that a simple anatomical-based law can be used to estimate local wall stress, demonstrating the importance of FE analyses.

Atrial fibrillation source area probability mapping using electrogram patterns of multipole catheters

BACKGROUND

Catheter ablation therapy involving isolation of pulmonary veins (PVs) from the left atrium is performed to terminate atrial fibrillation (AF). Unfortunately, standalone PV isolation procedure has shown to be a suboptimal success with AF continuation or recurrence. One reason, especially in patients with persistent or high-burden paroxysmal AF, is known to be due to the formation of repeating-pattern AF sources with a meandering core inside the atria. However, there is a need for accurate mapping and localization of these sources during catheter ablation.

METHODS

A novel AF source area probability (ASAP) mapping algorithm was developed and evaluated in 2D and 3D atrial simulated tissues with various arrhythmia scenarios and a retrospective study with three cases of clinical human AF. The ASAP mapping analyzes the electrograms collected from a multipole diagnostic catheter that is commonly used during catheter ablation procedure to intelligently sample the atria and delineate the trajectory path of a meandering repeating-pattern AF source. ASAP starts by placing the diagnostic catheter at an arbitrary location in the atria. It analyzes the recorded bipolar electrograms to build an ASAP map over the atrium anatomy and suggests an optimal location for the subsequent catheter location. ASAP then determines from the constructed ASAP map if an AF source has been delineated. If so, the catheter navigation is stopped and the algorithm provides the area of the AF source. Otherwise, the catheter is navigated to the suggested location, and the process is continued until an AF-source area is delineated.

RESULTS

ASAP delineated the AF source in over 95% of the simulated human AF cases within less than eight catheter placements regardless of the initial catheter placement. The success of ASAP in the clinical AF was confirmed by the ablation outcomes and the electrogram patterns at the delineated area.

CONCLUSION

Our analysis indicates the potential of the ASAP mapping to provide accurate information about the area of the meandering repeating-pattern AF sources as AF ablation targets for effective AF termination. Our algorithm could improve the success of AF catheter ablation therapy by locating and subsequently targeting patient-specific and repeating-pattern AF sources inside the atria.

Sensitivity and Uncertainty Analysis of Two Human Atrial Cardiac Cell Models Using Gaussian Process Emulators

Biophysically detailed cardiac cell models reconstruct the action potential and calcium dynamics of cardiac myocytes. They aim to capture the biophysics of current flow through ion channels, pumps, and exchangers in the cell membrane, and are highly detailed. However, the relationship between model parameters and model outputs is difficult to establish because the models are both complex and non-linear. The consequences of uncertainty and variability in model parameters are therefore difficult to determine without undertaking large numbers of model evaluations. The aim of the present study was to demonstrate how sensitivity and uncertainty analysis using Gaussian process emulators can be used for a systematic and quantitive analysis of biophysically detailed cardiac cell models. We selected the Courtemanche and Maleckar models of the human atrial action potential for analysis because these models describe a similar set of currents, with different formulations. In our approach Gaussian processes emulate the main features of the action potential and calcium transient. The emulators were trained with a set of design data comprising samples from parameter space and corresponding model outputs, initially obtained from 300 model evaluations. Variance based sensitivity indices were calculated using the emulators, and first order and total effect indices were calculated for each combination of parameter and output. The differences between the first order and total effect indices indicated that the effect of interactions between parameters was small. A second set of emulators were then trained using a new set of design data with a subset of the model parameters with a sensitivity index of more than 0.1 (10%). This second stage analysis enabled comparison of mechanisms in the two models. The second stage sensitivity indices enabled the relationship between the L-type Ca 2+ current and the action potential plateau to be quantified in each model. Our quantitative analysis predicted that changes in maximum conductance of the ultra-rapid K + channel I Kur would have opposite effects on action potential duration in the two models, and this prediction was confirmed by additional simulations. This study has demonstrated that Gaussian process emulators are an effective tool for sensitivity and uncertainty analysis of biophysically detailed cardiac cell models.

Phase-amplitude reduction far beyond the weakly perturbed paradigm

While phase reduction is a well-established technique for the analysis of perturbed limit cycle oscillators, practical application requires perturbations to be sufficiently weak thereby limiting its utility in many situations. Here, a general strategy is developed for constructing a set of phase-amplitude reduced equations that is valid to arbitrary orders of accuracy in the amplitude coordinates. This reduction framework can be used to investigate the behavior of oscillatory dynamical systems far beyond the weakly perturbed paradigm. Additionally, a patchwork phase-amplitude reduction method is suggested that is useful when exceedingly large magnitude perturbations are considered. This patchwork method incorporates the high-accuracy phase-amplitude reductions of multiple nearby periodic orbits that result from modifications to nominal parameters. The proposed method of high-accuracy phase-amplitude reduction can be readily implemented numerically and examples are provided where reductions are computed up to fourteenth order accuracy.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Locating Atrial Fibrillation Rotor and Focal Sources Using Iterative Navigation of Multipole Diagnostic Catheters

Purpose

Multi-polar diagnostic catheters are used to construct the 3D electro-anatomic mapping of the atrium during atrial fibrillation (AF) ablation procedures; however, it remains unclear how to use the electrograms recorded by these catheters to locate AF-driving sites known as focal and rotor source types. The purpose of this study is to present the first algorithm to iteratively navigate a circular multi-polar catheter to locate AF focal and rotor sources without the need to map the entire atria.

Methods

Starting from an initial location, the algorithm, which was blinded to the location and type of the AF source, iteratively advanced a Lasso catheter based on its electrogram characteristics. The algorithm stopped the catheter when it located of an AF source and identified the type. The efficiency of the algorithm is validated using a set of simulated focal and rotor-driven arrhythmias in fibrotic human 2D and 3D atrial tissue.

Results

Our study shows the feasibility of locating AF sources with a success rate of greater than 95.25% within average 7.56 ± 2.28 placements independently of the initial position of the catheter and the source type.

Conclusions

The algorithm could play a critical role in clinical electrophysiology laboratories for mapping patient-specific ablation of AF sources located outside the pulmonary veins and improving the procedure success.

Electronic supplementary material

The online version of this article (10.1007/s13239-019-00414-5) contains supplementary material, which is available to authorized users.

Phase reduction and phase-based optimal control for biological systems: a tutorial

A powerful technique for the analysis of nonlinear oscillators is the rigorous reduction to phase models, with a single variable describing the phase of the oscillation with respect to some reference state. An analog to phase reduction has recently been proposed for systems with a stable fixed point, and phase reduction for periodic orbits has recently been extended to take into account transverse directions and higher-order terms. This tutorial gives a unified treatment of such phase reduction techniques and illustrates their use through mathematical and biological examples. It also covers the use of phase reduction for designing control algorithms which optimally change properties of the system, such as the phase of the oscillation. The control techniques are illustrated for example neural and cardiac systems.

Iterative navigation of multipole diagnostic catheters to locate repeating-pattern atrial fibrillation drivers

Introduction

Targeting repeating‐pattern atrial fibrillation (AF) sources (reentry or focal drivers) can help in patient‐specific ablation therapy for AF; however, the development of reliable and accurate tools for locating such sources remains a major challenge. We describe iterative catheter navigation (ICAN) algorithm to locate AF drivers using a conventional circular Lasso catheter.

Methods and Results

At each step, the algorithm analyzes 10 bipolar electrograms recoded at a given catheter location and the history of previous catheter movements to determine if the source is inside the catheter loop. If not, it calculates new coordinates and selects a new position for the catheter. The process continues until a source is located. The algorithm was evaluated in a computer model of atrial tissue with various degrees of fibrosis under a broad range of arrhythmia scenarios. The latter included slow and fast reentry, macroreentry, figure‐of‐eight reentry, and fibrillatory conduction. Depending on the initial distance of the catheter from the source and scenario, it took about 3 to 16 steps to localize an AF source. In 94% of cases, the identified location was within 4 mm from the source, independently of the initial position of the catheter. The algorithm worked equally well in the presence of patchy fibrosis, low‐voltage areas, fragmented electrograms, and dominant‐frequency gradients.

Conclusions

AF repeating‐pattern sources can be localized using circular catheters without the need to map the entire tissue. The proposed algorithm has the potential to become a useful tool for patient‐specific ablation of AF sources located outside the pulmonary veins.

Local activation time sampling density for atrial tachycardia contact mapping: how much is enough?

Aims

Local activation time (LAT) mapping forms the cornerstone of atrial tachycardia diagnosis. Although anatomic and positional accuracy of electroanatomic mapping (EAM) systems have been validated, the effect of electrode sampling density on LAT map reconstruction is not known. Here, we study the effect of chamber geometry and activation complexity on optimal LAT sampling density using a combined in silico and in vivo approach.

Methods and results

In vivo 21 atrial tachycardia maps were studied in three groups: (1) focal activation, (2) macro-re-entry, and (3) localized re-entry. In silico activation was simulated on a 4×4cm atrial monolayer, sampled randomly at 0.25–10 points/cm² and used to re-interpolate LAT maps. Activation patterns were studied in the geometrically simple porcine right atrium (RA) and complex human left atrium (LA). Activation complexity was introduced into the porcine RA by incomplete inter-caval linear ablation. In all cases, optimal sampling density was defined as the highest density resulting in minimal further error reduction in the re-interpolated maps. Optimal sampling densities for LA tachycardias were 0.67 ± 0.17 points/cm² (focal activation), 1.05 ± 0.32 points/cm² (macro-re-entry) and 1.23 ± 0.26 points/cm² (localized re-entry), P = 0.0031. Increasing activation complexity was associated with increased optimal sampling density both in silico (focal activation 1.09 ± 0.14 points/cm²; re-entry 1.44 ± 0.49 points/cm²; spiral-wave 1.50 ± 0.34 points/cm², P < 0.0001) and in vivo (porcine RA pre-ablation 0.45 ± 0.13 vs. post-ablation 0.78 ± 0.17 points/cm², P = 0.0008). Increasing chamber geometry was also associated with increased optimal sampling density (0.61 ± 0.22 points/cm² vs. 1.0 ± 0.34 points/cm², P = 0.0015).

Conclusion

Optimal sampling densities can be identified to maximize diagnostic yield of LAT maps. Greater sampling density is required to correctly reveal complex activation and represent activation across complex geometries. Overall, the optimal sampling density for LAT map interpolation defined in this study was ∼1.0–1.5 points/cm².

Wavelength and Fibrosis Affect Phase Singularity Locations During Atrial Fibrillation

The mechanisms underlying atrial fibrillation (AF), the most common sustained cardiac rhythm disturbance, remain elusive. Atrial fibrosis plays an important role in the development of AF and rotor dynamics. Both electrical wavelength (WL) and the degree of atrial fibrosis change as AF progresses. However, their combined effect on rotor core location remains unknown. The aim of this study was to analyze the effects of WL change on rotor core location in both fibrotic and non-fibrotic atria. Three patient specific fibrosis distributions (total fibrosis content: 16.6, 22.8, and 19.2%) obtained from clinical imaging data of persistent AF patients were incorporated in a bilayer atrial computational model. Fibrotic effects were modeled as myocyte-fibroblast coupling + conductivity remodeling; structural remodeling; ionic current changes + conductivity remodeling; and combinations of these methods. To change WL, action potential duration (APD) was varied from 120 to 240ms, representing the range of clinically observed AF cycle length, by modifying the inward rectifier potassium current (IK1) conductance between 80 and 140% of the original value. Phase singularities (PSs) were computed to identify rotor core locations. Our results show that IK1 conductance variation resulted in a decrease of APD and WL across the atria. For large WL in the absence of fibrosis, PSs anchored to regions with high APD gradient at the center of the left atrium (LA) anterior wall and near the junctions of the inferior pulmonary veins (PVs) with the LA. Decreasing the WL induced more PSs, whose distribution became less clustered. With fibrosis, PS locations depended on the fibrosis distribution and the fibrosis implementation method. The proportion of PSs in fibrotic areas and along the borders varied with both WL and fibrosis modeling method: for patient one, this was 4.2-14.9% as IK1 varied for the structural remodeling representation, but 12.3-88.4% using the combination of structural remodeling with myocyte-fibroblast coupling. The degree and distribution of fibrosis and the choice of implementation technique had a larger effect on PS locations than the WL variation. Thus, distinguishing the fibrotic mechanisms present in a patient is important for interpreting clinical fibrosis maps to create personalized models.

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.

Atrial Rotor Dynamics Under Complex Fractional Order Diffusion

The mechanisms of atrial fibrillation (AF) are a challenging research topic. The rotor hypothesis states that the AF is sustained by a reentrant wave that propagates around an unexcited core. Cardiac tissue heterogeneities, both structural and cellular, play an important role during fibrillatory dynamics, so that the ionic characteristics of the currents, their spatial distribution and their structural heterogeneity determine the meandering of the rotor. Several studies about rotor dynamics implement the standard diffusion equation. However, this mathematical scheme carries some limitations. It assumes the myocardium as a continuous medium, ignoring, therefore, its discrete and heterogeneous aspects. A computational model integrating both, electrical and structural properties could complement experimental and clinical results. A new mathematical model of the action potential propagation, based on complex fractional order derivatives is presented. The complex derivative order appears of considering the myocardium as discrete-scale invariant fractal. The main aim is to study the role of a myocardial, with fractal characteristics, on atrial fibrillatory dynamics. For this purpose, the degree of structural heterogeneity is described through derivatives of complex order γ = α + jβ. A set of variations for γ is tested. The real part α takes values ranging from 1.1 to 2 and the imaginary part β from 0 to 0.28. Under this scheme, the standard diffusion is recovered when α = 2 and β = 0. The effect of γ on the action potential propagation over an atrial strand is investigated. Rotors are generated in a 2D model of atrial tissue under electrical remodeling due to chronic AF. The results show that the degree of structural heterogeneity, given by γ, modulates the electrophysiological properties and the dynamics of rotor-type reentrant mechanisms. The spatial stability of the rotor and the area of its unexcited core are modulated. As the real part decreases and the imaginary part increases, simulating a higher structural heterogeneity, the vulnerable window to reentrant is increased, as the total meandering of the rotor tip. This in silico study suggests that structural heterogeneity, described by means of complex order derivatives, modulates the stability of rotors and that a wide range of rotor dynamics can be generated.

Simulation of Spiral Waves and Point Sources in Atrial Fibrillation with Application to Rotor Localization

Cardiac simulations play an important role in studies involving understanding and investigating the mechanisms of cardiac arrhythmias. Today, studies of arrhythmogenesis and maintenance are largely being performed by creating simulations of a particular arrhythmia with high accuracy comparable to the results of clinical experiments. Atrial fibrillation (AF), the most common arrhythmia in the United States and many other parts of the world, is one of the major field where simulation and modeling is largely used. AF simulations not only assist in understanding its mechanisms but also help to develop, evaluate and improve the computer algorithms used in electrophysiology (EP) systems for ablation therapies. In this paper, we begin with a brief overeview of some common techniques used in simulations to simulate two major AF mechanisms - spiral waves (or rotors) and point (or focal) sources. We particularly focus on 2D simulations using Nygren et al.'s mathematical model of human atrial cell. Then, we elucidate an application of the developed AF simulation to an algorithm designed for localizing AF rotors for improving current AF ablation therapies. Our simulation methods and results, along with the other discussions presented in this paper is aimed to provide engineers and professionals with a working-knowledge of application-specific simulations of spirals and foci.

Machine learning methods for locating re-entrant drivers from electrograms in a model of atrial fibrillation

Mapping resolution has recently been identified as a key limitation in successfully locating the drivers of atrial fibrillation (AF). Using a simple cellular automata model of AF, we demonstrate a method by which re-entrant drivers can be located quickly and accurately using a collection of indirect electrogram measurements. The method proposed employs simple, out-of-the-box machine learning algorithms to correlate characteristic electrogram gradients with the displacement of an electrogram recording from a re-entrant driver. Such a method is less sensitive to local fluctuations in electrical activity. As a result, the method successfully locates 95.4% of drivers in tissues containing a single driver, and 95.1% (92.6%) for the first (second) driver in tissues containing two drivers of AF. Additionally, we demonstrate how the technique can be applied to tissues with an arbitrary number of drivers. In its current form, the techniques presented are not refined enough for a clinical setting. However, the methods proposed offer a promising path for future investigations aimed at improving targeted ablation for AF.

Stochastic Pacing Inhibits Spatially Discordant Cardiac Alternans

Depressed heart rate variability is a well-established risk factor for sudden cardiac death in survivors of acute myocardial infarction and for those with congestive heart failure. Although measurements of heart rate variability provide a valuable prognostic tool, it is unclear whether reduced heart rate variability itself is proarrhythmic or if it simply correlates with the severity of autonomic nervous system dysfunction. In this work, we investigate a possible mechanism by which heart rate variability could protect against cardiac arrhythmia. Specifically, in numerical simulations, we observe an inverse relationship between the variance of stochastic pacing and the occurrence of spatially discordant alternans, an arrhythmia that is widely believed to facilitate the development of cardiac fibrillation. By analyzing the effects of conduction velocity restitution, cellular dynamics, electrotonic coupling, and stochastic pacing on the nodal dynamics of spatially discordant alternans, we provide intuition for this observed behavior and propose control strategies to inhibit discordant alternans.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model

Accumulation of intracellular Na⁺ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na⁺ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na⁺ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na⁺]_i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na⁺]_i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na⁺]_i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na⁺]_i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na⁺]_i may play complex roles in cellular and tissue-level cardiac dynamics.

Emergent dynamics of spatio-temporal chaos in a heterogeneous excitable medium

Self-organized activation patterns in excitable media such as spiral waves and spatio-temporal chaos underlie dangerous cardiac arrhythmias. While the interaction of single spiral waves with different types of heterogeneity has been studied extensively, the effect of heterogeneity on fully developed spatio-temporal chaos remains poorly understood. We investigate how the complexity and stability properties of spatio-temporal chaos in the Bär-Eiswirth model of excitable media depend on the heterogeneity of the underlying medium. We employ different measures characterizing the chaoticity of the system and find that the spatial arrangement of multiple discrete lower excitability regions has a strong impact on the complexity of the dynamics. Varying the number, shape, and spatial arrangement of the heterogeneities, we observe strong emergent effects ranging from increases in chaoticity to the complete cessation of chaos, contrasting the expectation from the homogeneous behavior. The implications of our findings for the development and treatment of arrhythmias in the heterogeneous cardiac muscle are discussed.

Comparison of Detailed and Simplified Models of Human Atrial Myocytes to Recapitulate Patient Specific Properties

Computer studies are often used to study mechanisms of cardiac arrhythmias, including atrial fibrillation (AF). A crucial component in these studies is the electrophysiological model that describes the membrane potential of myocytes. The models vary from detailed, describing numerous ion channels, to simplified, grouping ionic channels into a minimal set of variables. The parameters of these models, however, are determined across different experiments in varied species. Furthermore, a single set of parameters may not describe variations across patients, and models have rarely been shown to recapitulate critical features of AF in a given patient. In this study we develop physiologically accurate computational human atrial models by fitting parameters of a detailed and of a simplified model to clinical data for five patients undergoing ablation therapy. Parameters were simultaneously fitted to action potential (AP) morphology, action potential duration (APD) restitution and conduction velocity (CV) restitution curves in these patients. For both models, our fitting procedure generated parameter sets that accurately reproduced clinical data, but differed markedly from published sets and between patients, emphasizing the need for patient-specific adjustment. Both models produced two-dimensional spiral wave dynamics for that were similar for each patient. These results show that simplified, computationally efficient models are an attractive choice for simulations of human atrial electrophysiology in spatially extended domains. This study motivates the development and validation of patient-specific model-based mechanistic studies to target therapy.

Catheter simulator software tool to generate electrograms of any multi-polar diagnostic catheter from 3D atrial tissue

Simulations are excellent tools for assessing new therapeutic strategies and are often conducted before implementing new therapy options in a clinical practice. For patients suffering from a heart arrhythmia, the main source of information comes from an intracardiac catheter. One of the common catheters is a Lasso multi-pole diagnostic catheter, which is a catheter that has 20 electrodes in a circular pattern. In this paper, we developed algorithm and simulation software that allows the users to place a multi-pole catheter on the atrial endocardial surface and record electrograms. In 3D atrial tissue, the plane of principal curvature is determined using eigenvectors of catheter vertices, from where the normals are projected and registered to the surface using 3D geodesic distance. This tool provides a platform for performing customized virtual cardiac experiments.

Toward a More Efficient Implementation of Antifibrillation Pacing

We devise a methodology to determine an optimal pattern of inputs to synchronize firing patterns of cardiac cells which only requires the ability to measure action potential durations in individual cells. In numerical bidomain simulations, the resulting synchronizing inputs are shown to terminate spiral waves with a higher probability than comparable inputs that do not synchronize the cells as strongly. These results suggest that designing stimuli which promote synchronization in cardiac tissue could improve the success rate of defibrillation, and point towards novel strategies for optimizing antifibrillation pacing.

The Multi-Domain Fibroblast/Myocyte Coupling in the Cardiac Tissue: A Theoretical Study

Cardiac fibroblast proliferation and concomitant collagenous matrix accumulation (fibrosis) develop during multiple cardiac pathologies. Recent studies have demonstrated direct electrical coupling between myocytes and fibroblasts in vitro, and assessed the electrophysiological implications of such coupling. However, in the living tissues, such coupling has not been demonstrated, and only indirect coupling via the extracellular space is likely to exist. In this study we employed a multi-domain model to assess the modulation of the cardiac electrophysiological properties by neighboring fibroblasts assuming only indirect coupling. Numerical simulations in 1D and 2D human atrial models showed that extracellular coupling sustains a significant impact on conduction velocity (CV) and a less significant effect on the action potential duration. Both CV and the slope of the CV restitution increased with increasing fibroblast density. This effect was more substantial for lower extracellular conductance. In 2D, spiral waves exhibited reduced frequency with increasing fibroblast density, and the propensity of wavebreaks and complex dynamics at high pacing rates significantly increased.

Computational models of atrial cellular electrophysiology and calcium handling, and their role in atrial fibrillation

The complexity of the heart makes an intuitive understanding of the relative contribution of ion channels, transporters and signalling pathways to cardiac electrophysiology challenging. Computational modelling of cardiac cellular electrophysiology has proven useful to integrate experimental findings, extrapolate results obtained in expression systems or animal models to other systems, test quantitatively ideas based on experimental data and provide novel hypotheses that are experimentally testable. While the bulk of computational modelling has traditionally been directed towards ventricular bioelectricity, increasing recognition of the clinical importance of atrial arrhythmias, particularly atrial fibrillation, has led to widespread efforts to apply computational approaches to understanding atrial electrical function. The increasing availability of detailed, atrial-specific experimental data has stimulated the development of novel computational models of atrial-cellular electrophysiology and Ca(2+) handling. To date, more than 300 studies have employed mathematical simulations to enhance our understanding of atrial electrophysiology, arrhythmogenesis and therapeutic responses. Future modelling studies are likely to move beyond current whole-cell models by incorporating new data on subcellular architecture, macromolecular protein complexes, and localized ion-channel regulation by signalling pathways. At the same time, more integrative multicellular models that take into account regional electrophysiological and Ca(2+) handling properties, mechano-electrical feedback and/or autonomic regulation will be needed to investigate the mechanisms governing atrial arrhythmias. A combined experimental and computational approach is expected to provide the more comprehensive understanding of atrial arrhythmogenesis that is required to develop improved diagnostic and therapeutic options. Here, we review this rapidly expanding area, with a particular focus on Ca(2+) handling, and provide ideas about potential future directions.

Characterization of Electrograms from Multipolar Diagnostic Catheters during Atrial Fibrillation

Atrial fibrillation (AF) is the most common arrhythmia in USA with more than 2.3 million people affected annually. Catheter ablation procedure is a method for treatment of AF, which involves 3D electroanatomic mapping of the patient's left atrium (LA) by maneuvering a conventional multipolar diagnostic catheter (MPDC) along the LA endocardial surface after which pulmonary vein (PV) isolation is performed, thus eliminating the AF triggers originating from the PVs. However, it remains unclear how to effectively utilize the information provided by the MPDC to locate the AF-sustaining sites, known as sustained rotor-like activities (RotAs). In this study, we use computer modeling to investigate the variations in the characteristics of the MPDC electrograms, namely, total conduction delay (TCD) and average cycle length (CL), as the MPDC moves towards a RotA source. Subsequently, a study with a human subject was performed in order to verify the predictions of the simulation study. The conclusions from this study may be used to iteratively direct an MPDC towards RotA sources thus allowing the RotAs to be localized for customized and improved AF ablation.

Human atrial cell models to analyse haemodialysis-related effects on cardiac electrophysiology: work in progress

During haemodialysis (HD) sessions, patients undergo alterations in the extracellular environment, mostly concerning plasma electrolyte concentrations, pH, and volume, together with a modification of sympathovagal balance. All these changes affect cardiac electrophysiology, possibly leading to an increased arrhythmic risk. Computational modeling may help to investigate the impact of HD-related changes on atrial electrophysiology. However, many different human atrial action potential (AP) models are currently available, all validated only with the standard electrolyte concentrations used in experiments. Therefore, they may respond in different ways to the same environmental changes. After an overview on how the computational approach has been used in the past to investigate the effect of HD therapy on cardiac electrophysiology, the aim of this work has been to assess the current state of the art in human atrial AP models, with respect to the HD context. All the published human atrial AP models have been considered and tested for electrolytes, volume changes, and different acetylcholine concentrations. Most of them proved to be reliable for single modifications, but all of them showed some drawbacks. Therefore, there is room for a new human atrial AP model, hopefully able to physiologically reproduce all the HD-related effects. At the moment, work is still in progress in this specific field.

Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation

AIMS

Human atrial electrophysiology exhibits high inter-subject variability in both sinus rhythm (SR) and chronic atrial fibrillation (cAF) patients. Variability is however rarely investigated in experimental and theoretical electrophysiological studies, thus hampering the understanding of its underlying causes but also its implications in explaining differences in the response to disease and treatment. In our study, we aim at investigating the ability of populations of human atrial cell models to capture the inter-subject variability in action potential (AP) recorded in 363 patients both under SR and cAF conditions.

METHODS AND RESULTS

Human AP recordings in atrial trabeculae (n = 469) from SR and cAF patients were used to calibrate populations of computational SR and cAF atrial AP models. Three populations of over 2000 sampled models were generated, based on three different human atrial AP models. Experimental calibration selected populations of AP models yielding AP with morphology and duration in range with experimental recordings. Populations using the three original models can mimic variability in experimental AP in both SR and cAF, with median conductance values in SR for most ionic currents deviating less than 30% from their original peak values. All cAF populations show similar variations in G(K1), G(Kur) and G(to), consistent with AF-related remodeling as reported in experiments. In all SR and cAF model populations, inter-subject variability in I(K1) and I(NaK) underlies variability in APD90, variability in I(Kur), I(CaL) and I(NaK) modulates variability in APD50 and combined variability in Ito and I(Kur) determines variability in APD20. The large variability in human atrial AP triangulation is mostly determined by I(K1) and either I(NaK) or I(NaCa) depending on the model.

CONCLUSION

Experimentally-calibrated human atrial AP models populations mimic AP variability in SR and cAF patient recordings, and identify potential ionic determinants of inter-subject variability in human atrial AP duration and morphology in SR versus cAF.

The interrelations among stochastic pacing, stability, and memory in the heart

Low pacing variability in the heart has been clinically reported as a risk factor for lethal cardiac arrhythmias and arrhythmic death. In ia previous simulation study, we demonstrated that stochastic pacing sustains an antiarrhythmic effect by moderating the slope of the action potential duration (APD) restitution curve, by reducing the propensity of APD alternans, converting discordant to concordant alternans, and ultimately preventing wavebreaks. However, the dynamic mechanisms relating pacing stochasticity to tissue stability are not yet known. In this work, we develop a mathematical framework to describe the APD signal using an autoregressive stochastic model, and we establish the interrelations between stochastic pacing, cardiac memory, and cardiac stability, as manifested by the degree of APD alternans. Employing stability analysis tools, we show that increased stochasticity in the ventricular tissue activation sequence works to lower the maximal absolute eigenvalues of the stochastic model, thereby contributing to increased stability. We also show that the memory coefficients of the autoregressive model are modulated by pacing stochasticity in a nonlinear, biphasic way, so that for exceedingly high levels of pacing stochasticity, the antiarrhythmic effect is hampered by increasing APD variance. This work may contribute to establishment of an optimal antiarrhythmic pacing protocol in a future study.

Mathematical approaches to understanding and imaging atrial fibrillation: significance for mechanisms and management

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans. The mechanisms that govern AF initiation and persistence are highly complex, of dynamic nature, and involve interactions across multiple temporal and spatial scales in the atria. This article aims to review the mathematical modeling and computer simulation approaches to understanding AF mechanisms and aiding in its management. Various atrial modeling approaches are presented, with descriptions of the methodological basis and advancements in both lower-dimensional and realistic geometry models. A review of the most significant mechanistic insights made by atrial simulations is provided. The article showcases the contributions that atrial modeling and simulation have made not only to our understanding of the pathophysiology of atrial arrhythmias, but also to the development of AF management approaches. A summary of the future developments envisioned for the field of atrial simulation and modeling is also presented. The review contends that computational models of the atria assembled with data from clinical imaging modalities that incorporate electrophysiological and structural remodeling could become a first line of screening for new AF therapies and approaches, new diagnostic developments, and new methods for arrhythmia prevention.

Stochastic cardiac pacing increases ventricular electrical stability--a computational study

The ventricular tissue is activated in a stochastic rather than in a deterministic rhythm due to the inherent heart rate variability (HRV). Low HRV is a known predictor for arrhythmia events and traditionally is attributed to autonomic nervous system tone damage. Yet, there is no model that directly assesses the antiarrhythmic effect of pacing stochasticity per se. One-dimensional (1D) and two-dimensional (2D) human ventricular tissues were modeled, and both deterministic and stochastic pacing protocols were applied. Action potential duration restitution (APDR) and conduction velocity restitution (CVR) curves were generated and analyzed, and the propensity and characteristics of action potential duration (APD) alternans were investigated. In the 1D model, pacing stochasticity was found to sustain a moderating effect on the APDR curve by reducing its slope, rendering the tissue less arrhythmogenic. Moreover, stochasticity was found to be a significant antagonist to the development of concordant APD alternans. These effects were generally amplified with increased variability in the pacing cycle intervals. In addition, in the 2D tissue configuration, stochastic pacing exerted a protective antiarrhythmic effect by reducing the spatial APD heterogeneity and converting discordant APD alternans to concordant ones. These results suggest that high cardiac pacing stochasticity is likely to reduce the risk of cardiac arrhythmias in patients.

A quantitative comparison of the behavior of human ventricular cardiac electrophysiology models in tissue

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

Arrhythmogenic substrate for atrial fibrillation: insights from an integrative computational model of pulmonary veins

Mechanisms underlying the genesis of re-entrant substrate for atrial fibrillation (AF) in the pulmonary veins (PVs) and left atrium (LA) are not well understood. We develop a biophysically detailed computational model for the PVs and surrounding LA tissue. The model integrates canine PV and LA single cell electrophysiology with the respective 3D tissue geometry and fiber orientation reconstructed from micro-CT data. The model simulations demonstrate that a combination of tissue anisotropy and electrical heterogeneity between the PVs and LA causes a break-down of normal electrical excitation wave-fronts. This leads to the generation of a high-frequency re-entrant source near the PV sleeves. Evidence of such sources have been seen clinically in AF patients. In summary, our modeling results provide new insights into the arrhythmogenic mechanisms of re-entrant excitation waves underlying AF.

Proarrhythmic effect of blocking the small conductance calcium activated potassium channel in isolated canine left atrium

BACKGROUND

Small conductance calcium activated potassium (SKCa) channels are voltage insensitive and are activated by intracellular calcium. Genome-wide association studies revealed that a variant of SKCa is associated with lone atrial fibrillation in humans. Roles of SKCa in atrial arrhythmias remain unclear.

OBJECTIVE

To determine roles of SKCa in atrial arrhythmias.

METHODS

Optical mapping using the isolated canine left atrium was performed. The optical action potential duration (APD) and induction of arrhythmia were evaluated before and after the addition of specific SKCa blockers-apamin or UCL-1684.

RESULTS

SKCa blockade significantly increased APD₈₀ (188 ± 19 ms vs 147 ± 11 ms; P<.001). The pacing cycle length thresholds to induce 2:2 alternans, and wave breaks were prolonged by SKCa blockade. Increased APD heterogeneity was observed after the SKCa blockade, as measured by the difference between the maximum and the minimum APD (39 ± 4 ms vs 26 ± 5 ms; P<.05), by standard deviation (12.43 ± 2.36 ms vs 7.49 ± 1.47 ms; P<.001), or by coefficient of variation (6.68% ± 0.97% vs 4.90% ± 0.84%; P<.05). No arrhythmia was induced at baseline by an S1-S2 protocol. After SKCa blockade, 4 of 6 atria developed arrhythmia.

CONCLUSIONS

SKCa blockade promotes arrhythmia and prolongs the pacing cycle length threshold of 2:2 alternans and wave breaks in the canine left atrium. The proarrhythmic effect could be attributed to increased APD heterogeneity in the canine left atrium. This study provides supportive evidence of genome-wide association studies showing association of KCNN3 and lone atrial fibrillation.

Simulation of biatrial conduction via different pathways during sinus rhythm with a detailed human atrial model

In order to better understand biatrial conduction, investigate various conduction pathways, and compare the differences between isotropic and anisotropic conductions in human atria, we present a simulation study of biatrial conduction with known/assumed conduction pathways using a recently developed human atrial model. In addition to known pathways: (1) Bachmann's bundle (BB), (2) limbus of fossa ovalis (LFO), and (3) coronary sinus (CS), we also hypothesize that there exist two fast conduction bundles that connect the crista terminalis (CT), LFO, and CS. Our simulation demonstrates that use of these fast conduction bundles results in a conduction pattern consistent with experimental data. The comparison of isotropic and anisotropoic conductions in the BB case showed that the atrial working muscles had small effect on conduction time and conduction speed, although the conductivities assigned in anisotropic conduction were two to four times higher than the isotropic conduction. In conclusion, we suggest that the hypothesized intercaval bundles play a significant role in the biatrial conduction and that myofiber orientation has larger effects on the conduction system than the atrial working muscles. This study presents readers with new insights into human atrial conduction.

Benchmarking electrophysiological models of human atrial myocytes

Mathematical modeling of cardiac electrophysiology is an insightful method to investigate the underlying mechanisms responsible for arrhythmias such as atrial fibrillation (AF). In past years, five models of human atrial electrophysiology with different formulations of ionic currents, and consequently diverging properties, have been published. The aim of this work is to give an overview of strengths and weaknesses of these models depending on the purpose and the general requirements of simulations. Therefore, these models were systematically benchmarked with respect to general mathematical properties and their ability to reproduce certain electrophysiological phenomena, such as action potential (AP) alternans. To assess the models' ability to replicate modified properties of human myocytes and tissue in cardiac disease, electrical remodeling in chronic atrial fibrillation (cAF) was chosen as test case. The healthy and remodeled model variants were compared with experimental results in single-cell, 1D and 2D tissue simulations to investigate AP and restitution properties, as well as the initiation of reentrant circuits.

Application of micro-computed tomography with iodine staining to cardiac imaging, segmentation, and computational model development

Micro-computed tomography (micro-CT) has been widely used to generate high-resolution 3-D tissue images from small animals nondestructively, especially for mineralized skeletal tissues. However, its application to the analysis of soft cardiovascular tissues has been limited by poor inter-tissue contrast. Recent ex vivo studies have shown that contrast between muscular and connective tissue in micro-CT images can be enhanced by staining with iodine. In the present study, we apply this novel technique for imaging of cardiovascular structures in canine hearts. We optimize the method to obtain high-resolution X-ray micro-CT images of the canine atria and its distinctive regions-including the Bachmann's bundle, atrioventricular node, pulmonary arteries and veins-with clear inter-tissue contrast. The imaging results are used to reconstruct and segment the detailed 3-D geometry of the atria. Structure tensor analysis shows that the arrangement of atrial fibers can also be characterized using the enhanced micro-CT images, as iodine preferentially accumulates within the muscular fibers rather than in connective tissues. This novel technique can be particularly useful in nondestructive imaging of 3-D cardiac architectures from large animals and humans, due to the combination of relatively high speed ( ~ 1 h/per scan of the large canine heart) and high voxel resolution (36 μm) provided. In summary, contrast micro-CT facilitates fast and nondestructive imaging and segmenting of detailed 3-D cardiovascular geometries, as well as measuring fiber orientation, which are crucial in constructing biophysically detailed computational cardiac models.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

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.

Nonlinear dynamics in cardiology

The dynamics of many cardiac arrhythmias, as well as the nature of transitions between different heart rhythms, have long been considered evidence of nonlinear phenomena playing a direct role in cardiac arrhythmogenesis. In most types of cardiac disease, the pathology develops slowly and gradually, often over many years. In contrast, arrhythmias often occur suddenly. In nonlinear systems, sudden changes in qualitative dynamics can, counterintuitively, result from a gradual change in a system parameter-this is known as a bifurcation. Here, we review how nonlinearities in cardiac electrophysiology influence normal and abnormal rhythms and how bifurcations change the dynamics. In particular, we focus on the many recent developments in computational modeling at the cellular level that are focused on intracellular calcium dynamics. We discuss two areas where recent experimental and modeling work has suggested the importance of nonlinearities in calcium dynamics: repolarization alternans and pacemaker cell automaticity.

Exploiting mathematical models to illuminate electrophysiological variability between individuals

Across individuals within a population, several levels of variability are observed, from the differential expression of ion channels at the molecular level, to the various action potential morphologies observed at the cellular level, to divergent responses to drugs at the organismal level. However, the limited ability of experiments to probe complex interactions between components has hitherto hindered our understanding of the factors that cause a range of behaviours within a population. Variability is a challenging issue that is encountered in all physiological disciplines, but recent work suggests that novel methods for analysing mathematical models can assist in illuminating its causes. In this review, we discuss mathematical modelling studies in cardiac electrophysiology and neuroscience that have enhanced our understanding of variability in a number of key areas. Specifically, we discuss parameter sensitivity analysis techniques that may be applied to generate quantitative predictions based on considering behaviours within a population of models, thereby providing novel insight into variability. Our discussion focuses on four issues that have benefited from the utilization of these methods: (1) the comparison of different electrophysiological models of cardiac myocytes, (2) the determination of the individual contributions of different molecular changes in complex disease phenotypes, (3) the identification of the factors responsible for the variable response to drugs, and (4) the constraining of free parameters in electrophysiological models of heart cells. Together, the studies that we discuss suggest that rigorous analyses of mathematical models can generate quantitative predictions regarding how molecular-level variations contribute to functional differences between experimental samples. These strategies may be applicable not just in cardiac electrophysiology, but in a wide range of disciplines.

The Na+/K+ pump is an important modulator of refractoriness and rotor dynamics in human atrial tissue

Pharmacological treatment of atrial fibrillation (AF) exhibits limited efficacy. Further developments require a comprehensive characterization of ionic modulators of electrophysiology in human atria. Our aim is to systematically investigate the relative importance of ionic properties in modulating excitability, refractoriness, and rotor dynamics in human atria before and after AF-related electrical remodeling (AFER). Computer simulations of single cell and tissue atrial electrophysiology were conducted using two human atrial action potential (AP) models. Changes in AP, refractory period (RP), conduction velocity (CV), and rotor dynamics caused by alterations in key properties of all atrial ionic currents were characterized before and after AFER. Results show that the investigated human atrial electrophysiological properties are primarily modulated by maximal value of Na(+)/K(+) pump current (G(NaK)) as well as conductances of inward rectifier potassium current (G(K1)) and fast inward sodium current (G(Na)). G(NaK) plays a fundamental role through both electrogenic and homeostatic modulation of AP duration (APD), APD restitution, RP, and reentrant dominant frequency (DF). G(K1) controls DF through modulation of AP, APD restitution, RP, and CV. G(Na) is key in determining DF through alteration of CV and RP, particularly in AFER. Changes in ionic currents have qualitatively similar effects in control and AFER, but effects are smaller in AFER. The systematic analysis conducted in this study unravels the important role of the Na(+)/K(+) pump current in determining human atrial electrophysiology.

Effects of electrical and structural remodeling on atrial fibrillation maintenance: a simulation study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength--electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely.

Heart rate variability effect on the myocyte action potential duration restitution: insights from switched systems theory

The physiological heart rate presents a stochastic behavior known as heart rate variability (HRV). In this framework the influence of HRV on the action potential duration (APD) of the atrial myocyte is analyzed in a computer model. We have found that introducing HRV into the myocyte action potential model decreases the APD of the extra beat S2 in an S1-S2 protocol compared to constant heart rate. A possible theoretical explanation for this is also presented and is derived from switched systems theory. It is suggested to consider the myocyte action potential phase 4 and phase 2 as two operation modes of a switching system and analyze the stability of switching between them. Since random switching is known to have a stabilization effect on a switching system, this might explain why HRV has a stabilization effect on the myocyte APD restitution. Implications of this finding include reduced system stability for conditions with low HRV. A possible application for this phenomenon regards artificial pacemakers, where a preset added HRV is predicted to reduce susceptibility to arrhythmias.

Teaching cardiac electrophysiology modeling to undergraduate students: laboratory exercises and GPU programming for the study of arrhythmias and spiral wave dynamics

As part of a 3-wk intersession workshop funded by a National Science Foundation Expeditions in Computing award, 15 undergraduate students from the City University of New York(1) collaborated on a study aimed at characterizing the voltage dynamics and arrhythmogenic behavior of cardiac cells for a broad range of physiologically relevant conditions using an in silico model. The primary goal of the workshop was to cultivate student interest in computational modeling and analysis of complex systems by introducing them through lectures and laboratory activities to current research in cardiac modeling and by engaging them in a hands-on research experience. The success of the workshop lay in the exposure of the students to active researchers and experts in their fields, the use of hands-on activities to communicate important concepts, active engagement of the students in research, and explanations of the significance of results as the students generated them. The workshop content addressed how spiral waves of electrical activity are initiated in the heart and how different parameter values affect the dynamics of these reentrant waves. Spiral waves are clinically associated with tachycardia, when the waves remain stable, and with fibrillation, when the waves exhibit breakup. All in silico experiments were conducted by simulating a mathematical model of cardiac cells on graphics processing units instead of the standard central processing units of desktop computers. This approach decreased the run time for each simulation to almost real time, thereby allowing the students to quickly analyze and characterize the simulated arrhythmias. Results from these simulations, as well as some of the background and methodology taught during the workshop, is presented in this article along with the programming code and the explanations of simulation results in an effort to allow other teachers and students to perform their own demonstrations, simulations, and studies.

Verification of cardiac tissue electrophysiology simulators using an N-version benchmark

Ongoing developments in cardiac modelling have resulted, in particular, in the development of advanced and increasingly complex computational frameworks for simulating cardiac tissue electrophysiology. The goal of these simulations is often to represent the detailed physiology and pathologies of the heart using codes that exploit the computational potential of high-performance computing architectures. These developments have rapidly progressed the simulation capacity of cardiac virtual physiological human style models; however, they have also made it increasingly challenging to verify that a given code provides a faithful representation of the purported governing equations and corresponding solution techniques. This study provides the first cardiac tissue electrophysiology simulation benchmark to allow these codes to be verified. The benchmark was successfully evaluated on 11 simulation platforms to generate a consensus gold-standard converged solution. The benchmark definition in combination with the gold-standard solution can now be used to verify new simulation codes and numerical methods in the future.