Data-Driven Uncertainty Quantification for Cardiac Electrophysiological Models: Impact of Physiological Variability on Action Potential and Spiral Wave Dynamics

Abstract

Computational modeling of cardiac electrophysiology (EP) has recently transitioned from a scientific research tool to clinical applications. To ensure reliability of clinical or regulatory decisions made using cardiac EP models, it is vital to evaluate the uncertainty in model predictions. Model predictions are uncertain because there is typically substantial uncertainty in model input parameters, due to measurement error or natural variability. While there has been much recent uncertainty quantification (UQ) research for cardiac EP models, all previous work has been limited by either: (i) considering uncertainty in only a subset of the full set of parameters; and/or (ii) assigning arbitrary variation to parameters (e.g., ±10 or 50% around mean value) rather than basing the parameter uncertainty on experimental data. In our recent work we overcame the first limitation by performing UQ and sensitivity analysis using a novel canine action potential model, allowing all parameters to be uncertain, but with arbitrary variation. Here, we address the second limitation by extending our previous work to use data-driven estimates of parameter uncertainty. Overall, we estimated uncertainty due to population variability in all parameters in five currents active during repolarization: inward potassium rectifier, transient outward potassium, L-type calcium, rapidly and slowly activating delayed potassium rectifier; 25 parameters in total (all model parameters except fast sodium current parameters). A variety of methods was used to estimate the variability in these parameters. We then propagated the uncertainties through the model to determine their impact on predictions of action potential shape, action potential duration (APD) prolongation due to drug block, and spiral wave dynamics. Parameter uncertainty had a significant effect on model predictions, especially L-type calcium current parameters. Correlation between physiological parameters was determined to play a role in physiological realism of action potentials. Surprisingly, even model outputs that were relative differences, specifically drug-induced APD prolongation, were heavily impacted by the underlying uncertainty. This is the first data-driven end-to-end UQ analysis in cardiac EP accounting for uncertainty in the vast majority of parameters, including first in tissue, and demonstrates how future UQ could be used to ensure model-based decisions are robust to all underlying parameter uncertainties.

Fiber Organization Has Little Effect on Electrical Activation Patterns During Focal Arrhythmias in the Left Atrium

“…One of the key obstacles towards the development of constructing patient-specific models is not only speed in simulations [11] but adequate quantification, verification and information for such models [29], [30]. Detailed models have a very high-dimensional parameter space; moreover, most of the parameters cannot be measured on an individual basis and need to be postulated based on some generic values or recomputed libraries [31], [32].…”

Neural Network Differential Equations For Ion Channel Modelling

“…We have shown two approaches for doing this: the first one uses a neural network to fully model the right-hand side of the ODEs; the second one uses a neural network to model only the missing dynamics of the model—discrepancy between a model and the true system. Assessing the model discrepancy in ion channel kinetics is vital to constructing accurate action potential models (Mirams et al, 2016 ; Clayton et al, 2020 ; Pathmanathan et al, 2020 ), but most studies assume correct ion channel kinetics models when fitting maximum conductances of different current types in an action potential model (Kaur et al, 2014 ; Groenendaal et al, 2015 ; Johnstone et al, 2016 ; Lei et al, 2017 ; Pouranbarani et al, 2019 ). Previous studies attempted to use different machine learning techniques and statistical methods to model the differences between the mechanistic model and the data.…”

Anatomical Model of Rat Ventricles to Study Cardiac Arrhythmias under Infarction Injury

“…Thus, we think that the instabilities which we observe are related to cellular (EAD) dynamics and not to conduction velocity restitution. The dependence of spiral waves on cellular properties has been demonstrated by Pathmanathan and co-authors in canine models of myocardial tissue [44]. They studied the effects of parameter uncertainties in a cellular model on spiral wave dynamics in 2D tissue.…”