ABSTRACT
Time series of spatially-extended two-dimensional recordings are the cornerstone of basic and clinical cardiac electrophysiology. The data source may be either multipolar catheters, multi-electrode arrays, optical mapping with the help of voltage and calcium-sensitive fluorescent dyes, or the output of simulation studies. The resulting data cubes (usually two spatial and one temporal dimension) are shared either as movie files or, after additional processing, various graphs and tables. However, such data products can only convey a limited view of the data. It will be beneficial if the data consumers can interactively process the data, explore different processing options and change its visualization. This paper presents the Unified Electrophysiology Mapping Framework (Unimapper) to facilitate the exchange of electrophysiology data. Its pedigree includes a Java-based optical mapping application. The core of Unimapper is a website and a collection of JavaScript utility functions for data import and visualization (including multi-channel visualization for simultaneous voltage/calcium mapping), basic image processing (e.g., smoothing), basic signal processing (e.g., signal detrending), and advanced processing (e.g., phase calculation using the Hilbert transform). Additionally, Unimapper can optionally use graphics processing units (GPUs) for computationally intensive operations. The Unimapper ecosystem also includes utility libraries for commonly used scientific programming languages (MATLAB, Python, and Julia) that allow the data producers to convert images and recorded signals into a standard format readable by Unimapper. Unimapper can act as a nexus to share electrophysiology data - whether recorded experimentally, clinically or generated by simulation - and enhance communication and collaboration among researchers.
ABSTRACT
Finding appropriate values for parameters in mathematical models of cardiac cells is a challenging task. Here, we show that it is possible to obtain good parameterizations in as little as 30-40 s when as many as 27 parameters are fit simultaneously using a genetic algorithm and two flexible phenomenological models of cardiac action potentials. We demonstrate how our implementation works by considering cases of "model recovery" in which we attempt to find parameter values that match model-derived action potential data from several cycle lengths. We assess performance by evaluating the parameter values obtained, action potentials at fit and non-fit cycle lengths, and bifurcation plots for fidelity to the truth as well as consistency across different runs of the algorithm. We also fit the models to action potentials recorded experimentally using microelectrodes and analyze performance. We find that our implementation can efficiently obtain model parameterizations that are in good agreement with the dynamics exhibited by the underlying systems that are included in the fitting process. However, the parameter values obtained in good parameterizations can exhibit a significant amount of variability, raising issues of parameter identifiability and sensitivity. Along similar lines, we also find that the two models differ in terms of the ease of obtaining parameterizations that reproduce model dynamics accurately, most likely reflecting different levels of parameter identifiability for the two models.
