3D virtual human atria: A computational platform for studying clinical atrial fibrillation

Lessons from computer simulations of ablation of atrial fibrillation

This paper reviews the simulations of catheter ablation in computer models of the atria, from the first attempts to the most recent anatomical models. It describes how postulated substrates of atrial fibrillation can be incorporated into mathematical models, how modelling studies can be designed to test ablation strategies, what their current trade-offs and limitations are, and what clinically relevant lessons can be learnt from these simulations. Drawing a parallel between clinical and modelling studies, six ablation targets are considered: pulmonary vein isolation, linear ablation, ectopic foci, complex fractionated atrial electrogram, rotors and ganglionated plexi. The examples presented for each ablation target illustrate a major advantage of computer models, the ability to identify why a therapy is successful or not in a given atrial fibrillation substrate. The integration of pathophysiological data to create detailed models of arrhythmogenic substrates is expected to solidify the understanding of ablation mechanisms and to provide theoretical arguments supporting substrate-specific ablation strategies.

Improving cardiomyocyte model fidelity and utility via dynamic electrophysiology protocols and optimization algorithms

Mathematical models of cardiac electrophysiology are instrumental in determining mechanisms of cardiac arrhythmias. However, the foundation of a realistic multiscale heart model is only as strong as the underlying cell model. While there have been myriad advances in the improvement of cellular-level models, the identification of model parameters, such as ion channel conductances and rate constants, remains a challenging problem. The primary limitations to this process include: (1) such parameters are usually estimated from data recorded using standard electrophysiology voltage-clamp protocols that have not been developed with model building in mind, and (2) model parameters are typically tuned manually to subjectively match a desired output. Over the last decade, methods aimed at overcoming these disadvantages have emerged. These approaches include the use of optimization or fitting tools for parameter estimation and incorporating more extensive data for output matching. Here, we review recent advances in parameter estimation for cardiomyocyte models, focusing on the use of more complex electrophysiology protocols and global search heuristics. We also discuss future applications of such parameter identification, including development of cell-specific and patient-specific mathematical models to investigate arrhythmia mechanisms and predict therapy strategies.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model

Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias.

Effects of Persistent Atrial Fibrillation-Induced Electrical Remodeling on Atrial Electro-Mechanics - Insights from a 3D Model of the Human Atria

Aims

Atrial stunning, a loss of atrial mechanical contraction, can occur following a successful cardioversion. It is hypothesized that persistent atrial fibrillation-induced electrical remodeling (AFER) on atrial electrophysiology may be responsible for such impaired atrial mechanics. This simulation study aimed to investigate the effects of AFER on atrial electro-mechanics.

Methods and Results

A 3D electromechanical model of the human atria was developed to investigate the effects of AFER on atrial electro-mechanics. Simulations were carried out in 3 conditions for 4 states: (i) the control condition, representing the normal tissue (state 1) and the tissue 2–3 months after cardioversion (state 2) when the atrial tissue recovers its electrophysiological properties after completion of reverse electrophysiological remodelling; (ii) AFER-SR condition for AF-remodeled tissue with normal sinus rhythm (SR) (state 3); and (iii) AFER-AF condition for AF-remodeled tissue with re-entrant excitation waves (state 4). Our results indicate that at the cellular level, AFER (states 3 & 4) abbreviated action potentials and reduced the Ca²⁺ content in the sarcoplasmic reticulum, resulting in a reduced amplitude of the intracellular Ca²⁺ transient leading to decreased cell active force and cell shortening as compared to the control condition (states 1 & 2). Consequently at the whole organ level, atrial contraction in AFER-SR condition (state 3) was dramatically reduced. In the AFER-AF condition (state 4) atrial contraction was almost abolished.

Conclusions

This study provides novel insights into understanding atrial electro-mechanics illustrating that AFER impairs atrial contraction due to reduced intracellular Ca²⁺ transients.

Detailed Anatomical and Electrophysiological Models of Human Atria and Torso for the Simulation of Atrial Activation

Atrial arrhythmias, and specifically atrial fibrillation (AF), induce rapid and irregular activation patterns that appear on the torso surface as abnormal P-waves in electrocardiograms and body surface potential maps (BSPM). In recent years both P-waves and the BSPM have been used to identify the mechanisms underlying AF, such as localizing ectopic foci or high-frequency rotors. However, the relationship between the activation of the different areas of the atria and the characteristics of the BSPM and P-wave signals are still far from being completely understood. In this work we developed a multi-scale framework, which combines a highly-detailed 3D atrial model and a torso model to study the relationship between atrial activation and surface signals in sinus rhythm. Using this multi scale model, it was revealed that the best places for recording P-waves are the frontal upper right and the frontal and rear left quadrants of the torso. Our results also suggest that only nine regions (of the twenty-one structures in which the atrial surface was divided) make a significant contribution to the BSPM and determine the main P-wave characteristics.

A 2D Electromechanical Model of Human Atrial Tissue Using the Discrete Element Method

Cardiac tissue is a syncytium of coupled cells with pronounced intrinsic discrete nature. Previous models of cardiac electromechanics often ignore such discrete properties and treat cardiac tissue as a continuous medium, which has fundamental limitations. In the present study, we introduce a 2D electromechanical model for human atrial tissue based on the discrete element method (DEM). In the model, single-cell dynamics are governed by strongly coupling the electrophysiological model of Courtemanche et al. to the myofilament model of Rice et al. with two-way feedbacks. Each cell is treated as a viscoelastic body, which is physically represented by a clump of nine particles. Cell aggregations are arranged so that the anisotropic nature of cardiac tissue due to fibre orientations can be modelled. Each cell is electrically coupled to neighbouring cells, allowing excitation waves to propagate through the tissue. Cell-to-cell mechanical interactions are modelled using a linear contact bond model in DEM. By coupling cardiac electrophysiology with mechanics via the intracellular Ca²⁺ concentration, the DEM model successfully simulates the conduction of cardiac electrical waves and the tissue's corresponding mechanical contractions. The developed DEM model is numerically stable and provides a powerful method for studying the electromechanical coupling problem in the heart.

Optimization of catheter ablation of atrial fibrillation: insights gained from clinically-derived computer models

Atrial fibrillation (AF) is the most common heart rhythm disturbance, and its treatment is an increasing economic burden on the health care system. Despite recent intense clinical, experimental and basic research activity, the treatment of AF with current antiarrhythmic drugs and catheter/surgical therapies remains limited. Radiofrequency catheter ablation (RFCA) is widely used to treat patients with AF. Current clinical ablation strategies are largely based on atrial anatomy and/or substrate detected using different approaches, and they vary from one clinical center to another. The nature of clinical ablation leads to ambiguity regarding the optimal patient personalization of the therapy partly due to the fact that each empirical configuration of ablation lines made in a patient is irreversible during one ablation procedure. To investigate optimized ablation lesion line sets, in silico experimentation is an ideal solution. 3D computer models give us a unique advantage to plan and assess the effectiveness of different ablation strategies before and during RFCA. Reliability of in silico assessment is ensured by inclusion of accurate 3D atrial geometry, realistic fiber orientation, accurate fibrosis distribution and cellular kinetics; however, most of this detailed information in the current computer models is extrapolated from animal models and not from the human heart. The predictive power of computer models will increase as they are validated with human experimental and clinical data. To make the most from a computer model, one needs to develop 3D computer models based on the same functionally and structurally mapped intact human atria with high spatial resolution. The purpose of this review paper is to summarize recent developments in clinically-derived computer models and the clinical insights they provide for catheter ablation.

A new algorithm to diagnose atrial ectopic origin from multi lead ECG systems--insights from 3D virtual human atria and torso

Rapid atrial arrhythmias such as atrial fibrillation (AF) predispose to ventricular arrhythmias, sudden cardiac death and stroke. Identifying the origin of atrial ectopic activity from the electrocardiogram (ECG) can help to diagnose the early onset of AF in a cost-effective manner. The complex and rapid atrial electrical activity during AF makes it difficult to obtain detailed information on atrial activation using the standard 12-lead ECG alone. Compared to conventional 12-lead ECG, more detailed ECG lead configurations may provide further information about spatio-temporal dynamics of the body surface potential (BSP) during atrial excitation. We apply a recently developed 3D human atrial model to simulate electrical activity during normal sinus rhythm and ectopic pacing. The atrial model is placed into a newly developed torso model which considers the presence of the lungs, liver and spinal cord. A boundary element method is used to compute the BSP resulting from atrial excitation. Elements of the torso mesh corresponding to the locations of the placement of the electrodes in the standard 12-lead and a more detailed 64-lead ECG configuration were selected. The ectopic focal activity was simulated at various origins across all the different regions of the atria. Simulated BSP maps during normal atrial excitation (i.e. sinoatrial node excitation) were compared to those observed experimentally (obtained from the 64-lead ECG system), showing a strong agreement between the evolution in time of the simulated and experimental data in the P-wave morphology of the ECG and dipole evolution. An algorithm to obtain the location of the stimulus from a 64-lead ECG system was developed. The algorithm presented had a success rate of 93%, meaning that it correctly identified the origin of atrial focus in 75/80 simulations, and involved a general approach relevant to any multi-lead ECG system. This represents a significant improvement over previously developed algorithms.

Three-dimensional computer model of the right atrium including the sinoatrial and atrioventricular nodes predicts classical nodal behaviours

The aim of the study was to develop a three-dimensional (3D) anatomically-detailed model of the rabbit right atrium containing the sinoatrial and atrioventricular nodes to study the electrophysiology of the nodes. A model was generated based on 3D images of a rabbit heart (atria and part of ventricles), obtained using high-resolution magnetic resonance imaging. Segmentation was carried out semi-manually. A 3D right atrium array model (∼3.16 million elements), including eighteen objects, was constructed. For description of cellular electrophysiology, the Rogers-modified FitzHugh-Nagumo model was further modified to allow control of the major characteristics of the action potential with relatively low computational resource requirements. Model parameters were chosen to simulate the action potentials in the sinoatrial node, atrial muscle, inferior nodal extension and penetrating bundle. The block zone was simulated as passive tissue. The sinoatrial node, crista terminalis, main branch and roof bundle were considered as anisotropic. We have simulated normal and abnormal electrophysiology of the two nodes. In accordance with experimental findings: (i) during sinus rhythm, conduction occurs down the interatrial septum and into the atrioventricular node via the fast pathway (conduction down the crista terminalis and into the atrioventricular node via the slow pathway is slower); (ii) during atrial fibrillation, the sinoatrial node is protected from overdrive by its long refractory period; and (iii) during atrial fibrillation, the atrioventricular node reduces the frequency of action potentials reaching the ventricles. The model is able to simulate ventricular echo beats. In summary, a 3D anatomical model of the right atrium containing the cardiac conduction system is able to simulate a wide range of classical nodal behaviours.

Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model

Aims

Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and sustenance of wave sources in a model of the PV tissue.

Methods and results

A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced electrical remodelling are introduced based on published experimental data, as changes of ion channel currents (I_CaL, I_K1, I_to, and I_Kur), the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5–12 ms) and anisotropy (conduction velocities of 0.75–1.25 and 0.21–0.31 m/s along and transverse to atrial fibres) can results in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wavelength (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulnerability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (I_K1 and I_Kur blockers) can also increase the heterogeneity (AP gradients of 26–27 ms) and the likelihood of wave breaks.

Conclusion

Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrangement enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological interventions.

Computational modeling of inhibition of voltage-gated Ca channels: identification of different effects on uterine and cardiac action potentials

The uterus and heart share the important physiological feature whereby contractile activation of the muscle tissue is regulated by the generation of periodic, spontaneous electrical action potentials (APs). Preterm birth arising from premature uterine contractions is a major complication of pregnancy and there remains a need to pursue avenues of research that facilitate the use of drugs, tocolytics, to limit these inappropriate contractions without deleterious actions on cardiac electrical excitation. A novel approach is to make use of mathematical models of uterine and cardiac APs, which incorporate many ionic currents contributing to the AP forms, and test the cell-specific responses to interventions. We have used three such models-of uterine smooth muscle cells (USMC), cardiac sinoatrial node cells (SAN), and ventricular cells-to investigate the relative effects of reducing two important voltage-gated Ca currents-the L-type (ICaL) and T-type (ICaT) Ca currents. Reduction of ICaL (10%) alone, or ICaT (40%) alone, blunted USMC APs with little effect on ventricular APs and only mild effects on SAN activity. Larger reductions in either current further attenuated the USMC APs but with also greater effects on SAN APs. Encouragingly, a combination of ICaL and ICaT reduction did blunt USMC APs as intended with little detriment to APs of either cardiac cell type. Subsequent overlapping maps of ICaL and ICaT inhibition profiles from each model revealed a range of combined reductions of ICaL and ICaT over which an appreciable diminution of USMC APs could be achieved with no deleterious action on cardiac SAN or ventricular APs. This novel approach illustrates the potential for computational biology to inform us of possible uterine and cardiac cell-specific mechanisms. Incorporating such computational approaches in future studies directed at designing new, or repurposing existing, tocolytics will be beneficial for establishing a desired uterine specificity of action.

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.

Towards computational modelling of the human foetal electrocardiogram: normal sinus rhythm and congenital heart block

AIMS

We aim to engineer a computational model of propagation during normal sinus rhythm in the foetal human heart, by modifying models for adult cardiac tissue to match foetal electrocardiogram (fECG) characteristics. The model will be partially validated by fECG data, and applied to explore possible mechanisms of arrhythmogenesis in the foetal heart.

METHODS AND RESULTS

Foetal electrocardiograms have been recorded during pregnancy, with P- and T-waves, and the QRS complex, identified by averaging and signal processing. Intervals of the fECG are extracted and used to modify currently available human adult cardiomyocyte models. RR intervals inform models of the pacemaking cells by constraining their rate, the QT interval and its rate dependence constrain models of ventricular cells, and the width of the P-wave, the QR and PR intervals constrain propagation times, conduction velocities, and intercellular coupling. These cell models are coupled into a one-dimensional (1D) model of propagation during normal sinus rhythm in the human foetal heart. We constructed a modular, heterogeneous 1D model for propagation in the foetal heart, and predicted the effects of reduction in L-type Ca(++) current. These include bradycardia and atrioventricular conduction blocks. These may account quantitatively for congenital heart block produced by positive IgG antibodies.

CONCLUSION

The fECG can be interpreted mechanistically and quantitatively by using a simple computational model for propagation. After further validation, by clinical recordings of the fECG and the electrophysiological experiments on foetal cardiac cells and tissues, the model may be used to predict the effects of maternally administered pharmaceuticals on the fECG.

A novel computational sheep atria model for the study of atrial fibrillation

Sheep are often used as animal models for experimental studies into the underlying mechanisms of cardiac arrhythmias. Previous studies have shown that biophysically detailed computer models of the heart provide a powerful alternative to experimental animal models for underpinning such mechanisms. In this study, we have developed a family of mathematical models for the electrical action potentials of various sheep atrial cell types. The developed cell models were then incorporated into a three-dimensional anatomical model of the sheep atria, which was recently reconstructed and segmented based on anatomical features within different regions. This created a novel biophysically detailed computational model of the three-dimensional sheep atria. Using the model, we then investigated the mechanisms by which paroxysmal rapid focal activity in the pulmonary veins can transit to sustained atrial fibrillation. It was found that the anisotropic property of the atria arising from the fibre structure plays an important role in facilitating the development of fibrillatory atrial excitation waves, and the electrical heterogeneity plays an important role in its initiation.

Antenatal architecture and activity of the human heart

We construct the components for a family of computational models of the electrophysiology of the human foetal heart from 60 days gestational age (DGA) to full term. This requires both cell excitation models that reconstruct the myocyte action potentials, and datasets of cardiac geometry and architecture. Fast low-angle shot and diffusion tensor magnetic resonance imaging (DT-MRI) of foetal hearts provides cardiac geometry with voxel resolution of approximately 100 µm. DT-MRI measures the relative diffusion of protons and provides a measure of the average intravoxel myocyte orientation, and the orientation of any higher order orthotropic organization of the tissue. Such orthotropic organization in the adult mammalian heart has been identified with myocardial sheets and cleavage planes between them. During gestation, the architecture of the human ventricular wall changes from being irregular and isotropic at 100 DGA to an anisotropic and orthotropic architecture by 140 DGA, when it has the smooth, approximately 120° transmural change in myocyte orientation that is characteristic of the adult mammalian ventricle. The DT obtained from DT-MRI provides the conductivity tensor that determines the spread of potential within computational models of cardiac tissue electrophysiology. The foetal electrocardiogram (fECG) can be recorded from approximately 60 DGA, and RR, PR and QT intervals between the P, R, Q and T waves of the fECG can be extracted by averaging from approximately 90 DGA. The RR intervals provide a measure of the pacemaker rate, the QT intervals an index of ventricular action potential duration, and its rate-dependence, and so these intervals constrain and inform models of cell electrophysiology. The parameters of models of adult human sinostrial node and ventricular cells that are based on adult cell electrophysiology and tissue molecular mapping have been modified to construct preliminary models of foetal cell electrophysiology, which reproduce these intervals from fECG recordings. The PR and QR intervals provide an index of conduction times, and hence propagation velocities (approx. 1-10 cm s(-1), increasing during gestation) and so inform models of tissue electrophysiology. Although the developing foetal heart is small and the cells are weakly coupled, it can support potentially lethal re-entrant arrhythmia.

Heterogeneous and anisotropic integrative model of pulmonary veins: computational study of arrhythmogenic substrate for atrial fibrillation

Mechanisms underlying the genesis of re-entrant substrate for the most common cardiac arrhythmia, atrial fibrillation (AF), are not well understood. In this study, we develop a multi-scale three-dimensional computational model that integrates cellular electrophysiology of the left atrium (LA) and pulmonary veins (PVs) with the respective tissue geometry and fibre orientation. The latter is reconstructed in unique detail from high-resolution (approx. 70 μm) contrast micro-computed tomography data. The model is used to explore the mechanisms of re-entry initiation and sustenance in the PV region, regarded as the primary source of high-frequency electrical activity in AF. Simulations of the three-dimensional model demonstrate that an initial break-down of normal electrical excitation wave-fronts can be caused by the electrical heterogeneity between the PVs and LA. High tissue anisotropy is then responsible for the slow conduction and generation of a re-entrant circuit near the PVs. Evidence of such circuits has been seen clinically in AF patients. Our computational study suggests that primarily the combination of electrical heterogeneity and conduction anisotropy between the PVs and LA tissues leads to the generation of a high-frequency (approx. 10 Hz) re-entrant source near the PV sleeves, thus providing new insights into the arrhythmogenic mechanisms of excitation waves underlying AF.

Rotor termination is critically dependent on kinetic properties of I kur inhibitors in an in silico model of chronic atrial fibrillation

Inhibition of the atrial ultra-rapid delayed rectifier potassium current (I_Kur) represents a promising therapeutic strategy in the therapy of atrial fibrillation. However, experimental and clinical data on the antiarrhythmic efficacy remain controversial. We tested the hypothesis that antiarrhythmic effects of I_Kur inhibitors are dependent on kinetic properties of channel blockade. A mathematical description of I_Kur blockade was introduced into Courtemanche-Ramirez-Nattel models of normal and remodeled atrial electrophysiology. Effects of five model compounds with different kinetic properties were analyzed. Although a reduction of dominant frequencies could be observed in two dimensional tissue simulations for all compounds, a reduction of spiral wave activity could be only be detected in two cases. We found that an increase of the percent area of refractory tissue due to a prolongation of the wavelength seems to be particularly important. By automatic tracking of spiral tip movement we find that increased refractoriness resulted in rotor extinction caused by an increased spiral-tip meandering. We show that antiarrhythmic effects of I_Kur inhibitors are dependent on kinetic properties of blockade. We find that an increase of the percent area of refractory tissue is the underlying mechanism for an increased spiral-tip meandering, resulting in the extinction of re-entrant circuits.

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.

Development of biophysically detailed electrophysiological models for pacemaking and non-pacemaking human pulmonary vein cardiomyocytes

Ectopic foci originating from the pulmonary veins (PVs) have been suggested as the underlying cause for generating atrial arrhythmias that include atrial fibrillation (AF). Recent experimental findings indicate two types of PV cells: pacemaking and non-pacemaking. In this study, we have developed two mathematical models for human PV cardiomyocytes with and without pacemaking activities. The models were reconstructed by modifying an existing model of the human right atrium to incorporate extant experimental data on the electrical differences between the two cell types. Differences in their action potential (AP) profiles and automaticity were reproduced by the models, which can be attributed to the observed differences in the current densities of I(NCX), I(to), I(Na) and I(Ca-L), as well as the difference in the channel kinetics of I(Ca-L) and inclusion of the I(f) and I(Ca-T) currents in the pacemaking cells. The developed models provide a useful tool suitable for studying the substrates for generating AF.

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.

Structure, function and clinical relevance of the cardiac conduction system, including the atrioventricular ring and outflow tract tissues

It is now over 100years since the discovery of the cardiac conduction system, consisting of three main parts, the sinus node, the atrioventricular node and the His-Purkinje system. The system is vital for the initiation and coordination of the heartbeat. Over the last decade, immense strides have been made in our understanding of the cardiac conduction system and these recent developments are reviewed here. It has been shown that the system has a unique embryological origin, distinct from that of the working myocardium, and is more extensive than originally thought with additional structures: atrioventricular rings, a third node (so called retroaortic node) and pulmonary and aortic sleeves. It has been shown that the expression of ion channels, intracellular Ca(2+)-handling proteins and gap junction channels in the system is specialised (different from that in the ordinary working myocardium), but appropriate to explain the functioning of the system, although there is continued debate concerning the ionic basis of pacemaking. We are beginning to understand the mechanisms (fibrosis and remodelling of ion channels and related proteins) responsible for dysfunction of the system (bradycardia, heart block and bundle branch block) associated with atrial fibrillation and heart failure and even athletic training. Equally, we are beginning to appreciate how naturally occurring mutations in ion channels cause congenital cardiac conduction system dysfunction. Finally, current therapies, the status of a new therapeutic strategy (use of a specific heart rate lowering drug) and a potential new therapeutic strategy (biopacemaking) are reviewed.

In-silico modeling of atrial repolarization in normal and atrial fibrillation remodeled state

Atrial fibrillation (AF) is the most common cardiac arrhythmia, and the total number of AF patients is constantly increasing. The mechanisms leading to and sustaining AF are not completely understood yet. Heterogeneities in atrial electrophysiology seem to play an important role in this context. Although some heterogeneities have been used in in-silico human atrial modeling studies, they have not been thoroughly investigated. In this study, the original electrophysiological (EP) models of Courtemanche et al., Nygren et al. and Maleckar et al. were adjusted to reproduce action potentials in 13 atrial regions. The parameter sets were validated against experimental action potential duration data and ECG data from patients with AV block. The use of the heterogeneous EP model led to a more synchronized repolarization sequence in a variety of 3D atrial anatomical models. Combination of the heterogeneous EP model with a model of persistent AF-remodeled electrophysiology led to a drastic change in cell electrophysiology. Simulated Ta-waves were significantly shorter under the remodeling. The heterogeneities in cell electrophysiology explain the previously observed Ta-wave effects. The results mark an important step toward the reliable simulation of the atrial repolarization sequence, give a deeper understanding of the mechanism of atrial repolarization and enable further clinical investigations.

Optimisation of ionic models to fit tissue action potentials: application to 3D atrial modelling

A 3D model of atrial electrical activity has been developed with spatially heterogeneous electrophysiological properties. The atrial geometry, reconstructed from the male Visible Human dataset, included gross anatomical features such as the central and peripheral sinoatrial node (SAN), intra-atrial connections, pulmonary veins, inferior and superior vena cava, and the coronary sinus. Membrane potentials of myocytes from spontaneously active or electrically paced in vitro rabbit cardiac tissue preparations were recorded using intracellular glass microelectrodes. Action potentials of central and peripheral SAN, right and left atrial, and pulmonary vein myocytes were each fitted using a generic ionic model having three phenomenological ionic current components: one time-dependent inward, one time-dependent outward, and one leakage current. To bridge the gap between the single-cell ionic models and the gross electrical behaviour of the 3D whole-atrial model, a simplified 2D tissue disc with heterogeneous regions was optimised to arrive at parameters for each cell type under electrotonic load. Parameters were then incorporated into the 3D atrial model, which as a result exhibited a spontaneously active SAN able to rhythmically excite the atria. The tissue-based optimisation of ionic models and the modelling process outlined are generic and applicable to image-based computer reconstruction and simulation of excitable tissue.

Drift laws for spiral waves on curved anisotropic surfaces

Rotating spiral waves organize spatial patterns in chemical, physical, and biological excitable systems. Factors affecting their dynamics, such as spatiotemporal drift, are of great interest for particular applications. Here, we propose a quantitative description for spiral wave dynamics on curved surfaces which shows that for a wide class of systems, including the Belousov-Zhabotinsky reaction and anisotropic cardiac tissue, the Ricci curvature scalar of the surface is the main determinant of spiral wave drift. The theory provides explicit equations for spiral wave drift direction, drift velocity, and the period of rotation. Depending on the parameters, the drift can be directed to the regions of either maximal or minimal Ricci scalar curvature, which was verified by direct numerical simulations.

Three-dimensional impulse propagation in myocardium: arrhythmogenic mechanisms at the tissue level

Impulse propagation in the heart depends on the excitability of individual cardiomyocytes, impulse transmission between adjacent myocytes, and the 3-dimensional arrangement of those cells. Here, we review the role of each of these factors in normal and aberrant cardiac electric activation, with particular emphasis on the effects of 3-dimensional myocyte architecture at the tissue scale. The analysis draws on findings from in vivo and in vitro experiments, as well as biophysically based computer models that have been used to integrate and interpret these experimental data. It indicates that discontinuous arrangement of myocytes and extracellular connective tissue at the tissue scale can give rise to current source-to-sink mismatch, spatiotemporal distribution of refractoriness, and rate-sensitive electric instability, which contribute to the initiation and maintenance of reentrant cardiac arrhythmia. This exacerbates the risk of rhythm disturbance associated with heart disease. We conclude that structure-based, multiscale computer models that incorporate accurate information about local cellular electric activity provide a powerful platform for investigating the basis of reentrant cardiac arrhythmia. However, it is important that these models capture key features of structure and related electric function at the tissue scale.

Pro-arrhythmogenic effects of atrial fibrillation-induced electrical remodelling: insights from the three-dimensional virtual human atria

Chronic atrial fibrillation (AF) is associated with structural and electrical remodelling in the atria, which are associated with a high recurrence of AF. Through biophysically detailed computer modelling, this study investigated mechanisms by which AF-induced electrical remodelling promotes and perpetuates AF. A family of Courtemanche–Ramirez–Nattel variant models of human atrial cell action potentials (APs), taking into account of intrinsic atrial electrophysiological properties, was modified to incorporate various experimental data sets on AF-induced changes of major ionic channel currents (I_CaL, I_Kur, I_to, I_K1, I_Ks, I_NaCa) and on intracellular Ca²⁺ handling. The single cell models for control and AF-remodelled conditions were incorporated into multicellular three-dimensional (3D) atrial tissue models. Effects of the AF-induced electrical remodelling were quantified as the changes of AP profile, AP duration (APD) and its dispersion across the atria, and the vulnerability of atrial tissue to the initiation of re-entry. The dynamic behaviour of re-entrant excitation waves in the 3D models was characterised. In our simulations, AF-induced electrical remodelling abbreviated atrial APD non-uniformly across the atria; this resulted in relatively short APDs co-existing with marked regional differences in the APD at junctions of the crista terminalis/pectinate muscle, pulmonary veins/left atrium. As a result, the measured tissue vulnerability to re-entry initiation at these tissue junctions was increased. The AF-induced electrical remodelling also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained re-entry. Under the AF-remodelled condition, re-entrant scroll waves in the 3D model degenerated into persistent and erratic wavelets, leading to fibrillation. In conclusion, realistic 3D atrial tissue models indicate that AF-induced electrical remodelling produces regionally heterogeneous and shortened APD; these respectively facilitate initiation and maintenance of re-entrant excitation waves.

A three-dimensional finite element model of human atrial anatomy: new methods for cubic Hermite meshes with extraordinary vertices

High-order cubic Hermite finite elements have been valuable in modeling cardiac geometry, fiber orientations, biomechanics, and electrophysiology, but their use in solving three-dimensional problems has been limited to ventricular models with simple topologies. Here, we utilized a subdivision surface scheme and derived a generalization of the "local-to-global" derivative mapping scheme of cubic Hermite finite elements to construct bicubic and tricubic Hermite models of the human atria with extraordinary vertices from computed tomography images of a patient with atrial fibrillation. To an accuracy of 0.6 mm, we were able to capture the left atrial geometry with only 142 bicubic Hermite finite elements, and the right atrial geometry with only 90. The left and right atrial bicubic Hermite meshes were G1 continuous everywhere except in the one-neighborhood of extraordinary vertices, where the mean dot products of normals at adjacent elements were 0.928 and 0.925. We also constructed two biatrial tricubic Hermite models and defined fiber orientation fields in agreement with diagrammatic data from the literature using only 42 angle parameters. The meshes all have good quality metrics, uniform element sizes, and elements with aspect ratios near unity, and are shared with the public. These new methods will allow for more compact and efficient patient-specific models of human atrial and whole heart physiology.

A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship

The most common sustained cardiac arrhythmias in humans are atrial tachyarrhythmias, mainly atrial fibrillation. Areas of complex fractionated atrial electrograms and high dominant frequency have been proposed as critical regions for maintaining atrial fibrillation; however, there is a paucity of data on the relationship between the characteristics of electrograms and the propagation pattern underlying them. In this study, a realistic 3D computer model of the human atria has been developed to investigate this relationship. The model includes a realistic geometry with fiber orientation, anisotropic conductivity and electrophysiological heterogeneity. We simulated different tachyarrhythmic episodes applying both transient and continuous ectopic activity. Electrograms and their dominant frequency and organization index values were calculated over the entire atrial surface. Our simulations show electrograms with simple potentials, with little or no cycle length variations, narrow frequency peaks and high organization index values during stable and regular activity as the observed in atrial flutter, atrial tachycardia (except in areas of conduction block) and in areas closer to ectopic activity during focal atrial fibrillation. By contrast, cycle length variations and polymorphic electrograms with single, double and fragmented potentials were observed in areas of irregular and unstable activity during atrial fibrillation episodes. Our results also show: (1) electrograms with potentials without negative deflection related to spiral or curved wavefronts that pass over the recording point and move away, (2) potentials with a much greater proportion of positive deflection than negative in areas of wave collisions, (3) double potentials related with wave fragmentations or blocking lines and (4) fragmented electrograms associated with pivot points. Our model is the first human atrial model with realistic fiber orientation used to investigate the relationship between different atrial arrhythmic propagation patterns and the electrograms observed at more than 43000 points on the atrial surface.

A bilayer representation of the human atria

Atrial fibrillation is the most commonly encountered clinical arrhythmia. Despite recent advances in treatment by catheter ablation, its origin is still incompletely understood and it may be difficult to treat. Computer modelling offers an attractive complement to experiment. Simulations of fibrillation, however, are computationally demanding since the phenomenon requires long periods of observation. Because the atria are thin walled structures, they are often modelled as surfaces. However, this may not always be appropriate as the crista terminalis and pectinate muscles are discrete fibrous structures lying on the endocardium and cannot be incorporated into the surface. In the left atrium, there are essentially two layers with an abrupt change in fibre orientation between them. We propose a double layer method, using shell elements to incorporate wall thickness, where fibre direction is independent in each layer and layers are electrically linked. Starting from human multi-detector CT (MDCT) images, we extracted surfaces for the atria and manually added a coronary sinus. Propagation of electrical activity was modelled with the monodomain equation. Results indicate that major features are retained while reducing computation cost considerably. Meshes based on the two layer approach will facilitate studies of AF.

Personalization of atrial anatomy and electrophysiology as a basis for clinical modeling of radio-frequency ablation of atrial fibrillation

Multiscale cardiac modeling has made great advances over the last decade. Highly detailed atrial models were created and used for the investigation of initiation and perpetuation of atrial fibrillation. The next challenge is the use of personalized atrial models in clinical practice. In this study, a framework of simple and robust tools is presented, which enables the generation and validation of patient-specific anatomical and electrophysiological atrial models. Introduction of rule-based atrial fiber orientation produced a realistic excitation sequence and a better correlation to the measured electrocardiograms. Personalization of the global conduction velocity lead to a precise match of the measured P-wave duration. The use of a virtual cohort of nine patient and volunteer models averaged out possible model-specific errors. Intra-atrial excitation conduction was personalized manually from left atrial local activation time maps. Inclusion of LE-MRI data into the simulations revealed possible gaps in ablation lesions. A fast marching level set approach to compute atrial depolarization was extended to incorporate anisotropy and conduction velocity heterogeneities and reproduced the monodomain solution. The presented chain of tools is an important step towards the use of atrial models for the patient-specific AF diagnosis and ablation therapy planing.

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.

Human atrial fibrillation: insights from computational electrophysiological models

Computational electrophysiology has proven useful to investigate the mechanisms of cardiac arrhythmias at various spatial scales, from isolated myocytes to the whole heart. This article reviews how mathematical modeling has aided our understanding of human atrial myocyte electrophysiology to study the contribution of structural and electrical remodeling to human atrial fibrillation. Potential new avenues of investigation and model development are suggested.

Towards personalized clinical in-silico modeling of atrial anatomy and electrophysiology

Computational atrial models aid the understanding of pathological mechanisms and therapeutic measures in basic research. The use of biophysical models in a clinical environment requires methods to personalize the anatomy and electrophysiology (EP). Strategies for the automation of model generation and for evaluation are needed. In this manuscript, the current efforts of clinical atrial modeling in the euHeart project are summarized within the context of recent publications in this field. Model-based segmentation methods allow for the automatic generation of ready-to-simulate patient-specific anatomical models. EP models can be adapted to patient groups based on a-priori knowledge and to the individual without significant further data acquisition. ECG and intracardiac data build the basis for excitation personalization. Information from late enhancement (LE) MRI can be used to evaluate the success of radio-frequency ablation (RFA) procedures and interactive virtual atria pave the way for RFA planning. Atrial modeling is currently in a transition from the sole use in basic research to future clinical applications. The proposed methods build the framework for model-based diagnosis and therapy evaluation and planning. Complex models allow to understand biophysical mechanisms and enable the development of simplified models for clinical applications.

Mechanisms of human atrial fibrillation initiation: clinical and computational studies of repolarization restitution and activation latency

BACKGROUND

Mechanisms of atrial fibrillation (AF) initiation are incompletely understood. We hypothesized that rate-dependent changes (restitution) in action potential duration (APD) and activation latency are central targets for clinical interventions that induce AF. We tested this hypothesis using clinical experiments and computer models.

METHODS AND RESULTS

In 50 patients (20 persistent, 23 paroxysmal AF, 7 controls), we used monophasic action potential catheters to define left atrial APD restitution, activation latency, and AF incidence from premature extrastimuli. Isoproterenol (n=14), adenosine (n=10), or rapid pacing (n=36) was then initiated to determine impact on these parameters. Compared with baseline in AF patients, isoproterenol and rapid pacing decreased activation latency (64±14 versus 31±13 versus 24±14 ms; P<0.05), steepened maximum APD restitution slope (0.8±0.7 versus 1.7±0.5 versus 1.1±0.5; P<0.05), and increased AF incidence (12% versus 64% versus 84%; P<0.05). Conversely, adenosine shortened APD (P<0.05), yet increased activation latency (86±27 ms; P=0.002) so that maximum APD restitution slope did not steepen (1.0±0.5; P=NS), and AF incidence was unchanged (10%; P=NS). In controls, no intervention steepened APD restitution or initiated AF. Computational modeling revealed that isoproterenol steepened APD restitution by increased L-type calcium current and decreased activation latency via enhanced rapid delayed potassium reactifier current inactivation, whereas rapid pacing steepened APD restitution via increased cardiac inward potassium rectifier current.

CONCLUSIONS

Steep APD restitution is a common pathway for AF initiation by isoproterenol and tachycardia via reduced activation latency that enables engagement of steep APD restitution at rapid rates. Modeling suggests that AF initiation from each intervention uses distinct ionic mechanisms. This insight may help design interventions to prevent AF.

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.

Structure-function relationship in the sinus and atrioventricular nodes

Recently published optical mapping studies of larger mammals, including humans, have identified functionally discrete sinoatrial exit pathways of activation. This is in line with earlier mapping studies of the dog and the human but in contrast with findings in the mouse and the rabbit, wherein a propagation wave front pattern of activation has been described. It underpins the complex three-dimensional (3D) organization of the cardiac pacemaking and conduction system in larger species, wherein sinoatrial and atrioventricular nodal physiologies both demonstrate identifiable activation pathways, which coincide with anatomic landmarks and histologic architecture, so that in addition to muscle fiber orientation and cell coupling, these intrinsic factors act to determine excitation pathways. This complex 3D organization increases the effect of source-to-sink mismatch both by greater variability in the space constant of tissue and by the 3D projection of this effect in all directions. Mathematical modeling provides a means to study these interactions, and newer models should incorporate these additional factors and their effect into the 3D structure of large mammal physiology.

Computational modeling of the human atrial anatomy and electrophysiology

This review article gives a comprehensive survey of the progress made in computational modeling of the human atria during the last 10 years. Modeling the anatomy has emerged from simple "peanut"-like structures to very detailed models including atrial wall and fiber direction. Electrophysiological models started with just two cellular models in 1998. Today, five models exist considering e.g. details of intracellular compartments and atrial heterogeneity. On the pathological side, modeling atrial remodeling and fibrotic tissue are the other important aspects. The bridge to data that are measured in the catheter laboratory and on the body surface (ECG) is under construction. Every measurement can be used either for model personalization or for validation. Potential clinical applications are briefly outlined and future research perspectives are suggested.

Correlation between P-wave morphology and origin of atrial focal tachycardia--insights from realistic models of the human atria and torso

Atrial arrhythmias resulting from abnormally rapid focal activity in the atria may be reflected in an altered P-wave morphology (PWM) in the ECG. Although clinically important, detailed relationships between PWM and origins of atrial focal excitations have not been established. To study such relationships, we developed computational models of the human atria and torso. The model simulation results were used to evaluate an extant clinical algorithm for locating the origin of atrial focal points from the ECG. The simulations showed that the algorithm was practical and could predict the atrial focal locations with 85% accuracy. We proposed a further refinement of the algorithm to distinguish between focal locations within the large atrial bundles.