Accelerating Cardiac Bidomain Simulations Using Graphics Processing Units

Anatomically realistic and biophysically detailed multiscale computer models of the heart are playing an increasingly important role in advancing our understanding of integrated cardiac function in health and disease. Such detailed simulations, however, are computationally vastly demanding, which is a limiting factor for a wider adoption of in-silico modeling. While current trends in high-performance computing (HPC) hardware promise to alleviate this problem, exploiting the potential of such architectures remains challenging since strongly scalable algorithms are necessitated to reduce execution times. Alternatively, acceleration technologies such as graphics processing units (GPUs) are being considered. While the potential of GPUs has been demonstrated in various applications, benefits in the context of bidomain simulations where large sparse linear systems have to be solved in parallel with advanced numerical techniques are less clear. In this study, the feasibility of multi-GPU bidomain simulations is demonstrated by running strong scalability benchmarks using a state-of-the-art model of rabbit ventricles. The model is spatially discretized using the finite element methods (FEM) on fully unstructured grids. The GPU code is directly derived from a large pre-existing code, the Cardiac Arrhythmia Research Package (CARP), with very minor perturbation of the code base. Overall, bidomain simulations were sped up by a factor of 11.8 to 16.3 in benchmarks running on 6-20 GPUs compared to the same number of CPU cores. To match the fastest GPU simulation which engaged 20 GPUs, 476 CPU cores were required on a national supercomputing facility

Predicting Atrial Fibrillation Recurrence by Combining Population Data and Virtual Cohorts of Patient-Specific Left Atrial Models.

BACKGROUND: Current ablation therapy for atrial fibrillation is suboptimal, and long-term response is challenging to predict. Clinical trials identify bedside properties that provide only modest prediction of long-term response in populations, while patient-specific models in small cohorts primarily explain acute response to ablation. We aimed to predict long-term atrial fibrillation recurrence after ablation in large cohorts, by using machine learning to complement biophysical simulations by encoding more interindividual variability. METHODS: Patient-specific models were constructed for 100 atrial fibrillation patients (43 paroxysmal, 41 persistent, and 16 long-standing persistent), undergoing first ablation. Patients were followed for 1 year using ambulatory ECG monitoring. Each patient-specific biophysical model combined differing fibrosis patterns, fiber orientation maps, electrical properties, and ablation patterns to capture uncertainty in atrial properties and to test the ability of the tissue to sustain fibrillation. These simulation stress tests of different model variants were postprocessed to calculate atrial fibrillation simulation metrics. Machine learning classifiers were trained to predict atrial fibrillation recurrence using features from the patient history, imaging, and atrial fibrillation simulation metrics. RESULTS: We performed 1100 atrial fibrillation ablation simulations across 100 patient-specific models. Models based on simulation stress tests alone showed a maximum accuracy of 0.63 for predicting long-term fibrillation recurrence. Classifiers trained to history, imaging, and simulation stress tests (average 10-fold cross-validation area under the curve, 0.85±0.09; recall, 0.80±0.13; precision, 0.74±0.13) outperformed those trained to history and imaging (area under the curve, 0.66±0.17) or history alone (area under the curve, 0.61±0.14). CONCLUSION: A novel computational pipeline accurately predicted long-term atrial fibrillation recurrence in individual patients by combining outcome data with patient-specific acute simulation response. This technique could help to personalize selection for atrial fibrillation ablation

Electrophysiology Model for a Human Heart with Ischemic Scar and Realistic Purkinje Network

The role of Purkinje fibres in the onset of arrhythmias is controversial and computer simulations may shed light on possible arrhythmic mechanisms involving the Purkinje fibres. However, few computational modelling studies currently include a detailed Purkinje network as part of the model. We present a coupled Purkinje-myocardium electrophysiology model that includes an explicit model for the ischemic scar plus a detailed Purkinje network, and compare simulated activation times to those obtained by electro-anatomical mapping in vivo during sinus rhythm pacing. The results illustrate the importance of using sufficiently dense Purkinje networks in patient-specific studies to capture correctly the myocardial early activation that may be influenced by surviving Purkinje fibres in the infarct region

Constructing bilayer and volumetric atrial models at scale.

To enable large in silico trials and personalized model predictions on clinical timescales, it is imperative that models can be constructed quickly and reproducibly. First, we aimed to overcome the challenges of constructing cardiac models at scale through developing a robust, open-source pipeline for bilayer and volumetric atrial models. Second, we aimed to investigate the effects of fibres, fibrosis and model representation on fibrillatory dynamics. To construct bilayer and volumetric models, we extended our previously developed coordinate system to incorporate transmurality, atrial regions and fibres (rule-based or data driven diffusion tensor magnetic resonance imaging (MRI)). We created a cohort of 1000 biatrial bilayer and volumetric models derived from computed tomography (CT) data, as well as models from MRI, and electroanatomical mapping. Fibrillatory dynamics diverged between bilayer and volumetric simulations across the CT cohort (correlation coefficient for phase singularity maps: left atrial (LA) 0.27 ± 0.19, right atrial (RA) 0.41 ± 0.14). Adding fibrotic remodelling stabilized re-entries and reduced the impact of model type (LA: 0.52 ± 0.20, RA: 0.36 ± 0.18). The choice of fibre field has a small effect on paced activation data (less than 12 ms), but a larger effect on fibrillatory dynamics. Overall, we developed an open-source user-friendly pipeline for generating atrial models from imaging or electroanatomical mapping data enabling in silico clinical trials at scale (https://github.com/pcmlab/atrialmtk)

A Regional Reduction in Ito and IKACh in the Murine Posterior Left Atrial Myocardium Is Associated with Action Potential Prolongation and Increased Ectopic Activity.

BACKGROUND: The left atrial posterior wall (LAPW) is potentially an important area for the development and maintenance of atrial fibrillation. We assessed whether there are regional electrical differences throughout the murine left atrial myocardium that could underlie regional differences in arrhythmia susceptibility. METHODS: We used high-resolution optical mapping and sharp microelectrode recordings to quantify regional differences in electrical activation and repolarisation within the intact, superfused murine left atrium and quantified regional ion channel mRNA expression by Taqman Low Density Array. We also performed selected cellular electrophysiology experiments to validate regional differences in ion channel function. RESULTS: Spontaneous ectopic activity was observed during sustained 1Hz pacing in 10/19 intact LA and this was abolished following resection of LAPW (0/19 resected LA, P<0.001). The source of the ectopic activity was the LAPW myocardium, distinct from the pulmonary vein sleeve and LAA, determined by optical mapping. Overall, LAPW action potentials (APs) were ca. 40% longer than the LAA and this region displayed more APD heterogeneity. mRNA expression of Kcna4, Kcnj3 and Kcnj5 was lower in the LAPW myocardium than in the LAA. Cardiomyocytes isolated from the LAPW had decreased Ito and a reduced IKACh current density at both positive and negative test potentials. CONCLUSIONS: The murine LAPW myocardium has a different electrical phenotype and ion channel mRNA expression profile compared with other regions of the LA, and this is associated with increased ectopic activity. If similar regional electrical differences are present in the human LA, then the LAPW may be a potential future target for treatment of atrial fibrillation

Spatiotemporally Controlled Cardiac Conduction Block Using High-Frequency Electrical Stimulation

Background: Methods for the electrical inhibition of cardiac excitation have long been sought to control excitability and conduction, but to date remain largely impractical. High-amplitude alternating current (AC) stimulation has been known to extend cardiac action potentials (APs), and has been recently exploited to terminate reentrant arrhythmias by producing reversible conduction blocks. Yet, low-amplitude currents at similar frequencies have been shown to entrain cardiac tissues by generation of repetitive APs, leading in some cases to ventricular fibrillation and hemodynamic collapse in vivo. Therefore, an inhibition method that does not lead to entrainment – irrespective of the stimulation amplitude (bound to fluctuate in an in vivo setting) – is highly desirable. Methodology/Principal Findings: We investigated the effects of broader amplitude and frequency ranges on the inhibitory effects of extracellular AC stimulation on HL-1 cardiomyocytes cultured on microelectrode arrays, using both sinusoidal and square waveforms. Our results indicate that, at sufficiently high frequencies, cardiac tissue exhibits a binary response to stimulus amplitude with either prolonged APs or no effect, thereby effectively avoiding the risks of entrainment by repetitive firing observed at lower frequencies. We further demonstrate the ability to precisely define reversible local conduction blocks in beating cultures without influencing the propagation activity in non-blocked areas. The conduction blocks were spatiotemporally controlled by electrode geometry and stimuli duration, respectively, and sustainable for long durations (300 s). Conclusion/Significance: Inhibition of cardiac excitation induced by high-frequency AC stimulation exhibits a binary response to amplitude above a threshold frequency, enabling the generation of reversible conduction blocks without the risks of entrainment. This inhibition method could yield novel approaches for arrhythmia modeling in vitro, as well as safer and more efficacious tools for in vivo cardiac mapping and radio-frequency ablation guidance applications

Electrophysiological and Structural Remodeling in Heart Failure Modulate Arrhythmogenesis. 1D Simulation Study

Background: Heart failure is a final common pathway or descriptor for various cardiac pathologies. It is associated with sudden cardiac death, which is frequently caused by ventricular arrhythmias. Electrophysiological remodeling, intercellular uncoupling, fibrosis and autonomic imbalance have been identified as major arrhythmogenic factors in heart failure etiology and progression. Objective: In this study we investigate in silico the role of electrophysiological and structural heart failure remodeling on the modulation of key elements of the arrhythmogenic substrate, i.e., electrophysiological gradients and abnormal impulse propagation. Methods: Two different mathematical models of the human ventricular action potential were used to formulate models of the failing ventricular myocyte. This provided the basis for simulations of the electrical activity within a transmural ventricular strand. Our main goal was to elucidate the roles of electrophysiological and structural remodeling in setting the stage for malignant life-threatening arrhythmias. Results: Simulation results illustrate how the presence of M cells and heterogeneous electrophysiological remodeling in the human failing ventricle modulate the dispersion of action potential duration and repolarization time. Specifically, selective heterogeneous remodeling of expression levels for the Na+ /Ca2+ exchanger and SERCA pump decrease these heterogeneities. In contrast, fibroblast proliferation and cellular uncoupling both strongly increase repolarization heterogeneities. Conduction velocity and the safety factor for conduction are also reduced by the progressive structural remodeling during heart failure. Conclusion: An extensive literature now establishes that in human ventricle, as heart failure progresses, gradients for repolarization are changed significantly by protein specific electrophysiological remodeling (either homogeneous or heterogeneous). Our simulations illustrate and provide new insights into this. Furthermore, enhanced fibrosis in failing hearts, as well as reduced intercellular coupling, combine to increase electrophysiological gradients and reduce electrical propagation. In combination these changes set the stage for arrhythmias.This work was partially supported by (i) the "VI Plan Nacional de Investigacion Cientifica, Desarrollo e Innovacion Tecnologica" from the Ministerio de Economia y Competitividad of Spain (grant number TIN2012-37546-C03-01) and the European Commission (European Regional Development Funds - ERDF - FEDER), (ii) the Direccion General de Politica Cientifica de la Generalitat Valenciana (grant number GV/2013/119), and (iii) Programa Prometeo (PROMETEO/2012/030) de la Conselleria d'Educacio Formacio I Ocupacio, Generalitat Valenciana. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Gómez García, JF.; Cardona, K.; Romero Pérez, L.; Ferrero De Loma-Osorio, JM.; Trénor Gomis, BA. (2014). Electrophysiological and Structural Remodeling in Heart Failure Modulate Arrhythmogenesis. 1D Simulation Study. PLoS ONE. 9(9). https://doi.org/10.1371/journal.pone.0106602S9

Protective Role of False Tendon in Subjects with Left Bundle Branch Block: A Virtual Population Study.

False tendons (FTs) are fibrous or fibromuscular bands that can be found in both the normal and abnormal human heart in various anatomical forms depending on their attachment points, tissue types, and geometrical properties. While FTs are widely considered to affect the function of the heart, their specific roles remain largely unclear and unexplored. In this paper, we present an in silico study of the ventricular activation time of the human heart in the presence of FTs. This study presents the first computational model of the human heart that includes a FT, Purkinje network, and papillary muscles. Based on this model, we perform simulations to investigate the effect of different types of FTs on hearts with the electrical conduction abnormality of a left bundle branch block (LBBB). We employ a virtual population of 70 human hearts derived from a statistical atlas, and run a total of 560 simulations to assess ventricular activation time with different FT configurations. The obtained results indicate that, in the presence of a LBBB, the FT reduces the total activation time that is abnormally augmented due to a branch block, to such an extent that surgical implant of cardiac resynchronisation devices might not be recommended by international guidelines. Specifically, the simulation results show that FTs reduce the QRS duration at least 10 ms in 80% of hearts, and up to 45 ms for FTs connecting to the ventricular free wall, suggesting a significant reduction of cardiovascular mortality risk. In further simulation studies we show the reduction in the QRS duration is more sensitive to the shape of the heart then the size of the heart or the exact location of the FT. Finally, the model suggests that FTs may contribute to reducing the activation time difference between the left and right ventricles from 12 ms to 4 ms. We conclude that FTs may provide an alternative conduction pathway that compensates for the propagation delay caused by the LBBB. Further investigation is needed to quantify the clinical impact of FTs on cardiovascular mortality risk