Verification of cardiac tissue electrophysiology simulators using an N-version benchmark

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.

Effects of cardiac growth on electrical dyssynchrony in the single ventricle patient

Single ventricle patients, including those with hypoplastic left heart syndrome (HLHS), typically undergo three palliative heart surgeries culminating in the Fontan procedure. HLHS is associated with high rates of morbidity and mortality, and many patients develop arrhythmias, electrical dyssynchrony, and eventually ventricular failure. However, the correlation between ventricular enlargement and electrical dysfunction in HLHS physiology remains poorly understood. Here we characterize the relationship between growth and electrophysiology in HLHS using computational modeling. We integrate a personalized finite element model, a volumetric growth model, and a personalized electrophysiology model to perform controlled in silico experiments. We show that right ventricle enlargement negatively affects QRS duration and interventricular dyssynchrony. Conversely, left ventricle enlargement can partially compensate for this dyssynchrony. These findings have potential implications on our understanding of the origins of electrical dyssynchrony and, ultimately, the treatment of HLHS patients.

Computational Modeling and Simulation of the Fibrotic Human Atria

Patient-specific modeling of atrial electrical activity enables the execution of simulations that can provide mechanistic insights and provide novel solutions to vexing clinical problems. The geometry and fibrotic remodeling of the heart can be reconstructed from clinical-grade medical scans and used to inform personalized models with detail incorporated at the cell- and tissue-scale to represent changes in image-identified diseased regions. Here, we provide a rubric for the reconstruction of realistic atrial models from pre-segmented 3D renderings of the left atrium with fibrotic tissue regions delineated, which are the output from clinical-grade systems for quantifying fibrosis. We then provide a roadmap for using those models to carry out patient-specific characterization of the fibrotic substrate in terms of its potential to harbor reentrant drivers via cardiac electrophysiology simulations.

Constructing bilayer and volumetric atrial models at scale

To enable large in silico trials and personalized model predictions on clinical timescales, it is imperative that models can be constructed quickly and reproducibly. First, we aimed to overcome the challenges of constructing cardiac models at scale through developing a robust, open-source pipeline for bilayer and volumetric atrial models. Second, we aimed to investigate the effects of fibres, fibrosis and model representation on fibrillatory dynamics. To construct bilayer and volumetric models, we extended our previously developed coordinate system to incorporate transmurality, atrial regions and fibres (rule-based or data driven diffusion tensor magnetic resonance imaging (MRI)). We created a cohort of 1000 biatrial bilayer and volumetric models derived from computed tomography (CT) data, as well as models from MRI, and electroanatomical mapping. Fibrillatory dynamics diverged between bilayer and volumetric simulations across the CT cohort (correlation coefficient for phase singularity maps: left atrial (LA) 0.27 ± 0.19, right atrial (RA) 0.41 ± 0.14). Adding fibrotic remodelling stabilized re-entries and reduced the impact of model type (LA: 0.52 ± 0.20, RA: 0.36 ± 0.18). The choice of fibre field has a small effect on paced activation data (less than 12 ms), but a larger effect on fibrillatory dynamics. Overall, we developed an open-source user-friendly pipeline for generating atrial models from imaging or electroanatomical mapping data enabling in silico clinical trials at scale (https://github.com/pcmlab/atrialmtk).

The simplified Kirchhoff network model (SKNM): a cell-based reaction-diffusion model of excitable tissue

Cell-based models of excitable tissues offer the advantage of cell-level precision, which cannot be achieved using traditional homogenized electrophysiological models. However, this enhanced accuracy comes at the cost of increased computational demands, necessitating the development of efficient cell-based models. The widely-accepted bidomain model serves as the standard in computational cardiac electrophysiology, and under certain anisotropy ratio conditions, it is well known that it can be reduced to the simpler monodomain model. Recently, the Kirchhoff Network Model (KNM) was developed as a cell-based counterpart to the bidomain model. In this paper, we aim to demonstrate that KNM can be simplified using the same steps employed to derive the monodomain model from the bidomain model. We present the cell-based Simplified Kirchhoff Network Model (SKNM), which produces results closely aligned with those of KNM while requiring significantly less computational resources.

AugmentA: Patient-specific augmented atrial model generation tool

Digital twins of patients' hearts are a promising tool to assess arrhythmia vulnerability and to personalize therapy. However, the process of building personalized computational models can be challenging and requires a high level of human interaction. We propose a patient-specific Augmented Atria generation pipeline (AugmentA) as a highly automated framework which, starting from clinical geometrical data, provides ready-to-use atrial personalized computational models. AugmentA identifies and labels atrial orifices using only one reference point per atrium. If the user chooses to fit a statistical shape model to the input geometry, it is first rigidly aligned with the given mean shape before a non-rigid fitting procedure is applied. AugmentA automatically generates the fiber orientation and finds local conduction velocities by minimizing the error between the simulated and clinical local activation time (LAT) map. The pipeline was tested on a cohort of 29 patients on both segmented magnetic resonance images (MRI) and electroanatomical maps of the left atrium. Moreover, the pipeline was applied to a bi-atrial volumetric mesh derived from MRI. The pipeline robustly integrated fiber orientation and anatomical region annotations in 38.4 ± 5.7 s. In conclusion, AugmentA offers an automated and comprehensive pipeline delivering atrial digital twins from clinical data in procedural time.

Evaluation of an open-source pipeline to create patient-specific left atrial models: A reproducibility study

This work presents an open-source software pipeline to create patient-specific left atrial models with fibre orientations and a fibrDEFAULTosis map, suitable for electrophysiology simulations, and quantifies the intra and inter observer reproducibility of the model creation. The semi-automatic pipeline takes as input a contrast enhanced magnetic resonance angiogram, and a late gadolinium enhanced (LGE) contrast magnetic resonance (CMR). Five operators were allocated 20 cases each from a set of 50 CMR datasets to create a total of 100 models to evaluate inter and intra-operator variability. Each output model consisted of: (1) a labelled surface mesh open at the pulmonary veins and mitral valve, (2) fibre orientations mapped from a diffusion tensor MRI (DTMRI) human atlas, (3) fibrosis map extracted from the LGE-CMR scan, and (4) simulation of local activation time (LAT) and phase singularity (PS) mapping. Reproducibility in our pipeline was evaluated by comparing agreement in shape of the output meshes, fibrosis distribution in the left atrial body, and fibre orientations. Reproducibility in simulations outputs was evaluated in the LAT maps by comparing the total activation times, and the mean conduction velocity (CV). PS maps were compared with the structural similarity index measure (SSIM). The users processed in total 60 cases for inter and 40 cases for intra-operator variability. Our workflow allows a single model to be created in 16.72 ± 12.25 min. Similarity was measured with shape, percentage of fibres oriented in the same direction, and intra-class correlation coefficient (ICC) for the fibrosis calculation. Shape differed noticeably only with users' selection of the mitral valve and the length of the pulmonary veins from the ostia to the distal end; fibrosis agreement was high, with ICC of 0.909 (inter) and 0.999 (intra); fibre orientation agreement was high with 60.63% (inter) and 71.77% (intra). The LAT showed good agreement, where the median ± IQR of the absolute difference of the total activation times was 2.02 ± 2.45 ms for inter, and 1.37 ± 2.45 ms for intra. Also, the average ± sd of the mean CV difference was -0.00404 ± 0.0155 m/s for inter, and 0.0021 ± 0.0115 m/s for intra. Finally, the PS maps showed a moderately good agreement in SSIM for inter and intra, where the mean ± sd SSIM for inter and intra were 0.648 ± 0.21 and 0.608 ± 0.15, respectively. Although we found notable differences in the models, as a consequence of user input, our tests show that the uncertainty caused by both inter and intra-operator variability is comparable with uncertainty due to estimated fibres, and image resolution accuracy of segmentation tools.

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.

An inner-outer subcycling algorithm for parallel cardiac electrophysiology simulations

This paper explores cardiac electrophysiological simulations of the monodomain equations and introduces a novel subcycling time integration algorithm to exploit the structure of the ionic model. The aim of this work is to improve upon the efficiency of parallel cardiac monodomain simulations by using our subcycling algorithm in the computation of the ionic model to handle the local sharp changes of the solution. This will reduce the turnaround time for the simulation of basic cardiac electrical function on both idealized and patient-specific geometry. Numerical experiments show that the proposed approach is accurate and also has close to linear parallel scalability on a computer with more than 1000 processor cores. Ultimately, the reduction in simulation time can be beneficial in clinical applications, where multiple simulations are often required to tune a model to match clinical measurements.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

A simple and efficient adaptive time stepping technique for low-order operator splitting schemes applied to cardiac electrophysiology

We present a simple, yet efficient adaptive time stepping scheme for cardiac electrophysiology (EP) simulations based on standard operator splitting techniques. The general idea is to exploit the relation between the splitting error and the reaction's magnitude-found in a previous one-dimensional analytical study by Spiteri and Ziaratgahi-to construct the new time step controller for three-dimensional problems. Accordingly, we propose to control the time step length of the operator splitting scheme as a function of the reaction magnitude, in addition to the common approach of adapting the reaction time step. This conforms with observations in numerical experiments supporting the need for a significantly smaller time step length during depolarization than during repolarization. The proposed scheme is compared with classical proportional-integral-differential controllers using state-of-the-art error estimators, which are also presented in details as they have not been previously applied in the context of cardiac EP with operator splitting techniques. Benchmarks show that choosing the time step as a sigmoidal function of the reaction magnitude is highly efficient and full cardiac cycles can be computed with precision even in a realistic biventricular setup. The proposed scheme outperforms common adaptive time stepping techniques, while depending on fewer tuning parameters.

The mechanisms of potassium loss in acute myocardial ischemia: New insights from computational simulations

Acute myocardial ischemia induces hyperkalemia (accumulation of extracellular potassium), a major perpetrator of lethal reentrant ventricular arrhythmias. Despite considerable experimental efforts to explain this pathology in the last decades, the intimate mechanisms behind hyperkalemia remain partially unknown. In order to investigate these mechanisms, we developed a novel computational model of acute myocardial ischemia which couples a) an electrophysiologically detailed human cardiomyocyte model that incorporates modifications to account for ischemia-induced changes in transmembrane currents, with b) a model of cardiac tissue and extracellular K ⁺ transport. The resulting model is able to reproduce and explain the triphasic time course of extracellular K ⁺ concentration within the ischemic zone, with values of [ K + ] o close to 14 mmol/L in the central ischemic zone after 30 min. In addition, the formation of a [ K + ] o border zone of approximately 1.2 cm 15 min after the onset of ischemia is predicted by the model. Our results indicate that the primary rising phase of [ K + ] o is mainly due to the imbalance between K ⁺ efflux, that increases slightly, and K ⁺ influx, that follows a reduction of the NaK pump activity by more than 50%. The onset of the plateau phase is caused by the appearance of electrical alternans (a novel mechanism identified by the model), which cause an abrupt reduction in the K ⁺ efflux. After the plateau, the secondary rising phase of [ K + ] o is caused by a subsequent imbalance between the K ⁺ influx, which continues to decrease slowly, and the K ⁺ efflux, which remains almost constant. Further, the study shows that the modulation of these mechanisms by the electrotonic coupling is the main responsible for the formation of the ischemic border zone in tissue, with K ⁺ transport playing only a minor role. Finally, the results of the model indicate that the injury current established between the healthy and the altered tissue is not sufficient to depolarize non-ischemic cells within the healthy tissue.

Material property alterations for phenotypes of heart failure with preserved ejection fraction: A numerical study of subject-specific porcine models

A substantial proportion of heart failure patients have a preserved left ventricular (LV) ejection fraction (HFpEF). This condition carries a high burden of morbidity and mortality and has limited therapeutic options. left ventricular pressure overload leads to an increase in myocardial collagen content, causing left ventricular stiffening that contributes to the development of heart failure patients have a preserved left ventricular ejection fraction. Although several heart failure patients have a preserved left ventricular ejection fraction models have been developed in recent years to aid the investigation of mechanical alterations, none has investigated different phenotypes of the disease and evaluated the alterations in material properties. In this study, two similar healthy swine were subjected to progressive and prolonged pressure overload to induce diastolic heart failure characteristics, providing a preclinical model of heart failure patients have a preserved left ventricular ejection fraction. Cardiac magnetic resonance imaging (cMRI) scans and intracardiac pressures were recorded before and after induction. In both healthy and disease states, a corresponding finite element (FE) cardiac model was developed via mesh morphing of the Living Heart Porcine model. The material properties were derived by calibrating to its passive and active behavior. The change in the passive behavior was predominantly isotropic when comparing the geometries before and after induction. Myocardial thickening allowed for a steady transition in the passive properties while maintaining tissue incompressibility. This study highlights the importance of hypertrophy as an initial compensatory response and might also pave the way for assessing disease severity.

Rule-based definition of muscle bundles in patient-specific models of the left atrium

Atrial fibrillation (AF) is the most common arrhythmia encountered clinically, and as the population ages, its prevalence is increasing. Although the CHA2DS2- VASc score is the most used risk-stratification system for stroke risk in AF, it lacks personalization. Patient-specific computer models of the atria can facilitate personalized risk assessment and treatment planning. However, a challenge faced in creating such models is the complexity of the atrial muscle arrangement and its influence on the atrial fiber architecture. This work proposes a semi-automated rule-based algorithm to generate the local fiber orientation in the left atrium (LA). We use the solutions of several harmonic equations to decompose the LA anatomy into subregions. Solution gradients define a two-layer fiber field in each subregion. The robustness of our approach is demonstrated by recreating the fiber orientation on nine models of the LA obtained from AF patients who underwent WATCHMAN device implantation. This cohort of patients encompasses a variety of morphology variants of the left atrium, both in terms of the left atrial appendages (LAAs) and the number of pulmonary veins (PVs). We test the fiber construction algorithm by performing electrophysiology (EP) simulations. Furthermore, this study is the first to compare its results with other rule-based algorithms for the LA fiber architecture definition available in the literature. This analysis suggests that a multi-layer fiber architecture is important to capture complex electrical activation patterns. A notable advantage of our approach is the ability to reconstruct the main LA fiber bundles in a variety of morphologies while solving for a small number of harmonic fields, leading to a comparatively straightforward and reproducible approach.

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

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

Patient-specific models are computational models that have been personalized using data from a patient. After decades of research, recent computational, data science and healthcare advances have opened the door to the fulfilment of the enormous potential of such models, from truly personalized medicine to efficient and cost-effective testing of new medical products. However, reliability (credibility) of patient-specific models is key to their success, and there are currently no general guidelines for evaluating credibility of patient-specific models. Here, we consider how frameworks and model evaluation activities that have been developed for generic (not patient-specific) computational models, can be extended to patient specific models. We achieve this through a detailed analysis of the activities required to evaluate cardiac electrophysiological models, chosen as an exemplar field due to its maturity and the complexity of such models. This is the first paper on the topic of reliability of patient-specific models and will help pave the way to reliable and trusted patient-specific modeling across healthcare applications.

Balancing conduction velocity error in cardiac electrophysiology using a modified quadrature approach

Conduction velocity error is often the main culprit behind the need for very fine spatial discretizations and high computational effort in cardiac electrophysiology problems. In light of this, a novel approach for simulating an accurate conduction velocity in coarse meshes with linear elements is suggested based on a modified quadrature approach. In this approach, the quadrature points are placed at arbitrary offsets of the isoparametric coordinates. A numerical study illustrates the dependence of the conduction velocity on the spatial discretization and the conductivity when using different quadrature rules and calculation approaches. Additionally, examples using the modified quadrature in coarse meshes for wave propagation demonstrate the improved accuracy of the conduction velocity with this method. This novel approach possesses great potential in reducing the computational effort required but remains limited to specific linear elements and experiences a reduction in accuracy for irregular meshes and heterogeneous conductivities. Further research can focus on developing an adaptive quadrature and extending the approach to other element formulations in order to make the approach more generally applicable.

Arrhythmogenic influence of mutations in a myocyte-based computational model of the pulmonary vein sleeve

In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte-to-myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI's cell-based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI's improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.

An Automata-Based Cardiac Electrophysiology Simulator to Assess Arrhythmia Inducibility

Personalized cardiac electrophysiology simulations have demonstrated great potential to study cardiac arrhythmias and help in therapy planning of radio-frequency ablation. Its application to analyze vulnerability to ventricular tachycardia and sudden cardiac death in infarcted patients has been recently explored. However, the detailed multi-scale biophysical simulations used in these studies are very demanding in terms of memory and computational resources, which prevents their clinical translation. In this work, we present a fast phenomenological system based on cellular automata (CA) to simulate personalized cardiac electrophysiology. The system is trained on biophysical simulations to reproduce cellular and tissue dynamics in healthy and pathological conditions, including action potential restitution, conduction velocity restitution and cell safety factor. We show that a full ventricular simulation can be performed in the order of seconds, emulate the results of a biophysical simulation and reproduce a patient’s ventricular tachycardia in a model that includes a heterogeneous scar region. The system could be used to study the risk of arrhythmia in infarcted patients for a large number of scenarios.

An Integrated Workflow for Building Digital Twins of Cardiac Electromechanics—A Multi-Fidelity Approach for Personalising Active Mechanics

Personalised computer models of cardiac function, referred to as cardiac digital twins, are envisioned to play an important role in clinical precision therapies of cardiovascular diseases. A major obstacle hampering clinical translation involves the significant computational costs involved in the personalisation of biophysically detailed mechanistic models that require the identification of high-dimensional parameter vectors. An important aspect to identify in electromechanics (EM) models are active mechanics parameters that govern cardiac contraction and relaxation. In this study, we present a novel, fully automated, and efficient approach for personalising biophysically detailed active mechanics models using a two-step multi-fidelity solution. In the first step, active mechanical behaviour in a given 3D EM model is represented by a purely phenomenological, low-fidelity model, which is personalised at the organ scale by calibration to clinical cavity pressure data. Then, in the second step, median traces of nodal cellular active stress, intracellular calcium concentration, and fibre stretch are generated and utilised to personalise the desired high-fidelity model at the cellular scale using a 0D model of cardiac EM. Our novel approach was tested on a cohort of seven human left ventricular (LV) EM models, created from patients treated for aortic coarctation (CoA). Goodness of fit, computational cost, and robustness of the algorithm against uncertainty in the clinical data and variations of initial guesses were evaluated. We demonstrate that our multi-fidelity approach facilitates the personalisation of a biophysically detailed active stress model within only a few (2 to 4) expensive 3D organ-scale simulations—a computational effort compatible with clinical model applications.

Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation

Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.

Predicting Atrial Fibrillation Recurrence by Combining Population Data and Virtual Cohorts of Patient-Specific Left Atrial Models

BACKGROUND

Current ablation therapy for atrial fibrillation is suboptimal, and long-term response is challenging to predict. Clinical trials identify bedside properties that provide only modest prediction of long-term response in populations, while patient-specific models in small cohorts primarily explain acute response to ablation. We aimed to predict long-term atrial fibrillation recurrence after ablation in large cohorts, by using machine learning to complement biophysical simulations by encoding more interindividual variability.

METHODS

Patient-specific models were constructed for 100 atrial fibrillation patients (43 paroxysmal, 41 persistent, and 16 long-standing persistent), undergoing first ablation. Patients were followed for 1 year using ambulatory ECG monitoring. Each patient-specific biophysical model combined differing fibrosis patterns, fiber orientation maps, electrical properties, and ablation patterns to capture uncertainty in atrial properties and to test the ability of the tissue to sustain fibrillation. These simulation stress tests of different model variants were postprocessed to calculate atrial fibrillation simulation metrics. Machine learning classifiers were trained to predict atrial fibrillation recurrence using features from the patient history, imaging, and atrial fibrillation simulation metrics.

RESULTS

We performed 1100 atrial fibrillation ablation simulations across 100 patient-specific models. Models based on simulation stress tests alone showed a maximum accuracy of 0.63 for predicting long-term fibrillation recurrence. Classifiers trained to history, imaging, and simulation stress tests (average 10-fold cross-validation area under the curve, 0.85±0.09; recall, 0.80±0.13; precision, 0.74±0.13) outperformed those trained to history and imaging (area under the curve, 0.66±0.17) or history alone (area under the curve, 0.61±0.14).

CONCLUSION

A novel computational pipeline accurately predicted long-term atrial fibrillation recurrence in individual patients by combining outcome data with patient-specific acute simulation response. This technique could help to personalize selection for atrial fibrillation ablation.

Deriving the Bidomain Model of Cardiac Electrophysiology From a Cell-Based Model; Properties and Comparisons

The bidomain model is considered to be the gold standard for numerical simulation of the electrophysiology of cardiac tissue. The model provides important insights into the conduction properties of the electrochemical wave traversing the cardiac muscle in every heartbeat. However, in normal resolution, the model represents the average over a large number of cardiomyocytes, and more accurate models based on representations of all individual cells have therefore been introduced in order to gain insight into the conduction properties close to the myocytes. The more accurate model considered here is referred to as the EMI model since both the extracellular space (E), the cell membrane (M) and the intracellular space (I) are explicitly represented in the model. Here, we show that the bidomain model can be derived from the cell-based EMI model and we thus reveal the close relation between the two models, and obtain an indication of the error introduced in the approximation. Also, we present numerical simulations comparing the results of the two models and thereby highlight both similarities and differences between the models. We observe that the deviations between the solutions of the models become larger for larger cell sizes. Furthermore, we observe that the bidomain model provides solutions that are very similar to the EMI model when conductive properties of the tissue are in the normal range, but large deviations are present when the resistance between cardiomyocytes is increased.

Seven challenges in the multiscale modeling of multicellular tissues

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

A Fully-Coupled Electro-Mechanical Whole-Heart Computational Model: Influence of Cardiac Contraction on the ECG

The ECG is one of the most commonly used non-invasive tools to gain insights into the electrical functioning of the heart. It has been crucial as a foundation in the creation and validation of in silico models describing the underlying electrophysiological processes. However, so far, the contraction of the heart and its influences on the ECG have mainly been overlooked in in silico models. As the heart contracts and moves, so do the electrical sources within the heart responsible for the signal on the body surface, thus potentially altering the ECG. To illuminate these aspects, we developed a human 4-chamber electro-mechanically coupled whole heart in silico model and embedded it within a torso model. Our model faithfully reproduces measured 12-lead ECG traces, circulatory characteristics, as well as physiological ventricular rotation and atrioventricular valve plane displacement. We compare our dynamic model to three non-deforming ones in terms of standard clinically used ECG leads (Einthoven and Wilson) and body surface potential maps (BSPM). The non-deforming models consider the heart at its ventricular end-diastatic, end-diastolic and end-systolic states. The standard leads show negligible differences during P-Wave and QRS-Complex, yet during T-Wave the leads closest to the heart show prominent differences in amplitude. When looking at the BSPM, there are no notable differences during the P-Wave, but effects of cardiac motion can be observed already during the QRS-Complex, increasing further during the T-Wave. We conclude that for the modeling of activation (P-Wave/QRS-Complex), the associated effort of simulating a complete electro-mechanical approach is not worth the computational cost. But when looking at ventricular repolarization (T-Wave) in standard leads as well as BSPM, there are areas where the signal can be influenced by cardiac motion of the heart to an extent that should not be ignored.

Multiscale Modeling of the Mitochondrial Origin of Cardiac Reentrant and Fibrillatory Arrhythmias

While mitochondrial dysfunction has been implicated in the pathogenesis of cardiac arrhythmias, how the abnormality occurring at the organelle level escalates to influence the rhythm of the heart remains incompletely understood. This is due, in part, to the complexity of the interactions formed by cardiac electrical, mechanical, and metabolic subsystems at various spatiotemporal scales that is difficult to fully comprehend solely with experiments. Computational models have emerged as a powerful tool to explore complicated and highly dynamic biological systems such as the heart, alone or in combination with experimental measurements. Here, we describe a strategy of integrating computer simulations with optical mapping of cardiomyocyte monolayers to examine how regional mitochondrial dysfunction elicits abnormal electrical activity, such as rebound and spiral waves, leading to reentry and fibrillation in cardiac tissue. We anticipate that this advanced modeling technology will enable new insights into the mechanisms by which changes in subcellular organelles can impact organ function.

A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation

Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.

A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation

Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.

From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology

Computational modeling has contributed significantly to present understanding of cardiac electrophysiology including cardiac conduction, excitation-contraction coupling, and the effects and side-effects of drugs. However, the accuracy of in silico analysis of electrochemical wave dynamics in cardiac tissue is limited by the homogenization procedure (spatial averaging) intrinsic to standard continuum models of conduction. Averaged models cannot resolve the intricate dynamics in the vicinity of individual cardiomyocytes simply because the myocytes are not present in these models. Here we demonstrate how recently developed mathematical models based on representing every myocyte can significantly increase the accuracy, and thus the utility of modeling electrophysiological function and dysfunction in collections of coupled cardiomyocytes. The present gold standard of numerical simulation for cardiac electrophysiology is based on the bidomain model. In the bidomain model, the extracellular (E) space, the cell membrane (M) and the intracellular (I) space are all assumed to be present everywhere in the tissue. Consequently, it is impossible to study biophysical processes taking place close to individual myocytes. The bidomain model represents the tissue by averaging over several hundred myocytes and this inherently limits the accuracy of the model. In our alternative approach both E, M, and I are represented in the model which is therefore referred to as the EMI model. The EMI model approach allows for detailed analysis of the biophysical processes going on in functionally important spaces very close to individual myocytes, although at the cost of significantly increased CPU-requirements.

Etiology-Specific Remodeling in Ventricular Tissue of Heart Failure Patients and Its Implications for Computational Modeling of Electrical Conduction

With an estimated 64.3 million cases worldwide, heart failure (HF) imposes an enormous burden on healthcare systems. Sudden death from arrhythmia is the major cause of mortality in HF patients. Computational modeling of the failing heart provides insights into mechanisms of arrhythmogenesis, risk stratification of patients, and clinical treatment. However, the lack of a clinically informed approach to model cardiac tissues in HF hinders progress in developing patient-specific strategies. Here, we provide a microscopy-based foundation for modeling conduction in HF tissues. We acquired 2D images of left ventricular tissues from HF patients (n = 16) and donors (n = 5). The composition and heterogeneity of fibrosis were quantified at a sub-micrometer resolution over an area of 1 mm². From the images, we constructed computational bidomain models of tissue electrophysiology. We computed local upstroke velocities of the membrane voltage and anisotropic conduction velocities (CV). The non-myocyte volume fraction was higher in HF than donors (39.68 ± 14.23 vs. 22.09 ± 2.72%, p < 0.01), and higher in ischemic (IC) than nonischemic (NIC) cardiomyopathy (47.2 ± 16.18 vs. 32.16 ± 6.55%, p < 0.05). The heterogeneity of fibrosis within each subject was highest for IC (27.1 ± 6.03%) and lowest for donors (7.47 ± 1.37%) with NIC (15.69 ± 5.76%) in between. K-means clustering of this heterogeneity discriminated IC and NIC with an accuracy of 81.25%. The heterogeneity in CV increased from donor to NIC to IC tissues. CV decreased with increasing fibrosis for longitudinal (R ² = 0.28, p < 0.05) and transverse conduction (R ² = 0.46, p < 0.01). The tilt angle of the CV vectors increased 2.1° for longitudinal and 0.91° for transverse conduction per 1% increase in fibrosis. Our study suggests that conduction fundamentally differs in the two etiologies due to the characteristics of fibrosis. Our study highlights the importance of the etiology-specific modeling of HF tissues and integration of medical history into electrophysiology models for personalized risk stratification and treatment planning.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.

Sex Differences in Drug-Induced Arrhythmogenesis

The electrical activity in the heart varies significantly between men and women and results in a sex-specific response to drugs. Recent evidence suggests that women are more than twice as likely as men to develop drug-induced arrhythmia with potentially fatal consequences. Yet, the sex-specific differences in drug-induced arrhythmogenesis remain poorly understood. Here we integrate multiscale modeling and machine learning to gain mechanistic insight into the sex-specific origin of drug-induced cardiac arrhythmia at differing drug concentrations. To quantify critical drug concentrations in male and female hearts, we identify the most important ion channels that trigger male and female arrhythmogenesis, and create and train a sex-specific multi-fidelity arrhythmogenic risk classifier. Our study reveals that sex differences in ion channel activity, tissue conductivity, and heart dimensions trigger longer QT-intervals in women than in men. We quantify the critical drug concentration for dofetilide, a high risk drug, to be seven times lower for women than for men. Our results emphasize the importance of including sex as an independent biological variable in risk assessment during drug development. Acknowledging and understanding sex differences in drug safety evaluation is critical when developing novel therapeutic treatments on a personalized basis. The general trends of this study have significant implications on the development of safe and efficacious new drugs and the prescription of existing drugs in combination with other drugs.

Cycle length statistics during human atrial fibrillation reveal refractory properties of the underlying substrate: a combined in silico and clinical test of concept study

AIMS

The treatment of atrial fibrillation beyond pulmonary vein isolation has remained an unsolved challenge. Targeting regions identified by different substrate mapping approaches for ablation resulted in ambiguous outcomes. With the effective refractory period being a fundamental prerequisite for the maintenance of fibrillatory conduction, this study aims at estimating the effective refractory period with clinically available measurements.

METHODS AND RESULTS

A set of 240 simulations in a spherical model of the left atrium with varying model initialization, combination of cellular refractory properties, and size of a region of lowered effective refractory period was implemented to analyse the capabilities and limitations of cycle length mapping. The minimum observed cycle length and the 25% quantile were compared to the underlying effective refractory period. The density of phase singularities was used as a measure for the complexity of the excitation pattern. Finally, we employed the method in a clinical test of concept including five patients. Areas of lowered effective refractory period could be distinguished from their surroundings in simulated scenarios with successfully induced multi-wavelet re-entry. Larger areas and higher gradients in effective refractory period as well as complex activation patterns favour the method. The 25% quantile of cycle lengths in patients with persistent atrial fibrillation was found to range from 85 to 190 ms.

CONCLUSION

Cycle length mapping is capable of highlighting regions of pathologic refractory properties. In combination with complementary substrate mapping approaches, the method fosters confidence to enhance the treatment of atrial fibrillation beyond pulmonary vein isolation particularly in patients with complex activation patterns.

Characterizing the arrhythmogenic substrate in personalized models of atrial fibrillation: sensitivity to mesh resolution and pacing protocol in AF models

AIMS

Computationally guided persistent atrial fibrillation (PsAF) ablation has emerged as an alternative to conventional treatment planning. To make this approach scalable, computational cost and the time required to conduct simulations must be minimized while maintaining predictive accuracy. Here, we assess the sensitivity of the process to finite-element mesh resolution. We also compare methods for pacing site distribution used to evaluate inducibility arrhythmia sustained by re-entrant drivers (RDs).

METHODS AND RESULTS

Simulations were conducted in low- and high-resolution models (average edge lengths: 400/350 µm) reconstructed from PsAF patients' late gadolinium enhancement magnetic resonance imaging scans. Pacing was simulated from 80 sites to assess RD inducibility. When pacing from the same site led to different outcomes in low-/high-resolution models, we characterized divergence dynamics by analysing dissimilarity index over time. Pacing site selection schemes prioritizing even spatial distribution and proximity to fibrotic tissue were evaluated. There were no RD sites observed in low-resolution models but not high-resolution models, or vice versa. Dissimilarity index analysis suggested that differences in simulation outcome arising from differences in discretization were the result of isolated conduction block incidents in one model but not the other; this never led to RD sites unique to one mesh resolution. Pacing site selection based on fibrosis proximity led to the best observed trade-off between number of stimulation locations and predictive accuracy.

CONCLUSION

Simulations conducted in meshes with 400 µm average edge length and ∼40 pacing sites proximal to fibrosis are sufficient to reveal the most comprehensive possible list of RD sites, given feasibility constraints.

Computer Modeling of the Heart for ECG Interpretation—A Review

Computer modeling of the electrophysiology of the heart has undergone significant progress. A healthy heart can be modeled starting from the ion channels via the spread of a depolarization wave on a realistic geometry of the human heart up to the potentials on the body surface and the ECG. Research is advancing regarding modeling diseases of the heart. This article reviews progress in calculating and analyzing the corresponding electrocardiogram (ECG) from simulated depolarization and repolarization waves. First, we describe modeling of the P-wave, the QRS complex and the T-wave of a healthy heart. Then, both the modeling and the corresponding ECGs of several important diseases and arrhythmias are delineated: ischemia and infarction, ectopic beats and extrasystoles, ventricular tachycardia, bundle branch blocks, atrial tachycardia, flutter and fibrillation, genetic diseases and channelopathies, imbalance of electrolytes and drug-induced changes. Finally, we outline the potential impact of computer modeling on ECG interpretation. Computer modeling can contribute to a better comprehension of the relation between features in the ECG and the underlying cardiac condition and disease. It can pave the way for a quantitative analysis of the ECG and can support the cardiologist in identifying events or non-invasively localizing diseased areas. Finally, it can deliver very large databases of reliably labeled ECGs as training data for machine learning.

Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals With a Hybrid in silico and in vivo Dataset

In patients with atrial fibrillation, intracardiac electrogram signal amplitude is known to decrease with increased structural tissue remodeling, referred to as fibrosis. In addition to the isolation of the pulmonary veins, fibrotic sites are considered a suitable target for catheter ablation. However, it remains an open challenge to find fibrotic areas and to differentiate their density and transmurality. This study aims to identify the volume fraction and transmurality of fibrosis in the atrial substrate. Simulated cardiac electrograms, combined with a generalized model of clinical noise, reproduce clinically measured signals. Our hybrid dataset approach combines in silico and clinical electrograms to train a decision tree classifier to characterize the fibrotic atrial substrate. This approach captures different in vivo dynamics of the electrical propagation reflected on healthy electrogram morphology and synergistically combines it with synthetic fibrotic electrograms from in silico experiments. The machine learning algorithm was tested on five patients and compared against clinical voltage maps as a proof of concept, distinguishing non-fibrotic from fibrotic tissue and characterizing the patient's fibrotic tissue in terms of density and transmurality. The proposed approach can be used to overcome a single voltage cut-off value to identify fibrotic tissue and guide ablation targeting fibrotic areas.

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.

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.

Computational modeling identifies embolic stroke of undetermined source patients with potential arrhythmic substrate

Cardiac magnetic resonance imaging (MRI) has revealed fibrosis in embolic stroke of undetermined source (ESUS) patients comparable to levels seen in atrial fibrillation (AFib). We used computational modeling to understand the absence of arrhythmia in ESUS despite the presence of putatively pro-arrhythmic fibrosis. MRI-based atrial models were reconstructed for 45 ESUS and 45 AFib patients. The fibrotic substrate's arrhythmogenic capacity in each patient was assessed computationally. Reentrant drivers were induced in 24/45 (53%) ESUS and 22/45 (49%) AFib models. Inducible models had more fibrosis (16.7 ± 5.45%) than non-inducible models (11.07 ± 3.61%; p<0.0001); however, inducible subsets of ESUS and AFib models had similar fibrosis levels (p=0.90), meaning that the intrinsic pro-arrhythmic substrate properties of fibrosis in ESUS and AFib are indistinguishable. This suggests that some ESUS patients have latent pre-clinical fibrotic substrate that could be a future source of arrhythmogenicity. Thus, our work prompts the hypothesis that ESUS patients with fibrotic atria are spared from AFib due to an absence of arrhythmia triggers.

A numerical study on the effects of spatial and temporal discretization in cardiac electrophysiology

Millions of degrees of freedom are often required to accurately represent the electrophysiology of the myocardium due to the presence of discretization effects. This study seeks to explore the influence of temporal and spatial discretization on the simulation of cardiac electrophysiology in conjunction with changes in modeling choices. Several finite element analyses are performed to examine how discretization affects solution time, conduction velocity and electrical excitation. Discretization effects are considered along with changes in the electrophysiology model and solution approach. Two action potential models are considered: the Aliev-Panfilov model and the ten Tusscher-Noble-Noble-Panfilov model. The solution approaches consist of two time integration schemes and different treatments for solving the local system of ordinary differential equations. The efficiency and stability of the calculation approaches are demonstrated to be dependent on the action potential model. The dependency of the conduction velocity on the element size and time step is shown to be different for changes in material parameters. Finally, the discrepancies between the wave propagation in coarse and fine meshes are analyzed based on the temporal evolution of the transmembrane potential at a node and its neighboring Gauss points. Insight obtained from this study can be used to suggest new methods to improve the efficiency of simulations in cardiac electrophysiology.

Parameter variations in personalized electrophysiological models of human heart ventricles

The objectives of this study were to evaluate the accuracy of personalized numerical simulations of the electrical activity in human ventricles by comparing simulated electrocardiograms (ECGs) with real patients’ ECGs and analyzing the sensitivity of the model output to variations in the model parameters. We used standard 12-lead ECGs and up to 224 unipolar body-surface ECGs to record three patients with cardiac resynchronization therapy devices and three patients with focal ventricular tachycardia. Patient-tailored geometrical models of the ventricles, atria, large vessels, liver, and spine were created using computed tomography data. Ten cases of focal ventricular activation were simulated using the bidomain model and the TNNP 2006 cellular model. The population-based values of electrical conductivities and other model parameters were used for accuracy analysis, and their variations were used for sensitivity analysis. The mean correlation coefficient between the simulated and real ECGs varied significantly (from r = 0.29 to r = 0.86) among the simulated cases. A strong mean correlation (r > 0.7) was found in eight of the ten model cases. The accuracy of the ECG simulation varied widely in the same patient depending on the localization of the excitation origin. The sensitivity analysis revealed that variations in the anisotropy ratio, blood conductivity, and cellular apicobasal heterogeneity had the strongest influence on transmembrane potential, while variation in lung conductivity had the greatest influence on body-surface ECGs. Futhermore, the anisotropy ratio predominantly affected the latest activation time and repolarization time dispersion, while the cellular apicobasal heterogeneity mainly affected the dispersion of action potential duration, and variation in lung conductivity mainly led to changes in the amplitudes of ECGs and cardiac electrograms. We also found that the effects of certain parameter variations had specific regional patterns on the cardiac and body surfaces. These observations are useful for further developing personalized cardiac models.

Data integration for the numerical simulation of cardiac electrophysiology

The increasing availability of extensive and accurate clinical data is rapidly shaping cardiovascular care by improving the understanding of physiological and pathological mechanisms of the cardiovascular system and opening new frontiers in designing therapies and interventions. In this direction, mathematical and numerical models provide a complementary relevant tool, able not only to reproduce patient-specific clinical indicators but also to predict and explore unseen scenarios. With this goal, clinical data are processed and provided as inputs to the mathematical model, which quantitatively describes the physical processes that occur in the cardiac tissue. In this paper, the process of integration of clinical data and mathematical models is discussed. Some challenges and contributions in the field of cardiac electrophysiology are reported.

Precision medicine in human heart modeling : Perspectives, challenges, and opportunities

Precision medicine is a new frontier in healthcare that uses scientific methods to customize medical treatment to the individual genes, anatomy, physiology, and lifestyle of each person. In cardiovascular health, precision medicine has emerged as a promising paradigm to enable cost-effective solutions that improve quality of life and reduce mortality rates. However, the exact role in precision medicine for human heart modeling has not yet been fully explored. Here, we discuss the challenges and opportunities for personalized human heart simulations, from diagnosis to device design, treatment planning, and prognosis. With a view toward personalization, we map out the history of anatomic, physical, and constitutive human heart models throughout the past three decades. We illustrate recent human heart modeling in electrophysiology, cardiac mechanics, and fluid dynamics and highlight clinically relevant applications of these models for drug development, pacing lead failure, heart failure, ventricular assist devices, edge-to-edge repair, and annuloplasty. With a view toward translational medicine, we provide a clinical perspective on virtual imaging trials and a regulatory perspective on medical device innovation. We show that precision medicine in human heart modeling does not necessarily require a fully personalized, high-resolution whole heart model with an entire personalized medical history. Instead, we advocate for creating personalized models out of population-based libraries with geometric, biological, physical, and clinical information by morphing between clinical data and medical histories from cohorts of patients using machine learning. We anticipate that this perspective will shape the path toward introducing human heart simulations into precision medicine with the ultimate goals to facilitate clinical decision making, guide treatment planning, and accelerate device design.

Action potential propagation and block in a model of atrial tissue with myocyte-fibroblast coupling

The electrical coupling between myocytes and fibroblasts and the spacial distribution of fibroblasts within myocardial tissues are significant factors in triggering and sustaining cardiac arrhythmias, but their roles are poorly understood. This article describes both direct numerical simulations and an asymptotic theory of propagation and block of electrical excitation in a model of atrial tissue with myocyte-fibroblast coupling. In particular, three idealized fibroblast distributions are introduced: uniform distribution, fibroblast barrier and myocyte strait-all believed to be constituent blocks of realistic fibroblast distributions. Primary action potential biomarkers including conduction velocity, peak potential and triangulation index are estimated from direct simulations in all cases. Propagation block is found to occur at certain critical values of the parameters defining each idealized fibroblast distribution, and these critical values are accurately determined. An asymptotic theory proposed earlier is extended and applied to the case of a uniform fibroblast distribution. Biomarker values are obtained from hybrid analytical-numerical solutions of coupled fast-time and slow-time periodic boundary value problems and compare well to direct numerical simulations. The boundary of absolute refractoriness is determined solely by the fast-time problem and is found to depend on the values of the myocyte potential and on the slow inactivation variable of the sodium current ahead of the propagating pulse. In turn, these quantities are estimated from the slow-time problem using a regular perturbation expansion to find the steady state of the coupled myocyte-fibroblast kinetics. The asymptotic theory gives a simple analytical expression that captures with remarkable accuracy the block of propagation in the presence of fibroblasts.

The 'Digital Twin' to enable the vision of precision cardiology

Providing therapies tailored to each patient is the vision of precision medicine, enabled by the increasing ability to capture extensive data about individual patients. In this position paper, we argue that the second enabling pillar towards this vision is the increasing power of computers and algorithms to learn, reason, and build the 'digital twin' of a patient. Computational models are boosting the capacity to draw diagnosis and prognosis, and future treatments will be tailored not only to current health status and data, but also to an accurate projection of the pathways to restore health by model predictions. The early steps of the digital twin in the area of cardiovascular medicine are reviewed in this article, together with a discussion of the challenges and opportunities ahead. We emphasize the synergies between mechanistic and statistical models in accelerating cardiovascular research and enabling the vision of precision medicine.

On the Role of Ionic Modeling on the Signature of Cardiac Arrhythmias for Healthy and Diseased Hearts

Computational cardiology is rapidly becoming the gold standard for innovative medical treatments and device development. Despite a worldwide effort in mathematical and computational modeling research, the complexity and intrinsic multiscale nature of the heart still limit our predictability power raising the question of the optimal modeling choice for large-scale whole-heart numerical investigations. We propose an extended numerical analysis among two different electrophysiological modeling approaches: a simplified phenomenological one and a detailed biophysical one. To achieve this, we considered three-dimensional healthy and infarcted swine heart geometries. Heterogeneous electrophysiological properties, fine-tuned DT-MRI -based anisotropy features, and non-conductive ischemic regions were included in a custom-built finite element code. We provide a quantitative comparison of the electrical behaviors during steady pacing and sustained ventricular fibrillation for healthy and diseased cases analyzing cardiac arrhythmias dynamics. Action potential duration (APD) restitution distributions, vortex filament counting, and pseudo-electrocardiography (ECG) signals were numerically quantified, introducing a novel statistical description of restitution patterns and ventricular fibrillation sustainability. Computational cost and scalability associated with the two modeling choices suggests that ventricular fibrillation signatures are mainly controlled by anatomy and structural parameters, rather than by regional restitution properties. Finally, we discuss limitations and translational perspectives of the different modeling approaches in view of large-scale whole-heart in silico studies.

Sensitivity analysis of an electrophysiology model for the left ventricle

Modelling the cardiac electrophysiology entails dealing with the uncertainties related to the input parameters such as the heart geometry and the electrical conductivities of the tissues, thus calling for an uncertainty quantification (UQ) of the results. Since the chambers of the heart have different shapes and tissues, in order to make the problem affordable, here we focus on the left ventricle with the aim of identifying which of the uncertain inputs mostly affect its electrophysiology. In a first phase, the uncertainty of the input parameters is evaluated using data available from the literature and the output quantities of interest (QoIs) of the problem are defined. According to the polynomial chaos expansion, a training dataset is then created by sampling the parameter space using a quasi-Monte Carlo method whereas a smaller independent dataset is used for the validation of the resulting metamodel. The latter is exploited to run a global sensitivity analysis with nonlinear variance-based indices and thus reduce the input parameter space accordingly. Thereafter, the uncertainty probability distribution of the QoIs are evaluated using a direct UQ strategy on a larger dataset and the results discussed in the light of the medical knowledge.

Using cardiac ionic cell models to interpret clinical data

For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

Stabilized hybrid discontinuous Galerkin finite element method for the cardiac monodomain equation

Numerical methods for solving the cardiac electrophysiology model, which describes the electrical activity in the heart, are proposed. The model problem consists of a nonlinear reaction-diffusion partial differential equation coupled to systems of ordinary differential equations that describes electrochemical reactions in cardiac cells. The proposed methods combine an operator splitting technique for the reaction-diffusion equation with primal hybrid methods for spatial discretization considering continuous or discontinuous approximations for the Lagrange multiplier. A static condensation is adopted to form a reduced global system in terms of the multiplier only. Convergence studies exhibit optimal rates of convergence and numerical experiments show that the proposed schemes can be more efficient than standard numerical techniques commonly used in this context when preconditioned iterative methods are used for the solution of linear systems.

Electrophysiological significance of the interatrial conduction including cavo-tricuspid isthmus during atrial fibrillation

KEY POINTS

The interatrial conduction, including Bachmann's bundle, the posterior septal conduction, the anterior septal conduction, and the cavo-tricuspid isthmus, contributes to the maintenance mechanisms of atrial fibrillation in a 3D biatrial model. The interatrial conduction ablation including a cavo-tricuspid isthmus ablation significantly affects the wave dynamics of atrial fibrillation (AF) and facilitates the AF termination or atrial tachycardia conversion of the AF after the circumferential pulmonary vein isolation. Additional cavo-tricuspid isthmus ablation after the circumferential pulmonary vein isolation improves long-term rhythm outcome after clinical AF catheter ablation.

ABSTRACT

Although it is known that atrial fibrillation (AF) is mainly a left atrial (LA) disease, the role of the right atrium (RA) and interatrial conduction (IAC), including the cavo-tricuspid isthmus (CTI), has not been clearly defined. We tested AF wave dynamics with or without IAC in computational modelling and the rhythm outcome of AF catheter ablation (AFCA) including CTI ablation in clinical cohort data. We evaluated the dominant frequency (DF) in 3D biatrial AF simulations integrated with 3D-computed tomograms obtained from 10 patients. The IAC was implemented at Bachmann's bundle, posterior septum and the CTI. After virtual circumferential PV isolation (CPVI), we disconnected IACs one by one, and observed the wave dynamics. We compared the long-term rhythm outcome after CPVI alone and additional CTI ablation in 846 patients with AFCA. LA-DF was higher than RA-DF in AF (P < 0.001). After CPVI, the DF decreased significantly by additional IAC ablation (P = 0.003), especially in the LA (P = 0.016). The amount of DF reduction (P = 0.020) and rates of AF termination (P < 0.001) or AT conversion (P = 0.021) were significantly higher after IAC ablations including CTI than those without. In clinical AFCA, the AF recurrence rate was significantly lower in patients with additional CTI ablation than CPVI alone during 25 ± 20 months' follow-up (hazard ratio 0.60 [0.46-0.79], P < 0.001, Log rank P < 0.001). IAC contributes to the maintenance mechanism of AF, and IAC including CTI ablation affects AF wave dynamics, facilitating AF termination in 3D biatrial modelling. Additional CTI ablation after CPVI improves the long-term rhythm outcome in clinical AFCA, potentially in a paroxysmal type with accompanying atrial flutter, or atrial dimension close to normal.