Arrhythmogenic vulnerability of re-entrant pathways in post-infarct ventricular tachycardia assessed by advanced computational modelling.

AIMS: Substrate assessment of scar-mediated ventricular tachycardia (VT) is frequently performed using late gadolinium enhancement (LGE) images. Although this provides structural information about critical pathways through the scar, assessing the vulnerability of these pathways for sustaining VT is not possible with imaging alone.This study evaluated the performance of a novel automated re-entrant pathway finding algorithm to non-invasively predict VT circuit and inducibility. METHODS: Twenty post-infarct VT-ablation patients were included for retrospective analysis. Commercially available software (ADAS3D left ventricular) was used to generate scar maps from 2D-LGE images using the default 40-60 pixel-signal-intensity (PSI) threshold. In addition, algorithm sensitivity for altered thresholds was explored using PSI 45-55, 35-65, and 30-70. Simulations were performed on the Virtual Induction and Treatment of Arrhythmias (VITA) framework to identify potential sites of block and assess their vulnerability depending on the automatically computed round-trip-time (RTT). Metrics, indicative of substrate complexity, were correlated with VT-recurrence during follow-up. RESULTS: Total VTs (85 ± 43 vs. 42 ± 27) and unique VTs (9 ± 4 vs. 5 ± 4) were significantly higher in patients with- compared to patients without recurrence, and were predictive of recurrence with area under the curve of 0.820 and 0.770, respectively. VITA was robust to scar threshold variations with no significant impact on total and unique VTs, and mean RTT between the four models. Simulation metrics derived from PSI 45-55 model had the highest number of parameters predictive for post-ablation VT-recurrence. CONCLUSION: Advanced computational metrics can non-invasively and robustly assess VT substrate complexity, which may aid personalized clinical planning and decision-making in the treatment of post-infarction VT.

Personalization of electro-mechanical models of the pressure-overloaded left ventricle: fitting of Windkessel-type afterload models.

Computer models of left ventricular (LV) electro-mechanics (EM) show promise as a tool for assessing the impact of increased afterload upon LV performance. However, the identification of unique afterload model parameters and the personalization of EM LV models remains challenging due to significant clinical input uncertainties. Here, we personalized a virtual cohort of N = 17 EM LV models under pressure overload conditions. A global-local optimizer was developed to uniquely identify parameters of a three-element Windkessel (Wk3) afterload model. The sensitivity of Wk3 parameters to input uncertainty and of the EM LV model to Wk3 parameter uncertainty was analysed. The optimizer uniquely identified Wk3 parameters, and outputs of the personalized EM LV models showed close agreement with clinical data in all cases. Sensitivity analysis revealed a strong dependence of Wk3 parameters on input uncertainty. However, this had limited impact on outputs of EM LV models. A unique identification of Wk3 parameters from clinical data appears feasible, but it is sensitive to input uncertainty, thus depending on accurate invasive measurements. By contrast, the EM LV model outputs were less sensitive, with errors of less than 8.14% for input data errors of 10%, which is within the bounds of clinical data uncertainty. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Patient-specific modeling of left ventricular electromechanics as a driver for haemodynamic analysis.

AIMS: Models of blood flow in the left ventricle (LV) and aorta are an important tool for analysing the interplay between LV deformation and flow patterns. Typically, image-based kinematic models describing endocardial motion are used as an input to blood flow simulations. While such models are suitable for analysing the hemodynamic status quo, they are limited in predicting the response to interventions that alter afterload conditions. Mechano-fluidic models using biophysically detailed electromechanical (EM) models have the potential to overcome this limitation, but are more costly to build and compute. We report our recent advancements in developing an automated workflow for the creation of such CFD ready kinematic models to serve as drivers of blood flow simulations. METHODS AND RESULTS: EM models of the LV and aortic root were created for four pediatric patients treated for either aortic coarctation or aortic valve disease. Using MRI, ECG and invasive pressure recordings, anatomy as well as electrophysiological, mechanical and circulatory model components were personalized. RESULTS: The implemented modeling pipeline was highly automated and allowed model construction and execution of simulations of a patient's heartbeat within 1 day. All models reproduced clinical data with acceptable accuracy. CONCLUSION: Using the developed modeling workflow, the use of EM LV models as driver of fluid flow simulations is becoming feasible. While EM models are costly to construct, they constitute an important and nontrivial step towards fully coupled electro-mechano-fluidic (EMF) models and show promise as a tool for predicting the response to interventions which affect afterload conditions.

pyCEPS: A cross-platform electroanatomic mapping data to computational model conversion platform for the calibration of digital twin models of cardiac electrophysiology.

BACKGROUND AND OBJECTIVE: Data from electro-anatomical mapping (EAM) systems are playing an increasingly important role in computational modeling studies for the patient-specific calibration of digital twin models. However, data exported from commercial EAM systems are challenging to access and parse. Converting to data formats that are easily amenable to be viewed and analyzed with commonly used cardiac simulation software tools such as openCARP remains challenging. We therefore developed an open-source platform, pyCEPS, for parsing and converting clinical EAM data conveniently to standard formats widely adopted within the cardiac modeling community. METHODS AND RESULTS: pyCEPS is an open-source Python-based platform providing the following functions: (i) access and interrogate the EAM data exported from clinical mapping systems; (ii) efficient browsing of EAM data to preview mapping procedures, electrograms (EGMs), and electro-cardiograms (ECGs); (iii) conversion to modeling formats according to the openCARP standard, to be amenable to analysis with standard tools and advanced workflows as used for in silico EAM data. Documentation and training material to facilitate access to this complementary research tool for new users is provided. We describe the technological underpinnings and demonstrate the capabilities of pyCEPS first, and showcase its use in an exemplary modeling application where we use clinical imaging data to build a patient-specific anatomical model. CONCLUSION: With pyCEPS we offer an open-source framework for accessing EAM data, and converting these to cardiac modeling standard formats. pyCEPS provides the core functionality needed to integrate EAM data in cardiac modeling research. We detail how pyCEPS could be integrated into model calibration workflows facilitating the calibration of a computational model based on EAM data.

ForCEPSS-A framework for cardiac electrophysiology simulations standardization.

BACKGROUND AND OBJECTIVE: Simulation of cardiac electrophysiology (CEP) is an important research tool that is increasingly being adopted in industrial and clinical applications. Typical workflows for CEP simulation consist of a sequence of processing stages starting with building an anatomical model and then calibrating its electrophysiological properties to match observable data. While the calibration stages are common and generalizable, most CEP studies re-implement these steps in complex and highly variable workflows. This lack of standardization renders the execution of computational CEP studies in an efficient, robust, and reproducible manner a significant challenge. Here, we propose ForCEPSS as an efficient and robust, yet flexible, software framework for standardizing CEP simulation studies. METHODS AND RESULTS: Key processing stages of CEP simulation studies are identified and implemented in a standardized workflow that builds on openCARP1 Plank et al. (2021) and the Python-based carputils2 framework. Stages include (i) the definition and initialization of action potential phenotypes, (ii) the tissue scale calibration of conduction properties, (iii) the functional initialization to approximate a limit cycle corresponding to the dynamic reference state according to an experimental protocol, and, (iv) the execution of the CEP study where the electrophysiological response to a perturbation of the limit cycle is probed. As an exemplar application, we employ ForCEPSS to prepare a CEP study according to the Virtual Arrhythmia Risk Prediction protocol used for investigating the arrhythmogenic risk of developing infarct-related ventricular tachycardia (VT) in ischemic cardiomyopathy patients. We demonstrate that ForCEPSS enables a fully automated execution of all stages of this complex protocol. CONCLUSION: ForCEPSS offers a novel comprehensive, standardized, and automated CEP simulation workflow. The high degree of automation accelerates the execution of CEP simulation studies, reduces errors, improves robustness, and makes CEP studies reproducible. Verification of simulation studies within the CEP modeling community is thus possible. As such, ForCEPSS makes an important contribution towards increasing transparency, standardization, and reproducibility of in silico CEP experiments.

Hybrid Neural State-Space Modeling for Supervised and Unsupervised Electrocardiographic Imaging.

State-space modeling (SSM) provides a general framework for many image reconstruction tasks. Error in a priori physiological knowledge of the imaging physics, can bring incorrectness to solutions. Modern deep-learning approaches show great promise but lack interpretability and rely on large amounts of labeled data. In this paper, we present a novel hybrid SSM framework for electrocardiographic imaging (ECGI) to leverage the advantage of state-space formulations in data-driven learning. We first leverage the physics-based forward operator to supervise the learning. We then introduce neural modeling of the transition function and the associated Bayesian filtering strategy. We applied the hybrid SSM framework to reconstruct electrical activity on the heart surface from body-surface potentials. In unsupervised settings of both in-silico and in-vivo data without cardiac electrical activity as the ground truth to supervise the learning, we demonstrated improved ECGI performances of the hybrid SSM framework trained from a small number of ECG observations in comparison to the fixed SSM. We further demonstrated that, when in-silico simulation data becomes available, mixed supervised and unsupervised training of the hybrid SSM achieved a further 40.6% and 45.6% improvements, respectively, in comparison to traditional ECGI baselines and supervised data-driven ECGI baselines for localizing the origin of ventricular activations in real data.

Predicting postinfarct ventricular tachycardia by integrating cardiac MRI and advanced computational reentrant pathway analysis.

BACKGROUND: Implantable cardiac defibrillator (ICD) implantation can protect against sudden cardiac death after myocardial infarction. However, improved risk stratification for device requirement is still needed. OBJECTIVE: The purpose of this study was to improve assessment of postinfarct ventricular electropathology and prediction of appropriate ICD therapy by combining late gadolinium enhancement (LGE) and advanced computational modeling. METHODS: ADAS 3D LV (ADAS LV Medical, Barcelona, Spain) and custom-made software were used to generate 3-dimensional patient-specific ventricular models in a prospective cohort of patients with a myocardial infarction (N = 40) having undergone LGE imaging before ICD implantation. Corridor metrics and 3-dimensional surface features were computed from LGE images. The Virtual Induction and Treatment of Arrhythmias (VITA) framework was applied to patient-specific models to comprehensively probe the vulnerability of the scar substrate to sustaining reentrant circuits. Imaging and VITA metrics, related to the numbers of induced ventricular tachycardias and their corresponding round trip times (RTTs), were compared with ICD therapy during follow-up. RESULTS: Patients with an event (n = 17) had a larger interface between healthy myocardium and scar and higher VITA metrics. Cox regression analysis demonstrated a significant independent association with an event: interface (hazard ratio [HR] 2.79; 95% confidence interval [CI] 1.44-5.44; P < .01), unique ventricular tachycardias (HR 1.67; 95% CI 1.04-2.68; P = .03), mean RTT (HR 2.14; 95% CI 1.11-4.12; P = .02), and maximum RTT (HR 2.13; 95% CI 1.19-3.81; P = .01). CONCLUSION: A detailed quantitative analysis of LGE-based scar maps, combined with advanced computational modeling, can accurately predict ICD therapy and could facilitate the early identification of high-risk patients in addition to left ventricular ejection fraction.

Frequency Domain Analysis of Endocardial Electrograms for Detection of Nontransmural Myocardial Fibrosis in Nonischemic Cardiomyopathy.

BACKGROUND: Voltage mapping in nonischemic cardiomyopathy can fail to identify midmyocardial substrate for ventricular arrhythmias, an important cause of ablation failure. OBJECTIVES: The aim of this study was to assess whether frequency domain analysis of endocardial left ventricular electrograms (EGMs) can better predict the presence of midmyocardial fibrosis (MMF) compared with voltage amplitude. METHODS: Nonischemic cardiomyopathy patients undergoing ventricular tachycardia ablation with registered preprocedural cardiac computed tomography and late iodine enhancement were included. Presence of fibrosis at each EGM site was assessed. Bipolar and unipolar EGMs were transformed to the frequency domain using multitaper spectral analysis. Singular value decomposition of the EGM frequency spectrum was used within a supervised machine learning process to select features to predict the presence of MMF and compare against predictions using voltage amplitude. RESULTS: Thirteen patients were included (median age 57 years [IQR: 28-73 years], median ejection fraction 40% [IQR: 15%-57%]). A total of 6,015 EGM pairs were processed: 2,459 EGM pairs in MMF areas and 3,556 EGM pairs in non-MMF areas. Supervised classifiers were trained with stratified k-fold cross-validation within patients. The distribution of mean area under the curve metrics using frequency features, f, was significantly greater than voltage feature area under the curve metrics, v, (mean f = 0.841 [95% CI: 0.789-0.884] vs mean v = 0.591 [95% CI: 0.530-0.658]; P < 0.001), indicating that frequency-trained classifiers better predicted the presence of MMF. CONCLUSIONS: These data indicate the promising discriminatory value of endocardial EGM frequency content in the assessment of concealed myocardial substrate. Further studies are needed to investigate the importance of the specific frequency features identified.

Predicting arrhythmia recurrence following catheter ablation for ventricular tachycardia using late gadolinium enhancement magnetic resonance imaging: Implications of varying scar ranges.

BACKGROUND: Thresholding-based analysis of late gadolinium enhancement cardiac magnetic resonance (LGE-CMR) can create scar maps and identify corridors that might provide a reentrant substrate for ventricular tachycardia (VT). Current recommendations use a full-width-at-half-maximum approach, effectively classifying areas with a pixel signal intensity (PSI) >40% as border zone (BZ) and >60% as core. OBJECTIVE: The purpose of this study was to investigate the impact of 4 different threshold settings on scar and corridor quantification and to correlate this with postablation VT recurrence. METHODS: Twenty-seven patients with ischemic cardiomyopathy who had undergone catheter ablation for VT were included for retrospective analysis. LGE-CMR images were analyzed using ADAS3D LV. Scar maps were created for 4 PSI thresholds (40-60, 35-65, 30-70, and 45-55), and the extent of variation in BZ and core, as well as the number and weight of conduction corridors, were quantified. Three-dimensional representations were reconstructed from exported segmentations and used to quantify the surface area between healthy myocardium and scar (BZ + core), and between BZ and core. RESULTS: A wider PSI threshold was associated with an increase in BZ mass and decrease in scar (P <.001). No significant differences were observed for the total number of corridors and their mass with increasing PSI threshold. The best correlation in predicting arrhythmia recurrence was observed for PSI 45-55 (area under the curve 0.807; P = .001). CONCLUSION: Varying PSI has a significant impact on quantification of LGE-CMR parameters and may have incremental clinical value in predicting arrhythmia recurrence. Further prospective investigation is warranted to clarify the functional implications of these findings for LGE-CMR-guided ventricular ablation.

A personalized real-time virtual model of whole heart electrophysiology.

Computer models capable of representing the intrinsic personal electrophysiology (EP) of the heart in silico are termed virtual heart technologies. When anatomy and EP are tailored to individual patients within the model, such technologies are promising clinical and industrial tools. Regardless of their vast potential, few virtual technologies simulating the entire organ-scale EP of all four-chambers of the heart have been reported and widespread clinical use is limited due to high computational costs and difficulty in validation. We thus report on the development of a novel virtual technology representing the electrophysiology of all four-chambers of the heart aiming to overcome these limitations. In our previous work, a model of ventricular EP embedded in a torso was constructed from clinical magnetic resonance image (MRI) data and personalized according to the measured 12 lead electrocardiogram (ECG) of a single subject under normal sinus rhythm. This model is then expanded upon to include whole heart EP and a detailed representation of the His-Purkinje system (HPS). To test the capacities of the personalized virtual heart technology to replicate standard clinical morphological ECG features under such conditions, bundle branch blocks within both the right and the left ventricles under two different conduction velocity settings are modeled alongside sinus rhythm. To ensure clinical viability, model generation was completely automated and simulations were performed using an efficient real-time cardiac EP simulator. Close correspondence between the measured and simulated 12 lead ECG was observed under normal sinus conditions and all simulated bundle branch blocks manifested relevant clinical morphological features.

Automated Framework for the Inclusion of a His-Purkinje System in Cardiac Digital Twins of Ventricular Electrophysiology.

Personalized models of cardiac electrophysiology (EP) that match clinical observation with high fidelity, referred to as cardiac digital twins (CDTs), show promise as a tool for tailoring cardiac precision therapies. Building CDTs of cardiac EP relies on the ability of models to replicate the ventricular activation sequence under a broad range of conditions. Of pivotal importance is the His-Purkinje system (HPS) within the ventricles. Workflows for the generation and incorporation of HPS models are needed for use in cardiac digital twinning pipelines that aim to minimize the misfit between model predictions and clinical data such as the 12 lead electrocardiogram (ECG). We thus develop an automated two stage approach for HPS personalization. A fascicular-based model is first introduced that modulates the endocardial Purkinje network. Only emergent features of sites of earliest activation within the ventricular myocardium and a fast-conducting sub-endocardial layer are accounted for. It is then replaced by a topologically realistic Purkinje-based representation of the HPS. Feasibility of the approach is demonstrated. Equivalence between both HPS model representations is investigated by comparing activation patterns and 12 lead ECGs under both sinus rhythm and right-ventricular apical pacing. Predominant ECG morphology is preserved by both HPS models under sinus conditions, but elucidates differences during pacing.

The openCARP simulation environment for cardiac electrophysiology.

BACKGROUND AND OBJECTIVE: Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community. METHODS AND RESULTS: openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential. CONCLUSION: As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.

A Framework for the generation of digital twins of cardiac electrophysiology from clinical 12-leads ECGs.

Cardiac digital twins (Cardiac Digital Twin (CDT)s) of human electrophysiology (Electrophysiology (EP)) are digital replicas of patient hearts derived from clinical data that match like-for-like all available clinical observations. Due to their inherent predictive potential, CDTs show high promise as a complementary modality aiding in clinical decision making and also in the cost-effective, safe and ethical testing of novel EP device therapies. However, current workflows for both the anatomical and functional twinning phases within CDT generation, referring to the inference of model anatomy and parameters from clinical data, are not sufficiently efficient, robust and accurate for advanced clinical and industrial applications. Our study addresses three primary limitations impeding the routine generation of high-fidelity CDTs by introducing; a comprehensive parameter vector encapsulating all factors relating to the ventricular EP; an abstract reference frame within the model allowing the unattended manipulation of model parameter fields; a novel fast-forward electrocardiogram (Electrocardiogram (ECG)) model for efficient and bio-physically-detailed simulation required for parameter inference. A novel workflow for the generation of CDTs is then introduced as an initial proof of concept. Anatomical twinning was performed within a reasonable time compatible with clinical workflows (<4h) for 12 subjects from clinically-attained magnetic resonance images. After assessment of the underlying fast forward ECG model against a gold standard bidomain ECG model, functional twinning of optimal parameters according to a clinically-attained 12 lead ECG was then performed using a forward Saltelli sampling approach for a single subject. The achieved results in terms of efficiency and fidelity demonstrate that our workflow is well-suited and viable for generating biophysically-detailed CDTs at scale.

Automating image-based mesh generation and manipulation tasks in cardiac modeling workflows using Meshtool.

Advanced cardiac modeling studies rely on the ability to generate and functionalize personalized in silico models from tomographic multi-label image stacks. Eventually, this is used for building virtual cohorts that capture the variability in size, shape, and morphology of individual hearts. Typical modeling workflows involve a multitude of interactive mesh manipulation steps, rendering model generation expensive. Meshtool is software specifically designed for automating all complex mesh manipulation tasks emerging in such workflows by implementing algorithms for tasks describable as operations on label fields and/or geometric features. We illustrate how Meshtool increases efficiency and reduces costs by offering an automatable, high performance mesh manipulation toolbox.

Double-kissing culotte technique for coronary bifurcation stenting.

AIMS: The aim of this study was to assess whether the culotte technique could be improved by an additional kissing dilation prior to main branch (MB) stenting. METHODS AND RESULTS: Double-kissing (DK) culotte was compared to the culotte and DK-crush techniques in a bench model (n=24). Results were evaluated for stent apposition, luminal opening and flow dynamics. The total procedure duration of DK-culotte was 18.3±3.4 minutes, significantly lower than for DK-crush (24.3±5.7 min; p=0.015), but similar to culotte (21.6±5.9 min, p=0.104). In DK-culotte the overall rate of moderate (200-500 µm) and significant (>500 µm) malapposition was 2.1±1.9% and 0.4±0.2%, similar as compared to culotte (3.7±3.8%, p=0.459 and 1.0±1.0%, p=0.517, respectively), and lower as compared to DK-crush (8.1±2.5%, p<0.001 and 3.7±5.3%, p=0.002, respectively). The lower malapposition rate of DK-culotte as compared to DK-crush was due to less moderate and significant malapposition in the proximal MB (0.0±0.0% vs 14.0±7.6%, p<0.001 and 0.0±0.0% vs 4.2±9.1%, p=0.026, respectively). Micro-computed tomography did not show a difference in luminal opening at the proximal MB, distal MB or SB. There was no difference either in the maximum shear rate or in areas of high shear or recirculation. CONCLUSIONS: Bench test data suggest that the DK approach facilitates the culotte technique. The clinical validity and relevance remain to be confirmed in a larger in vivo population.

Simulating ventricular systolic motion in a four-chamber heart model with spatially varying robin boundary conditions to model the effect of the pericardium.

The pericardium affects cardiac motion by limiting epicardial displacement normal to the surface. In computational studies, it is important for the model to replicate realistic motion, as this affects the physiological fidelity of the model. Previous computational studies showed that accounting for the effect of the pericardium allows for a more realistic motion simulation. In this study, we describe the mechanism through which the pericardium causes improved cardiac motion. We simulated electrical activation and contraction of the ventricles on a four-chamber heart in the presence and absence of the effect of the pericardium. We simulated the mechanical constraints imposed by the pericardium by applying normal Robin boundary conditions on the ventricular epicardium. We defined a regional scaling of normal springs stiffness based on image-derived motion from CT images. The presence of the pericardium reduced the error between simulated and image-derived end-systolic configurations from 12.8±4.1 mm to 5.7±2.5 mm. First, the pericardium prevents the ventricles from spherising during isovolumic contraction, reducing the outward motion of the free walls normal to the surface and the upwards motion of the apex. Second, by restricting the inward motion of the free and apical walls of the ventricles the pericardium increases atrioventricular plane displacement by four folds during ejection. Our results provide a mechanistic explanation of the importance of the pericardium in physiological simulations of electromechanical cardiac function.

A publicly available virtual cohort of four-chamber heart meshes for cardiac electro-mechanics simulations.

Computational models of the heart are increasingly being used in the development of devices, patient diagnosis and therapy guidance. While software techniques have been developed for simulating single hearts, there remain significant challenges in simulating cohorts of virtual hearts from multiple patients. To facilitate the development of new simulation and model analysis techniques by groups without direct access to medical data, image analysis techniques and meshing tools, we have created the first publicly available virtual cohort of twenty-four four-chamber hearts. Our cohort was built from heart failure patients, age 67±14 years. We segmented four-chamber heart geometries from end-diastolic (ED) CT images and generated linear tetrahedral meshes with an average edge length of 1.1±0.2mm. Ventricular fibres were added in the ventricles with a rule-based method with an orientation of -60° and 80° at the epicardium and endocardium, respectively. We additionally refined the meshes to an average edge length of 0.39±0.10mm to show that all given meshes can be resampled to achieve an arbitrary desired resolution. We ran simulations for ventricular electrical activation and free mechanical contraction on all 1.1mm-resolution meshes to ensure that our meshes are suitable for electro-mechanical simulations. Simulations for electrical activation resulted in a total activation time of 149±16ms. Free mechanical contractions gave an average left ventricular (LV) and right ventricular (RV) ejection fraction (EF) of 35±1% and 30±2%, respectively, and a LV and RV stroke volume (SV) of 95±28mL and 65±11mL, respectively. By making the cohort publicly available, we hope to facilitate large cohort computational studies and to promote the development of cardiac computational electro-mechanics for clinical applications.

Image-based models of cardiac structure with applications in arrhythmia and defibrillation studies.

The objective of this article is to present a set of methods for constructing realistic computational models of cardiac structure from high-resolution structural and diffusion tensor magnetic resonance images and to demonstrate the applicability of the models in simulation studies. The structural image is segmented to identify various regions such as normal myocardium, ventricles, and infarct. A finite element mesh is generated from the processed structural data, and fiber orientations are assigned to the elements. The Purkinje system, when visible, is modeled using linear elements that interconnect a set of manually identified points. The methods were applied to construct 2 different models; and 2 simulation studies, which demonstrate the applicability of the models in the analysis of arrhythmia and defibrillation, were performed. The models represent cardiac structure with unprecedented detail for simulation studies.

Universal ventricular coordinates: A generic framework for describing position within the heart and transferring data.

Being able to map a particular set of cardiac ventricles to a generic topologically equivalent representation has many applications, including facilitating comparison of different hearts, as well as mapping quantities and structures of interest between them. In this paper we describe Universal Ventricular Coordinates (UVC), which can be used to describe position within any biventricular heart. UVC comprise four unique coordinates that we have chosen to be intuitive, well defined, and relevant for physiological descriptions. We describe how to determine these coordinates for any volumetric mesh by illustrating how to properly assign boundary conditions and utilize solutions to Laplace's equation. Using UVC, we transferred scalar, vector, and tensor data between four unstructured ventricular meshes from three different species. Performing the mappings was very fast, on the order of a few minutes, since mesh nodes were searched in a KD tree. Distance errors in mapping mesh nodes back and forth between meshes were less than the size of an element. Analytically derived fiber directions were also mapped across meshes and compared, showing  < 5° difference over most of the ventricles. The ability to transfer gradients was also demonstrated. Topologically variable structures, like papillary muscles, required further definition outside of the UVC framework. In conclusion, UVC can aid in transferring many types of data between different biventricular geometries.

Assessment of wall stresses and mechanical heart power in the left ventricle: Finite element modeling versus Laplace analysis.

INTRODUCTION: Stenotic aortic valve disease (AS) causes pressure overload of the left ventricle (LV) that may trigger adverse remodeling and precipitate progression towards heart failure (HF). As myocardial energetics can be impaired during AS, LV wall stresses and biomechanical power provide a complementary view of LV performance that may aide in better assessing the state of disease. OBJECTIVES: Using a high-resolution electro-mechanical (EM) in silico model of the LV as a reference, we evaluated clinically feasible Laplace-based methods for assessing global LV wall stresses and biomechanical power. METHODS: We used N = 4 in silico finite element (FE) EM models of LV and aorta of patients suffering from AS. All models were personalized with clinical data under pretreatment conditions. Left ventricle wall stresses and biomechanical power were computed accurately from FE kinematic data and compared with Laplace-based estimation methods, which were applied to the same FE model data. RESULTS AND CONCLUSION: Laplace estimates of LV wall stress are able to provide a rough approximation of global mean stress in the circumferential-longitudinal plane of the LV. However, according to FE results, spatial heterogeneity of stresses in the LV wall is significant, leading to major discrepancies between local stresses and global mean stress. Assessment of mechanical power with Laplace methods is feasible, but these are inferior in accuracy compared with FE models. The accurate assessment of stress and power density distribution in the LV wall is only feasible based on patient-specific FE modeling.