Characterizing Conduction Channels in Postinfarction Patients Using a Personalized Virtual Heart

Computational Medicine: What Electrophysiologists Should Know to Stay Ahead of the Curve

PURPOSE OF REVIEW

Technology drives the field of cardiac electrophysiology. Recent computational advances will bring exciting changes. To stay ahead of the curve, we recommend electrophysiologists develop a robust appreciation for novel computational techniques, including deterministic, statistical, and hybrid models.

RECENT FINDINGS

In clinical applications, deterministic models use biophysically detailed simulations to offer patient-specific insights. Statistical techniques like machine learning and artificial intelligence recognize patterns in data. Emerging clinical tools are exploring avenues to combine all the above methodologies. We review three ways that computational medicine will aid electrophysiologists by: (1) improving personalized risk assessments, (2) weighing treatment options, and (3) guiding ablation procedures. Leveraging clinical data that are often readily available, computational models will offer valuable insights to improve arrhythmia patient care. As emerging tools promote personalized medicine, physicians must continue to critically evaluate technology-driven tools they consider using to ensure their appropriate implementation.

Cardiac MRI Oversampling in Heart Digital Twins Improves Preprocedure Ventricular Tachycardia Identification in Postinfarction Patients

BACKGROUND

Ventricular tachycardia (VT), which can lead to sudden cardiac death, occurs frequently in patients after myocardial infarction. Radiofrequency catheter ablation (RFA) is a modestly effective treatment of VT, but it has limitations and risks. Cardiac magnetic resonance (CMR)-based heart digital twins have emerged as a useful tool for identifying VT circuits for RFA treatment planning. However, the CMR resolution used to reconstruct these digital twins may impact VT circuit predictions, leading to incorrect RFA treatment planning.

OBJECTIVES

This study sought to predict RFA targets in the arrhythmogenic substrate using heart digital twins reconstructed from both clinical and high-resolution 2-dimensional CMR datasets and compare the predictions.

METHODS

High-resolution (1.35 × 1.35 × 3 mm), or oversampled resolution (Ov-Res), short-axis late gadolinium-enhanced CMR was acquired by combining 2 subsequent clinical resolution (Clin-Res) (1.35 × 1.35 × 6 mm) short-axis late gadolinium-enhanced CMR scans from 6 post-myocardial infarction patients undergoing VT ablation and used to reconstruct a total of 3 digital twins (1 Ov-Res, 2 Clin-Res) for each patient. Rapid pacing was used to assess VT circuits and identify the optimal ablation targets in each digital twin. VT circuits predicted by the digital twins were compared with intraprocedural electroanatomic mapping data and used to identify emergent VT.

RESULTS

The Ov-Res digital twins reduced partial volume effects and better predicted unique VT circuits compared with the Clin-Res digital twins (66.6% vs 54.5%; P < 0.01). Only the Ov-Res digital twin successfully identified emergent VT after a failed initial ablation.

CONCLUSIONS

Digital twin infarct geometry and VT circuit predictions depend on the magnetic resonance resolution. Ov-Res digital twins better predict VT circuits and emergent VT, which may improve RFA outcomes.

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Studying the Influence of Finite Element Mesh Size on the Accuracy of Ventricular Tachycardia Simulation

Background

Ventricular tachycardia (VT) is a life-threatening heart condition commonly seen in patients with myocardial infarction (MI). Although personalized computational modeling has been used to understand VT and its treatment noninvasively, this approach can be computationally intensive and time consuming. Therefore, finding a balance between mesh size and computational efficiency is important. This study aimed to find an optimal mesh resolution that minimizes the need for computational resources while maintaining numerical accuracy and to investigate the effect of mesh resolution variation on the simulation results.

Methods

We constructed ventricular models from contrast-enhanced magnetic resonance imaging data from six patients with MI. We created seven different models for each patient, with average edge lengths ranging from 315 to 645 µm using commercial software, Mimics. Programmed electrical stimulation was used to assess VT inducibility from 19 sites in each heart model.

Results

The simulation results in the slab model with adaptive tetrahedral mesh (same as in the patient-specific model) showed that the absolute and relative differences in conduction velocity (CV) were 6.1 cm/s and 7.8% between average mesh sizes of 142 and 600 µm, respectively. However, the simulation results in the six patient-specific models showed that average mesh sizes with 350 µm yielded over 85% accuracy for clinically relevant VT. Although average mesh sizes of 417 and 478 µm could also achieve approximately 80% accuracy for clinically relevant VT, the percentage of incorrectly predicted VTs increases. When conductivity was modified to match the CV in the model with the finest mesh size, the overall ratio of positively predicted VT increased.

Conclusions

The proposed personalized heart model could achieve an optimal balance between simulation time and VT prediction accuracy when discretized with adaptive tetrahedral meshes with an average edge length about 350 µm.

Predicting ventricular tachycardia circuits in patients with arrhythmogenic right ventricular cardiomyopathy using genotype-specific heart digital twins

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a genetic cardiac disease that leads to ventricular tachycardia (VT), a life-threatening heart rhythm disorder. Treating ARVC remains challenging due to the complex underlying arrhythmogenic mechanisms, which involve structural and electrophysiological (EP) remodeling. Here, we developed a novel genotype-specific heart digital twin (Geno-DT) approach to investigate the role of pathophysiological remodeling in sustaining VT reentrant circuits and to predict the VT circuits in ARVC patients of different genotypes. This approach integrates the patient's disease-induced structural remodeling reconstructed from contrast-enhanced magnetic-resonance imaging and genotype-specific cellular EP properties. In our retrospective study of 16 ARVC patients with two genotypes: plakophilin-2 (PKP2, n = 8) and gene-elusive (GE, n = 8), we found that Geno-DT accurately and non-invasively predicted the VT circuit locations for both genotypes (with 100%, 94%, 96% sensitivity, specificity, and accuracy for GE patient group, and 86%, 90%, 89% sensitivity, specificity, and accuracy for PKP2 patient group), when compared to VT circuit locations identified during clinical EP studies. Moreover, our results revealed that the underlying VT mechanisms differ among ARVC genotypes. We determined that in GE patients, fibrotic remodeling is the primary contributor to VT circuits, while in PKP2 patients, slowed conduction velocity and altered restitution properties of cardiac tissue, in addition to the structural substrate, are directly responsible for the formation of VT circuits. Our novel Geno-DT approach has the potential to augment therapeutic precision in the clinical setting and lead to more personalized treatment strategies in ARVC.

Predicting Ventricular Tachycardia Circuits in Patients with Arrhythmogenic Right Ventricular Cardiomyopathy using Genotype-specific Heart Digital Twins

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a genetic cardiac disease that leads to ventricular tachycardia (VT), a life-threatening heart rhythm disorder. Treating ARVC remains challenging due to the complex underlying arrhythmogenic mechanisms, which involve structural and electrophysiological (EP) remodeling. Here, we developed a novel genotype-specific heart digital twin (Geno-DT) approach to investigate the role of pathophysiological remodeling in sustaining VT reentrant circuits and to predict the VT circuits in ARVC patients of different genotypes. This approach integrates the patient's disease-induced structural remodeling reconstructed from contrast-enhanced magnetic-resonance imaging and genotype-specific cellular EP properties. In our retrospective study of 16 ARVC patients with two genotypes: plakophilin-2 (PKP2, n = 8) and gene-elusive (GE, n = 8), we found that Geno-DT accurately and non-invasively predicted the VT circuit locations for both genotypes (with 100%, 94%, 96% sensitivity, specificity, and accuracy for GE patient group, and 86%, 90%, 89% sensitivity, specificity, and accuracy for PKP2 patient group), when compared to VT circuit locations identified during clinical EP studies. Moreover, our results revealed that the underlying VT mechanisms differ among ARVC genotypes. We determined that in GE patients, fibrotic remodeling is the primary contributor to VT circuits, while in PKP2 patients, slowed conduction velocity and altered restitution properties of cardiac tissue, in addition to the structural substrate, are directly responsible for the formation of VT circuits. Our novel Geno-DT approach has the potential to augment therapeutic precision in the clinical setting and lead to more personalized treatment strategies in ARVC.

Fat infiltration in the infarcted heart as a paradigm for ventricular arrhythmias

Infiltrating adipose tissue (inFAT) has been recently found to co-localize with scar in infarcted hearts and may contribute to ventricular arrhythmias (VAs), a life-threatening heart rhythm disorder. However, the contribution of inFAT to VA has not been well-established. We investigated the role of inFAT versus scar in VA through a combined prospective clinical and mechanistic computational study. Using personalized computational heart models and comparing the results from simulations of VA dynamics with measured electrophysiological abnormalities during the clinical procedure, we demonstrate that inFAT, rather than scar, is a primary driver of arrhythmogenic propensity and is frequently present in critical regions of the VA circuit. We determined that, within the VA circuitry, inFAT, as opposed to scar, is primarily responsible for conduction slowing in critical sites, mechanistically promoting VA. Our findings implicate inFAT as a dominant player in infarct-related VA, challenging existing paradigms and opening the door for unexplored anti-arrhythmic strategies.

An automated near-real time computational method for induction and treatment of scar-related ventricular tachycardias

Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.

How to use pace mapping for ventricular tachycardia ablation in post-infarct patients

We aim to describe the technical aspects of pace mapping (PM), as well as the two typical patterns of pacing correlation maps during ventricular tachycardia (VT) ablation. The first main pattern is focal, with a gradual and eccentric decrease of the QRS correlation from the area with the best PM correlation. This focal pattern may be associated with two clinical situations: (1) with some endocardial points showing a good correlation compared to VT morphology: true endocardial exit of VT or endocardial breakthrough of either an intramural or an epicardial circuit; (2) without any endocardial points showing a good correlation compared to VT morphology: the VT may originate from the other ventricle, but the presence of an intramural or an epicardial circuit should be considered in patients with a structural heart disease. The second pattern is the presence of PM points exhibiting a good correlation close to other PM points showing a poor correlation compared to VT morphology: this abrupt change in paced QRS morphology over a short distance indicates divergence of activation wavefronts between these sites and suggests the presence of a slow conduction channel: the VT isthmus.

An Automata-Based Cardiac Electrophysiology Simulator to Assess Arrhythmia Inducibility

Personalized cardiac electrophysiology simulations have demonstrated great potential to study cardiac arrhythmias and help in therapy planning of radio-frequency ablation. Its application to analyze vulnerability to ventricular tachycardia and sudden cardiac death in infarcted patients has been recently explored. However, the detailed multi-scale biophysical simulations used in these studies are very demanding in terms of memory and computational resources, which prevents their clinical translation. In this work, we present a fast phenomenological system based on cellular automata (CA) to simulate personalized cardiac electrophysiology. The system is trained on biophysical simulations to reproduce cellular and tissue dynamics in healthy and pathological conditions, including action potential restitution, conduction velocity restitution and cell safety factor. We show that a full ventricular simulation can be performed in the order of seconds, emulate the results of a biophysical simulation and reproduce a patient’s ventricular tachycardia in a model that includes a heterogeneous scar region. The system could be used to study the risk of arrhythmia in infarcted patients for a large number of scenarios.

Simulated late gadolinium enhanced cardiac magnetic resonance imaging dataset from mechanical XCAT phantom including a myocardial infarct

The late enhanced magnetic resonance image dataset in this article is simulated using a mechanistic cardiac phantom that includes an myocardial infarct. Settings of the image simulation pipeline are adjusted such that high- and low-resolution images, with and without slice alignment artifacts, are simulated. Our article on the influence of image artifacts on image-based models of the cardiac electrophysiology is based on this data (Kruithof et al., 2021). This dataset provides image-analysis researchers a reference to perform validation of their methods using the included high-resolution ground truth image, a resource that is often unavailable clinically.

Preliminary Study: Learning the Impact of Simulation Time on Reentry Location and Morphology Induced by Personalized Cardiac Modeling

Personalized cardiac modeling is widely used for studying the mechanisms of cardiac arrythmias. Due to the high demanding of computational resource of modeling, the arrhythmias induced in the models are usually simulated for just a few seconds. In clinic, it is common that arrhythmias last for more than several minutes and the morphologies of reentries are not always stable, so it is not clear that whether the simulation of arrythmias for just a few seconds is long enough to match the arrhythmias detected in patients. This study aimed to observe how long simulation of the induced arrhythmias in the personalized cardiac models is sufficient to match the arrhythmias detected in patients. A total of 5 contrast enhanced MRI datasets of patient hearts with myocardial infarction were used in this study. Then, a classification method based on Gaussian mixture model was used to detect the infarct tissue. For each reentry, 3 s and 10 s were simulated. The characteristics of each reentry simulated for different duration were studied. Reentries were induced in all 5 ventricular models and sustained reentries were induced at 39 stimulation sites in the model. By analyzing the simulation results, we found that 41% of the sustained reentries in the 3 s simulation group terminated in the longer simulation groups (10 s). The second finding in our simulation was that only 23.1% of the sustained reentries in the 3 s simulation did not change location and morphology in the extended 10 s simulation. The third finding was that 35.9% reentries were stable in the 3 s simulation and should be extended for the simulation time. The fourth finding was that the simulation results in 10 s simulation matched better with the clinical measurements than the 3 s simulation. It was shown that 10 s simulation was sufficient to make simulation results stable. The findings of this study not only improve the simulation accuracy, but also reduce the unnecessary simulation time to achieve the optimal use of computer resources to improve the simulation efficiency and shorten the simulation time to meet the time node requirements of clinical operation on patients.

Whole-heart ventricular arrhythmia modeling moving forward: Mechanistic insights and translational applications

Ventricular arrhythmias are the primary cause of sudden cardiac death and one of the leading causes of mortality worldwide. Whole-heart computational modeling offers a unique approach for studying ventricular arrhythmias, offering vast potential for developing both a mechanistic understanding of ventricular arrhythmias and clinical applications for treatment. In this review, the fundamentals of whole-heart ventricular modeling and current methods of personalizing models using clinical data are presented. From this foundation, the authors summarize recent advances in whole-heart ventricular arrhythmia modeling. Efforts in gaining mechanistic insights into ventricular arrhythmias are discussed, in addition to other applications of models such as the assessment of novel therapeutics. The review emphasizes the unique benefits of computational modeling that allow for insights that are not obtainable by contemporary experimental or clinical means. Additionally, the clinical impact of modeling is explored, demonstrating how patient care is influenced by the information gained from ventricular arrhythmia models. The authors conclude with future perspectives about the direction of whole-heart ventricular arrhythmia modeling, outlining how advances in neural network methodologies hold the potential to reduce computational expense and permit for efficient whole-heart modeling.

Influence of image artifacts on image-based computer simulations of the cardiac electrophysiology

Myocardial infarct patients have an increased risk of scar-based ventricular tachycardia. Late gadolinium enhanced magnetic resonance (MR) imaging provides the geometric extent of myocardial infarct. Computational electrophysiological models based on such images can provide a personalized prediction of the patient's tachycardia risk. In this work, the effect of respiratory slice alignment image artifacts on image-based electrophysiological simulations is investigated in two series of models. For the first series, a clinical MR image is used in which slice translations are applied to artificially induce and correct for slice misalignment. For the second series, computer simulated MR images with and without slice misalignments are created using a mechanistic anatomical phantom of the torso. From those images, personalized models are created in which electrical stimuli are applied in an attempt to induce tachycardia. The response of slice-aligned and slice-misaligned models to different interval stimuli is used to assess tachycardia risk. The presented results indicate that slice misalignments affect image-based simulation outcomes. The extent to which the assessed risk is affected is found to depend upon the geometry of the infarct area. The number of unidirectional block tachycardias varied from 1 to 3 inducible patterns depending on slice misalignment severity and, along with it, the number of tachycardia inducing stimuli locations varied from 2 to 4 from 6 different locations. For tachycardias sustained by conducting channels through the scar core, no new patterns are induced by altering the slice alignment in the corresponding image. However, it affected the assessed risk as tachycardia inducing stimuli locations varied from 1 to 5 from the 6 stimuli locations. In addition, if the conducting channel is not maintained in the image due to slice misalignments, the channel-dependent tachycardia is not inducible anymore in the image-based model.

Impact of augmented-reality improvement in ablation catheter navigation as assessed by virtual-heart simulations of ventricular tachycardia ablation

BACKGROUND

Recently, an augmented reality (AR) solution allows the physician to place the ablation catheter at the designated lesion site more accurately during cardiac electrophysiology studies. The improvement in navigation accuracy may positively affect ventricular tachycardia (VT) ablation termination, however assessment of this in the clinic would be difficult. Novel personalized virtual heart technology enables non-invasive identification of optimal lesion targets for infarct-related VT. This study aims to evaluate the potential impact of such catheter navigation accuracy improvement in virtual VT ablations.

METHODS

2 MRI-based virtual hearts with 2 in silico induced VTs (VT 1, VT 2) were included. VTs were terminated with virtual "ground truth" endocardial ablation lesions. 106 navigation error values that were previously assessed in a clinical study evaluating the improvement of ablation catheter navigation accuracy guided with AR (53 with, 53 without) were used to displace the "ground truth" ablation targets. The corresponding ablations were simulated based on these errors and VT termination for each simulation was assessed.

RESULTS

In 54 VT 1 ablation simulations, smaller error with AR significantly resulted in more VT termination (25) compared to the error without AR (16) (P < 0.01). In 52 VT 2 ablation simulations, no significant difference was observed from error with (11) and without AR (13) (P = 0.58). The substrate characteristic may impact the effect of improved accuracy to an improved VT termination.

CONCLUSION

Virtual heart shows that the increased catheter navigation accuracy provided by AR guidance can affect the VT termination.

Fully Automatic Scar Segmentation for Late Gadolinium Enhancement MRI Images in Left Ventricle with Myocardial Infarction

Numerous methods have been published to segment the infarct tissue in the left ventricle, most of them either need manual work, post-processing, or suffer from poor reproducibility. We proposed an automatic segmentation method for segmenting the infarct tissue in left ventricle with myocardial infarction. Cardiac images of a total of 60 diseased hearts (55 human hearts and 5 porcine hearts) were used in this study. The epicardial and endocardial boundaries of the ventricles in every 2D slice of the cardiac magnetic resonance with late gadolinium enhancement images were manually segmented. The subsequent pipeline of infarct tissue segmentation is fully automatic. The segmentation results with the automatic algorithm proposed in this paper were compared to the consensus ground truth. The median of Dice overlap between our automatic method and the consensus ground truth is 0.79. We also compared the automatic method with the consensus ground truth using different image sources from different centers with different scan parameters and different scan machines. The results showed that the Dice overlap with the public dataset was 0.83, and the overall Dice overlap was 0.79. The results show that our method is robust with respect to different MRI image sources, which were scanned by different centers with different image collection parameters. The segmentation accuracy we obtained is comparable to or better than that of the conventional semi-automatic methods. Our segmentation method may be useful for processing large amount of dataset in clinic.

Role of cardiac imaging in patients undergoing catheter ablation of ventricular tachycardia

Ventricular tachycardia is a major health issue in patients with structural heart disease (SHD). Implantable cardioverter defibrillator (ICD) therapy has significantly reduced the risk of sudden cardiac death (SCD) in such patients, but on the other hand, it has led to frequent ICD shocks as an emerging problem, being associated with poor quality of life, frequent hospitalizations and increased mortality. Myocardial scar plays a central role in the genesis and maintenance of re-entrant arrhythmias, as the coexistence of surviving myocardial fibres within fibrotic tissue leads to the formation of slow conduction pathways and to a dispersion of activation and refractoriness that constitutes the milieu for ventricular tachycardia circuits. Catheter ablation has repeatedly proven to be well tolerated and highly effective in treating VT and in the last two decades has benefited from continuous efforts to determine ventricular tachycardia mechanisms by integration with a wide range of invasive and noninvasive imaging techniques such as intracardiac echocardiography, cardiac magnetic resonance, multidetector computed tomography and nuclear imaging. Cardiovascular imaging has become a fundamental aid in planning and guiding catheter ablation procedures by integrating structural and electrophysiological information, enabling the ventricular tachycardia arrhythmogenic substrate to be characterized and effective ablation targets to be identified with increasing precision, and allowing the development of new ablation strategies with improved outcomes. In this review, we provide an overview of the role of cardiac imaging in patients undergoing catheter ablation of ventricular tachycardia.

How personalized heart modeling can help treatment of lethal arrhythmias: A focus on ventricular tachycardia ablation strategies in post-infarction patients

Precision Cardiology is a targeted strategy for cardiovascular disease prevention and treatment that accounts for individual variability. Computational heart modeling is one of the novel approaches that have been developed under the umbrella of Precision Cardiology. Personalized computational modeling of patient hearts has made strides in the development of models that incorporate the individual geometry and structure of the heart as well as other patient-specific information. Of these developments, one of the potentially most impactful is the research aimed at noninvasively predicting the targets of ablation of lethal arrhythmia, ventricular tachycardia (VT), using patient-specific models. The approach has been successfully applied to patients with ischemic cardiomyopathy in proof-of-concept studies. The goal of this paper is to review the strategies for computational VT ablation guidance in ischemic cardiomyopathy patients, from model developments to the intricacies of the actual clinical application. To provide context in describing the road these computational modeling applications have undertaken, we first review the state of the art in VT ablation in the clinic, emphasizing the benefits that personalized computational prediction of ablation targets could bring to the clinical electrophysiology practice. This article is characterized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Translational, Genomic, and Systems Medicine > Translational Medicine.