Computational generation of the Purkinje network driven by clinical measurements: the case of pathological propagations

Enhanced optimization-based method for the generation of patient-specific models of Purkinje networks

Cardiac Purkinje networks are a fundamental part of the conduction system and are known to initiate a variety of cardiac arrhythmias. However, patient-specific modeling of Purkinje networks remains a challenge due to their high morphological complexity. This work presents a novel method based on optimization principles for the generation of Purkinje networks that combines geometric and activation accuracy in branch size, bifurcation angles, and Purkinje-ventricular-junction activation times. Three biventricular meshes with increasing levels of complexity are used to evaluate the performance of our approach. Purkinje-tissue coupled monodomain simulations are executed to evaluate the generated networks in a realistic scenario using the most recent Purkinje/ventricular human cellular models and physiological values for the Purkinje-ventricular-junction characteristic delay. The results demonstrate that the new method can generate patient-specific Purkinje networks with controlled morphological metrics and specified local activation times at the Purkinje-ventricular junctions.

Computational electrophysiology of the coronary sinus branches based on electro-anatomical mapping for the prediction of the latest activated region

This work dealt with the assessment of a computational tool to estimate the electrical activation in the left ventricle focusing on the latest electrically activated segment (LEAS) in patients with left bundle branch block and possible myocardial fibrosis. We considered the Eikonal-diffusion equation and to recover the electrical activation maps in the myocardium. The model was calibrated by using activation times acquired in the coronary sinus (CS) branches or in the CS solely with an electroanatomic mapping system (EAMS) during cardiac resynchronization therapy (CRT). We applied our computational tool to ten patients founding an excellent accordance with EAMS measures; in particular, the error for LEAS location was less than 4 mm. We also calibrated our model using only information in the CS, still obtaining an excellent agreement with the measured LEAS. The proposed tool was able to accurately reproduce the electrical activation maps and in particular LEAS location in the CS branches, with an almost real-time computational effort, regardless of the presence of myocardial fibrosis, even when information only at CS was used to calibrate the model. This could be useful in the clinical practice since LEAS is often used as a target site for the left lead placement during CRT.

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.

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 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.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

An Inverse Eikonal Method for Identifying Ventricular Activation Sequences from Epicardial Activation Maps

A key mechanism controlling cardiac function is the electrical activation sequence of the heart's main pumping chambers termed the ventricles. As such, personalization of the ventricular activation sequences is of pivotal importance for the clinical utility of computational models of cardiac electrophysiology. However, a direct observation of the activation sequence throughout the ventricular volume is virtually impossible. In this study, we report on a novel method for identification of activation sequences from activation maps measured at the outer surface of the heart termed the epicardium. Conceptually, the method attempts to identify the key factors governing the ventricular activation sequence - the timing of earliest activation sites (EAS) and the velocity tensor field within the ventricular walls - from sparse and noisy activation maps sampled from the epicardial surface and fits an Eikonal model to the observations. Regularization methods are first investigated to overcome the severe ill-posedness of the inverse problem in a simplified 2D example. These methods are then employed in an anatomically accurate biventricular model with two realistic activation models of varying complexity - a simplified trifascicular model (3F) and a topologically realistic model of the His-Purkinje system (HPS). Using epicardial activation maps at full resolution, we first demonstrate that reconstructing the volumetric activation sequence is, in principle, feasible under the assumption of known location of EAS and later evaluate robustness of the method against noise and reduced spatial resolution of observations. Our results suggest that the FIMIN algorithm is able to robustly recover the full 3D activation sequence using epicardial activation maps at a spatial resolution achievable with current mapping systems and in the presence of noise. Comparing the accuracy achieved in the reconstructed activation maps with clinical data uncertainties suggests that the FIMIN method may be suitable for the patient- specific parameterization of activation models.

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.

Automatic estimation of Purkinje-Myocardial junction hot-spots from noisy endocardial samples: A simulation study

The reconstruction of the ventricular cardiac conduction system (CCS) from patient-specific data is a challenging problem. High-resolution imaging techniques have allowed only the segmentation of proximal sections of the CCS from images acquired ex vivo. In this paper, we present an algorithm to estimate the location of a set of Purkinje-myocardial junctions (PMJs) from electro-anatomical maps, as those acquired during radio-frequency ablation procedures. The method requires a mesh representing the myocardium with local activation time measurements on a subset of nodes. We calculate the backwards propagation of the electrical signal from the measurement points to all the points in the mesh to define a set of candidate PMJs that is iteratively refined. The algorithm has been tested on several Purkinje network configurations, with simulated activation maps, subject to different error amplitudes. The results show that the method is able to build a set of PMJs that explain the observed activation map for different synthetic CCS configurations. In the tests, the average error in the predicted activation time is below the amplitude of the error applied to the data.

Numerical approximation of the electromechanical coupling in the left ventricle with inclusion of the Purkinje network

In this work, we consider the numerical approximation of the electromechanical coupling in the left ventricle with inclusion of the Purkinje network. The mathematical model couples the 3D elastodynamics and bidomain equations for the electrophysiology in the myocardium with the 1D monodomain equation in the Purkinje network. For the numerical solution of the coupled problem, we consider a fixed-point iterative algorithm that enables a partitioned solution of the myocardium and Purkinje network problems. Different levels of myocardium-Purkinje network splitting are considered and analyzed. The results are compared with those obtained using standard strategies proposed in the literature to trigger the electrical activation. Finally, we present a numerical study that, although performed in an idealized computational domain, features all the physiological issues that characterize a heartbeat simulation, including the initiation of the signal in the Purkinje network and the systolic and diastolic phases.

Inverse estimation of cardiac activation times via gradient-based optimization

Computational modeling may provide a quantitative framework for integrating multiscale data to gain insight into mechanisms of heart disease, identify and test pharmacological and electrical therapy and interventions, and support clinical decisions. Patient-specific computational cardiac models can help guide such procedures, and cardiac inverse modeling is a promising alternative to adequately personalize these models. Indeed, full cardiac inverse modeling is currently becoming computationally feasible; however, fundamental work to assess the feasibility of emerging techniques is still needed. In this study, we use a partial differential equation-constrained optimal control approach to numerically investigate the identifiability of an initial activation sequence from synthetic (partial) observations of the extracellular potential using the bidomain approximation and 2D representations of cardiac tissue. Our results demonstrate that activation times and duration of several stimuli can be recovered even with high levels of noise, that it is sufficient to sample the observations at the electrocardiogram-relevant sampling frequency of 1 kHz, and that spatial resolutions that are coarser than the standard in electrophysiological simulations can be used. The optimization of activation times is still effective when synthetic data are generated with a different cell membrane kinetics model than optimized for. The findings thus indicate that the presented approach has potential for finding activation sequences from clinical data modalities, as an extension to existing cardiac imaging approaches.

Generating Purkinje networks in the human heart

The Purkinje network is an integral part of the excitation system in the human heart. Yet, to date, there is no in vivo imaging technique to accurately reconstruct its geometry and structure. Computational modeling of the Purkinje network is increasingly recognized as an alternative strategy to visualize, simulate, and understand the role of the Purkinje system. However, most computational models either have to be generated manually, or fail to smoothly cover the irregular surfaces inside the left and right ventricles. Here we present a new algorithm to reliably create robust Purkinje networks within the human heart. We made the source code of this algorithm freely available online. Using Monte Carlo simulations, we demonstrate that the fractal tree algorithm with our new projection method generates denser and more compact Purkinje networks than previous approaches on irregular surfaces. Under similar conditions, our algorithm generates a network with 1219±61 branches, three times more than a conventional algorithm with 419±107 branches. With a coverage of 11±3mm, the surface density of our new Purkije network is twice as dense as the conventional network with 22±7mm. To demonstrate the importance of a dense Purkinje network in cardiac electrophysiology, we simulated three cases of excitation: with our new Purkinje network, with left-sided Purkinje network, and without Purkinje network. Simulations with our new Purkinje network predicted more realistic activation sequences and activation times than simulations without. Six-lead electrocardiograms of the three case studies agreed with the clinical electrocardiograms under physiological conditions, under pathological conditions of right bundle branch block, and under pathological conditions of trifascicular block. Taken together, our results underpin the importance of the Purkinje network in realistic human heart simulations. Human heart modeling has the potential to support the design of personalized strategies for single- or bi-ventricular pacing, radiofrequency ablation, and cardiac defibrillation with the common goal to restore a normal heart rhythm.

Estimation of Purkinje trees from electro-anatomical mapping of the left ventricle using minimal cost geodesics

The electrical activation of the heart is a complex physiological process that is essential for the understanding of several cardiac dysfunctions, such as ventricular tachycardia (VT). Nowadays, patient-specific activation times on ventricular chambers can be estimated from electro-anatomical maps, providing crucial information to clinicians for guiding cardiac radio-frequency ablation treatment. However, some relevant electrical pathways such as those of the Purkinje system are very difficult to interpret from these maps due to sparsity of data and the limited spatial resolution of the system. We present here a novel method to estimate these fast electrical pathways from the local activations maps (LATs) obtained from electro-anatomical maps. The location of Purkinje-myocardial junctions (PMJs) is estimated considering them as critical points of a distance map defined by the activation maps, and then minimal cost geodesic paths are computed on the ventricular surface between the detected junctions. Experiments to validate the proposed method have been carried out in simplified and realistic simulated data, showing good performance on recovering the main characteristics of simulated Purkinje networks (e.g. PMJs). A feasibility study with real cases of fascicular VT was also performed, showing promising results.

Three-dimensional cardiac computational modelling: methods, features and applications

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

Patient-specific generation of the Purkinje network driven by clinical measurements of a normal propagation

The propagation of the electrical signal in the Purkinje network is the starting point for the activation of the ventricular muscular cells leading to the contraction of the ventricle. In the computational models, describing the electrical activity of the ventricle is therefore important to account for the Purkinje fibers. Until now, the inclusion of such fibers has been obtained either by using surrogates such as space-dependent conduction properties or by generating a network based on an a priori anatomical knowledge. The aim of this work was to propose a new method for the generation of the Purkinje network using clinical measures of the activation times on the endocardium related to a normal electrical propagation, allowing to generate a patient-specific network. The measures were acquired by means of the EnSite NavX system. This system allows to measure for each point of the ventricular endocardium the time at which the activation front, that spreads through the ventricle, has reached the subjacent muscle. We compared the accuracy of the proposed method with the one of other strategies proposed so far in the literature for three subjects with a normal electrical propagation. The results showed that with our method we were able to reduce the absolute errors, intended as the difference between the measured and the computed data, by a factor in the range 9-25 %, with respect to the best of the other strategies. This highlighted the reliability of the proposed method and the importance of including a patient-specific Purkinje network in computational models.