Systematic reduction of a detailed atrial myocyte model

Cardiac arrhythmias are a major health concern and often involve poorly understood mechanisms. Mathematical modeling is able to provide insights into these mechanisms which might result in better treatment options. A key element of this modeling is a description of the electrophysiological properties of cardiac cells. A number of electrophysiological models have been developed, ranging from highly detailed and complex models, containing numerous parameters and variables, to simplified models in which variables and parameters no longer directly correspond to electrophysiological quantities. In this study, we present a systematic reduction of the complexity of the detailed model of Koivumaki et al. using the recently developed manifold boundary approximation method. We reduce the original model, containing 42 variables and 37 parameters, to a model with only 11 variables and 5 parameters and show that this reduced model can accurately reproduce the action potential shape and restitution curve of the original model. The reduced model contains only five currents and all variables and parameters can be directly linked to electrophysiological quantities. Due to its reduction in complexity, simulation times of our model are decreased more than three-fold. Furthermore, fitting the reduced model to clinical data is much more efficient, a potentially important step towards patient-specific modeling.