A stochastic mathematical model for coupling the cytosolic and sarcoplasmic calcium movements in diseased cardiac myocytes

Several computational studies have been undertaken to explore the Ca2+-induced Ca2+ release (CICR) events in cardiac myocytes and along with experimental studies it has given us invaluable insight into the mechanism of CICR from spark/blink initiation to termination and regulation, and their interplay under normal and pathological conditions. The computational modelling of this mechanism has mainly been investigated using coupled differential equations (DEs). However, there is a lack of computational investigation into (1) how the different formulation of coupled DEs capture the Ca2+ movement in the cytosol and sarcoplasmic reticulum (SR), (2) the buffer and dye inclusion in both compartments, and (3) the effect of buffer and dye properties on the calcium behaviour. This work is set out to explore (1) the effect of different coupled formulation of DEs on spark/blink occurrence, (2) the inclusion of improved sarcoplasmic buffering properties, and (3) the effects of cytosolic and sarcoplasmic dye and buffer properties on Ca2+ movement. The simulation results show large discrepancies between different formulations of the governing equations. Additionally, extension of the model to include sarcoplasmic buffering properties show normalised fluorescent dye profiles to be in good agreement with experimental and amongst its one- and two-dimensional representations.