A computational model of a human single sinoatrial node cell

For the investigation of the spontaneous rhythmical activity response in the application of cardiac neuromodulation, we formulated a human sinoatrial node (SAN) cell model. With the aim of decreasing elevated heart rate (HR), we want to establish a hardware-in-the-loop system including this model for the analysis of optimal stimulation patterns of the neurostimulation system. Base model structures are adopted from rabbit SAN cell models available in literature and conveyed with Hodgkin-Huxley-type model equations describing the complex time and voltage dependent activation and deactivation processes of individual ion channels. The resulting model consists of 15 currents which are currently known to be responsible for the generation of the membrane action potential (AP). The model reproduces AP frequencies equivalent to those measured in isolated human SAN cells with a resulting HR of 71.8 bpm. Model validation via simulation of the inhibitory effect of ivabradine showed accordance with experimental results obtained in human studies. Furthermore, we could validate the model in regard to its HR effects upon parasympathetic stimulation with results obtained in a human trial study.