A rule-based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts

Harnessing 12-lead ECG and MRI data to personalise repolarisation profiles in cardiac digital twin models for enhanced virtual drug testing

Cardiac digital twins are computational tools capturing key functional and anatomical characteristics of patient hearts for investigating disease phenotypes and predicting responses to therapy. When paired with large-scale computational resources and large clinical datasets, digital twin technology can enable virtual clinical trials on virtual cohorts to fast-track therapy development. Here, we present an open-source automated pipeline for personalising ventricular electrophysiological function based on routinely acquired magnetic resonance imaging (MRI) data and the standard 12-lead electrocardiogram (ECG). Using MRI-based anatomical models, a sequential Monte-Carlo approximate Bayesian computational inference method is extended to infer electrical activation and repolarisation characteristics from the ECG. Fast simulations are conducted with a reaction-Eikonal model, including the Purkinje network and biophysically-detailed subcellular ionic current dynamics for repolarisation. For each patient, parameter uncertainty is represented by inferring an envelope of plausible ventricular models rather than a single one, which means that parameter uncertainty can be propagated to therapy evaluation. Furthermore, we have developed techniques for translating from reaction-Eikonal to monodomain simulations, which allows more realistic simulations of cardiac electrophysiology. The pipeline is demonstrated in three healthy subjects, where our inferred pseudo-diffusion reaction-Eikonal models reproduced the patient's ECG with a median Pearson's correlation coefficient of 0.9, and then translated to monodomain simulations with a median correlation coefficient of 0.84 across all subjects. We then demonstrate our digital twins for virtual evaluation of Dofetilide with uncertainty quantification. These evaluations using our cardiac digital twins reproduced dose-dependent QTc and T peak to T end prolongations that are in keeping with large population drug response data. The methodologies for cardiac digital twinning presented here are a step towards personalised virtual therapy testing and can be scaled to generate virtual populations for clinical trials to fast-track therapy evaluation. The tools developed for this paper are open-source, documented, and made publicly available.

Fast and accurate prediction of drug induced proarrhythmic risk with sex specific cardiac emulators

In silico trials for drug safety assessment require many high-fidelity 3D cardiac simulations to predict drug-induced QT interval prolongation, which is often computationally prohibitive. To streamline this process, we developed sex-specific emulators for a fast prediction of QT interval, trained on a dataset of 900 simulations. Our results show significant differences between 3D and 0D single-cell models as risk levels increase, underscoring the ability of 3D modeling to capture more complex cardiac responses. The emulators demonstrated an average error of 4% compared to simulations, allowing for efficient global sensitivity analysis and fast replication of in silico clinical trials. This approach enables rapid, multi-dose drug testing on standard hardware, addressing critical industry challenges around trial design, assay variability, and cost-effective safety evaluations. By integrating these emulators into drug development, we can improve preclinical reliability and advance the practical application of digital twins in biomedicine.

A Matlab Toolbox for cardiac electrophysiology simulations on patient-specific geometries

In this paper, we present CardioMat, a Matlab toolbox for cardiac electrophysiology simulation based on patient-specific anatomies. The strength of CardioMat is the easy and fast construction of electrophysiology cardiac digital twins from segmented anatomical images in a general-purpose software such as Matlab. CardioMat implements a quasi-automatic pipeline that guides the user toward the construction of anatomically detailed cardiac electrophysiology models. Importantly, the CardioMat framework includes the generation of physiologically plausible fiber orientation and Purkinje networks. The main novelty of our framework is its ability to handle voxel-based geometries as produced by segmentation procedures directly, without the need for an unstructured mesh. Indeed, the CardioMat monodomain solver uses a smoothed boundary approach and runs completely on GPU for fast simulations. We employed CardioMat in different application scenarios to show its potentialities and provide preliminary assessment of the feasibility, diagnostic performance, and accuracy of the toolbox. In particular, we showed that CardioMat simulations derived from post-infarction patients hold high sensitivity, specificity, predictive value, and accuracy for localization of deceleration zones in sinus rhythm.

A scalable framework for learning the geometry-dependent solution operators of partial differential equations

Solving partial differential equations (PDEs) using numerical methods is a ubiquitous task in engineering and medicine. However, the computational costs can be prohibitively high when many-query evaluations of PDE solutions on multiple geometries are needed. Here we aim to address this challenge by introducing Diffeomorphic Mapping Operator Learning (DIMON), a generic artificial intelligence framework that learns geometry-dependent solution operators of different types of PDE on a variety of geometries. We present several examples to demonstrate the performance, efficiency and scalability of the framework in learning both static and time-dependent PDEs on parameterized and non-parameterized domains; these include solving the Laplace equations, reaction-diffusion equations and a system of multiscale PDEs that characterize the electrical propagation on thousands of personalized heart digital twins. DIMON can reduce the computational costs of solution approximations on multiple geometries from hours to seconds with substantially less computational resources.

ECG analysis of ventricular fibrillation dynamics reflects ischaemic progression subject to variability in patient anatomy and electrode location

Background

Ventricular fibrillation (VF) is the deadliest arrhythmia, often caused by myocardial ischaemia. VF patients require urgent intervention planned quickly and non-invasively. However, the accuracy with which electrocardiographic (ECG) markers reflect the underlying arrhythmic substrate is unknown.

Methods

We analysed how ECG metrics reflect the fibrillatory dynamics of electrical excitation and ischaemic substrate. For this, we developed a human-based computational modelling and simulation framework for the quantification of ECG metrics, namely, frequency, slope, and amplitude spectrum area (AMSA) during VF in acute ischaemia for several electrode configurations. Simulations reproduced experimental and clinical findings in 21 scenarios presenting variability in the location and transmural extent of regional ischaemia, and severity of ischaemia in the remote myocardium secondary to VF.

Results

Regional acute myocardial ischaemia facilitated re-entries, potentially breaking up into VF. Ischaemia in the remote myocardium modulated fibrillation dynamics. Cases presenting a mildly ischaemic remote myocardium yielded sustained VF, enabled by the high proliferation of phase singularities (PS, 11-22) causing remarkably disorganised activation patterns. Conversely, global acute ischaemia induced stable rotors (3-12 PS). Changes in frequency and morphology of the ECG during VF reproduced clinical findings but did not show a direct correlation with the underlying wave dynamics. AMSA allowed the precise stratification of VF according to ischaemic severity in the remote myocardium (healthy: 23.62-24.45 mV Hz; mild ischaemia: 10.58-21.47 mV Hz; moderate ischaemia: 4.82-11.12 mV Hz). Within the context of clinical reference values, apex-anterior and apex-posterior electrode configurations were the most discriminatory in stratifying VF based on the underlying ischaemic substrate.

Conclusion

This in silico study provides further insights into non-invasive patient-specific strategies for assessing acute ventricular arrhythmias. The use of reliable ECG markers to characterise VF is critical for developing tailored resuscitation strategies.

In silico evaluation of cell therapy in acute versus chronic infarction: role of automaticity, heterogeneity and Purkinje in human

Human-based modelling and simulation offer an ideal testbed for novel medical therapies to guide experimental and clinical studies. Myocardial infarction (MI) is a common cause of heart failure and mortality, for which novel therapies are urgently needed. Although cell therapy offers promise, electrophysiological heterogeneity raises pro-arrhythmic safety concerns, where underlying complex spatio-temporal dynamics cannot be investigated experimentally. Here, after demonstrating credibility of the modelling and simulation framework, we investigate cell therapy in acute versus chronic MI and the role of cell heterogeneity, scar size and the Purkinje system. Simulations agreed with experimental and clinical recordings from ionic to ECG dynamics in acute and chronic infarction. Following cell delivery, spontaneous beats were facilitated by heterogeneity in cell populations, chronic MI due to tissue depolarisation and slow sinus rhythm. Subsequent re-entrant arrhythmias occurred, in some instances with Purkinje involvement and their susceptibility was enhanced by impaired Purkinje-myocardium coupling, large scars and acute infarction. We conclude that homogeneity in injected ventricular-like cell populations minimises their spontaneous beating, which is enhanced by chronic MI, whereas a healthy Purkinje-myocardium coupling is key to prevent subsequent re-entrant arrhythmias, particularly for large scars.

Evaluation of a deep learning-enabled automated computational heart modelling workflow for personalized assessment of ventricular arrhythmias

Personalized, image-based computational heart modelling is a powerful technology that can be used to improve patient-specific arrhythmia risk stratification and ventricular tachycardia (VT) ablation targeting. However, most state-of-the-art methods still require manual interactions by expert users. The goal of this study is to evaluate the feasibility of an automated, deep learning-based workflow for reconstructing personalized computational electrophysiological heart models to guide patient-specific treatment of VT. Contrast-enhanced computed tomography (CE-CT) images with expert ventricular myocardium segmentations were acquired from 111 patients across five cohorts from three different institutions. A deep convolutional neural network (CNN) for segmenting left ventricular myocardium from CE-CT was developed, trained and evaluated. From both CNN-based and expert segmentations in a subset of patients, personalized electrophysiological heart models were reconstructed and rapid pacing was used to induce VTs. CNN-based and expert segmentations were more concordant in the middle myocardium than in the heart's base or apex. Wavefront propagation during pacing was similar between CNN-based and original heart models. Between most sets of heart models, VT inducibility was the same, the number of induced VTs was strongly correlated, and VT circuits co-localized. Our results demonstrate that personalized computational heart models reconstructed from deep learning-based segmentations even with a small training set size can predict similar VT inducibility and circuit locations as those from expertly-derived heart models. Hence, a user-independent, automated framework for simulating arrhythmias in personalized heart models could feasibly be used in clinical settings to aid VT risk stratification and guide VT ablation therapy. KEY POINTS: Personalized electrophysiological heart modelling can aid in patient-specific ventricular tachycardia (VT) risk stratification and VT ablation targeting. Current state-of-the-art, image-based heart models for VT prediction require expert-dependent, manual interactions that may not be accessible across clinical settings. In this study, we develop an automated, deep learning-based workflow for reconstructing personalized heart models capable of simulating arrhythmias and compare its predictions with that of expert-generated heart models. The number and location of VTs was similar between heart models generated from the deep learning-based workflow and expert-generated heart models. These results demonstrate the feasibility of using an automated computational heart modelling workflow to aid in VT therapeutics and has implications for generalizing personalized computational heart technology to a broad range of clinical centres.

Left bundle branch pacing guide by uninterrupted recording of intrinsic filtered and unfiltered intracardiac electrograms

BACKGROUND

Currently, the interrupted recording technique is commonly used to perform left bundle branch (LBB) pacing (LBBP). However, this method requires repeated testing to confirm that the LBB is captured and perforations are avoided. An automated solution may make this repetitive work easier.

CASE SUMMARY

LBBP was performed using an uninterrupted recording technique in an 86-year-old woman. Lead position and LBB capture was confirmed by the characteristics of the intrinsic filtered and unfiltered intracardiac electrograms.

CONCLUSION

Continuous mapping and recording technique may help achieve more accurate positioning of LBBP lead in the ventricular septum.

A universal biventricular coordinate system incorporating valve annuli: Validation in congenital heart disease

Universal coordinate systems have been proposed to facilitate anatomic registration between three-dimensional images, data and models of the ventricles of the heart. However, current universal ventricular coordinate systems do not account for the outflow tracts and valve annuli where the anatomy is complex. Here we propose an extension to the 'Cobiveco' biventricular coordinate system that also accounts for the intervalvular bridges of the base and provides a tool for anatomically consistent registration between widely varying biventricular shapes. CobivecoX uses a novel algorithm to separate intervalvular bridges and assign new coordinates, including an inflow-outflow coordinate, to describe local positions in these regions uniquely and consistently. Anatomic consistency of registration was validated using curated three-dimensional biventricular shape models derived from cardiac MRI measurements in normal hearts and hearts from patients with congenital heart diseases. This new method allows the advantages of universal cardiac coordinates to be used for three-dimensional ventricular imaging data and models that include the left and right ventricular outflow tracts and valve annuli.

A transferable in-silico augmented ischemic model for virtual myocardial perfusion imaging and myocardial infarction detection

This paper proposes an innovative approach to generate a generalized myocardial ischemia database by modeling the virtual electrophysiology of the heart and the 12-lead electrocardiography projected by the in-silico model can serve as a ready-to-use database for automatic myocardial infarction/ischemia (MI) localization and classification. Although the virtual heart can be created by an established technique combining the cell model with personalized heart geometry to observe the spatial propagation of depolarization and repolarization waves, we developed a strategy based on the clinical pathophysiology of MI to generate a heterogeneous database with a generic heart while maintaining clinical relevance and reduced computational complexity. First, the virtual heart is simplified into 11 regions that match the types and locations, which can be diagnosed by 12-lead ECG; the major arteries were divided into 3-5 segments from the upstream to the downstream based on the general anatomy. Second, the stenosis or infarction of the major or minor coronary artery branches can cause different perfusion drops and infarct sizes. We simulated the ischemic sites in different branches of the arteries by meandering the infarction location to elaborate on possible ECG representations, which alters the infraction's size and changes the transmembrane potential (TMP) of the myocytes associated with different levels of perfusion drop. A total of 8190 different case combinations of cardiac potentials with ischemia and MI were simulated, and the corresponding ECGs were generated by forward calculations. Finally, we trained and validated our in-silico database with a sparse representation classification (SRC) and tested the transferability of the model on the real-world Physikalisch Technische Bundesanstalt (PTB) database. The overall accuracies for localizing the MI region on the PTB data achieved 0.86, which is only 2% drop compared to that derived from the simulated database (0.88). In summary, we have shown a proof-of-concept for transferring an in-silico model to real-world database to compensate for insufficient data.

Virtual clinical QT exposure-response studies - A translational computational approach

BACKGROUND AND PURPOSE

A recent paradigm shift in proarrhythmic risk assessment suggests that the integration of clinical, non-clinical, and computational evidence can be used to reach a comprehensive understanding of the proarrhythmic potential of drug candidates. While current computational methodologies focus on predicting the incidence of proarrhythmic events after drug administration, the objective of this study is to predict concentration-response relationships of QTc as a clinical endpoint.

EXPERIMENTAL APPROACH

Full heart computational models reproducing human cardiac populations were created to predict the concentration-response relationship of changes in the QT interval as recommended for clinical trials. The concentration-response relationship of the QT-interval prolongation obtained from the computational cardiac population was compared against the relationship from clinical trial data for a set of well-characterized compounds: moxifloxacin, dofetilide, verapamil, and ondansetron.

KEY RESULTS

Computationally derived concentration-response relationships of QT interval changes for three of the four drugs had slopes within the confidence interval of clinical trials (dofetilide, moxifloxacin and verapamil) when compared to placebo-corrected concentration-ΔQT and concentration-ΔQT regressions. Moxifloxacin showed a higher intercept, outside the confidence interval of the clinical data, demonstrating that in this example, the standard linear regression does not appropriately capture the concentration-response results at very low concentrations. The concentrations corresponding to a mean QTc prolongation of 10 ms were consistently lower in the computational model than in clinical data. The critical concentration varied within an approximate ratio of 0.5 (moxifloxacin and ondansetron) and 1 times (dofetilide, verapamil) the critical concentration observed in human clinical trials. Notably, no other in silico methodology can approximate the human critical concentration values for a QT interval prolongation of 10 ms.

CONCLUSION AND IMPLICATIONS

Computational concentration-response modelling of a virtual population of high-resolution, 3-dimensional cardiac models can provide comparable information to clinical data and could be used to complement pre-clinical and clinical safety packages. It provides access to an unlimited exposure range to support trial design and can improve the understanding of pre-clinical-clinical translation.

HPC Framework for Performing in Silico Trials Using a 3D Virtual Human Cardiac Population as Means to Assess Drug-Induced Arrhythmic Risk

Following the 3 R's principles of animal research-replacement, reduction, and refinement-a high-performance computational framework was produced to generate a platform to perform human cardiac in-silico clinical trials as means to assess the pro-arrhythmic risk after the administrations of one or combination of two potentially arrhythmic drugs. The drugs assessed in this study were hydroxychloroquine and azithromycin. The framework employs electrophysiology simulations on high-resolution three-dimensional, biventricular human heart anatomies including phenotypic variabilities, so as to determine if differential QT-prolongation responds to drugs as observed clinically. These simulations also reproduce sex-specific ionic channel characteristics. The derived changes in the pseudo-electrocardiograms, calcium concentrations, as well as activation patterns within 3D geometries were evaluated for signs of induced arrhythmia. The virtual subjects could be evaluated at two different cycle lengths: at a normal heart rate and at a heart rate associated with stress as means to analyze the proarrhythmic risks after the administrations of hydroxychloroquine and azithromycin. Additionally, a series of experiments performed on reanimated swine hearts utilizing Visible Heart® methodologies in a four-chamber working heart model were performed to verify the arrhythmic behaviors observed in the in silico trials.The obtained results indicated similar pro-arrhythmic risk assessments within the virtual population as compared to published clinical trials (21% clinical risk vs 21.8% in silico trial risk). Evidence of transmurally heterogeneous action potential prolongations after providing a large dose of hydroxychloroquine was found as the observed mechanisms for elicited arrhythmias, both in the in vitro and the in silico models. The proposed workflow for in silico clinical drug cardiotoxicity trials allows for reproducing the complex behavior of cardiac electrophysiology in a varied population, in a matter of a few days as compared to the months or years it requires for most in vivo human clinical trials. Importantly, our results provided evidence of the common phenotype variants that produce distinct drug-induced arrhythmogenic outcomes.

Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region

BACKGROUND

While the microstructure of the left ventricle (LV) has been largely described, only a few studies investigated the right ventricular insertion point (RVIP). It was accepted that the aggregate cardiomyocytes organization was much more complex due to the intersection of the ventricular cavities but a precise structural characterization in the human heart was lacking even if clinical phenotypes related to right ventricular wall stress or arrhythmia were observed in this region.

METHODS

MRI-derived anatomical imaging (150 µm3) and diffusion tensor imaging (600 µm3) were performed in large mammalian whole hearts (human: N = 5, sheep: N = 5). Fractional anisotropy, aggregate cardiomyocytes orientations and tractography were compared within both species. Aggregate cardiomyocytes orientation on one ex-vivo sheep whole heart was then computed using structure tensor imaging (STI) from 21 µm isotropic acquisition acquired with micro computed tomography (MicroCT) imaging. Macroscopic and histological examination were performed. Lastly, experimental cardiomyocytes orientation distribution was then compared to the usual rule-based model using electrophysiological (EP) modeling. Electrical activity was modeled with the monodomain formulation.

RESULTS

The RVIP at the level of the inferior ventricular septum presented a unique arrangement of aggregate cardiomyocytes. An abrupt, mid-myocardial change in cardiomyocytes orientation was observed, delimiting a triangle-shaped region, present in both sheep and human hearts. FA's histogram distribution (mean ± std: 0.29 ± 0.06) of the identified region as well as the main dimension (22.2 mm ± 5.6 mm) was found homogeneous across samples and species. Averaged volume is 0.34 cm3 ± 0.15 cm3. Both local activation time (LAT) and morphology of pseudo-ECGs were strongly impacted with delayed LAT and change in peak-to-peak amplitude in the simulated wedge model.

CONCLUSION

The study was the first to describe the 3D cardiomyocytes architecture of the basal inferoseptal left ventricle region in human hearts and identify the presence of a well-organized aggregate cardiomyocytes arrangement and cardiac structural discontinuities. The results might offer a better appreciation of clinical phenotypes like RVIP-late gadolinium enhancement or uncommon idiopathic ventricular arrhythmias (VA) originating from this region.

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

A comprehensive mathematical model for cardiac perfusion

The aim of this paper is to introduce a new mathematical model that simulates myocardial blood perfusion that accounts for multiscale and multiphysics features. Our model incorporates cardiac electrophysiology, active and passive mechanics, hemodynamics, valve modeling, and a multicompartment Darcy model of perfusion. We consider a fully coupled electromechanical model of the left heart that provides input for a fully coupled Navier-Stokes-Darcy Model for myocardial perfusion. The fluid dynamics problem is modeled in a left heart geometry that includes large epicardial coronaries, while the multicompartment Darcy model is set in a biventricular myocardium. Using a realistic and detailed cardiac geometry, our simulations demonstrate the biophysical fidelity of our model in describing cardiac perfusion. Specifically, we successfully validate the model reliability by comparing in-silico coronary flow rates and average myocardial blood flow with clinically established values ranges reported in relevant literature. Additionally, we investigate the impact of a regurgitant aortic valve on myocardial perfusion, and our results indicate a reduction in myocardial perfusion due to blood flow taken away by the left ventricle during diastole. To the best of our knowledge, our work represents the first instance where electromechanics, hemodynamics, and perfusion are integrated into a single computational framework.

Transmural fibre orientations based on Laplace-Dirichlet-Rule-Based-Methods and their influence on human heart simulations

It is well known that the orthotropic tissue structure decisively influences the mechanical and electrical properties of the heart. Numerous approaches to compute the orthotropic tissue structure in computational heart models have been developed in the past decades. In this study, we investigate to what extent different Laplace-Dirichlet-Rule-Based-Methods (LDRBMs) influence the local orthotropic tissue structure and thus the electromechanical behaviour of the subsequent cardiac simulation. In detail, we are utilising three Laplace-Dirichlet-Rule-Based-Methods and compare: (i) the local myofibre orientation; (ii) important global characteristics (ejection fraction, peak pressure, apex shortening, myocardial volume reduction, fractional wall thickening); (iii) local characteristics (active fibre stress, fibre strain). We observe that the orthotropic tissue structures for the three LDRBMs show significant differences in the local myofibre orientation. The global characteristics myocardial volume reduction and peak pressure are rather insensitive to a change in local myofibre orientation, while the ejection fraction is moderately influenced by the different LDRBMs. Moreover, the apical shortening and fractional wall thickening exhibit a sensitive behaviour to a change in the local myofibre orientation. The highest sensitivity can be observed for the local characteristics.

Assessment of Risk for Ventricular Tachycardia based on Extensive Electrophysiology Simulations

Patients that have suffered a myocardial infarction are at high risk of developing ventricular tachycardia. Patient stratification is often determined by characterization of the underlying myocardial substrate by cardiac imaging methods. In this study, we show that computer modeling of cardiac electrophysiology based on personalized fast 3D simulations can help to assess patient risk to arrhythmia. We perform a large simulation study on 21 patient digital twins and reproduce successfully the clinical outcomes. In addition, we provide the sites which are prone to sustain ventricular tachycardias, i.e, onset sites around the scar region, and validate if they colocalize with exit sites from slow conduction channels.Clinical relevance- Fast electrophysiological simulations can provide advanced patient stratification indices and predict arrhythmic susceptibility to suffer from ventricular tachycardia in patients that have suffered a myocardial infarction.

GPU accelerated digital twins of the human heart open new routes for cardiovascular research

The recruitment of patients for rare or complex cardiovascular diseases is a bottleneck for clinical trials and digital twins of the human heart have recently been proposed as a viable alternative. In this paper we present an unprecedented cardiovascular computer model which, relying on the latest GPU-acceleration technologies, replicates the full multi-physics dynamics of the human heart within a few hours per heartbeat. This opens the way to extensive simulation campaigns to study the response of synthetic cohorts of patients to cardiovascular disorders, novel prosthetic devices or surgical procedures. As a proof-of-concept we show the results obtained for left bundle branch block disorder and the subsequent cardiac resynchronization obtained by pacemaker implantation. The in-silico results closely match those obtained in clinical practice, confirming the reliability of the method. This innovative approach makes possible a systematic use of digital twins in cardiovascular research, thus reducing the need of real patients with their economical and ethical implications. This study is a major step towards in-silico clinical trials in the era of digital medicine.

Fiber Organization Has Little Effect on Electrical Activation Patterns During Focal Arrhythmias in the Left Atrium

Over the past two decades there has been a steady trend towards the development of realistic models of cardiac conduction with increasing levels of detail. However, making models more realistic complicates their personalization and use in clinical practice due to limited availability of tissue and cellular scale data. One such limitation is obtaining information about myocardial fiber organization in the clinical setting. In this study, we investigated a chimeric model of the left atrium utilizing clinically derived patient-specific atrial geometry and a realistic, yet foreign for a given patient fiber organization. We discovered that even significant variability of fiber organization had a relatively small effect on the spatio-temporal activation pattern during regular pacing. For a given pacing site, the activation maps were very similar across all fiber organizations tested.

lifex-fiber: an open tool for myofibers generation in cardiac computational models

BACKGROUND

Modeling the whole cardiac function involves the solution of several complex multi-physics and multi-scale models that are highly computationally demanding, which call for simpler yet accurate, high-performance computational tools. Despite the efforts made by several research groups, no software for whole-heart fully-coupled cardiac simulations in the scientific community has reached full maturity yet.

RESULTS

In this work we present [Formula: see text]-fiber, an innovative tool for the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, which are the essential building blocks for modeling the electrophysiological, mechanical and electromechanical cardiac function, from single-chamber to whole-heart simulations. [Formula: see text]-fiber is the first publicly released module for cardiac simulations based on [Formula: see text], an open-source, high-performance Finite Element solver for multi-physics, multi-scale and multi-domain problems developed in the framework of the iHEART project, which aims at making in silico experiments easily reproducible and accessible to a wide community of users, including those with a background in medicine or bio-engineering.

CONCLUSIONS

The tool presented in this document is intended to provide the scientific community with a computational tool that incorporates general state of the art models and solvers for simulating the cardiac function within a high-performance framework that exposes a user- and developer-friendly interface. This report comes with an extensive technical and mathematical documentation to welcome new users to the core structure of [Formula: see text]-fiber and to provide them with a possible approach to include the generated cardiac fibers into more sophisticated computational pipelines. In the near future, more modules will be successively published either as pre-compiled binaries for x86-64 Linux systems or as open source software.

A mathematical model that integrates cardiac electrophysiology, mechanics, and fluid dynamics: Application to the human left heart

We propose a mathematical and numerical model for the simulation of the heart function that couples cardiac electrophysiology, active and passive mechanics and hemodynamics, and includes reduced models for cardiac valves and the circulatory system. Our model accounts for the major feedback effects among the different processes that characterize the heart function, including electro-mechanical and mechano-electrical feedback as well as force-strain and force-velocity relationships. Moreover, it provides a three-dimensional representation of both the cardiac muscle and the hemodynamics, coupled in a fluid-structure interaction (FSI) model. By leveraging the multiphysics nature of the problem, we discretize it in time with a segregated electrophysiology-force generation-FSI approach, allowing for efficiency and flexibility in the numerical solution. We employ a monolithic approach for the numerical discretization of the FSI problem. We use finite elements for the spatial discretization of partial differential equations. We carry out a numerical simulation on a realistic human left heart model, obtaining results that are qualitatively and quantitatively in agreement with physiological ranges and medical images.

The role of pulse timing in cardiac defibrillation

Life-threatening cardiac arrhythmias require immediate defibrillation. For state-of-the-art shock treatments, a high field strength is required to achieve a sufficient success rate for terminating the complex spiral wave (rotor) dynamics underlying cardiac fibrillation. However, such high energy shocks have many adverse side effects due to the large electric currents applied. In this study, we show, using 2D simulations based on the Fenton-Karma model, that also pulses of relatively low energy may terminate the chaotic activity if applied at the right moment in time. In our simplified model for defibrillation, complex spiral waves are terminated by local perturbations corresponding to conductance heterogeneities acting as virtual electrodes in the presence of an external electric field. We demonstrate that time series of the success rate for low energy shocks exhibit pronounced peaks which correspond to short intervals in time during which perturbations aiming at terminating the chaotic fibrillation state are (much) more successful. Thus, the low energy shock regime, although yielding very low temporal average success rates, exhibits moments in time for which success rates are significantly higher than the average value shown in dose-response curves. This feature might be exploited in future defibrillation protocols for achieving high termination success rates with low or medium pulse energies.

Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O'Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studies.

A meshless fragile points method for rule-based definition of myocardial fiber orientation

BACKGROUND AND OBJECTIVE

Rule-based methods are commonly used to estimate the arrangement of myocardial fibers by solving the Laplace problem with appropriate Dirichlet boundary conditions. Existing algorithms are using the Finite Element Method (FEM) to solve the Laplace-Dirichlet problem. However, meshless methods are under development for cardiac electrophysiology simulation. The objective of this work is to propose a meshless rule based method for the determination of myocardial fiber arrangement without requiring a mesh discretization as it is required by FEM.

METHODS

The proposed method employs the Fragile Points Method (FPM) for the solution of the Laplace-Dirichlet problem. FPM uses simple discontinuous trial functions and single-point exact integration for linear trial functions that set it as a promising alternative to the Finite Element Method. We derive the FPM formulation of the Laplace-Dirichlet and we estimate ventricular and atrial fiber arrangements according to rules based on histology findings for four different geometries. The obtained fiber arrangements from FPM are compared with the ones obtained from FEM by calculating the angle between the fiber vector fields of the two methods for three different directions (i.e., longitudinal, sheet, transverse).

RESULTS

The fiber arrangements that were generated with FPM were in close agreement with the generated arrangements from FEM for all three directions. The mean angle difference between the FPM and FEM vector fields were lower than 0.030^∘ for the ventricular fiber arrangements and lower than 0.036^∘ for the atrial fiber arrangements.

DISCUSSION

The proposed meshless rule-based method was proven to generate myocardial fiber arrangements with very close agreement with FEM while alleviates the requirement for a mesh of the latter. This is of great value for cardiac electrophysiology solvers that are based on meshless methods since they require a well defined myocardial fiber arrangement to simulate accurately the propagation of electrical signals in the heart. Combining a meshless solution for both the determination of the fibers and the electrical signal propagation can allow for solution that do not require the definition of a mesh. To our knowledge, this work is the first one to propose a meshless rule-based method for myocardial fiber arrangement determination.

An Atlas-Based Analysis of Biventricular Mechanics in Tetralogy of Fallot

The current study proposes an efficient strategy for exploiting the statistical power of cardiac atlases to investigate whether clinically significant variations in ventricular shape are sufficient to explain corresponding differences in ventricular wall motion directly, or if they are indirect markers of altered myocardial mechanical properties. This study was conducted in a cohort of patients with repaired tetralogy of Fallot (rTOF) that face long-term right ventricular (RV) and/or left ventricular (LV) dysfunction as a consequence of adverse remodeling. Features of biventricular end-diastolic (ED) shape associated with RV apical dilation, LV dilation, RV basal bulging, and LV conicity correlated with components of systolic wall motion (SWM) that contribute most to differences in global systolic function. A finite element analysis of systolic biventricular mechanics was employed to assess the effect of perturbations in these ED shape modes on corresponding components of SWM. Perturbations to ED shape modes and myocardial contractility explained observed variation in SWM to varying degrees. In some cases, shape markers were partial determinants of systolic function and, in other cases, they were indirect markers for altered myocardial mechanical properties. Patients with rTOF may benefit from an atlas-based analysis of biventricular mechanics to improve prognosis and gain mechanistic insight into underlying myocardial pathophysiology.

Training machine learning models with synthetic data improves the prediction of ventricular origin in outflow tract ventricular arrhythmias

In order to determine the site of origin (SOO) in outflow tract ventricular arrhythmias (OTVAs) before an ablation procedure, several algorithms based on manual identification of electrocardiogram (ECG) features, have been developed. However, the reported accuracy decreases when tested with different datasets. Machine learning algorithms can automatize the process and improve generalization, but their performance is hampered by the lack of large enough OTVA databases. We propose the use of detailed electrophysiological simulations of OTVAs to train a machine learning classification model to predict the ventricular origin of the SOO of ectopic beats. We generated a synthetic database of 12-lead ECGs (2,496 signals) by running multiple simulations from the most typical OTVA SOO in 16 patient-specific geometries. Two types of input data were considered in the classification, raw and feature ECG signals. From the simulated raw 12-lead ECG, we analyzed the contribution of each lead in the predictions, keeping the best ones for the training process. For feature-based analysis, we used entropy-based methods to rank the obtained features. A cross-validation process was included to evaluate the machine learning model. Following, two clinical OTVA databases from different hospitals, including ECGs from 365 patients, were used as test-sets to assess the generalization of the proposed approach. The results show that V2 was the best lead for classification. Prediction of the SOO in OTVA, using both raw signals or features for classification, presented high accuracy values (>0.96). Generalization of the network trained on simulated data was good for both patient datasets (accuracy of 0.86 and 0.84, respectively) and presented better values than using exclusively real ECGs for classification (accuracy of 0.84 and 0.76 for each dataset). The use of simulated ECG data for training machine learning-based classification algorithms is critical to obtain good SOO predictions in OTVA compared to real data alone. The fast implementation and generalization of the proposed methodology may contribute towards its application to a clinical routine.

Meshless Electrophysiological Modeling of Cardiac Resynchronization Therapy—Benchmark Analysis with Finite-Element Methods in Experimental Data

Computational models of cardiac electrophysiology are promising tools for reducing the rates of non-response patients suitable for cardiac resynchronization therapy (CRT) by optimizing electrode placement. The majority of computational models in the literature are mesh-based, primarily using the finite element method (FEM). The generation of patient-specific cardiac meshes has traditionally been a tedious task requiring manual intervention and hindering the modeling of a large number of cases. Meshless models can be a valid alternative due to their mesh quality independence. The organization of challenges such as the CRT-EPiggy19, providing unique experimental data as open access, enables benchmarking analysis of different cardiac computational modeling solutions with quantitative metrics. We present a benchmark analysis of a meshless-based method with finite-element methods for the prediction of cardiac electrical patterns in CRT, based on a subset of the CRT-EPiggy19 dataset. A data assimilation strategy was designed to personalize the most relevant parameters of the electrophysiological simulations and identify the optimal CRT lead configuration. The simulation results obtained with the meshless model were equivalent to FEM, with the most relevant aspect for accurate CRT predictions being the parameter personalization strategy (e.g., regional conduction velocity distribution, including the Purkinje system and CRT lead distribution).

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

A Transmural Path Model Improves the Definition of the Orthotropic Tissue Structure in Heart Simulations

In the past decades, the structure of the heart, human as well as other species, has been explored in a detailed way, e.g., via histological studies or diffusion tensor magnetic resonance imaging. Nevertheless, the assignment of the characteristic orthotropic structure in a patient-specific finite element model remains a challenging task. Various types of rule-based models, which define the local fiber and sheet orientation depending on the transmural depth, have been developed. However, the correct assessment of the transmural depth is not trivial. Its accuracy has a substantial influence on the overall mechanical and electrical properties in rule-based models. The main purpose of this study is the development of a finite element-based approach to accurately determine the transmural depth on a general unstructured grid. Instead of directly using the solution of the Laplace problem as the transmural depth, we make use of a well-established model for the assessment of the transmural thickness. It is based on two hyperbolic first-order partial differential equations for the definition of a transmural path, whereby the transmural thickness is defined as the arc length of this path. Subsequently, the transmural depth is determined based on the position on the transmural path. Originally, the partial differential equations were solved via finite differences on structured grids. In order to circumvent the need of two grids and mapping between the structured (to determine the transmural depth) and unstructured (electromechanical heart simulation) grids, we solve the equations directly on the same unstructured tetrahedral mesh. We propose a finite-element-based discontinuous Galerkin approach. Based on the accurate transmural depth, we assign the local material orientation of the orthotropic tissue structure in a usual fashion. We show that this approach leads to a more accurate definition of the transmural depth. Furthermore, for the left ventricle, we propose functions for the transmural fiber and sheet orientation by fitting them to literature-based diffusion tensor magnetic resonance imaging data. The proposed functions provide a distinct improvement compared to existing rules from the literature.

Effect of anisotropy in myocardial electrical conductivity on lesion characteristics during radiofrequency cardiac ablation: a numerical study

BACKGROUND

Traditional computer simulation studies of radiofrequency catheter ablation (RFCA) usually neglect the anisotropy in myocardial electrical conductivity (MEC), which is likely an essential factor in governing the ablation outcome. Here, a numerical study of lesion characteristics during RFCA based on an anatomy-based model incorporating fiber orientation was performed to investigate the anisotropy in MEC.

METHODS

A three-dimensional thorax model including atria, blood, connective tissue, muscle, fat, and skin was constructed. The myocardial fiber was established through a rule-based method (RBM) based on the anatomical structure of the heart. The anisotropic MEC were 0.40 and 0.28 S m^-1 in longitudinal and transverse directions, respectively. The ablation result was compared with the isotropic scenario where the isotropic MEC was the average of the anisotropic conductivities as 0.34 S m^-1.

RESULTS

The complexity of fiber architecture varied with that of the local anatomical structure. At RF power of 20 W for 30 s, the tissue temperature and lesion volume were reduced by 2.8 ± 0.1% and 6.9 ± 0.5%, respectively, under anisotropic MEC around the ostium of the pulmonary vein and left atrial appendage. Those for the posterior wall and roof of the left atrium, and the inside of the superior vena cava were 1.9 ± 0.3% and 5.6 ± 1.2%, respectively.

CONCLUSIONS

Anisotropy in MEC has a greater reduction effect on lesion volume than on tissue temperature during RFCA; this effect tends to be restrained at positions with more uniform fiber distributions and can be enhanced where significant variation in fiber architecture occurred.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.

Analyzing the Role of Repolarization Gradients in Post-infarct Ventricular Tachycardia Dynamics Using Patient-Specific Computational Heart Models

Aims: Disease-induced repolarization heterogeneity in infarcted myocardium contributes to VT arrhythmogenesis but how apicobasal and transmural (AB-TM) repolarization gradients additionally affect post-infarct VT dynamics is unknown. The goal of this study is to assess how AB-TM repolarization gradients impact post-infarct VT dynamics using patient-specific heart models.

Method: 3D late gadolinium-enhanced cardiac magnetic resonance images were acquired from seven post-infarct patients. Models representing the patient-specific scar and infarct border zone distributions were reconstructed without (baseline) and with repolarization gradients along both the AB-TM axes. AB only and TM only models were created to assess the effects of each ventricular gradient on VT dynamics. VTs were induced in all models via rapid pacing.

Results: Ten baseline VTs were induced. VT inducibility in AB-TM models was not significantly different from baseline (p>0.05). Reentry pathways in AB-TM models were different than baseline pathways due to alterations in the location of conduction block (p<0.05). VT exit sites in AB-TM models were different than baseline VT exit sites (p<0.05). VT inducibility of AB only and TM only models were not significantly different than that of baseline (p>0.05) or AB-TM models (p>0.05). Reentry pathways and VT exit sites in AB only and TM only models were different than in baseline (p<0.05). Lastly, repolarization gradients uncovered multiple VT morphologies with different reentrant pathways and exit sites within the same structural, conducting channels.

Conclusion: VT inducibility was not impacted by the addition of AB-TM repolarization gradients, but the VT reentrant pathway and exit sites were greatly affected due to modulation of conduction block. Thus, during ablation procedures, physiological and pharmacological factors that impact the ventricular repolarization gradient might need to be considered when targeting the VTs.

An Implementation of Patient-Specific Biventricular Mechanics Simulations With a Deep Learning and Computational Pipeline

Parameterised patient-specific models of the heart enable quantitative analysis of cardiac function as well as estimation of regional stress and intrinsic tissue stiffness. However, the development of personalised models and subsequent simulations have often required lengthy manual setup, from image labelling through to generating the finite element model and assigning boundary conditions. Recently, rapid patient-specific finite element modelling has been made possible through the use of machine learning techniques. In this paper, utilising multiple neural networks for image labelling and detection of valve landmarks, together with streamlined data integration, a pipeline for generating patient-specific biventricular models is applied to clinically-acquired data from a diverse cohort of individuals, including hypertrophic and dilated cardiomyopathy patients and healthy volunteers. Valve motion from tracked landmarks as well as cavity volumes measured from labelled images are used to drive realistic motion and estimate passive tissue stiffness values. The neural networks are shown to accurately label cardiac regions and features for these diverse morphologies. Furthermore, differences in global intrinsic parameters, such as tissue anisotropy and normalised active tension, between groups illustrate respective underlying changes in tissue composition and/or structure as a result of pathology. This study shows the successful application of a generic pipeline for biventricular modelling, incorporating artificial intelligence solutions, within a diverse cohort.

A finite element model of the cardiac ventricles with coupled circulation: Biventricular mesh generation with hexahedral elements, airbags and a functional mockup interface to the circulation

INTRODUCTION

Finite element (FE) mechanics models of the heart are becoming more sophisticated. However, there is lack of consensus about optimal element type and coupling of FE models to the circulation. We describe biventricular (left (LV) and right (RV) ventricles) FE mechanics model creation using hexahedral elements, airbags and a functional mockup interface (FMI) to lumped-parameter models of the circulation.

METHODS

Cardiac MRI (CMR) was performed in two healthy volunteers and a single patient with ischemic heart disease (IHD). CMR images were segmented and surfaced, meshing with hexahedral elements was performed with a "thin butterfly with septum" topology. LV and RV inflow and outflow airbags were coupled to lumped-parameter circulation models with an FMI interface. Pulmonary constriction (PAC) and vena cava occlusion (VCO) were simulated and end-systolic pressure-volume relations (ESPVR) were calculated.

RESULTS

Mesh construction was prompt with representative contouring and mesh adjustment requiring 32 and 26 min Respectively. The numbers of elements ranged from 4104 to 5184 with a representative Jacobian of 1.0026 ± 0.4531. Agreement between CMR-based surfaces and mesh was excellent with root-mean-squared error of 0.589 ± 0.321 mm. The LV ESPVR slope was 3.37 ± 0.09 in volunteers but 2.74 in the IHD patient. The effect of PAC and VCO on LV ESPVR was consistent with ventricular interaction (p = 0.0286).

CONCLUSION

Successful co-simulation using a biventricular FE mechanics model with hexahedral elements, airbags and an FMI interface to lumped-parameter model of the circulation was demonstrated. Future studies will include comparison of element type and study of cardiovascular pathologies and device therapies.

Estimation of Personalized Minimal Purkinje Systems From Human Electro-Anatomical Maps

The Purkinje system is a heart structure responsible for transmitting electrical impulses through the ventricles in a fast and coordinated way to trigger mechanical contraction. Estimating a patient-specific compatible Purkinje Network from an electro-anatomical map is a challenging task, that could help to improve models for electrophysiology simulations or provide aid in therapy planning, such as radiofrequency ablation. In this study, we present a methodology to inversely estimate a Purkinje network from a patient's electro-anatomical map. First, we carry out a simulation study to assess the accuracy of the method for different synthetic Purkinje network morphologies and myocardial junction densities. Second, we estimate the Purkinje network from a set of 28 electro-anatomical maps from patients, obtaining an optimal conduction velocity in the Purkinje network of 1.95 ± 0.25 m/s, together with the location of their Purkinje-myocardial junctions, and Purkinje network structure. Our results showed an average local activation time error of 6.8±2.2 ms in the endocardium. Finally, using the personalized Purkinje network, we obtained correlations higher than 0.85 between simulated and clinical 12-lead ECGs.

A dual adaptive explicit time integration algorithm for efficiently solving the cardiac monodomain equation

The monodomain model is widely used in in-silico cardiology to describe excitation propagation in the myocardium. Frequently, operator splitting is used to decouple the stiff reaction term and the diffusion term in the monodomain model so that they can be solved separately. Commonly, the diffusion term is solved implicitly with a large time step while the reaction term is solved by using an explicit method with adaptive time stepping. In this work, we propose a fully explicit method for the solution of the decoupled monodomain model. In contrast to semi-implicit methods, fully explicit methods present lower memory footprint and higher scalability. However, such methods are only conditionally stable. We overcome the conditional stability limitation by proposing a dual adaptive explicit method in which adaptive time integration is applied for the solution of both the reaction and diffusion terms. We perform a set of numerical examples where cardiac propagation is simulated under physiological and pathophysiological conditions. Results show that the proposed method presents preserved accuracy and improved computational efficiency as compared to standard operator splitting-based methods.

Multiphysics computational modelling of the cardiac ventricles

Development of cardiac multiphysics models has progressed significantly over the decades and simulations combining multiple physics interactions have become increasingly common. In this review, we summarise the progress in this field focusing on various approaches of integrating ventricular structures. electrophysiological properties, myocardial mechanics, as well as incorporating blood hemodynamics and the circulatory system. Common coupling approaches are discussed and compared, including the advantages and shortcomings of each. Currently used strategies for patient-specific implementations are highlighted and potential future improvements considered.

Inference of ventricular activation properties from non-invasive electrocardiography

The realisation of precision cardiology requires novel techniques for the non-invasive characterisation of individual patients’ cardiac function to inform therapeutic and diagnostic decision-making. Both electrocardiography and imaging are used for the clinical diagnosis of cardiac disease. The integration of multi-modal datasets through advanced computational methods could enable the development of the cardiac ‘digital twin’, a comprehensive virtual tool that mechanistically reveals a patient's heart condition from clinical data and simulates treatment outcomes. The adoption of cardiac digital twins requires the non-invasive efficient personalisation of the electrophysiological properties in cardiac models. This study develops new computational techniques to estimate key ventricular activation properties for individual subjects by exploiting the synergy between non-invasive electrocardiography, cardiac magnetic resonance (CMR) imaging and modelling and simulation. More precisely, we present an efficient sequential Monte Carlo approximate Bayesian computation-based inference method, integrated with Eikonal simulations and torso-biventricular models constructed based on clinical CMR imaging. The method also includes a novel strategy to treat combined continuous (conduction speeds) and discrete (earliest activation sites) parameter spaces and an efficient dynamic time warping-based ECG comparison algorithm. We demonstrate results from our inference method on a cohort of twenty virtual subjects with cardiac ventricular myocardial-mass volumes ranging from 74 cm³ to 171 cm³ and considering low versus high resolution for the endocardial discretisation (which determines possible locations of the earliest activation sites). Results show that our method can successfully infer the ventricular activation properties in sinus rhythm from non-invasive epicardial activation time maps and ECG recordings, achieving higher accuracy for the endocardial speed and sheet (transmural) speed than for the fibre or sheet-normal directed speeds.

•

Estimation of the ventricular speeds and earliest activation sites from ECG and CMR.

•

Evaluation with twenty virtual subjects shows the effect of anatomical variability.

•

Bayesian-inspired simultaneous estimation of continuous and discrete parameters.

•

Efficient dynamic time warping-based comparison of electrocardiograms (ECG).

•

Changing fibre and sheet-normal speed does not affect healthy activation sequence.

Electro-Mechanical Whole-Heart Digital Twins: A Fully Coupled Multi-Physics Approach

Mathematical models of the human heart are evolving to become a cornerstone of precision medicine and support clinical decision making by providing a powerful tool to understand the mechanisms underlying pathophysiological conditions. In this study, we present a detailed mathematical description of a fully coupled multi-scale model of the human heart, including electrophysiology, mechanics, and a closed-loop model of circulation. State-of-the-art models based on human physiology are used to describe membrane kinetics, excitation-contraction coupling and active tension generation in the atria and the ventricles. Furthermore, we highlight ways to adapt this framework to patient specific measurements to build digital twins. The validity of the model is demonstrated through simulations on a personalized whole heart geometry based on magnetic resonance imaging data of a healthy volunteer. Additionally, the fully coupled model was employed to evaluate the effects of a typical atrial ablation scar on the cardiovascular system. With this work, we provide an adaptable multi-scale model that allows a comprehensive personalization from ion channels to the organ level enabling digital twin modeling.

A 3D personalized cardiac myocyte aggregate orientation model using MRI data-driven low-rank basis functions

Cardiac myocyte aggregate orientation has a strong impact on cardiac electrophysiology and mechanics. Studying the link between structural characteristics, strain, and stresses over the cardiac cycle and cardiac function requires a full volumetric representation of the microstructure. In this work, we exploit the structural similarity across hearts to extract a low-rank representation of predominant myocyte orientation in the left ventricle from high-resolution magnetic resonance ex-vivo cardiac diffusion tensor imaging (cDTI) in porcine hearts. We compared two reduction methods, Proper Generalized Decomposition combined with Singular Value Decomposition and Proper Orthogonal Decomposition. We demonstrate the existence of a general set of basis functions of aggregated myocyte orientation which defines a data-driven, personalizable, parametric model featuring higher flexibility than existing atlas and rule-based approaches. A more detailed representation of microstructure matching the available patient data can improve the accuracy of personalized computational models. Additionally, we approximate the myocyte orientation of one ex-vivo human heart and demonstrate the feasibility of transferring the basis functions to humans.

Volumetric parcellation of the cardiac right ventricle for regional geometric and functional assessment

3D echocardiography is an increasingly popular tool for assessing cardiac remodelling in the right ventricle (RV). It allows quantification of the cardiac chambers without any geometric assumptions, which is the main weakness of 2D echocardiography. However, regional quantification of geometry and function is limited by the lower spatial and temporal resolution and the scarcity of identifiable anatomical landmarks, especially within the ventricular cavity. We developed a technique for regionally assessing the volume of 3 relevant RV volumetric regions: apical, inlet and outflow. The proposed parcellation method is based on the geodesic distances to anatomical landmarks that are easily identifiable in the images: the apex and the tricuspid and pulmonary valves, each associated to a region. Based on these distances, we define a partition in the endocardium at end-diastole (ED). This partition is then interpolated to the blood cavity using the Laplace equation, which allows to compute regional volumes. For obtaining an end-systole (ES) partition, the endocardial partition is transported from ED to ES using a commercial image-based tracking software, and then the interpolation process is repeated. We assessed the intra- and inter-observer reproducibility using a 10-subjects dataset containing repeated quantifications of the same images, obtaining intra- and inter- observer errors (7-12% and 10-23% respectively). Finally, we propose a novel synthetic mesh generation algorithm that deforms a template mesh imposing a user-defined strain to a template mesh. We used this method to create a new dataset for involving distinct types of remodelling that were used to assess the sensitivity of the parcellation method to identify volume changes affecting different parts. We show that the parcellation method is adequate for capturing local circumferential and global circumferential and longitudinal RV remodelling, which are the most clinically relevant cases.

A microdevice platform for characterizing the effect of mechanical strain magnitudes on the maturation of iPSC-Cardiomyocytes

The use of human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) as an in vitro model of the heart is limited by their structurally and functionally immature phenotypes. During heart development, mechanical stimuli from in vivo microenvironments are known to regulate cardiomyocyte gene expression and maturation. Accordingly, protocols for culturing iPSC-CMs have recently incorporated mechanical or electromechanical stimulation to induce cellular maturation in vitro; however, the response of iPSC-CMs to different mechanical strain magnitudes is unknown, and existing techniques lack the capability to dynamically measure changes to iPSC-CM contractility in situ as maturation progresses. We developed a microdevice platform which applies cyclical strains of varying magnitudes (5%, 10%, 15% and 20%) to a monolayer of iPSC-CMs, coincidentally measuring contractile stress during mechanical stimulation using fluorescent nanobeads embedded in the microdevice's suspended membrane. Cyclic strain was found to induce circumferential cell alignment on the actuated membranes. In situ contractility measurements revealed that cyclic stimulation gradually increased cardiomyocyte contractility during a 10-day culture period. The contractile stress of iPSC-CM monolayers was found to increase with a higher strain magnitude and plateaued at 15% strain. Cardiomyocyte contractility positively correlated with the elongation of sarcomeres and an increased expression of β-myosin heavy chain (MYH7) in a strain magnitude-dependent manner, illustrating how mechanical stress can be optimized for the phenotypic and proteomic maturation of the cells. iPSC-CMs with improved maturity have the potential to create a more accurate heart model in vitro for applications in disease modeling and therapeutic discovery.

Integration of activation maps of epicardial veins in computational cardiac electrophysiology

In this work we address the issue of validating the monodomain equation used in combination with the Bueno-Orovio ionic model for the prediction of the activation times in cardiac electro-physiology of the left ventricle. To this aim, we consider four patients who suffered from Left Bundle Branch Block (LBBB). We use activation maps performed at the septum as input data for the model and maps at the epicardial veins for the validation. In particular, a first set (half) of the latter are used to estimate the conductivities of the patient and a second set (the remaining half) to compute the errors of the numerical simulations. We find an excellent agreement between measures and numerical results. Our validated computational tool could be used to accurately predict activation times at the epicardial veins with a short mapping, i.e. by using only a part (the most proximal) of the standard acquisition points, thus reducing the invasive procedure and exposure to radiation.

In silico pace-mapping: prediction of left vs. right outflow tract origin in idiopathic ventricular arrhythmias with patient-specific electrophysiological simulations

AIMS

A pre-operative non-invasive identification of the site of origin (SOO) of outflow tract ventricular arrhythmias (OTVAs) is important to properly plan radiofrequency ablation procedures. Although some algorithms based on electrocardiograms (ECGs) have been developed to predict left vs. right ventricular origins, their accuracy is still limited, especially in complex anatomies. The aim of this work is to use patient-specific electrophysiological simulations of the heart to predict the SOO in OTVA patients.

METHODS AND RESULTS

An in silico pace-mapping procedure was designed and used on 11 heart geometries, generating for each case simulated ECGs from 12 clinically plausible SOO. Subsequently, the simulated ECGs were compared with patient ECG data obtained during the clinical tachycardia using the 12-lead correlation coefficient (12-lead ρ). Left ventricle (LV) vs. right ventricle (RV) SOO was estimated by computing the LV/RV ratio for each patient, obtained by dividing the average 12-lead ρ value of the LV- and RV-SOO simulated ECGs, respectively. Simulated ECGs that had virtual sites close to the ablation points that stopped the arrhythmia presented higher correlation coefficients. The LV/RV ratio correctly predicted LV vs. RV SOO in 10/11 cases; 1.07 vs. 0.93 P < 0.05 for 12-lead ρ.

CONCLUSION

The obtained results demonstrate the potential of the developed in silico pace-mapping technique to complement standard ECG for the pre-operative planning of complex ventricular arrhythmias.