Abstract
Spatial–temporal calcium dynamics due to calcium release, buffering and re-uptaking plays a central role in studying excitation–contraction (E–C) coupling in both normal and diseased cardiac myocytes. In this paper, we employ a meshless method, namely, the local radial basis function collocation method (LRBFCM), to model such calcium behaviors by solving a nonlinear system of reaction–diffusion partial differential equations. In particular, a simplified structural unit containing a single transverse tubule (T-tubule) and its surrounding half sarcomeres is investigated using the meshless method. Numerical results are compared with those generated by finite element methods, showing the capability and efficiency of the LRBFCM in modeling calcium dynamics in ventricular myocytes. The single T-tubule model is also extended to the whole-cell scale with T-tubules excluded to demonstrate the scalability of the proposed meshless method in handling very large domains. The experiments have shown that the LRBFCM is suitable to multiscale modeling of calcium dynamics in ventricular myocytes with high accuracy and efficiency.