Corrigendum to "effects of fibrosis morphology on reentrant ventricular tachycardia inducibility and simulation fidelity in patient-derived models"

Computational models in cardiology

The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.

Optimal contrast-enhanced MRI image thresholding for accurate prediction of ventricular tachycardia using ex-vivo high resolution models

Patient specific models created from contrast-enhanced (i.e. late-gadolinium, LGE) MRI images can be used for prediction of reentry location and clinical ablation planning. However, there is still a need for direct and systematic comparison between characteristics of ventricular tachycardia (VT) morphologies predicted in computational models and those acquired in clinical or experimental protocols. In this study, we aimed to: 1) assess the differences in VT morphologies predicted by modeling and recorded in experiments in terms of patterns and location of reentries, earliest and latest activation sites, and cycle lengths; and 2) define the optimal range of infarct tissue threshold values which provide best match between simulation and experimental results. To achieve these goals, we utilized LGE-MRI images from 4 swine hearts with inducible monomorphic VT. The images were segmented to identify non-infarcted myocardium, semi viable gray zone (GZ), and core scar based on pixel intensity. Several models were reconstructed from each LGE-MRI scan, with voxels of intensity between that of non-infarcted myocardium and 20-50% of the maximum intensity (in 10% increments) in the infarct region classified as GZ. VT induction was simulated in each model. Our simulation results showed that using GZ intensity thresholds of 20% or 30% resulted in the best match of simulated propagation patterns and reentry locations with those from the experiment. Overall, we matched 70% (7/10) morphologies for all the hearts. Our simulation shows that MRI-based computational models of hearts with myocardial infarction can accurately reproduce the majority of experimentally recorded post-infarction VTs.

Predicting the risk of sudden cardiac death

Sudden cardiac death (SCD) is the result of a change of cardiac activity from normal (typically sinus) rhythm to a rhythm that does not pump adequate blood to the brain. The most common rhythms leading to SCD are ventricular tachycardia (VT) or ventricular fibrillation (VF). These result from an accelerated ventricular pacemaker or ventricular reentrant waves. Despite significant efforts to develop accurate predictors for the risk of SCD, current methods for risk stratification still need to be improved. In this article we briefly review current approaches to risk stratification. Then we discuss the mathematical basis for dynamical transitions (called bifurcations) that may lead to VT and VF. One mechanism for transition to VT or VF involves a perturbation by a premature ventricular complex (PVC) during sinus rhythm. We describe the main mechanisms of PVCs (reentry, independent pacemakers and abnormal depolarizations). An emerging approach to risk stratification for SCD involves the development of individualized dynamical models of a patient based on measured anatomy and physiology. Careful analysis and modelling of dynamics of ventricular arrhythmia on an individual basis will be essential in order to improve risk stratification for SCD and to lay a foundation for personalized (precision) medicine in cardiology.

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.

Accuracy of prediction of infarct-related arrhythmic circuits from image-based models reconstructed from low and high resolution MRI

Identification of optimal ablation sites in hearts with infarct-related ventricular tachycardia (VT) remains difficult to achieve with the current catheter-based mapping techniques. Limitations arise from the ambiguities in determining the reentrant pathways location(s). The goal of this study was to develop experimentally validated, individualized computer models of infarcted swine hearts, reconstructed from high-resolution ex-vivo MRI and to examine the accuracy of the reentrant circuit location prediction when models of the same hearts are instead reconstructed from low clinical-resolution MRI scans. To achieve this goal, we utilized retrospective data obtained from four pigs ~10 weeks post infarction that underwent VT induction via programmed stimulation and epicardial activation mapping via a multielectrode epicardial sock. After the experiment, high-resolution ex-vivo MRI with late gadolinium enhancement was acquired. The Hi-res images were downsampled into two lower resolutions (Med-res and Low-res) in order to replicate image quality obtainable in the clinic. The images were segmented and models were reconstructed from the three image stacks for each pig heart. VT induction similar to what was performed in the experiment was simulated. Results of the reconstructions showed that the geometry of the ventricles including the infarct could be accurately obtained from Med-res and Low-res images. Simulation results demonstrated that induced VTs in the Med-res and Low-res models were located close to those in Hi-res models. Importantly, all models, regardless of image resolution, accurately predicted the VT morphology and circuit location induced in the experiment. These results demonstrate that MRI-based computer models of hearts with ischemic cardiomyopathy could provide a unique opportunity to predict and analyze VT resulting for from specific infarct architecture, and thus may assist in clinical decisions to identify and ablate the reentrant circuit(s).