Abstract:

Generally, models of cardiac electrophysiology describe physiologic conditions in detail. However, other conditions, such as drug interactions or mutations of ion channels are of interest for research. Therefore, the simulated ion currents have to be fitted to measured voltage or patch clamp data. In this work, three different methods for the model parametrization were compared: one based on Powells algorithm implemented in a modular C++ framework and two optimization techniques realized in Matlab. The latter two approaches differed in solving the ordinary differential equations describing the channel gating. They can either be approximated numerically or solved analytically, since the transmembrane voltage is a piecewise constant function during the applied clamp protocol. All three methods were compared regarding computing time and quality of the fit using least squares. The modular C++ framework was slower than the numerical Matlab method, which took longer than the analytical one. The quality of the fit was similar for almost all analyzed methods. Therefore, the analytical method grants a fast and reliable solution for the calibration of ion current models for applications with constant membrane voltage, as e.g. in case of voltage or patch clamp data.