Computational models in cardiology

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

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

Identifying locations susceptible to micro-anatomical reentry using a spatial network representation of atrial fibre maps

Micro-anatomical reentry has been identified as a potential driver of atrial fibrillation (AF). In this paper, we introduce a novel computational method which aims to identify which atrial regions are most susceptible to micro-reentry. The approach, which considers the structural basis for micro-reentry only, is based on the premise that the accumulation of electrically insulating interstitial fibrosis can be modelled by simulating percolation-like phenomena on spatial networks. Our results suggest that at high coupling, where micro-reentry is rare, the micro-reentrant substrate is highly clustered in areas where the atrial walls are thin and have convex wall morphology, likely facilitating localised treatment via ablation. However, as transverse connections between fibres are removed, mimicking the accumulation of interstitial fibrosis, the substrate becomes less spatially clustered, and the bias to forming in thin, convex regions of the atria is reduced, possibly restricting the efficacy of localised ablation. Comparing our algorithm on image-based models with and without atrial fibre structure, we find that strong longitudinal fibre coupling can suppress the micro-reentrant substrate, whereas regions with disordered fibre orientations have an enhanced risk of micro-reentry. With further development, these methods may be useful for modelling the temporal development of the fibrotic substrate on an individualised basis.

A Novel In Silico Electromechanical Model of Human Ventricular Cardiomyocyte

Contractility has become one of the main readouts in computational and experimental studies on cardiomyocytes. Following this trend, we propose a novel mathematical model of human ventricular cardiomyocytes electromechanics, BPSLand, by coupling a recent human contractile element to the BPS2020 model of electrophysiology. BPSLand is the result of a hybrid optimization process and it reproduces all the electrophysiology experimental indices captured by its predecessor BPS2020, simultaneously enabling the simulation of realistic human active tension and its potential abnormalities. The transmural heterogeneity in both electrophysiology and contractility departments was simulated consistent with previous computational and in vitro studies. Furthermore, our model could capture delayed afterdepolarizations (DADs), early afterdepolarizations (EADs), and contraction abnormalities in terms of aftercontractions triggered by either drug action or special pacing modes. Finally, we further validated the mechanical results of the model against previous experimental and in silico studies, e.g., the contractility dependence on pacing rate. Adding a new level of applicability to the normative models of human cardiomyocytes, BPSLand represents a robust, fully-human in silico model with promising capabilities for translational cardiology.

The Cell Physiome: What Do We Need in a Computational Physiology Framework for Predicting Single-Cell Biology?

Modern biology and biomedicine are undergoing a big data explosion, needing advanced computational algorithms to extract mechanistic insights on the physiological state of living cells. We present the motivation for the Cell Physiome project: a framework and approach for creating, sharing, and using biophysics-based computational models of single-cell physiology. Using examples in calcium signaling, bioenergetics, and endosomal trafficking, we highlight the need for spatially detailed, biophysics-based computational models to uncover new mechanisms underlying cell biology. We review progress and challenges to date toward creating cell physiome models. We then introduce bond graphs as an efficient way to create cell physiome models that integrate chemical, mechanical, electromagnetic, and thermal processes while maintaining mass and energy balance. Bond graphs enhance modularization and reusability of computational models of cells at scale. We conclude with a look forward at steps that will help fully realize this exciting new field of mechanistic biomedical data science. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

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.

Pacemaking function of two simplified cell models

Simplified nonlinear models of biological cells are widely used in computational electrophysiology. The models reproduce qualitatively many of the characteristics of various organs, such as the heart, brain, and intestine. In contrast to complex cellular ion-channel models, the simplified models usually contain a small number of variables and parameters, which facilitates nonlinear analysis and reduces computational load. In this paper, we consider pacemaking variants of the Aliev-Panfilov and Corrado two-variable excitable cell models. We conducted a numerical simulation study of these models and investigated the main nonlinear dynamic features of both isolated cells and 1D coupled pacemaker-excitable systems. Simulations of the 2D sinoatrial node and 3D intestine tissue as application examples of combined pacemaker-excitable systems demonstrated results similar to obtained previously. The uniform formulation for the conventional excitable cell models and proposed pacemaker models allows a convenient and easy implementation for the construction of personalized physiological models, inverse tissue modeling, and development of real-time simulation systems for various organs that contain both pacemaker and excitable cells.

Modelling the interaction between stem cells derived cardiomyocytes patches and host myocardium to aid non-arrhythmic engineered heart tissue design

Application of epicardial patches constructed from human-induced pluripotent stem cell- derived cardiomyocytes (hiPSC-CMs) has been proposed as a long-term therapy to treat scarred hearts post myocardial infarction (MI). Understanding electrical interaction between engineered heart tissue patches (EHT) and host myocardium represents a key step toward a successful patch engraftment. EHT retain different electrical properties with respect to the host heart tissue due to the hiPSC-CMs immature phenotype, which may lead to increased arrhythmia risk. We developed a modelling framework to examine the influence of patch design on electrical activation at the engraftment site. We performed an in silico investigation of different patch design approaches to restore pre-MI activation properties and evaluated the associated arrhythmic risk. We developed an in silico cardiac electrophysiology model of a transmural cross section of host myocardium. The model featured an infarct region, an epicardial patch spanning the infarct region and a bath region. The patch is modelled as a layer of hiPSC-CM, combined with a layer of conductive polymer (CP). Tissue and patch geometrical dimensions and conductivities were incorporated through 10 modifiable model parameters. We validated our model against 4 independent experimental studies and showed that it can qualitatively reproduce their findings. We performed a global sensitivity analysis (GSA) to isolate the most important parameters, showing that the stimulus propagation is mainly governed by the scar depth, radius and conductivity when the scar is not transmural, and by the EHT patch conductivity when the scar is transmural. We assessed the relevance of small animal studies to humans by comparing simulations of rat, rabbit and human myocardium. We found that stimulus propagation paths and GSA sensitivity indices are consistent across species. We explored which EHT design variables have the potential to restore physiological propagation. Simulations predict that increasing EHT conductivity from 0.28 to 1-1.1 S/m recovered physiological activation in rat, rabbit and human. Finally, we assessed arrhythmia risk related to increasing EHT conductivity and tested increasing the EHT Na+ channel density as an alternative strategy to match healthy activation. Our results revealed a greater arrhythmia risk linked to increased EHT conductivity compared to increased Na+ channel density. We demonstrated that our modeling framework could capture the interaction between host and EHT patches observed in in vitro experiments. We showed that large (patch and tissue dimensions) and small (cardiac myocyte electrophysiology) scale differences between small animals and humans do not alter EHT patch effect on infarcted tissue. Our model revealed that only when the scar is transmural do EHT properties impact activation times and isolated the EHT conductivity as the main parameter influencing propagation. We predicted that restoring physiological activation by tuning EHT conductivity is possible but may promote arrhythmic behavior. Finally, our model suggests that acting on hiPSC-CMs low action potential upstroke velocity and lack of IK1 may restore pre-MI activation while not promoting arrhythmia.

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.

Weighted Time Warping Improves T-wave Morphology Markers Clinical Significance

Background: T-wave (TW) morphology indices based on time-warping (d_w) have shown significant cardiovascular risk stratification value. However, errors in the location of TW boundaries may impact their prognostic power. Our aim was to test the hypothesis that a weighted time-warping function (WF) would reduce the sensitivity of d_w to these errors and improve their clinical significance. Methods: The WFs were proportional to (i) the reference TW (T), and (ii) the absolute value of its derivative (D). The index d_w was recalculated using these WFs, and its performance was compared to the unweighted control case (C) in four different scenarios: 1) robustness against simulated TW boundaries location errors; 2) ability to retain physiological information in an electrophysiological cardiac model; 3) ability to monitor blood potassium concentration changes ([K⁺]) in 29 hemodialysis (HD) patients; 4) and the sudden cardiac death (SCD) risk stratification value of the TW morphology restitution (TMR) index, derived from d_w, in 651 chronic heart failure (CHF) patients. Results and Discussion: The WFs led to a reduced sensitivity (R) of d_w to TW boundary location errors as compared to C (median R=0.19 and 0.22 and 0.35 for T, D and C, respectively). They also preserved the physiological relationship between d_w and repolarization dispersion changes at ventricular level. No improvements in [K⁺] tracking were observed for the HD patients (Pearsons median correlation [r] between [K⁺] and d_w was 0.86r0.90 for T, D and C). In CHF patients, the SCD risk stratification value of TMR was improved by applying T (hazard ratio, HAR, of 2.80), followed by D (HAR=2.32) and C (HAR=2.23). Conclusions and Significance: The proposed WFs, with T showing the best performance, increased the robustness of time-warping based markers against TW location errors preserving their physiological information content and boosting their SCD risk stratification value. Results from this work support the use of T when deriving d_w for future clinical applications.

Functional and structural differences between skinned and intact muscle preparations

Myofilaments and their associated proteins, which together constitute the sarcomeres, provide the molecular-level basis for contractile function in all muscle types. In intact muscle, sarcomere-level contraction is strongly coupled to other cellular subsystems, in particular the sarcolemmal membrane. Skinned muscle preparations (where the sarcolemma has been removed or permeabilized) are an experimental system designed to probe contractile mechanisms independently of the sarcolemma. Over the last few decades, experiments performed using permeabilized preparations have been invaluable for clarifying the understanding of contractile mechanisms in both skeletal and cardiac muscle. Today, the technique is increasingly harnessed for preclinical and/or pharmacological studies that seek to understand how interventions will impact intact muscle contraction. In this context, intrinsic functional and structural differences between skinned and intact muscle pose a major interpretational challenge. This review first surveys measurements that highlight these differences in terms of the sarcomere structure, passive and active tension generation, and calcium dependence. We then highlight the main practical challenges and caveats faced by experimentalists seeking to emulate the physiological conditions of intact muscle. Gaining an awareness of these complexities is essential for putting experiments in due perspective.

Disease Progression of Hypertrophic Cardiomyopathy: Modeling Using Machine Learning

BACKGROUND

Cardiovascular disorders in general are responsible for 30% of deaths worldwide. Among them, hypertrophic cardiomyopathy (HCM) is a genetic cardiac disease that is present in about 1 of 500 young adults and can cause sudden cardiac death (SCD).

OBJECTIVE

Although the current state-of-the-art methods model the risk of SCD for patients, to the best of our knowledge, no methods are available for modeling the patient's clinical status up to 10 years ahead. In this paper, we propose a novel machine learning (ML)-based tool for predicting disease progression for patients diagnosed with HCM in terms of adverse remodeling of the heart during a 10-year period.

METHODS

The method consisted of 6 predictive regression models that independently predict future values of 6 clinical characteristics: left atrial size, left atrial volume, left ventricular ejection fraction, New York Heart Association functional classification, left ventricular internal diastolic diameter, and left ventricular internal systolic diameter. We supplemented each prediction with the explanation that is generated using the Shapley additive explanation method.

RESULTS

The final experiments showed that predictive error is lower on 5 of the 6 constructed models in comparison to experts (on average, by 0.34) or a consortium of experts (on average, by 0.22). The experiments revealed that semisupervised learning and the artificial data from virtual patients help improve predictive accuracies. The best-performing random forest model improved R² from 0.3 to 0.6.

CONCLUSIONS

By engaging medical experts to provide interpretation and validation of the results, we determined the models' favorable performance compared to the performance of experts for 5 of 6 targets.

Determining anatomical and electrophysiological detail requirements for computational ventricular models of porcine myocardial infarction

BACKGROUND

Computational models of the heart built from cardiac MRI and electrophysiology (EP) data have shown promise for predicting the risk of and ablation targets for myocardial infarction (MI) related ventricular tachycardia (VT), as well as to predict paced activation sequences in heart failure patients. However, most recent studies have relied on low resolution imaging data and little or no EP personalisation, which may affect the accuracy of model-based predictions.

OBJECTIVE

To investigate the impact of model anatomy, MI scar morphology, and EP personalisation strategies on paced activation sequences and VT inducibility to determine the level of detail required to make accurate model-based predictions.

METHODS

Imaging and EP data were acquired from a cohort of six pigs with experimentally induced MI. Computational models of ventricular anatomy, incorporating MI scar, were constructed including bi-ventricular or left ventricular (LV) only anatomy, and MI scar morphology with varying detail. Tissue conductivities and action potential duration (APD) were fitted to 12-lead ECG data using the QRS duration and the QT interval, respectively, in addition to corresponding literature parameters. Paced activation sequences and VT induction were simulated. Simulated paced activation and VT inducibility were compared between models and against experimental data.

RESULTS

Simulations predict that the level of model anatomical detail has little effect on simulated paced activation, with all model predictions comparing closely with invasive EP measurements. However, detailed scar morphology from high-resolution images, bi-ventricular anatomy, and personalized tissue conductivities are required to predict experimental VT outcome.

CONCLUSION

This study provides clear guidance for model generation based on clinical data. While a representing high level of anatomical and scar detail will require high-resolution image acquisition, EP personalisation based on 12-lead ECG can be readily incorporated into modelling pipelines, as such data is widely available.

Data-driven identification of dynamical models using adaptive parameter sets

This paper presents two data-driven model identification techniques for dynamical systems with fixed point attractors. Both strategies implement adaptive parameter update rules to limit truncation errors in the inferred dynamical models. The first strategy can be considered an extension of the dynamic mode decomposition with control (DMDc) algorithm. The second strategy uses a reduced order isostable coordinate basis that captures the behavior of the slowest decaying modes of the Koopman operator. The accuracy and robustness of both model identification algorithms is considered in a simple model with dynamics near a Hopf bifurcation. A more complicated model for nonlinear convective flow past an obstacle is also considered. In these examples, the proposed strategies outperform a collection of other commonly used data-driven model identification algorithms including Koopman model predictive control, Galerkin projection, and DMDc.

Synthetic PPG signal Generation to Improve Coronary Artery Disease Classification: Study with Physical Model of Cardiovascular System

This paper presents a novel approach of generating synthetic Photoplethysmogram (PPG) data using a physical model of the cardiovascular system to improve classifier performance with a combination of synthetic and real data. The physical model is an in-silico cardiac computational model, consisting of a four-chambered heart with electrophysiology, hemodynamic, and blood pressure autoregulation functionality. Starting with a small number of measured PPG data, the cardiac model is used to synthesize healthy as well as PPG time-series pertaining to coronary artery disease (CAD) by varying pathophysiological parameters. A Variational Autoencoder (VAE) structure is proposed to derive a statistical feature space for CAD classification. Results are presented in two perspectives namely, (i) using artificially reduced real disease data and (ii) using all the real disease data. In both cases, by augmenting with the synthetic data for training, the performance (sensitivity, specificity) of the classifier changes from (i) (0.65, 1) to (1, 0.9) and (ii) (1, 0.95) to (1,1). The proposed hybrid approach of combining physical modelling and statistical feature space selection generates realistic PPG data with pathophysiological interpretation and can outperform a baseline Generative Adversarial Network (GAN) architecture with a relatively small amount of real data for training. This proposed method could aid as a substitution technique for handling the problem of bulk data required for training machine learning algorithms for cardiac health-care applications.

Review of applications and challenges of quantitative systems pharmacology modeling and machine learning for heart failure

Quantitative systems pharmacology (QSP) is an important approach in pharmaceutical research and development that facilitates in silico generation of quantitative mechanistic hypotheses and enables in silico trials. As demonstrated by applications from numerous industry groups and interest from regulatory authorities, QSP is becoming an increasingly critical component in clinical drug development. With rapidly evolving computational tools and methods, QSP modeling has achieved important progress in pharmaceutical research and development, including for heart failure (HF). However, various challenges exist in the QSP modeling and clinical characterization of HF. Machine/deep learning (ML/DL) methods have had success in a wide variety of fields and disciplines. They provide data-driven approaches in HF diagnosis and modeling, and offer a novel strategy to inform QSP model development and calibration. The combination of ML/DL and QSP modeling becomes an emergent direction in the understanding of HF and clinical development new therapies. In this work, we review the current status and achievement in QSP and ML/DL for HF, and discuss remaining challenges and future perspectives in the field.

A computational analysis of atrial fibrillation effects on coronary perfusion across the different myocardial layers

Patients with atrial fibrillation (AF) may present ischemic chest pain in the absence of classical obstructive coronary disease. Among the possible causes, the direct hemodynamic effect exerted by the irregular arrhythmia has not been studied in detail. We performed a computational fluid dynamics analysis by means of a 1D-0D multiscale model of the entire human cardiovascular system, enriched by a detailed mathematical modeling of the coronary arteries and their downstream distal microcirculatory districts (subepicardial, midwall and subendocardial layers). Three mean ventricular rates were simulated (75, 100, 125 bpm) in both sinus rhythm (SR) and atrial fibrillation, and an inter-layer and inter-frequency analysis was conducted focusing on the ratio between mean beat-to-beat blood flow in AF compared to SR. Our results show that AF exerts direct hemodynamic consequences on the coronary microcirculation, causing a reduction in microvascular coronary flow particularly at higher ventricular rates; the most prominent reduction was seen in the subendocardial layers perfused by left coronary arteries (left anterior descending and left circumflex arteries).

A guide for assessment of myocardial stiffness in health and disease

Myocardial stiffness is an intrinsic property of the myocardium that influences both diastolic and systolic cardiac function. Myocardial stiffness represents the resistance of this tissue to being deformed and depends on intracellular components of the cardiomyocyte, particularly the cytoskeleton, and on extracellular components, such as collagen fibers. Myocardial disease is associated with changes in myocardial stiffness, and its assessment is a key diagnostic marker of acute or chronic pathological myocardial disease with the potential to guide therapeutic decision-making. In this Review, we appraise the different techniques that can be used to estimate myocardial stiffness, evaluate their advantages and disadvantages, and discuss potential clinical applications.

Molecular modeling in cardiovascular pharmacology: Current state of the art and perspectives

Molecular modeling in pharmacology is a promising emerging tool for exploring drug interactions with cellular components. Recent advances in molecular simulations, big data analysis, and artificial intelligence (AI) have opened new opportunities for rationalizing drug interactions with their pharmacological targets. Despite the obvious utility and increasing impact of computational approaches, their development is not progressing at the same speed in different fields of pharmacology. Here, we review current in silico techniques used in cardiovascular diseases (CVDs), cardiological drug discovery, and assessment of cardiotoxicity. In silico techniques are paving the way to a new era in cardiovascular medicine, but their use somewhat lags behind that in other fields.

Mechanical communication in fibrosis progression

Mechanical hallmarks of fibrotic microenvironments are both outcomes and causes of fibrosis progression. Understanding how cells sense and transmit mechanical cues in the interplay with extracellular matrix (ECM) and hemodynamic forces is a significant challenge. Recent advances highlight the evolvement of intracellular mechanotransduction pathways responding to ECM remodeling and abnormal hemodynamics (i.e., low and disturbed shear stress, pathological stretch, and increased pressure), which are prevalent biomechanical characteristics of fibrosis in multiple organs (e.g., liver, lung, and heart). Here, we envisage the mechanical communication in cell-ECM, cell-hemodynamics and cell-ECM-cell crosstalk (namely paratensile signaling) during fibrosis expansion. We also provide a comprehensive overview of in vitro and in silico engineering systems for disease modeling that will aid the identification and prediction of mechano-based therapeutic targets to ameliorate fibrosis progression.

What keeps us ticking? Sinoatrial node mechano-sensitivity: the grandfather clock of cardiac rhythm

The rhythmic and spontaneously generated electrical excitation that triggers the heartbeat originates in the sinoatrial node (SAN). SAN automaticity has been thoroughly investigated, which has uncovered fundamental mechanisms involved in cardiac pacemaking that are generally categorised into two interacting and overlapping systems: the 'membrane' and 'Ca²⁺ clock'. The principal focus of research has been on these two systems of oscillators, which have been studied primarily in single cells and isolated tissue, experimental preparations that do not consider mechanical factors present in the whole heart. SAN mechano-sensitivity has long been known to be a contributor to SAN pacemaking-both as a driver and regulator of automaticity-but its essential nature has been underappreciated. In this review, following a description of the traditional 'clocks' of SAN automaticity, we describe mechanisms of SAN mechano-sensitivity and its vital role for SAN function, making the argument that the 'mechanics oscillator' is, in fact, the 'grandfather clock' of cardiac rhythm.

Uncertainty Quantification of Regional Cardiac Tissue Properties in Arrhythmogenic Cardiomyopathy Using Adaptive Multiple Importance Sampling

Introduction: Computational models of the cardiovascular system are widely used to simulate cardiac (dys)function. Personalization of such models for patient-specific simulation of cardiac function remains challenging. Measurement uncertainty affects accuracy of parameter estimations. In this study, we present a methodology for patient-specific estimation and uncertainty quantification of parameters in the closed-loop CircAdapt model of the human heart and circulation using echocardiographic deformation imaging. Based on patient-specific estimated parameters we aim to reveal the mechanical substrate underlying deformation abnormalities in patients with arrhythmogenic cardiomyopathy (AC). Methods: We used adaptive multiple importance sampling to estimate the posterior distribution of regional myocardial tissue properties. This methodology is implemented in the CircAdapt cardiovascular modeling platform and applied to estimate active and passive tissue properties underlying regional deformation patterns, left ventricular volumes, and right ventricular diameter. First, we tested the accuracy of this method and its inter- and intraobserver variability using nine datasets obtained in AC patients. Second, we tested the trueness of the estimation using nine in silico generated virtual patient datasets representative for various stages of AC. Finally, we applied this method to two longitudinal series of echocardiograms of two pathogenic mutation carriers without established myocardial disease at baseline. Results: Tissue characteristics of virtual patients were accurately estimated with a highest density interval containing the true parameter value of 9% (95% CI [0-79]). Variances of estimated posterior distributions in patient data and virtual data were comparable, supporting the reliability of the patient estimations. Estimations were highly reproducible with an overlap in posterior distributions of 89.9% (95% CI [60.1-95.9]). Clinically measured deformation, ejection fraction, and end-diastolic volume were accurately simulated. In presence of worsening of deformation over time, estimated tissue properties also revealed functional deterioration. Conclusion: This method facilitates patient-specific simulation-based estimation of regional ventricular tissue properties from non-invasive imaging data, taking into account both measurement and model uncertainties. Two proof-of-principle case studies suggested that this cardiac digital twin technology enables quantitative monitoring of AC disease progression in early stages of disease.

Multiparameter optimization of nonuniform passive diffusion properties for creating coarse-grained equivalent models of cardiac propagation

The arrhythmogenic role of discrete cardiac propagation may be assessed by comparing discrete (fine-grained) and equivalent continuous (coarse-grained) models. We aim to develop an optimization algorithm for estimating the smooth conductivity field that best reproduces the diffusion properties of a given discrete model. Our algorithm iteratively adjusts local conductivity of the coarse-grained continuous model by simulating passive diffusion from white noise initial conditions during 3-10 ms and computing the root mean square error with respect to the discrete model. The coarse-grained conductivity field was interpolated from up to 300 evenly spaced control points. We derived an approximate formula for the gradient of the cost function that required (in two dimensions) only two additional simulations per iteration regardless of the number of estimated parameters. Conjugate gradient solver facilitated simultaneous optimization of multiple conductivity parameters. The method was tested in rectangular anisotropic tissues with uniform and nonuniform conductivity (slow regions with sinusoidal profile) and random diffuse fibrosis, as well as in a monolayer interconnected cable model of the left atrium with spatially-varying fibrosis density. Comparison of activation maps served as validation. The results showed that after convergence the errors in activation time were < 1 ms for rectangular geometries and 1-3 ms in the atrial model. Our approach based on the comparison of passive properties (<10 ms simulation) avoids performing active propagation simulations (>100 ms) at each iteration while reproducing activation maps, with possible applications to investigating the impact of microstructure on arrhythmias.

Spatial Changes in the Atrial Fibrillation Wave-Dynamics After Using Antiarrhythmic Drugs: A Computational Modeling Study

Background: We previously reported that a computational modeling-guided antiarrhythmic drug (AAD) test was feasible for evaluating multiple AADs in patients with atrial fibrillation (AF). We explored the anti-AF mechanisms of AADs and spatial change in the AF wave-dynamics by a realistic computational model. Methods: We used realistic computational modeling of 25 AF patients (68% male, 59.8 ± 9.8 years old, 32.0% paroxysmal AF) reflecting the anatomy, histology, and electrophysiology of the left atrium (LA) to characterize the effects of five AADs (amiodarone, sotalol, dronedarone, flecainide, and propafenone). We evaluated the spatial change in the AF wave-dynamics by measuring the mean dominant frequency (DF) and its coefficient of variation [dominant frequency-coefficient of variation (DF-COV)] in 10 segments of the LA. The mean DF and DF-COV were compared according to the pulmonary vein (PV) vs. extra-PV, maximal slope of the restitution curves (Smax), and defragmentation of AF. Results: The mean DF decreased after the administration of AADs in the dose dependent manner (p < 0.001). Under AADs, the DF was significantly lower (p < 0.001) and COV-DF higher (p = 0.003) in the PV than extra-PV region. The mean DF was significantly lower at a high Smax (≥1.4) than a lower Smax condition under AADs. During the episodes of AF defragmentation, the mean DF was lower (p < 0.001), but the COV-DF was higher (p < 0.001) than that in those without defragmentation. Conclusions: The DF reduction with AADs is predominant in the PVs and during a high Smax condition and causes AF termination or defragmentation during a lower DF and spatially unstable (higher DF-COV) condition.

Artificial intelligence in the diagnosis and management of arrhythmias

The field of cardiac electrophysiology (EP) had adopted simple artificial intelligence (AI) methodologies for decades. Recent renewed interest in deep learning techniques has opened new frontiers in electrocardiography analysis including signature identification of diseased states. Artificial intelligence advances coupled with simultaneous rapid growth in computational power, sensor technology, and availability of web-based platforms have seen the rapid growth of AI-aided applications and big data research. Changing lifestyles with an expansion of the concept of internet of things and advancements in telecommunication technology have opened doors to population-based detection of atrial fibrillation in ways, which were previously unimaginable. Artificial intelligence-aided advances in 3D cardiac imaging heralded the concept of virtual hearts and the simulation of cardiac arrhythmias. Robotics, completely non-invasive ablation therapy, and the concept of extended realities show promise to revolutionize the future of EP. In this review, we discuss the impact of AI and recent technological advances in all aspects of arrhythmia care.

Reconstruction of three-dimensional biventricular activation based on the 12-lead electrocardiogram via patient-specific modelling

Aims

Non-invasive imaging of electrical activation requires high-density body surface potential mapping. The nine electrodes of the 12-lead electrocardiogram (ECG) are insufficient for a reliable reconstruction with standard inverse methods. Patient-specific modelling may offer an alternative route to physiologically constraint the reconstruction. The aim of the study was to assess the feasibility of reconstructing the fully 3D electrical activation map of the ventricles from the 12-lead ECG and cardiovascular magnetic resonance (CMR).

Methods and results

Ventricular activation was estimated by iteratively optimizing the parameters (conduction velocity and sites of earliest activation) of a patient-specific model to fit the simulated to the recorded ECG. Chest and cardiac anatomy of 11 patients (QRS duration 126–180 ms, documented scar in two) were segmented from CMR images. Scar presence was assessed by magnetic resonance (MR) contrast enhancement. Activation sequences were modelled with a physiologically based propagation model and ECGs with lead field theory. Validation was performed by comparing reconstructed activation maps with those acquired by invasive electroanatomical mapping of coronary sinus/veins (CS) and right ventricular (RV) and left ventricular (LV) endocardium. The QRS complex was correctly reproduced by the model (Pearson’s correlation r = 0.923). Reconstructions accurately located the earliest and latest activated LV regions (median barycentre distance 8.2 mm, IQR 8.8 mm). Correlation of simulated with recorded activation time was very good at LV endocardium (r = 0.83) and good at CS (r = 0.68) and RV endocardium (r = 0.58).

Conclusion

Non-invasive assessment of biventricular 3D activation using the 12-lead ECG and MR imaging is feasible. Potential applications include patient-specific modelling and pre-/per-procedural evaluation of ventricular activation.

Atlas-based methods for efficient characterization of patient-specific ventricular activation patterns

AIMS

Ventricular activation patterns can aid clinical decision-making directly by providing spatial information on cardiac electrical activation or indirectly through derived clinical indices. The aim of this work was to derive an atlas of the major modes of variation of ventricular activation from model-predicted 3D bi-ventricular activation time distributions and to relate these modes to corresponding vectorcardiograms (VCGs). We investigated how the resulting dimensionality reduction can improve and accelerate the estimation of activation patterns from surface electrogram measurements.

METHODS AND RESULTS

Atlases of activation time (AT) and VCGs were derived using principal component analysis on a dataset of simulated electrophysiology simulations computed on eight patient-specific bi-ventricular geometries. The atlases provided significant dimensionality reduction, and the modes of variation in the two atlases described similar features. Utility of the atlases was assessed by resolving clinical waveforms against them and the VCG atlas was able to accurately reconstruct the patient VCGs with fewer than 10 modes. A sensitivity analysis between the two atlases was performed by calculating a compact Jacobian. Finally, VCGs generated by varying AT atlas modes were compared with clinical VCGs to estimate patient-specific activation maps, and the resulting errors between the clinical and atlas-based VCGs were less than those from more computationally expensive method.

CONCLUSION

Atlases of activation and VCGs represent a new method of identifying and relating the features of these high-dimensional signals that capture the major sources of variation between patients and may aid in identifying novel clinical indices of arrhythmia risk or therapeutic outcome.

Applications of artificial intelligence in cardiovascular imaging

Research into artificial intelligence (AI) has made tremendous progress over the past decade. In particular, the AI-powered analysis of images and signals has reached human-level performance in many applications owing to the efficiency of modern machine learning methods, in particular deep learning using convolutional neural networks. Research into the application of AI to medical imaging is now very active, especially in the field of cardiovascular imaging because of the challenges associated with acquiring and analysing images of this dynamic organ. In this Review, we discuss the clinical questions in cardiovascular imaging that AI can be used to address and the principal methodological AI approaches that have been developed to solve the related image analysis problems. Some approaches are purely data-driven and rely mainly on statistical associations, whereas others integrate anatomical and physiological information through additional statistical, geometric and biophysical models of the human heart. In a structured manner, we provide representative examples of each of these approaches, with particular attention to the underlying computational imaging challenges. Finally, we discuss the remaining limitations of AI approaches in cardiovascular imaging (such as generalizability and explainability) and how they can be overcome.

Causes of altered ventricular mechanics in hypertrophic cardiomyopathy: an in-silico study

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is typically caused by mutations in sarcomeric genes leading to cardiomyocyte disarray, replacement fibrosis, impaired contractility, and elevated filling pressures. These varying tissue properties are associated with certain strain patterns that may allow to establish a diagnosis by means of non-invasive imaging without the necessity of harmful myocardial biopsies or contrast agent application. With a numerical study, we aim to answer: how the variability in each of these mechanisms contributes to altered mechanics of the left ventricle (LV) and if the deformation obtained in in-silico experiments is comparable to values reported from clinical measurements.

METHODS

We conducted an in-silico sensitivity study on physiological and pathological mechanisms potentially underlying the clinical HCM phenotype. The deformation of the four-chamber heart models was simulated using a finite-element mechanical solver with a sliding boundary condition to mimic the tissue surrounding the heart. Furthermore, a closed-loop circulatory model delivered the pressure values acting on the endocardium. Deformation measures and mechanical behavior of the heart models were evaluated globally and regionally.

RESULTS

Hypertrophy of the LV affected the course of strain, strain rate, and wall thickening-the root-mean-squared difference of the wall thickening between control (mean thickness 10 mm) and hypertrophic geometries (17 mm) was >10%. A reduction of active force development by 40% led to less overall deformation: maximal radial strain reduced from 26 to 21%. A fivefold increase in tissue stiffness caused a more homogeneous distribution of the strain values among 17 heart segments. Fiber disarray led to minor changes in the circumferential and radial strain. A combination of pathological mechanisms led to reduced and slower deformation of the LV and halved the longitudinal shortening of the LA.

CONCLUSIONS

This study uses a computer model to determine the changes in LV deformation caused by pathological mechanisms that are presumed to underlay HCM. This knowledge can complement imaging-derived information to obtain a more accurate diagnosis of HCM.

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

Computationally efficient model of myocardial electromechanics for multiscale simulations

A model of myocardial electromechanics is suggested. It combines modified and simplified versions of previously published models of cardiac electrophysiology, excitation-contraction coupling, and mechanics. The mechano-calcium and mechano-electrical feedbacks, including the strain-dependence of the propagation velocity of the action potential, are also accounted for. The model reproduces changes in the twitch amplitude and Ca2+-transients upon changes in muscle strain including the slow response. The model also reproduces the Bowditch effect and changes in the twitch amplitude and duration upon changes in the interstimulus interval, including accelerated relaxation at high stimulation frequency. Special efforts were taken to reduce the stiffness of the differential equations of the model. As a result, the equations can be integrated numerically with a relatively high time step making the model suitable for multiscale simulation of the human heart and allowing one to study the impact of myocardial mechanics on arrhythmias.

Understanding PITX2-Dependent Atrial Fibrillation Mechanisms through Computational Models

Atrial fibrillation (AF) is a common arrhythmia. Better prevention and treatment of AF are needed to reduce AF-associated morbidity and mortality. Several major mechanisms cause AF in patients, including genetic predispositions to AF development. Genome-wide association studies have identified a number of genetic variants in association with AF populations, with the strongest hits clustering on chromosome 4q25, close to the gene for the homeobox transcription PITX2. Because of the inherent complexity of the human heart, experimental and basic research is insufficient for understanding the functional impacts of PITX2 variants on AF. Linking PITX2 properties to ion channels, cells, tissues, atriums and the whole heart, computational models provide a supplementary tool for achieving a quantitative understanding of the functional role of PITX2 in remodelling atrial structure and function to predispose to AF. It is hoped that computational approaches incorporating all we know about PITX2-related structural and electrical remodelling would provide better understanding into its proarrhythmic effects leading to development of improved anti-AF therapies. In the present review, we discuss advances in atrial modelling and focus on the mechanistic links between PITX2 and AF. Challenges in applying models for improving patient health are described, as well as a summary of future perspectives.

Atrial remodeling and atrial fibrillation recurrence after catheter ablation : Past, present, and future developments

The term "atrial remodeling" is used to describe the electrical, mechanical, and structural changes associated with the presence of an arrhythmogenic substrate for atrial fibrillation. Rhythm control therapy may slow down or even reverse progressive atrial remodeling. Atrial remodeling has long been recognized as an important predictor of clinical outcomes and therapeutic success, but recent advances have highlighted its clinical relevance and revealed the implications of specific anatomical changes such as atrial asymmetry or shape. This has opened the path to computational precision medicine that captures these data in detail and combines them with other factors, to provide patient-specific solutions. The goal of precision medicine lies in improving clinical outcomes, reducing costs, and avoiding unnecessary procedures. In this article, we review the history of atrial remodeling and we summarize the insights from our research on anatomical atrial remodeling and its association with rhythm outcomes after catheter ablation. Finally, we present recent advances in the field, reflecting the beginning of a new technological era that will enable us to improve patient care by personalized patient-specific medicine.

A space-fractional bidomain framework for cardiac electrophysiology: 1D alternans dynamics

Cardiac electrophysiology modeling deals with a complex network of excitable cells forming an intricate syncytium: the heart. The electrical activity of the heart shows recurrent spatial patterns of activation, known as cardiac alternans, featuring multiscale emerging behavior. On these grounds, we propose a novel mathematical formulation for cardiac electrophysiology modeling and simulation incorporating spatially non-local couplings within a physiological reaction-diffusion scenario. In particular, we formulate, a space-fractional electrophysiological framework, extending and generalizing similar works conducted for the monodomain model. We characterize one-dimensional excitation patterns by performing an extended numerical analysis encompassing a broad spectrum of space-fractional derivative powers and various intra- and extracellular conductivity combinations. Our numerical study demonstrates that (i) symmetric properties occur in the conductivity parameters' space following the proposed theoretical framework, (ii) the degree of non-local coupling affects the onset and evolution of discordant alternans dynamics, and (iii) the theoretical framework fully recovers classical formulations and is amenable for parametric tuning relying on experimental conduction velocity and action potential morphology.

Conductance heterogeneities induced by multistability in the dynamics of coupled cardiac gap junctions

In this paper, we study the propagation of the cardiac action potential in a one-dimensional fiber, where cells are electrically coupled through gap junctions (GJs). We consider gap junctional gate dynamics that depend on the intercellular potential. We find that different GJs in the tissue can end up in two different states: a low conducting state and a high conducting state. We first present evidence of the dynamical multistability that occurs by setting specific parameters of the GJ dynamics. Subsequently, we explain how the multistability is a direct consequence of the GJ stability problem by reducing the dynamical system's dimensions. The conductance dispersion usually occurs on a large time scale, i.e., thousands of heartbeats. The full cardiac model simulations are computationally demanding, and we derive a simplified model that allows for a reduction in the computational cost of four orders of magnitude. This simplified model reproduces nearly quantitatively the results provided by the original full model. We explain the discrepancies between the two models due to the simplified model's lack of spatial correlations. This simplified model provides a valuable tool to explore cardiac dynamics over very long time scales. That is highly relevant in studying diseases that develop on a large time scale compared to the basic heartbeat. As in the brain, plasticity and tissue remodeling are crucial parameters in determining the action potential wave propagation's stability.

Automated Localization of Focal Ventricular Tachycardia From Simulated Implanted Device Electrograms: A Combined Physics-AI Approach

Background: Focal ventricular tachycardia (VT) is a life-threating arrhythmia, responsible for high morbidity rates and sudden cardiac death (SCD). Radiofrequency ablation is the only curative therapy against incessant VT; however, its success is dependent on accurate localization of its source, which is highly invasive and time-consuming. Objective: The goal of our study is, as a proof of concept, to demonstrate the possibility of utilizing electrogram (EGM) recordings from cardiac implantable electronic devices (CIEDs). To achieve this, we utilize fast and accurate whole torso electrophysiological (EP) simulations in conjunction with convolutional neural networks (CNNs) to automate the localization of focal VTs using simulated EGMs. Materials and Methods: A highly detailed 3D torso model was used to simulate ∼4000 focal VTs, evenly distributed across the left ventricle (LV), utilizing a rapid reaction-eikonal environment. Solutions were subsequently combined with lead field computations on the torso to derive accurate electrocardiograms (ECGs) and EGM traces, which were used as inputs to CNNs to localize focal sources. We compared the localization performance of a previously developed CNN architecture (Cartesian probability-based) with our novel CNN algorithm utilizing universal ventricular coordinates (UVCs). Results: Implanted device EGMs successfully localized VT sources with localization error (8.74 mm) comparable to ECG-based localization (6.69 mm). Our novel UVC CNN architecture outperformed the existing Cartesian probability-based algorithm (errors = 4.06 mm and 8.07 mm for ECGs and EGMs, respectively). Overall, localization was relatively insensitive to noise and changes in body compositions; however, displacements in ECG electrodes and CIED leads caused performance to decrease (errors 16-25 mm). Conclusion: EGM recordings from implanted devices may be used to successfully, and robustly, localize focal VT sources, and aid ablation planning.

Inference of ventricular activation properties from non-invasive electrocardiography

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

•

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

•

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

•

Bayesian-inspired simultaneous estimation of continuous and discrete parameters.

•

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

•

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

In silico deceased donor intervention research: A potential accelerant for progress

Progress in deceased donor intervention research has been limited. Development of an in silico model of deceased donor physiology may elucidate potential therapeutic targets and provide an efficient mechanism for testing proposed deceased donor interventions. In this study, we report a preliminary in silico model of deceased kidney donor injury built, calibrated, and validated based on data from published animal and human studies. We demonstrate that the in silico model behaves like animal studies of brain death pathophysiology with respect to upstream markers of renal injury including hemodynamics, oxygenation, cytokines expression, and inflammation. Therapeutic hypothermia, a deceased donor intervention studied in human trials, is performed to demonstrate the model's ability to mimic an established clinical trial. Finally, future directions for developing this concept into a functional, clinically applicable model are discussed.

Percolation and tortuosity in heart-like cells

In the last several years, quite a few papers on the joint question of transport, tortuosity and percolation have appeared in the literature, dealing with passage of miscellaneous liquids or electrical currents in different media. However, these methods have not been applied to the passage of action potential in heart fibrosis (HF), which is crucial for problems of heart arrhythmia, especially of atrial tachycardia and fibrillation. In this work we address the HF problem from these aspects. A cellular automaton model is used to analyze percolation and transport of a distributed-fibrosis inflicted heart-like tissue. Although based on a rather simple mathematical model, it leads to several important outcomes: (1) It is shown that, for a single wave front (as the one emanated by the heart's sinus node), the percolation of heart-like matrices is exactly similar to the forest fire case. (2) It is shown that, on the average, the shape of the transport (a question not dealt with in relation to forest fire, and deals with the delay of action potential when passing a fibrotic tissue) behaves like a Gaussian. (3) Moreover, it is shown that close to the percolation threshold the parameters of this Gaussian behave in a critical way. From the physical point of view, these three results are an important contribution to the general percolation investigation. The relevance of our results to cardiological issues, specifically to the question of reentry initiation, are discussed and it is shown that: (A) Without an ectopic source and under a mere sinus node operation, no arrhythmia is generated, and (B) A sufficiently high refractory period could prevent some reentry mechanisms, even in partially fibrotic heart tissue.

Model Systems for Addressing Mechanism of Arrhythmogenesis in Cardiac Repair

PURPOSE OF REVIEW

Cardiac cell-based therapy represents a promising approach for cardiac repair. However, one of the main challenges is cardiac arrhythmias associated with stem cell transplantation. The current review summarizes the recent progress in model systems for addressing mechanisms of arrhythmogenesis in cardiac repair.

RECENT FINDINGS

Animal models have been extensively developed for mechanistic studies of cardiac arrhythmogenesis. Advances in human induced pluripotent stem cells (hiPSCs), patient-specific disease models, tissue engineering, and gene editing have greatly enhanced our ability to probe the mechanistic bases of cardiac arrhythmias. Additionally, recent development in multiscale computational studies and machine learning provides yet another powerful tool to quantitatively decipher the mechanisms of cardiac arrhythmias. Advancing efforts towards the integrations of experimental and computational studies are critical to gain insights into novel mitigation strategies for cardiac arrhythmias in cell-based therapy.

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.

Scaling digital twins from the artisanal to the industrial

Mathematical modeling and simulation are moving from being powerful development and analysis tools towards having increased roles in operational monitoring, control and decision support, in which models of specific entities are continually updated in the form of a digital twin. However, current digital twins are largely the result of bespoke technical solutions that are difficult to scale. We discuss two exemplar applications that motivate challenges and opportunities for scaling digital twins, and that underscore potential barriers to wider adoption of this technology.

Inflammation as a Risk Factor in Cardiotoxicity: An Important Consideration for Screening During Drug Development

Numerous commonly prescribed drugs, including antiarrhythmics, antihistamines, and antibiotics, carry a proarrhythmic risk and may induce dangerous arrhythmias, including the potentially fatal Torsades de Pointes. For this reason, cardiotoxicity testing has become essential in drug development and a required step in the approval of any medication for use in humans. Blockade of the hERG K⁺ channel and the consequent prolongation of the QT interval on the ECG have been considered the gold standard to predict the arrhythmogenic risk of drugs. In recent years, however, preclinical safety pharmacology has begun to adopt a more integrative approach that incorporates mathematical modeling and considers the effects of drugs on multiple ion channels. Despite these advances, early stage drug screening research only evaluates QT prolongation in experimental and computational models that represent healthy individuals. We suggest here that integrating disease modeling with cardiotoxicity testing can improve drug risk stratification by predicting how disease processes and additional comorbidities may influence the risks posed by specific drugs. In particular, chronic systemic inflammation, a condition associated with many diseases, affects heart function and can exacerbate medications’ cardiotoxic effects. We discuss emerging research implicating the role of inflammation in cardiac electrophysiology, and we offer a perspective on how in silico modeling of inflammation may lead to improved evaluation of the proarrhythmic risk of drugs at their early stage of development.

Transverse isotropic modelling of left-ventricle passive filling: Mechanical characterization for epicardial biomaterial manufacturing

Biomaterials applied to the epicardium have been studied intensively in recent years for different therapeutic purposes. Their mechanical influence on the heart, however, has not been clearly identified. Most biomaterials for epicardial applications are manufactured as membranes or cardiac patches that have isotropic geometry, which is not well suited to myocardial wall motion. Myocardial wall motion during systole and diastole produces a complex force in different directions. Membrane or cardiac patches that cannot adapt to these specific directions will exert an inappropriate force on the heart, at the risk of overly restricting or dilating it. Accurately characterizing the mechanical properties of the myocardial wall is thus essential, through analysis of muscle orientation and elasticity. In this study, we investigated the Hertz contact theory for characterizing cardiac tissue, using nanoindentation measurements to distinguish different patterns in the local myocardium. We then evaluated the predictive accuracy of this model using Finite Element Analysis (FEA) to mimic the diastolic phase of the heart. Our results, extracted from instrumented nanoindentation experiments in a liquid environment using five pig hearts, revealed variations in elasticity according to the local orientation of the myocardial tissue. In addition, applying the Finite Element Method (FEM) in our model based on transverse isotropy and local tissue orientation proved able to accurately simulate the passive filling of a left ventricle (LV) in a representative 3D geometry. Our model enables improved understanding of the underlying mechanical properties of the LV wall and can serve as a guide for designing and manufacturing biomedical material better adapted to the local epicardial tissue.

Polygonal surface processing and mesh generation tools for the numerical simulation of the cardiac function

In order to simulate the cardiac function for a patient-specific geometry, the generation of the computational mesh is crucially important. In practice, the input is typically a set of unprocessed polygonal surfaces coming either from a template geometry or from medical images. These surfaces need ad-hoc processing to be suitable for a volumetric mesh generation. In this work we propose a set of new algorithms and tools aiming to facilitate the mesh generation process. In particular, we focus on different aspects of a cardiac mesh generation pipeline: (1) specific polygonal surface processing for cardiac geometries, like connection of different heart chambers or segmentation outputs; (2) generation of accurate boundary tags; (3) definition of mesh-size functions dependent on relevant geometric quantities; (4) processing and connecting together several volumetric meshes. The new algorithms-implemented in the open-source software vmtk-can be combined with each other allowing the creation of personalized pipelines, that can be optimized for each cardiac geometry or for each aspect of the cardiac function to be modeled. Thanks to these features, the proposed tools can significantly speed-up the mesh generation process for a large range of cardiac applications, from single-chamber single-physics simulations to multi-chambers multi-physics simulations. We detail all the proposed algorithms motivating them in the cardiac context and we highlight their flexibility by showing different examples of cardiac mesh generation pipelines.