Credibility assessment of patient-specific computational modeling using patient-specific cardiac modeling as an exemplar.

Reliable and robust simulation of individual patients using patient-specific models (PSMs) is one of the next frontiers for modeling and simulation (M&S) in healthcare. PSMs, which form the basis of digital twins, can be employed as clinical tools to, for example, assess disease state, predict response to therapy, or optimize therapy. They may also be used to construct virtual cohorts of patients, for in silico evaluation of medical product safety and/or performance. Methods and frameworks have recently been proposed for evaluating the credibility of M&S in healthcare applications. However, such efforts have generally been motivated by models of medical devices or generic patient models; how best to evaluate the credibility of PSMs has largely been unexplored. The aim of this paper is to understand and demonstrate the credibility assessment process for PSMs using patient-specific cardiac electrophysiological (EP) modeling as an exemplar. We first review approaches used to generate cardiac PSMs and consider how verification, validation, and uncertainty quantification (VVUQ) apply to cardiac PSMs. Next, we execute two simulation studies using a publicly available virtual cohort of 24 patient-specific ventricular models, the first a multi-patient verification study, the second investigating the impact of uncertainty in personalized and non-personalized inputs in a virtual cohort. We then use the findings from our analyses to identify how important characteristics of PSMs can be considered when assessing credibility with the approach of the ASME V&V40 Standard, accounting for PSM concepts such as inter- and intra-user variability, multi-patient and "every-patient" error estimation, uncertainty quantification in personalized vs non-personalized inputs, clinical validation, and others. The results of this paper will be useful to developers of cardiac and other medical image based PSMs, when assessing PSM credibility.

High mean entropy calculated from cardiac MRI texture analysis is associated with antitachycardia pacing failure.

BACKGROUND: Antitachycardia pacing (ATP), which may avoid unnecessary implantable cardioverter-defibrillator (ICD) shocks, does not always terminate ventricular arrhythmias (VAs). Mean entropy calculated using cardiac magnetic resonance texture analysis (CMR-TA) has been shown to predict appropriate ICD therapy. We examined whether scar heterogeneity, quantified by mean entropy, is associated with ATP failure and explore potential mechanisms using computer modeling. METHODS: A subanalysis of 114 patients undergoing CMR-TA where the primary endpoint was delivery of appropriate ICD therapy (ATP or shock therapy) was performed. Patients receiving appropriate ICD therapy (n = 33) were dichotomized into "successful ATP" versus "shock therapy" groups. In silico computer modeling was used to explore underlying mechanisms. RESULTS: A total of 16 of 33 (48.5%) patients had successful ATP to terminate VA, and 17 of 33 (51.5%) patients required shock therapy. Mean entropy was significantly higher in the shock versus successful ATP group (6.1 ± 0.5 vs 5.5 ± 0.7, P = .037). Analysis of patients receiving ATP (n = 22) showed significantly higher mean entropy in the six of 22 patients that failed ATP (followed by rescue ICD shock) compared to 16 of 22 that had successful ATP (6.3 ± 0.7 vs 5.5 ± 0.7, P = .048). Computer modeling suggested inability of the paced wavefront in ATP to successfully propagate from the electrode site through patchy fibrosis as a possible mechanism of failed ATP. CONCLUSIONS: Our findings suggest lower scar heterogeneity (mean entropy) is associated with successful ATP, whereas higher scar heterogeneity is associated with more aggressive VAs unresponsive to ATP requiring shock therapy that may be due to inability of the paced wavefront to propagate through scar and terminate the VA circuit.

Evaluation of a real-time magnetic resonance imaging-guided electrophysiology system for structural and electrophysiological ventricular tachycardia substrate assessment.

AIMS: Potential advantages of real-time magnetic resonance imaging (MRI)-guided electrophysiology (MR-EP) include contemporaneous three-dimensional substrate assessment at the time of intervention, improved procedural guidance, and ablation lesion assessment. We evaluated a novel real-time MR-EP system to perform endocardial voltage mapping and assessment of delayed conduction in a porcine ischaemia-reperfusion model. METHODS AND RESULTS: Sites of low voltage and slow conduction identified using the system were registered and compared to regions of late gadolinium enhancement (LGE) on MRI. The Sorensen-Dice similarity coefficient (DSC) between LGE scar maps and voltage maps was computed on a nodal basis. A total of 445 electrograms were recorded in sinus rhythm (range: 30-186) using the MR-EP system including 138 electrograms from LGE regions. Pacing captured at 103 sites; 47 (45.6%) sites had a stimulus-to-QRS (S-QRS) delay of ≥40 ms. Using conventional (0.5-1.5 mV) bipolar voltage thresholds, the sensitivity and specificity of voltage mapping using the MR-EP system to identify MR-derived LGE was 57% and 96%, respectively. Voltage mapping had a better predictive ability in detecting LGE compared to S-QRS measurements using this system (area under curve: 0.907 vs. 0.840). Using an electrical threshold of 1.5 mV to define abnormal myocardium, the total DSC, scar DSC, and normal myocardium DSC between voltage maps and LGE scar maps was 79.0 ± 6.0%, 35.0 ± 10.1%, and 90.4 ± 8.6%, respectively. CONCLUSION: Low-voltage zones and regions of delayed conduction determined using a real-time MR-EP system are moderately associated with LGE areas identified on MRI.

Effect of scar and pacing location on repolarization in a porcine myocardial infarction model.

Background: The effect of chronic ischemic scar on repolarization is unclear, with conflicting results from human and animal studies. An improved understanding of electrical remodeling within scar and border zone tissue may enhance substrate-guided ablation techniques for treatment of ventricular tachycardia. Computational modeling studies have suggested increased dispersion of repolarization during epicardial, but not endocardial, left ventricular pacing, in close proximity to scar. However, the effect of endocardial pacing near scar in vivo is unknown. Objective: The purpose of this study was to investigate the effect of scar and pacing location on local repolarization in a porcine myocardial infarction model. Methods: Six model pigs underwent late gadolinium enhancement cardiac magnetic resonance (LGE-CMR) imaging followed by electroanatomic mapping of the left ventricular endocardium. LGE-CMR images were registered to the anatomic shell and scar defined by LGE. Activation recovery intervals (ARIs), a surrogate for action potential duration, and local ARI gradients were calculated from unipolar electrograms within areas of late gadolinium enhancement (aLGE) and healthy myocardium. Results: There was no significant difference between aLGE and healthy myocardium in mean ARI (304.20 ± 19.44 ms vs 300.59 ± 19.22 ms; P = .43), ARI heterogeneity (23.32 ± 11.43 ms vs 24.85 ± 12.99 ms; P = .54), or ARI gradients (6.18 ± 2.09 vs 5.66 ± 2.32 ms/mm; P = .39). Endocardial pacing distance from scar did not affect ARI gradients. Conclusion: Our findings suggest that changes in ARI are not an intrinsic property of surviving myocytes within scar, and endocardial pacing close to scar does not affect local repolarization.

3D Electrophysiological Modeling of Interstitial Fibrosis Networks and Their Role in Ventricular Arrhythmias in Non-Ischemic Cardiomyopathy.

OBJECTIVE: Interstitial fibrosis is a pathological expansion of the heart's inter-cellular collagen matrix. It is a potential complication of nonischemic cardiomyopathy (NICM), a class of diseases involving electrical and or mechanical dysfunction of cardiac tissue not caused by atherosclerosis. Patients with NICM and interstitial fibrosis often suffer from life threatening arrhythmias, which we aim to simulate in this study. METHODS: Our methodology builds on an efficient discrete finite element (DFE) method which allows for the representation of fibrosis as infinitesimal splits in a mesh. We update the DFE method with a local connectivity analysis which creates a consistent topology in the fibrosis network. This is particularly important in nonischemic disease due to the potential presence of large and contiguous fibrotic regions and therefore potentially complex fibrosis networks. RESULTS: In experiments with an image-based model, we demonstrate that our methodology is able to simulate reentrant electrical events associated with cardiac arrhythmias. These reentries depended crucially upon sufficient fibrosis density, which was marked by conduction slowing at high pacing rates. We also created a 2D test-case which demonstrated that fibrosis topologies can modulate transient conduction block, and thereby reentrant activations. CONCLUSION: Ventricular arrhythmias due to interstitial fibrosis in NICM can be efficiently simulated using our methods in medical image based geometries. Furthermore, fibrosis topology modulates transient conduction block, and should be accounted for in electrophysiological simulations with interstitial fibrosis. SIGNIFICANCE: Our study provides methodology which has the potential to predict arrhythmias and to optimize treatments non-invasively for nonischemic cardiomyopathies.

Computational Modeling for Cardiac Resynchronization Therapy.

Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.

Mind the Gap: A Semicontinuum Model for Discrete Electrical Propagation in Cardiac Tissue.

Electrical propagation in cardiac tissue is a discrete or discontinuous phenomenon that reflects the complexity of the anatomical structures and their organization in the heart, such as myocytes, gap junctions, microvessels, and extracellular matrix, just to name a few. Discrete models or microscopic and discontinuous models are, so far, the best options to accurately study how structural properties of cardiac tissue influence electrical propagation. These models are, however, inappropriate in the context of large scale simulations, which have been traditionally performed by the use of continuum and macroscopic models, such as the monodomain and the bidomain models. However, continuum models may fail to reproduce many important physiological and physiopathological aspects of cardiac electrophysiology, for instance, those related to slow conduction. In this study, we develop a new mathematical model that combines characteristics of both continuum and discrete models. The new model was evaluated in scenarios of low gap-junctional coupling, where slow conduction is observed, and was able to reproduce conduction block, increase of the maximum upstroke velocity and of the repolarization dispersion. None of these features can be captured by continuum models. In addition, the model overcomes a great disadvantage of discrete models, as it allows variation of the spatial resolution within a certain range.

An efficient finite element approach for modeling fibrotic clefts in the heart.

Advanced medical imaging technologies provide a wealth of information on cardiac anatomy and structure at a paracellular resolution, allowing to identify microstructural discontinuities which disrupt the intracellular matrix. Current state-of-the-art computer models built upon such datasets account for increasingly finer anatomical details, however, structural discontinuities at the paracellular level are typically discarded in the model generation process, owing to the significant costs which incur when using high resolutions for explicit representation. In this study, a novel discontinuous finite element (dFE) approach for discretizing the bidomain equations is presented, which accounts for fine-scale structures in a computer model without the need to increase spatial resolution. In the dFE method, this is achieved by imposing infinitely thin lines of electrical insulation along edges of finite elements which approximate the geometry of discontinuities in the intracellular matrix. Simulation results demonstrate that the dFE approach accounts for effects induced by microscopic size scale discontinuities, such as the formation of microscopic virtual electrodes, with vast computational savings as compared to high resolution continuous finite element models. Moreover, the method can be implemented in any standard continuous finite element code with minor effort.

Parametrization strategies for matching activation sequences in models of ventricular electrophysiology.

Driven by recent advances in medical imaging, image segmentation and numerical techniques computer models of ventricular electrophysiology account for increasingly finer levels of anatomical and biophysical detail. However, considering the large number of model parameters involved parametrization poses a major challenge. A minimum requirement in combined experimental and modeling studies which aim at making specific predictions on a case by case basis is to achieve good agreement in activation and repolarization sequences between model and experiment or patient data. In this study we propose basic techniques which aide in determining bidomain parameters to match ventricular activation sequences. Two specific aspects will be considered. First, conduction velocity in the ventricles is orthotropic and varies in space. An iterative parametrization algorithm is implemented which determines appropriate conductivities which yield prescribed velocities. Secondly, impulse propagation in the ventricles is initiated subendocardially at Purkinje-ventricular junctions, the terminal endings of Purkinje system (PS), and, thus, the PS plays a key role in determining the shape of activation wave fronts as reflected in the QRS complex of the electro-cardiogram (ECG). While ventricular models equipped with generic PS topologies match well with experimental observation in terms of epicardial breakthrough sites, predicted ECGs match poorly with known key ECG characteristics.

Automatic Parameterization Strategy for Cardiac Electrophysiology Simulations.

Driven by recent advances in medical imaging, image segmentation and numerical techniques, computer models of ventricular electrophysiology account for increasingly finer levels of anatomical and biophysical detail. However, considering the large number of model parameters involved parameterization poses a major challenge. A minimum requirement in combined experimental and modeling studies is to achieve good agreement in activation and repolarization sequences between model and experiment or patient data. In this study, we propose basic techniques which aid in determining bidomain parameters to match activation sequences. An iterative parameterization algorithm is implemented which determines appropriate bulk conductivities which yield prescribed velocities. In addition, a method is proposed for splitting the computed bulk conductivities into individual bidomain conductivities by prescribing anisotropy ratios.

A finite element approach for modeling micro-structural discontinuities in the heart.

The presence of connective tissue as well as interstitial clefts forms a natural barrier to the electrical propagation in the heart. At a microscopic scale, such uncoupling structures change the pattern of the electrical conduction from uniform towards complex and may play a role in the genesis of cardiac arrhythmias. The anatomical diversity of conduction structures and their topology at a microscopic size scale is overwhelming for experimental techniques. Mathematical models have been often employed to study the behavior of the electrical propagation at a sub-cellular level. However, very fine and computationally expensive meshes are required to capture all microscopic details found in the cardiac tissue. In this work, we present a numerical technique based on the finite element method which allows to reproduce the effects of microscopic conduction barriers caused by the presence of uncoupling structures without actually resolving these structures in a high resolution mesh, thereby reducing the computational costs significantly.

Limitations of the homogenized cardiac Monodomain model for the case of low gap junctional coupling.

The cardiac Monodomain model is a mathematical model extensively used in studies of propagation of bioelectric wavefronts in the heart. To be able to use the model for complex and large cardiac simulations, such as the case of whole heart and 3D simulations, some parameters of the model that are known to physiologically vary in space, such as the intracellular conductivity, are traditionally kept constant at effective values. These effective values can be obtained via a mathematical procedure called homogenization. In this work we revisit the classical homogenized monodomain formulation to evaluate its ability to reproduce the situation of low gap junctional coupling. This situation arises in many pathological conditions such as during ischemia. Our numerical results suggest some limitations of the homogenized cardiac Monodomain model under these conditions in terms of computed conduction velocity and Action Potential waverforms.