Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia

Simultaneous data assimilation and cardiac electrophysiology model correction using differentiable physics and deep learning

Modelling complex systems, like the human heart, has made great progress over the last decades. Patient-specific models, called 'digital twins', can aid in diagnosing arrhythmias and personalizing treatments. However, building highly accurate predictive heart models requires a delicate balance between mathematical complexity, parameterization from measurements and validation of predictions. Cardiac electrophysiology (EP) models range from complex biophysical models to simplified phenomenological models. Complex models are accurate but computationally intensive and challenging to parameterize, while simplified models are computationally efficient but less realistic. In this paper, we propose a hybrid approach by leveraging deep learning to complete a simplified cardiac model from data. Our novel framework has two components, decomposing the dynamics into a physics based and a data-driven term. This construction allows our framework to learn from data of different complexity, while simultaneously estimating model parameters. First, using in silico data, we demonstrate that this framework can reproduce the complex dynamics of cardiac transmembrane potential even in the presence of noise in the data. Second, using ex vivo optical data of action potentials (APs), we demonstrate that our framework can identify key physical parameters for anatomical zones with different electrical properties, as well as to reproduce the AP wave characteristics obtained from various pacing locations. Our physics-based data-driven approach may improve cardiac EP modelling by providing a robust biophysical tool for predictions.

Efficient Patient-Specific Simulations of Ventricular Tachycardia Based on Computed Tomography-Defined Wall Thickness Heterogeneity

BACKGROUND

Electrophysiological mapping of ventricular tachycardia (VT) is tedious and poorly reproducible. Substrate analysis on imaging cannot explicitly display VT circuits.

OBJECTIVES

This study sought to introduce a computed tomography-based model personalization approach, allowing for the simulation of postinfarction VT in a clinically compatible time frame.

METHODS

In 10 patients (age 65 ± 11 years, 9 male) referred for post-VT ablation, computed tomography-derived wall thickness maps were registered to 25 electroanatomical maps (sinus rhythm, paced, and VT). The relationship between wall thickness and electrophysiological characteristics (activation-recovery interval) was analyzed. Wall thickness was then employed to parameterize a fast and tractable organ-scale wave propagation model. Pacing protocols were simulated from multiple sites to test VT induction in silico. In silico VTs were compared to VT circuits mapped clinically.

RESULTS

Clinically, 6 different VTs could be induced with detailed maps in 9 patients. The proposed model allowed for fast simulation (median: 6 min/pacing site). Simulations of steady pacing (600 milliseconds) from 100 different sites/patient never triggered any arrhythmia. Applying S1-S2 or S1-S2-S3 induction schemes allowed for the induction of in silico VTs in the 9 of 10 patients who were clinically inducible. The patient who was not inducible clinically was also noninducible in silico. A total of 42 different VTs were simulated (4.2 ± 2 per patient). Six in silico VTs matched a VT circuit mapped clinically.

CONCLUSIONS

The proposed framework allows for personalized simulations in a matter of hours. In 6 of 9 patients, simulations show re-entrant patterns matching intracardiac recordings.

An automated near-real time computational method for induction and treatment of scar-related ventricular tachycardias

Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.

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.

Characterization of the Electrophysiologic Remodeling of Patients With Ischemic Cardiomyopathy by Clinical Measurements and Computer Simulations Coupled With Machine Learning

RATIONALE

Patients with ischemic cardiomyopathy (ICMP) are at high risk for malignant arrhythmias, largely due to electrophysiological remodeling of the non-infarcted myocardium. The electrophysiological properties of the non-infarcted myocardium of patients with ICMP remain largely unknown.

OBJECTIVES

To assess the pro-arrhythmic behavior of non-infarcted myocardium in ICMP patients and couple computational simulations with machine learning to establish a methodology for the development of disease-specific action potential models based on clinically measured action potential duration restitution (APDR) data.

METHODS AND RESULTS

We enrolled 22 patients undergoing left-sided ablation (10 ICMP) and compared APDRs between ICMP and structurally normal left ventricles (SNLVs). APDRs were clinically assessed with a decremental pacing protocol. Using genetic algorithms (GAs), we constructed populations of action potential models that incorporate the cohort-specific APDRs. The variability in the populations of ICMP and SNLV models was captured by clustering models based on their similarity using unsupervised machine learning. The pro-arrhythmic potential of ICMP and SNLV models was assessed in cell- and tissue-level simulations. Clinical measurements established that ICMP patients have a steeper APDR slope compared to SNLV (by 38%, p < 0.01). In cell-level simulations, APD alternans were induced in ICMP models at a longer cycle length compared to SNLV models (385-400 vs 355 ms). In tissue-level simulations, ICMP models were more susceptible for sustained functional re-entry compared to SNLV models.

CONCLUSION

Myocardial remodeling in ICMP patients is manifested as a steeper APDR compared to SNLV, which underlies the greater arrhythmogenic propensity in these patients, as demonstrated by cell- and tissue-level simulations using action potential models developed by GAs from clinical measurements. The methodology presented here captures the uncertainty inherent to GAs model development and provides a blueprint for use in future studies aimed at evaluating electrophysiological remodeling resulting from other cardiac diseases.

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.

Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning

We present a new, computationally efficient framework to perform forward uncertainty quantification (UQ) in cardiac electrophysiology. We consider the monodomain model to describe the electrical activity in the cardiac tissue, coupled with the Aliev-Panfilov model to characterize the ionic activity through the cell membrane. We address a complete forward UQ pipeline, including both: (i) a variance-based global sensitivity analysis for the selection of the most relevant input parameters, and (ii) a way to perform uncertainty propagation to investigate the impact of intra-subject variability on outputs of interest depending on the cardiac potential. Both tasks exploit stochastic sampling techniques, thus implying overwhelming computational costs because of the huge amount of queries to the high-fidelity, full-order computational model obtained by approximating the coupled monodomain/Aliev-Panfilov system through the finite element method. To mitigate this computational burden, we replace the full-order model with computationally inexpensive projection-based reduced-order models (ROMs) aimed at reducing the state-space dimensionality. Resulting approximation errors on the outputs of interest are finally taken into account through artificial neural network (ANN)-based models, enhancing the accuracy of the whole UQ pipeline. Numerical results show that the proposed physics-based ROMs outperform regression-based emulators relying on ANNs built with the same amount of training data, in terms of both numerical accuracy and overall computational efficiency.

A porous circulation model of the human brain for in silico clinical trials in ischaemic stroke

The advancement of ischaemic stroke treatment relies on resource-intensive experiments and clinical trials. In order to improve ischaemic stroke treatments, such as thrombolysis and thrombectomy, we target the development of computational tools for in silico trials which can partially replace these animal and human experiments with fast simulations. This study proposes a model that will serve as part of a predictive unit within an in silico clinical trial estimating patient outcome as a function of treatment. In particular, the present work aims at the development and evaluation of an organ-scale microcirculation model of the human brain for perfusion prediction. The model relies on a three-compartment porous continuum approach. Firstly, a fast and robust method is established to compute the anisotropic permeability tensors representing arterioles and venules. Secondly, vessel encoded arterial spin labelling magnetic resonance imaging and clustering are employed to create an anatomically accurate mapping between the microcirculation and large arteries by identifying superficial perfusion territories. Thirdly, the parameter space of the problem is reduced by analysing the governing equations and experimental data. Fourthly, a parameter optimization is conducted. Finally, simulations are performed with the tuned model to obtain perfusion maps corresponding to an open and an occluded (ischaemic stroke) scenario. The perfusion map in the occluded vessel scenario shows promising qualitative agreement with computed tomography images of a patient with ischaemic stroke caused by large vessel occlusion. The results highlight that in the case of vessel occlusion (i) identifying perfusion territories is essential to capture the location and extent of underperfused regions and (ii) anisotropic permeability tensors are required to give quantitatively realistic estimation of perfusion change. In the future, the model will be thoroughly validated against experiments.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

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.

Considering discrepancy when calibrating a mechanistic electrophysiology model

Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions-that of uncertainty in the model structure or the equations themselves. The difference between imperfect models and reality is termed model discrepancy, and we are often uncertain as to the size and consequences of this discrepancy. Here, we provide two examples of the consequences of discrepancy when calibrating models at the ion channel and action potential scales. Furthermore, we attempt to account for this discrepancy when calibrating and validating an ion channel model using different methods, based on modelling the discrepancy using Gaussian processes and autoregressive-moving-average models, then highlight the advantages and shortcomings of each approach. Finally, suggestions and lines of enquiry for future work are provided. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

An audit of uncertainty in multi-scale cardiac electrophysiology models

Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

How personalized heart modeling can help treatment of lethal arrhythmias: A focus on ventricular tachycardia ablation strategies in post-infarction patients

Precision Cardiology is a targeted strategy for cardiovascular disease prevention and treatment that accounts for individual variability. Computational heart modeling is one of the novel approaches that have been developed under the umbrella of Precision Cardiology. Personalized computational modeling of patient hearts has made strides in the development of models that incorporate the individual geometry and structure of the heart as well as other patient-specific information. Of these developments, one of the potentially most impactful is the research aimed at noninvasively predicting the targets of ablation of lethal arrhythmia, ventricular tachycardia (VT), using patient-specific models. The approach has been successfully applied to patients with ischemic cardiomyopathy in proof-of-concept studies. The goal of this paper is to review the strategies for computational VT ablation guidance in ischemic cardiomyopathy patients, from model developments to the intricacies of the actual clinical application. To provide context in describing the road these computational modeling applications have undertaken, we first review the state of the art in VT ablation in the clinic, emphasizing the benefits that personalized computational prediction of ablation targets could bring to the clinical electrophysiology practice. This article is characterized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Translational, Genomic, and Systems Medicine > Translational Medicine.

Fast personalized electrophysiological models from computed tomography images for ventricular tachycardia ablation planning

Aims

Clinical application of patient-specific cardiac computer models requires fast and robust processing pipelines that can be seamlessly integrated into clinical workflows. We aim at building such a pipeline from computed tomography (CT) images to personalized cardiac electrophysiology (EP) model. The simulation output could be useful in the context of post-infarct ventricular tachycardia (VT) radiofrequency ablation (RFA) planning for pre-operative targets prediction.

Methods and results

The support for model personalization is a patient-specific virtual three-dimensional heart obtained from CT images. Here, the scar is identified as thinning of the myocardial wall on automatically computed thickness maps. We then use an Eikonal model of wave front propagation with reduced velocity in the damaged areas. An image-based vessel enhancement algorithm can automatically identify VT isthmuses. The personalized model is used for virtual pacing. We obtained a very fast pipeline that enables simulations in only a few minutes. It is fully automated starting from the semi-automated image segmentation phase. The computational time frame is compatible with the construction of a virtual pacing tool. In this tool, onset points and an optional directional block could be interactively selected. The directional block is a simple way to model tissue refractoriness. Output activation maps are compared with EP data acquired pre-operatively. We show that this framework allows the reproduction of recorded re-entrant VT activation patterns.

Conclusion

Our simulation framework has an application in VT RFA intervention planning. It could be used to guide EP explorations and even predict ablation targets pre-operatively. This could reduce intervention duration and improve success rate.

Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia

In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.

Computational models in cardiology

The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.

Quantifying the uncertainty in model parameters using Gaussian process-based Markov chain Monte Carlo in cardiac electrophysiology

Model personalization requires the estimation of patient-specific tissue properties in the form of model parameters from indirect and sparse measurement data. Moreover, a low-dimensional representation of the parameter space is needed, which often has a limited ability to reveal the underlying tissue heterogeneity. As a result, significant uncertainty can be associated with the estimated values of the model parameters which, if left unquantified, will lead to unknown variability in model outputs that will hinder their reliable clinical adoption. Probabilistic estimation of model parameters, however, remains an unresolved challenge. Direct Markov Chain Monte Carlo (MCMC) sampling of the posterior distribution function (pdf) of the parameters is infeasible because it involves repeated evaluations of the computationally expensive simulation model. To accelerate this inference, one popular approach is to construct a computationally efficient surrogate and sample from this approximation. However, by sampling from an approximation, efficiency is gained at the expense of sampling accuracy. In this paper, we address this issue by integrating surrogate modeling of the posterior pdf into accelerating the Metropolis-Hastings (MH) sampling of the exact posterior pdf. It is achieved by two main components: (1) construction of a Gaussian process (GP) surrogate of the exact posterior pdf by actively selecting training points that allow for a good global approximation accuracy with a focus on the regions of high posterior probability; and (2) use of the GP surrogate to improve the proposal distribution in MH sampling, in order to improve the acceptance rate. The presented framework is evaluated in its estimation of the local tissue excitability of a cardiac electrophysiological model in both synthetic data experiments and real data experiments. In addition, the obtained posterior distributions of model parameters are interpreted in relation to the factors contributing to parameter uncertainty, including different low-dimensional representations of the parameter space, parameter non-identifiability, and parameter correlations.

Patient-Specific Cardiovascular Computational Modeling: Diversity of Personalization and Challenges

Patient-specific computer models have been developed representing a variety of aspects of the cardiovascular system spanning the disciplines of electrophysiology, electromechanics, solid mechanics, and fluid dynamics. These physiological mechanistic models predict macroscopic phenomena such as electrical impulse propagation and contraction throughout the entire heart as well as flow and pressure dynamics occurring in the ventricular chambers, aorta, and coronary arteries during each heartbeat. Such models have been used to study a variety of clinical scenarios including aortic aneurysms, coronary stenosis, cardiac valvular disease, left ventricular assist devices, cardiac resynchronization therapy, ablation therapy, and risk stratification. After decades of research, these models are beginning to be incorporated into clinical practice directly via marketed devices and indirectly by improving our understanding of the underlying mechanisms of health and disease within a clinical context.

Model-Based Feature Augmentation for Cardiac Ablation Target Learning From Images

GOAL

We present a model-based feature augmentation scheme to improve the performance of a learning algorithm for the detection of cardiac radio-frequency ablation (RFA) targets with respect to learning from images alone.

METHODS

Initially, we compute image features from delayed-enhanced magnetic resonance imaging (DE-MRI) to describe local tissue heterogeneities and feed them into a machine learning framework with uncertainty assessment for the identification of potential ablation targets. Next, we introduce the use of a patient-specific image-based model derived from DE-MRI coupled with the Mitchell-Schaeffer electrophysiology model and a dipole formulation for the simulation of intracardiac electrograms. Relevant features are extracted from these simulated signals which serve as a feature augmentation scheme for the learning algorithm. We assess the classifier's performance when using only image features and with model-based feature augmentation.

RESULTS

We obtained average classification scores of 97.2 % accuracy, 82.4 % sensitivity, and 95.0 % positive predictive value by using a model-based feature augmentation scheme. Preliminary results also show that training the algorithm on the closest patient from the database, instead of using all the patients, improves the classification results.

CONCLUSION

We presented a feature augmentation scheme based on biophysical cardiac electrophysiology modeling to increase the prediction scores of a machine learning framework for the RFA target prediction.

SIGNIFICANCE

The results derived from this study are a proof of concept that the use of model-based feature augmentation strengthens the performance of a purely image driven learning scheme for the prediction of cardiac ablation targets.

Using physiologically based models for clinical translation: predictive modelling, data interpretation or something in-between?

Heart disease continues to be a significant clinical problem in Western society. Predictive models and simulations that integrate physiological understanding with patient information derived from clinical data have huge potential to contribute to improving our understanding of both the progression and treatment of heart disease. In particular they provide the potential to improve patient selection and optimisation of cardiovascular interventions across a range of pathologies. Currently a significant proportion of this potential is still to be realised. In this paper we discuss the opportunities and challenges associated with this realisation. Reviewing the successful elements of model translation for biophysically based models and the emerging supporting technologies, we propose three distinct modes of clinical translation. Finally we outline the challenges ahead that will be fundamental to overcome if the ultimate goal of fully personalised clinical cardiac care is to be achieved.

Noninvasive Personalization of a Cardiac Electrophysiology Model From Body Surface Potential Mapping

GOAL

We use noninvasive data (body surface potential mapping, BSPM) to personalize the main parameters of a cardiac electrophysiological (EP) model for predicting the response to different pacing conditions.

METHODS

First, an efficient forward model is proposed, coupling the Mitchell-Schaeffer transmembrane potential model with a current dipole formulation. Then, we estimate the main parameters of the cardiac model: activation onset location and tissue conductivity. A large patient-specific database of simulated BSPM is generated, from which specific features are extracted to train a machine learning algorithm. The activation onset location is computed from a Kernel Ridge Regression and a second regression calibrates the global ventricular conductivity.

RESULTS

The evaluation of the results is done both on a benchmark dataset of a patient with premature ventricular contraction (PVC) and on five nonischaemic implanted cardiac resynchonization therapy (CRT) patients with a total of 21 different pacing conditions. Good personalization results were found in terms of the activation onset location for the PVC (mean distance error, MDE = 20.3 mm), for the pacing sites (MDE = 21.7 mm) and for the CRT patients (MDE = 24.6 mm). We tested the predictive power of the personalized model for biventricular pacing and showed that we could predict the new electrical activity patterns with a good accuracy in terms of BSPM signals.

CONCLUSION

We have personalized the cardiac EP model and predicted new patient-specific pacing conditions.

SIGNIFICANCE

This is an encouraging first step towards a noninvasive preoperative prediction of the response to different pacing conditions to assist clinicians for CRT patient selection and therapy planning.

ECG imaging of ventricular tachycardia: evaluation against simultaneous non-contact mapping and CMR-derived grey zone

ECG imaging is an emerging technology for the reconstruction of cardiac electric activity from non-invasively measured body surface potential maps. In this case report, we present the first evaluation of transmurally imaged activation times against endocardially reconstructed isochrones for a case of sustained monomorphic ventricular tachycardia (VT). Computer models of the thorax and whole heart were produced from MR images. A recently published approach was applied to facilitate electrode localization in the catheter laboratory, which allows for the acquisition of body surface potential maps while performing non-contact mapping for the reconstruction of local activation times. ECG imaging was then realized using Tikhonov regularization with spatio-temporal smoothing as proposed by Huiskamp and Greensite and further with the spline-based approach by Erem et al. Activation times were computed from transmurally reconstructed transmembrane voltages. The results showed good qualitative agreement between the non-invasively and invasively reconstructed activation times. Also, low amplitudes in the imaged transmembrane voltages were found to correlate with volumes of scar and grey zone in delayed gadolinium enhancement cardiac MR. The study underlines the ability of ECG imaging to produce activation times of ventricular electric activity-and to represent effects of scar tissue in the imaged transmembrane voltages.

Transcatheter Myocardial Needle Chemoablation During Real-Time Magnetic Resonance Imaging: A New Approach to Ablation Therapy for Rhythm Disorders

BACKGROUND

Radiofrequency ablation for ventricular arrhythmias is limited by inability to visualize tissue destruction, by reversible conduction block resulting from edema surrounding lesions, and by insufficient lesion depth. We hypothesized that transcatheter needle injection of caustic agents doped with gadolinium contrast under real-time magnetic resonance imaging (MRI) could achieve deep, targeted, and irreversible myocardial ablation, which would be immediately visible.

METHODS AND RESULTS

Under real-time MRI guidance, ethanol or acetic acid was injected into the myocardium of 8 swine using MRI-conspicuous needle catheters. Chemoablation lesions had identical geometry by in vivo and ex vivo MRI and histopathology, both immediately and after 12 (7-17) days. Ethanol caused stellate lesions with patchy areas of normal myocardium, whereas acetic acid caused homogeneous circumscribed lesions of irreversible necrosis. Ischemic cardiomyopathy was created in 10 additional swine by subselective transcoronary ethanol administration into noncontiguous territories. After 12 (8-15) days, real-time MRI-guided chemoablation-with 2 to 5 injections to create a linear lesion-successfully eliminated the isthmus and local abnormal voltage activities.

CONCLUSIONS

Real-time MRI-guided chemoablation with acetic acid enabled the intended arrhythmic substrate, whether deep or superficial, to be visualized immediately and ablated irreversibly. In an animal model of ischemic cardiomyopathy, obliteration of a conductive isthmus both anatomically and functionally and abolition of local abnormal voltage activities in areas of heterogeneous scar were feasible. This represents the first report of MRI-guided myocardial chemoablation, an approach that could improve the efficacy of arrhythmic substrate ablation in the thick ventricular myocardium.

Personalized Computer Simulation of Diastolic Function in Heart Failure

The search for a parameter representing left ventricular relaxation from non-invasive and invasive diagnostic tools has been extensive, since heart failure (HF) with preserved ejection fraction (HF-pEF) is a global health problem. We explore here the feasibility using patient-specific cardiac computer modeling to capture diastolic parameters in patients suffering from different degrees of systolic HF. Fifty eight patients with idiopathic dilated cardiomyopathy have undergone thorough clinical evaluation, including cardiac magnetic resonance imaging (MRI), heart catheterization, echocardiography, and cardiac biomarker assessment. A previously-introduced framework for creating multi-scale patient-specific cardiac models has been applied on all these patients. Novel parameters, such as global stiffness factor and maximum left ventricular active stress, representing cardiac active and passive tissue properties have been computed for all patients. Invasive pressure measurements from heart catheterization were then used to evaluate ventricular relaxation using the time constant of isovolumic relaxation Tau (τ). Parameters from heart catheterization and the multi-scale model have been evaluated and compared to patient clinical presentation. The model parameter global stiffness factor, representing diastolic passive tissue properties, is correlated significantly across the patient population with τ. This study shows that multi-modal cardiac models can successfully capture diastolic (dys) function, a prerequisite for future clinical trials on HF-pEF.

Computational Representations of Myocardial Infarct Scars and Implications for Arrhythmogenesis

Image-based computational modeling is becoming an increasingly used clinical tool to provide insight into the mechanisms of reentrant arrhythmias. In the context of ischemic heart disease, faithful representation of the electrophysiological properties of the infarct region within models is essential, due to the scars known for arrhythmic properties. Here, we review the different computational representations of the infarcted region, summarizing the experimental measurements upon which they are based. We then focus on the two most common representations of the scar core (complete insulator or electrically passive tissue) and perform simulations of electrical propagation around idealized infarct geometries. Our simulations highlight significant differences in action potential duration and focal effective refractory period (ERP) around the scar, driven by differences in electrotonic loading, depending on the choice of scar representation. Finally, a novel mechanism for arrhythmia induction, following a focal ectopic beat, is demonstrated, which relies on localized gradients in ERP directly caused by the electrotonic sink effects of the neighboring passive scar.

Uncertainty and variability in computational and mathematical models of cardiac physiology

KEY POINTS

Mathematical and computational models of cardiac physiology have been an integral component of cardiac electrophysiology since its inception, and are collectively known as the Cardiac Physiome. We identify and classify the numerous sources of variability and uncertainty in model formulation, parameters and other inputs that arise from both natural variation in experimental data and lack of knowledge. The impact of uncertainty on the outputs of Cardiac Physiome models is not well understood, and this limits their utility as clinical tools. We argue that incorporating variability and uncertainty should be a high priority for the future of the Cardiac Physiome. We suggest investigating the adoption of approaches developed in other areas of science and engineering while recognising unique challenges for the Cardiac Physiome; it is likely that novel methods will be necessary that require engagement with the mathematics and statistics community.

ABSTRACT

The Cardiac Physiome effort is one of the most mature and successful applications of mathematical and computational modelling for describing and advancing the understanding of physiology. After five decades of development, physiological cardiac models are poised to realise the promise of translational research via clinical applications such as drug development and patient-specific approaches as well as ablation, cardiac resynchronisation and contractility modulation therapies. For models to be included as a vital component of the decision process in safety-critical applications, rigorous assessment of model credibility will be required. This White Paper describes one aspect of this process by identifying and classifying sources of variability and uncertainty in models as well as their implications for the application and development of cardiac models. We stress the need to understand and quantify the sources of variability and uncertainty in model inputs, and the impact of model structure and complexity and their consequences for predictive model outputs. We propose that the future of the Cardiac Physiome should include a probabilistic approach to quantify the relationship of variability and uncertainty of model inputs and outputs.

Biophysical Modeling Predicts Ventricular Tachycardia Inducibility and Circuit Morphology: A Combined Clinical Validation and Computer Modeling Approach

INTRODUCTION

Computational modeling of cardiac arrhythmogenesis and arrhythmia maintenance has made a significant contribution to the understanding of the underlying mechanisms of arrhythmia. We hypothesized that a cardiac model using personalized electro-anatomical parameters could define the underlying ventricular tachycardia (VT) substrate and predict reentrant VT circuits. We used a combined modeling and clinical approach in order to validate the concept.

METHODS AND RESULTS

Non-contact electroanatomic mapping studies were performed in 7 patients (5 ischemics, 2 non-ischemics). Three ischemic cardiomyopathy patients underwent a clinical VT stimulation study. Anatomical information was obtained from cardiac magnetic resonance imaging (CMR) including high-resolution scar imaging. A simplified biophysical mono-domain action potential model personalized with the patients' anatomical and electrical information was used to perform in silico VT stimulation studies for comparison. The personalized in silico VT stimulations were able to predict VT inducibility as well as the macroscopic characteristics of the VT circuits in patients who had clinical VT stimulation studies. The patients with positive clinical VT stimulation studies had wider distribution of action potential duration restitution curve (APD-RC) slopes and APDs than the patient with a negative VT stimulation study. The exit points of reentrant VT circuits encompassed a higher percentage of the maximum APD-RC slope compared to the scar and non-scar areas, 32%, 4%, and 0.2%, respectively.

CONCLUSIONS

VT stimulation studies can be simulated in silico using a personalized biophysical cardiac model. Myocardial spatial heterogeneity of APD restitution properties and conductivity may help predict the location of crucial entry/exit points of reentrant VT circuits.

Personalized Models of Human Atrial Electrophysiology Derived From Endocardial Electrograms

OBJECTIVE

Computational models represent a novel framework for understanding the mechanisms behind atrial fibrillation (AF) and offer a pathway for personalizing and optimizing treatment. The characterization of local electrophysiological properties across the atria during procedures remains a challenge. The aim of this work is to characterize the regional properties of the human atrium from multielectrode catheter measurements.

METHODS

We propose a novel method that characterizes regional electrophysiology properties by fitting parameters of an ionic model to conduction velocity and effective refractory period restitution curves obtained by a s₁-s₂ pacing protocol applied through a multielectrode catheter. Using an in-silico dataset we demonstrate that the fitting method can constrain parameters with a mean error of 21.9 ± 16.1% and can replicate conduction velocity and effective refractory curves not used in the original fitting with a relative error of 4.4 ± 6.9%.

RESULTS

We demonstrate this parameter estimation approach on five clinical datasets recorded from AF patients. Recordings and parametrization took approx. 5 and 6 min, respectively. Models fitted restitution curves with an error of ~ 5% and identify a unique parameter set. Tissue properties were predicted using a two-dimensional atrial tissue sheet model. Spiral wave stability in each case was predicted using tissue simulations, identifying distinct stable (2/5), meandering and breaking up (2/5), and unstable self-terminating (1/5) spiral tip patterns for different cases.

CONCLUSION AND SIGNIFICANCE

We have developed and demonstrated a robust and rapid approach for personalizing local ionic models from a clinically tractable.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

Examining the Impact of Prior Models in Transmural Electrophysiological Imaging: A Hierarchical Multiple-Model Bayesian Approach

Noninvasive cardiac electrophysiological (EP) imaging aims to mathematically reconstruct the spatiotemporal dynamics of cardiac sources from body-surface electrocardiographic (ECG) data. This ill-posed problem is often regularized by a fixed constraining model. However, a fixed-model approach enforces the source distribution to follow a pre-assumed structure that does not always match the varying spatiotemporal distribution of actual sources. To understand the model-data relation and examine the impact of prior models, we present a multiple-model approach for volumetric cardiac EP imaging where multiple prior models are included and automatically picked by the available ECG data. Multiple models are incorporated as an Lp-norm prior for sources, where p is an unknown hyperparameter with a prior uniform distribution. To examine how different combinations of models may be favored by different measurement data, the posterior distribution of cardiac sources and hyperparameter p is calculated using a Markov Chain Monte Carlo (MCMC) technique. The importance of multiple-model prior was assessed in two sets of synthetic and real-data experiments, compared to fixed-model priors (using Laplace and Gaussian priors). The results showed that the posterior combination of models (the posterior distribution of p) as determined by the ECG data differed substantially when reconstructing sources with different sizes and structures. While the use of fixed models is best suited in situations where the prior assumption fits the actual source structures, the use of an automatically adaptive set of models may have the ability to better address model-data mismatch and to provide consistent performance in reconstructing sources with different properties.

A model model: a commentary on DiFrancesco and Noble (1985) 'A model of cardiac electrical activity incorporating ionic pumps and concentration changes'

This paper summarizes the advances made by the DiFrancesco and Noble (DFN) model of cardiac cellular electrophysiology, which was published in Philosophical Transactions B in 1985. This model was developed at a time when the introduction of new techniques and provision of experimental data had resulted in an explosion of knowledge about the cellular and biophysical properties of the heart. It advanced the cardiac modelling field from a period when computer models considered only the voltage-dependent channels in the surface membrane. In particular, it included a consideration of changes of both intra- and extracellular ionic concentrations. In this paper, we summarize the most important contributions of the DiFrancesco and Noble paper. We also describe how computer modelling has developed subsequently with the extension from the single cell to the whole heart as well as its use in understanding disease and predicting the effects of pharmaceutical interventions. This commentary was written to celebrate the 350th anniversary of the journal Philosophical Transactions of the Royal Society.

Clinical Diagnostic Biomarkers from the Personalization of Computational Models of Cardiac Physiology

Computational modelling of the heart is rapidly advancing to the point of clinical utility. However, the difficulty of parameterizing and validating models from clinical data indicates that the routine application of truly predictive models remains a significant challenge. We argue there is significant value in an intermediate step towards prediction. This step is the use of biophysically based models to extract clinically useful information from existing patient data. Specifically in this paper we review methodologies for applying modelling frameworks for this goal in the areas of quantifying cardiac anatomy, estimating myocardial stiffness and optimizing measurements of coronary perfusion. Using these indicative examples of the general overarching approach, we finally discuss the value, ongoing challenges and future potential for applying biophysically based modelling in the clinical context.

Towards Personalized Cardiology: Multi-Scale Modeling of the Failing Heart

BACKGROUND

Despite modern pharmacotherapy and advanced implantable cardiac devices, overall prognosis and quality of life of HF patients remain poor. This is in part due to insufficient patient stratification and lack of individualized therapy planning, resulting in less effective treatments and a significant number of non-responders.

METHODS AND RESULTS

State-of-the-art clinical phenotyping was acquired, including magnetic resonance imaging (MRI) and biomarker assessment. An individualized, multi-scale model of heart function covering cardiac anatomy, electrophysiology, biomechanics and hemodynamics was estimated using a robust framework. The model was computed on n=46 HF patients, showing for the first time that advanced multi-scale models can be fitted consistently on large cohorts. Novel multi-scale parameters derived from the model of all cases were analyzed and compared against clinical parameters, cardiac imaging, lab tests and survival scores to evaluate the explicative power of the model and its potential for better patient stratification. Model validation was pursued by comparing clinical parameters that were not used in the fitting process against model parameters.

CONCLUSION

This paper illustrates how advanced multi-scale models can complement cardiovascular imaging and how they could be applied in patient care. Based on obtained results, it becomes conceivable that, after thorough validation, such heart failure models could be applied for patient management and therapy planning in the future, as we illustrate in one patient of our cohort who received CRT-D implantation.

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

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

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

To properly describe the electrical activity of the left ventricle, it is necessary to model the Purkinje fibers, responsible for the fast and coordinate ventricular activation, and their interaction with the muscular propagation. The aim of this work is to propose a methodology for the generation of a patient-specific Purkinje network driven by clinical measurements of the activation times related to pathological propagations. In this case, one needs to consider a strongly coupled problem between the network and the muscle, where the feedback from the latter to the former cannot be neglected as in a normal propagation. We apply the proposed strategy to data acquired on three subjects, one of them suffering from muscular conduction problems owing to a scar and the other two with a muscular pre-excitation syndrome (Wolff-Parkinson-White). To assess the accuracy of the proposed method, we compare the results obtained by using the patient-specific Purkinje network generated by our strategy with the ones obtained by using a non-patient-specific network. The results show that the mean absolute errors in the activation time is reduced for all the cases, highlighting the importance of including a patient-specific Purkinje network in computational models.

Methodology for image-based reconstruction of ventricular geometry for patient-specific modeling of cardiac electrophysiology

Patient-specific modeling of ventricular electrophysiology requires an interpolated reconstruction of the 3-dimensional (3D) geometry of the patient ventricles from the low-resolution (Lo-res) clinical images. The goal of this study was to implement a processing pipeline for obtaining the interpolated reconstruction, and thoroughly evaluate the efficacy of this pipeline in comparison with alternative methods. The pipeline implemented here involves contouring the epi- and endocardial boundaries in Lo-res images, interpolating the contours using the variational implicit functions method, and merging the interpolation results to obtain the ventricular reconstruction. Five alternative interpolation methods, namely linear, cubic spline, spherical harmonics, cylindrical harmonics, and shape-based interpolation were implemented for comparison. In the thorough evaluation of the processing pipeline, Hi-res magnetic resonance (MR), computed tomography (CT), and diffusion tensor (DT) MR images from numerous hearts were used. Reconstructions obtained from the Hi-res images were compared with the reconstructions computed by each of the interpolation methods from a sparse sample of the Hi-res contours, which mimicked Lo-res clinical images. Qualitative and quantitative comparison of these ventricular geometry reconstructions showed that the variational implicit functions approach performed better than others. Additionally, the outcomes of electrophysiological simulations (sinus rhythm activation maps and pseudo-ECGs) conducted using models based on the various reconstructions were compared. These electrophysiological simulations demonstrated that our implementation of the variational implicit functions-based method had the best accuracy.

Data-driven estimation of cardiac electrical diffusivity from 12-lead ECG signals

Diagnosis and treatment of dilated cardiomyopathy (DCM) is challenging due to a large variety of causes and disease stages. Computational models of cardiac electrophysiology (EP) can be used to improve the assessment and prognosis of DCM, plan therapies and predict their outcome, but require personalization. In this work, we present a data-driven approach to estimate the electrical diffusivity parameter of an EP model from standard 12-lead electrocardiograms (ECG). An efficient forward model based on a mono-domain, phenomenological Lattice-Boltzmann model of cardiac EP, and a boundary element-based mapping of potentials to the body surface is employed. The electrical diffusivity of myocardium, left ventricle and right ventricle endocardium is then estimated using polynomial regression which takes as input the QRS duration and electrical axis. After validating the forward model, we computed 9500 EP simulations on 19 different DCM patients in just under three seconds each to learn the regression model. Using this database, we quantify the intrinsic uncertainty of electrical diffusion for given ECG features and show in a leave-one-patient-out cross-validation that the regression method is able to predict myocardium diffusion within the uncertainty range. Finally, our approach is tested on the 19 cases using their clinical ECG. 84% of them could be personalized using our method, yielding mean prediction errors of 18.7ms for the QRS duration and 6.5° for the electrical axis, both values being within clinical acceptability. By providing an estimate of diffusion parameters from readily available clinical data, our data-driven approach could therefore constitute a first calibration step toward a more complete personalization of cardiac EP.

Advances in modeling ventricular arrhythmias: from mechanisms to the clinic

Modern cardiovascular research has increasingly recognized that heart models and simulation can help interpret an array of experimental data and dissect important mechanisms and interrelationships, with developments rooted in the iterative interaction between modeling and experimentation. This article reviews the progress made in simulating cardiac electrical behavior at the level of the organ and, specifically, in the development of models of ventricular arrhythmias and fibrillation, as well as their termination (defibrillation). The ability to construct multiscale models of ventricular arrhythmias, representing integrative behavior from the molecule to the entire organ, has enabled mechanistic inquiry into the dynamics of ventricular arrhythmias in the diseased myocardium, in understanding drug-induced proarrhythmia, and in the development of new modalities for defibrillation, to name a few. In this article, we also review the initial use of ventricular models of arrhythmia in personalized diagnosis, treatment planning, and prevention of sudden cardiac death. Implementing individualized cardiac simulations at the patient bedside is poised to become one of the most thrilling examples of computational science and engineering approaches in translational medicine.

Bridging the gap between computation and clinical biology: validation of cable theory in humans

Introduction: Computerized simulations of cardiac activity have significantly contributed to our understanding of cardiac electrophysiology, but techniques of simulations based on patient-acquired data remain in their infancy. We sought to integrate data acquired from human electrophysiological studies into patient-specific models, and validated this approach by testing whether electrophysiological responses to sequential premature stimuli could be predicted in a quantitatively accurate manner.

Methods: Eleven patients with structurally normal hearts underwent electrophysiological studies. Semi-automated analysis was used to reconstruct activation and repolarization dynamics for each electrode. This S₂ extrastimuli data was used to inform individualized models of cardiac conduction, including a novel derivation of conduction velocity restitution. Activation dynamics of multiple premature extrastimuli were then predicted from this model and compared against measured patient data as well as data derived from the ten-Tusscher cell-ionic model.

Results: Activation dynamics following a premature S₃ were significantly different from those after an S₂. Patient specific models demonstrated accurate prediction of the S₃ activation wave, (Pearson's R² = 0.90, median error 4%). Examination of the modeled conduction dynamics allowed inferences into the spatial dispersion of activation delay. Further validation was performed against data from the ten-Tusscher cell-ionic model, with our model accurately recapitulating predictions of repolarization times (R² = 0.99).

Conclusions: Simulations based on clinically acquired data can be used to successfully predict complex activation patterns following sequential extrastimuli. Such modeling techniques may be useful as a method of incorporation of clinical data into predictive models.

Computational cardiology: how computer simulations could be used to develop new therapies and advance existing ones

This article reviews the latest developments in computational cardiology. It focuses on the contribution of cardiac modelling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modelling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Computational medicine: translating models to clinical care

Because of the inherent complexity of coupled nonlinear biological systems, the development of computational models is necessary for achieving a quantitative understanding of their structure and function in health and disease. Statistical learning is applied to high-dimensional biomolecular data to create models that describe relationships between molecules and networks. Multiscale modeling links networks to cells, organs, and organ systems. Computational approaches are used to characterize anatomic shape and its variations in health and disease. In each case, the purposes of modeling are to capture all that we know about disease and to develop improved therapies tailored to the needs of individuals. We discuss advances in computational medicine, with specific examples in the fields of cancer, diabetes, cardiology, and neurology. Advances in translating these computational methods to the clinic are described, as well as challenges in applying models for improving patient health.

Understanding the mechanisms amenable to CRT response: from pre-operative multimodal image data to patient-specific computational models

This manuscript describes our recent developments towards better understanding of the mechanisms amenable to cardiac resynchronization therapy response. We report the results from a full multimodal dataset corresponding to eight patients from the euHeart project. The datasets include echocardiography, MRI and electrophysiological studies. We investigate two aspects. The first one focuses on pre-operative multimodal image data. From 2D echocardiography and 3D tagged MRI images, we compute atlas based dyssynchrony indices. We complement these indices with presence and extent of scar tissue and correlate them with CRT response. The second one focuses on computational models. We use pre-operative imaging to generate a patient-specific computational model. We show results of a fully automatic personalized electromechanical simulation. By case-per-case discussion of the results, we highlight the potential and key issues of this multimodal pipeline for the understanding of the mechanisms of CRT response and a better patient selection.

Characterization and modeling of the peripheral cardiac conduction system

The development of biophysical models of the heart has the potential to get insights in the patho-physiology of the heart, which requires to accurately modeling anatomy and function. The electrical activation sequence of the ventricles depends strongly on the cardiac conduction system (CCS). Its morphology and function cannot be observed in vivo, and therefore data available come from histological studies. We present a review on data available of the peripheral CCS including new experiments. In order to build a realistic model of the CCS we designed a procedure to extract morphological characteristics of the CCS from stained calf tissue samples. A CCS model personalized with our measurements has been built using L-systems. The effect of key unknown parameters of the model in the electrical activation of the left ventricle has been analyzed. The CCS models generated share the main characteristics of observed stained Purkinje networks. The timing of the simulated electrical activation sequences were in the physiological range for CCS models that included enough density of PMJs. These results show that this approach is a potential methodology for collecting knowledge-domain data and build improved CCS models of the heart automatically.