Sensitivity and uncertainty analysis of two human atrial cardiac cell models using Gaussian process emulators

Cardiac cell models reconstruct the action potential and calcium dynamics of cardiac myocytes, and are becoming widely used research tools. These models are highly detailed, with many parameters in the equations that describe current flow through ion channels, pumps, and exchangers in the cell membrane, and so it is difficult to link changes in model inputs to model behaviours. The aim of the present study was to undertake sensitivity and uncertainty analysis of two models of the human atrial action potential. We used Gaussian processes to emulate the way that 11 features of the action potential and calcium transient produced by each model depended on a set of. The emulators were trained by maximising likelihood conditional on a set of design data, obtained from 300 model evaluations. For each model evaluation, the set of inputs was obtained from uniform distributions centred on the default values for each parameter, using latin-hypercube sampling. First order and total effect sensitivity indices were calculated for each combination of input and output. First order indices were well correlated with the square root of sensitivity indices obtained by partial least squares regression of the design data. The sensitivity indices highlighted a difference in the balance of inward and outward currents during the plateau phase of the action potential in each model, with the consequence that changes to one parameter can have opposite effects in the two models. Overall the interactions among inputs were not as important as the first order effects, indicating that model parameters tend to have independent effects on the model outputs. This study has shown that Gaussian process emulators are an effective tool for sensitivity and uncertainty analysis of cardiac cell models.Author summaryThe time course of the cardiac action potential is determined by the balance of inward and outward currents across the cell membrane, and these in turn depend on dynamic behaviour of ion channels, pumps and exchangers in the cell membrane. Cardiac cell models reconstruct the action potential by representing transmembrane current as a set of stiff and nonlinear ordinary differential equations. These models capture biophysical detail, but are complex and have large numbers of parameters, so cause and effect relationships are difficult to identify. In recent years there has been an increasing interest in uncertainty and variability in computational models, and a number of tools have been developed. In this study we have used one of these tools, Gaussian process emulators, to compare and contrast two models of the human atrial action potential. We obtained sensitivity indices based on the proportion of variance in a model output that is accounted for by variance in each of the model parameters. These sensitivity indices highlighted the model parameters that had the most influence on the model outputs, and provided a means to make a quantitative comparison between the models.