Digital twinning of cardiac electrophysiology for congenital heart disease.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in paediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and using rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in paediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.

Whole-heart electromechanical simulations using Latent Neural Ordinary Differential Equations.

Cardiac digital twins provide a physics and physiology informed framework to deliver personalized medicine. However, high-fidelity multi-scale cardiac models remain a barrier to adoption due to their extensive computational costs. Artificial Intelligence-based methods can make the creation of fast and accurate whole-heart digital twins feasible. We use Latent Neural Ordinary Differential Equations (LNODEs) to learn the pressure-volume dynamics of a heart failure patient. Our surrogate model is trained from 400 simulations while accounting for 43 parameters describing cell-to-organ cardiac electromechanics and cardiovascular hemodynamics. LNODEs provide a compact representation of the 3D-0D model in a latent space by means of an Artificial Neural Network that retains only 3 hidden layers with 13 neurons per layer and allows for numerical simulations of cardiac function on a single processor. We employ LNODEs to perform global sensitivity analysis and parameter estimation with uncertainty quantification in 3 hours of computations, still on a single processor.

A Modular Framework for Implicit 3D-0D Coupling in Cardiac Mechanics.

In numerical simulations of cardiac mechanics, coupling the heart to a model of the circulatory system is essential for capturing physiological cardiac behavior. A popular and efficient technique is to use an electrical circuit analogy, known as a lumped parameter network or zero-dimensional (0D) fluid model, to represent blood flow throughout the cardiovascular system. Due to the strong physical interaction between the heart and the blood circulation, developing accurate and efficient numerical coupling methods remains an active area of research. In this work, we present a modular framework for implicitly coupling three-dimensional (3D) finite element simulations of cardiac mechanics to 0D models of blood circulation. The framework is modular in that the circulation model can be modified independently of the 3D finite element solver, and vice versa. The numerical scheme builds upon a previous work that combines 3D blood flow models with 0D circulation models (3D fluid - 0D fluid). Here, we extend it to couple 3D cardiac tissue mechanics models with 0D circulation models (3D structure - 0D fluid), showing that both mathematical problems can be solved within a unified coupling scheme. The effectiveness, temporal convergence, and computational cost of the algorithm are assessed through multiple examples relevant to the cardiovascular modeling community. Importantly, in an idealized left ventricle example, we show that the coupled model yields physiological pressure-volume loops and naturally recapitulates the isovolumic contraction and relaxation phases of the cardiac cycle without any additional numerical techniques. Furthermore, we provide a new derivation of the scheme inspired by the Approximate Newton Method of Chan (1985), explaining how the proposed numerical scheme combines the stability of monolithic approaches with the modularity and flexibility of partitioned approaches.

Learning the intrinsic dynamics of spatio-temporal processes through Latent Dynamics Networks.

Predicting the evolution of systems with spatio-temporal dynamics in response to external stimuli is essential for scientific progress. Traditional equations-based approaches leverage first principles through the numerical approximation of differential equations, thus demanding extensive computational resources. In contrast, data-driven approaches leverage deep learning algorithms to describe system evolution in low-dimensional spaces. We introduce an architecture, termed Latent Dynamics Network, capable of uncovering low-dimensional intrinsic dynamics in potentially non-Markovian systems. Latent Dynamics Networks automatically discover a low-dimensional manifold while learning the system dynamics, eliminating the need for training an auto-encoder and avoiding operations in the high-dimensional space. They predict the evolution, even in time-extrapolation scenarios, of space-dependent fields without relying on predetermined grids, thus enabling weight-sharing across query-points. Lightweight and easy-to-train, Latent Dynamics Networks demonstrate superior accuracy (normalized error 5 times smaller) in highly-nonlinear problems with significantly fewer trainable parameters (more than 10 times fewer) compared to state-of-the-art methods.

Digital twinning of cardiac electrophysiology for congenital heart disease.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in pediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and utilizing rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in pediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations.

BACKGROUND: Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community. RESULTS: This work introduces [Formula: see text]-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. [Formula: see text]-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, [Formula: see text]-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within [Formula: see text]-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying [Formula: see text]-ep, along with comprehensive implementation details and instructions for users. [Formula: see text]-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of [Formula: see text]-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files. CONCLUSIONS: [Formula: see text]-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. [Formula: see text]-ep represents a valuable tool for conducting in silico patient-specific simulations.

Fast and robust parameter estimation with uncertainty quantification for the cardiac function.

BACKGROUND AND OBJECTIVES: Parameter estimation and uncertainty quantification are crucial in computational cardiology, as they enable the construction of digital twins that faithfully replicate the behavior of physical patients. Many model parameters regarding cardiac electromechanics and cardiovascular hemodynamics need to be robustly fitted by starting from a few, possibly non-invasive, noisy observations. Moreover, short execution times and a small amount of computational resources are required for the effective clinical translation. METHODS: In the framework of Bayesian statistics, we combine Maximum a Posteriori estimation and Hamiltonian Monte Carlo to find an approximation of model parameters and their posterior distributions. Fast simulations and minimal memory requirements are achieved by using an accurate and geometry-specific Artificial Neural Network surrogate model for the cardiac function, matrix-free methods, automatic differentiation and automatic vectorization. Furthermore, we account for the surrogate modeling error and measurement error. RESULTS: We perform three different in silico test cases, ranging from the ventricular function to the entire cardiocirculatory system, involving whole-heart mechanics, arterial and venous hemodynamics. By employing a single central processing unit on a standard laptop, we attain highly accurate estimations for all model parameters in short computational times. Furthermore, we obtain posterior distributions that contain the true values inside the 90% credibility regions. CONCLUSIONS: Many model parameters regarding the entire cardiovascular system can be fastly and robustly identified with minimal hardware requirements. This can be achieved when a small amount of non-invasive data is available and when high levels of signal-to-noise ratio are present in the quantities of interest. With these features, our approach meets the requirements for clinical exploitation, while being compliant with Green Computing practices.

The role of mechano-electric feedbacks and hemodynamic coupling in scar-related ventricular tachycardia.

Mechano-electric feedbacks (MEFs), which model how mechanical stimuli are transduced into electrical signals, have received sparse investigation by considering electromechanical simulations in simplified scenarios. In this paper, we study the effects of different MEFs modeling choices for myocardial deformation and nonselective stretch-activated channels (SACs) in the monodomain equation. We perform numerical simulations during ventricular tachycardia (VT) by employing a biophysically detailed and anatomically accurate 3D electromechanical model for the left ventricle (LV) coupled with a 0D closed-loop model of the cardiocirculatory system. We model the electromechanical substrate responsible for scar-related VT with a distribution of infarct and peri-infarct zones. Our mathematical framework takes into account the hemodynamic effects of VT due to myocardial impairment and allows for the classification of their hemodynamic nature, which can be either stable or unstable. By combining electrophysiological, mechanical and hemodynamic models, we observe that all MEFs may alter the propagation of the transmembrane potential. In particular, we notice that the presence of myocardial deformation in the monodomain equation may change the VT basis cycle length and the conduction velocity but do not affect the hemodynamic nature of the VT. Finally, nonselective SACs may affect VT stability, by possibly turning a hemodynamically stable VT into a hemodynamically unstable one.

Electromechanical modeling of human ventricles with ischemic cardiomyopathy: numerical simulations in sinus rhythm and under arrhythmia.

We developed a novel patient-specific computational model for the numerical simulation of ventricular electromechanics in patients with ischemic cardiomyopathy (ICM). This model reproduces the activity both in sinus rhythm (SR) and in ventricular tachycardia (VT). The presence of scars, grey zones and non-remodeled regions of the myocardium is accounted for by the introduction of a spatially heterogeneous coefficient in the 3D electromechanics model. This 3D electromechanics model is firstly coupled with a 2-element Windkessel afterload model to fit the pressure-volume (PV) loop of a patient-specific left ventricle (LV) with ICM in SR. Then, we employ the coupling with a 0D closed-loop circulation model to analyze a VT circuit over multiple heartbeats on the same LV. We highlight similarities and differences on the solutions obtained by the electrophysiology model and those of the electromechanics model, while considering different scenarios for the circulatory system. We observe that very different parametrizations of the circulation model induce the same hemodynamical considerations for the patient at hand. Specifically, we classify this VT as unstable. We conclude by stressing the importance of combining electrophysiological, mechanical and hemodynamical models to provide relevant clinical indicators in how arrhythmias evolve and can potentially lead to sudden cardiac death.

Presence of Left Atrial Fibrosis May Contribute to Aberrant Hemodynamics and Increased Risk of Stroke in Atrial Fibrillation Patients.

Atrial fibrillation (AF) patients are at high risk of stroke, with the left atrial appendage (LAA) found to be the most common site of clot formation. Presence of left atrial (LA) fibrosis has also been associated with higher stroke risk. However, the mechanisms for increased stroke risk in patients with atrial fibrotic remodeling are poorly understood. We sought to explore these mechanisms using fluid dynamic analysis and to test the hypothesis that the presence of LA fibrosis leads to aberrant hemodynamics in the LA, contributing to increased stroke risk in AF patients. We retrospectively collected late-gadolinium-enhanced MRI (LGE-MRI) images of eight AF patients (four persistent and four paroxysmal) and reconstructed their 3D LA surfaces. Personalized computational fluid dynamic simulations were performed, and hemodynamics at the LA wall were quantified by wall shear stress (WSS, friction of blood), oscillatory shear index (OSI, temporal directional change of WSS), endothelial cell activation potential (ECAP, ratio of OSI and WSS), and relative residence time (RRT, residence time of blood near the LA wall). For each case, these hemodynamic metrics were compared between fibrotic and non-fibrotic portions of the wall. Our results showed that WSS was lower, and OSI, ECAP, and RRT was higher in the fibrotic region as compared to the non-fibrotic region, with ECAP (p = 0.001) and RRT (p = 0.002) having significant differences. Case-wise analysis showed that these differences in hemodynamics were statistically significant for seven cases. Furthermore, patients with higher fibrotic burden were exposed to larger regions of high ECAP, which represents regions of low WSS and high OSI. Consistently, high ECAP in the vicinity of the fibrotic wall suggest that local blood flow was slow and oscillating that represents aberrant hemodynamic conditions, thus enabling prothrombotic conditions for circulating blood. AF patients with high LA fibrotic burden had more prothrombotic regions, providing more sites for potential clot formation, thus increasing their risk of stroke.