Biomaterial injection is a novel therapy to treat ischemic heart failure (HF) that has shown to reduce remodeling and restore cardiac function in recent preclinical studies. While the effect of biomaterial injection in reducing mechanical wall stress has been recently demonstrated, the influence of biomaterials on the electrical behavior of treated hearts has not been elucidated. In this work, we developed computational models of swine hearts to study the electrophysiological vulnerability associated with biomaterial injection therapy. The propagation of action potentials on realistic biventricular geometries was simulated by numerically solving the monodomain electrophysiology equations on anatomically-detailed models of normal, HF untreated, and HF treated hearts. Heart geometries were constructed from high-resolution magnetic resonance images (MRI) where the healthy, peri-infarcted, infarcted and gel regions were identified, and the orientation of cardiac fibers was informed from diffusion-tensor MRI. Regional restitution properties in each case were evaluated by constructing a probability density function of the action potential duration (APD) at different cycle lengths. A comparative analysis of the ventricular fibrillation (VF) dynamics for every heart was carried out by measuring the number of filaments formed after wave braking. Our results suggest that biomaterial injection therapy does not affect the regional dispersion of repolarization when comparing untreated and treated failing hearts. Further, we found that the treated failing heart is more prone to sustain VF than the normal heart, and is at least as susceptible to sustained VF as the untreated failing heart. Moreover, we show that the main features of VF dynamics in a treated failing heart are not affected by the level of electrical conductivity of the biogel injectates. This work represents a novel proof-of-concept study demonstrating the feasibility of computer simulations of the heart in understanding the arrhythmic behavior in novel therapies for HF.

Computational cardiology is rapidly becoming the gold standard for innovative medical treatments and device development. Despite a worldwide effort in mathematical and computational modeling research, the complexity and intrinsic multiscale nature of the heart still limit our predictability power raising the question of the optimal modeling choice for large-scale whole-heart numerical investigations. We propose an extended numerical analysis among two different electrophysiological modeling approaches: a simplified phenomenological one and a detailed biophysical one. To achieve this, we considered three-dimensional healthy and infarcted swine heart geometries. Heterogeneous electrophysiological properties, fine-tuned DT-MRI -based anisotropy features, and non-conductive ischemic regions were included in a custom-built finite element code. We provide a quantitative comparison of the electrical behaviors during steady pacing and sustained ventricular fibrillation for healthy and diseased cases analyzing cardiac arrhythmias dynamics. Action potential duration (APD) restitution distributions, vortex filament counting, and pseudo-electrocardiography (ECG) signals were numerically quantified, introducing a novel statistical description of restitution patterns and ventricular fibrillation sustainability. Computational cost and scalability associated with the two modeling choices suggests that ventricular fibrillation signatures are mainly controlled by anatomy and structural parameters, rather than by regional restitution properties. Finally, we discuss limitations and translational perspectives of the different modeling approaches in view of large-scale whole-heart in silico studies.