Estimation of stochastic behaviour in cardiac myocytes: I. Ca 2+ movements inside the cytosol and sarcoplasmic reticulum on curvilinear domains

Since the discovery of Ca sparks and their stochastic behaviour in cardiac myocytes, models have focused on the inclusion of stochasticity in their studies. While most models pay much attention to the stochastic modelling of cytosolic Ca concentration the coupling of Ca sparks and blinks in a stochastic model has not been explored fully. The cell morphology in in silico studies in the past is assumed to be Cartesian, spherical or cylindrical. The application on curvilinear grids can easily address certain restrictions posed by such grid set up and provide more realistic cell morphology. In this paper, we present a stochastic reaction-diffusion model that couples Ca sparks and blinks in realistic shapes of cells in curvilinear domains. Transformation of the model was performed to the curvilinear coordinate system. The set of equations is used to produce Ca waves initiated from sparks and blinks. A non-buffered and non-dyed version as well as a buffered and dyed version of these equations were studied in light of observing the dynamics on the two different systems. For comparison, results for both the Cartesian and curvilinear grids are provided. A successful demonstration of the application of curvilinear grids serving as basis for future developments.