Nagaiah C. and Rüdiger S.

We developed a 3D finite element simulator interface for the numerical simulation of stochastic and deterministic equations for single and multiple clusters of Ca^(2+) releasing channels. Mathematically, diffusion, reactions and membrane transport of calcium ions in cells are represented by a coupled system of reaction-diffusion equations. We adopted a hybrid algorithm to address the coupling of the deterministic reaction-diffusion equations for Ca^(2+) and Ca^(2+) buffers and Markovian dynamics of IP_3 R channel gating. Using highly unstructured meshes, our method bridges many orders of magnitude to represent accurately the Ca^(2+)distribution from the single channel to entire cells with multiple clusters of channels. To save computational time, a conforming finite element method is employed for the spatial discretization and adaptive and higher order linearly implicit methods, Rosenbrock type methods, are used for the time integration. This allows an efficient representation of inhomogeneous intra-cluster Ca^(2+) distribution at the nanometer scale even for whole-cell simulations with multiple clusters. Numerical results are demonstrated for different fine spatial resolution meshes as well as different higher order time integrators to insure the numerical convergence of schemes which we apply to study the long time behavior. The parallelization is shown to be essential by the numerical study of long time behavior of calcium concentration. We further present the parallel scalability of the deterministic equations for different arrangements of clusters. The main emphasis is on large scale and long time behavior of the studied equations that capture the detailed local dynamics as well as the temporal hierarchy of dynamical processes. Our approach thus extends our earlier simulations of release from single channels and clusters of channels and systematically integrates stochasticity on all scales of a cell’s calcium dynamics.

Intracellular calcium signals, Gillespie method, finite elements, parallel computations.

[1]ALFONSI A, CANCES E, TURINICI G, DI VENTURA B & HUISINGA W. 2004. Exact simulation of hybrid stochastic and determinis- tic models for biochemical systems. INRIA Rapport de Recherche, Themes NUM et BIO, 5435.

[2]BASTIAN P, BIRKEN K, LANG S, JOHANNSEN K, NEUß N, RENTZ-REICHERT H & WIENERS C. 1997. UG: A flexible software toolbox for solving partial differential equations. Computing and Visualization in Science, 1: 27-40.

[3]BASTIAN P, BLATT M, DEDNER A, ENGWER C, KLOFKORN R, KORNHUBER R, OHLBERGER M & SANDER O. 2008. A generic grid interface for parallel and adaptive scientific computing. Part II: implementation and tests in DUNE. Computing, 82(2):121-138, July.

[4] BERRIDGE MJ. 1990. Calcium oscillations. J. Biol. Chem.,265(17): 9583-9586.

[5] BERRIDGE MJ. 1993. Inositol trisphosphate and calcium signalling. Nature, 361: 315.

[6] BERRIDGE MJ, LIPP P & BOOTMAN MD. 2000. The versatility and universality of calcium signalling. Nature Rev. Mol. Cell Biol., 1: 11-22.

[7] BEZPROZVANNY I, WATRAS J & EHRLICH BE. 1991. Bell-shaped calcium-response curves of Ins(1,4,5)P3- and calcium-gated channels from endoplasmatic reticulum of cerebellum. Nature, 351: 751-754.

[8] DEKKER K & VERWER JG. 1984. Stability of Runge Kutta methods for stiff nonlinear differential equations. North Holland Elsevier Science Publishers.

[9] DEMURO A & PARKER I. 2008. Multi-dimensional resolution of elementary Ca2+ signals by simultaneous multi-focal imaging. Cell Calcium, 43(4): 367-374.

[10] DeYOUNG G & KEIZER J. 1992. A single-pool inositol 1,4,5 trisphosphate-receptor-based model for agonist-stimulated oscil- lations in Ca2+ concentration. Proc. Natl. Acad. Sci. USA, 89:9895-9899.

[11] GILLESPIE DT. 1977. Exact stochastic simulating of coupled chemical reactions. J. Phys. Chem., 81: 2340-2361.

[12] GOBBERT MK. 2008. Long-time simulations on high resolution meshes to model calcium waves in a heart cell. SIAM J. Sci. Com- put., 30: 2922-2947, October.

[13] GUSTAFSSON K, LUNDH M & SODERLIND G. 1988. A PI step-size control for the numerical solution of ordinary differential equations. BIT, 28(2): 270-287.

[14] HANHART AL, GOBBERT MK & IZU LT. 2004. A memory-efficient finite element method for systems of reaction-diffusion equations with non-smooth forcing. J. Comput. Appl. Math., 169(2): 431-458, August.

[15] LANG J. 2001. Adaptive Multilevel Solution of Nonlinear Parabolic PDE Systems, volume 16 of Lecture Notes in Computational Science and Engineering. Springer-Verlag, Berlin.

[16] LANG J & VERWER J. 2001. ROS3P - an accurate third-order Rosenbrock solver designed for parabolic problems. BIT, 41(4): 730-737.

[17] LANG J & TELEAGA D. 2008. Towards a fully space-time adaptive FEM for magnetoquasistatics. IEEE Transactions on Magnet- ics, 44(6): 1238-1241.

18] MARCHANT JS & PARKER I. 2001. Role of elementary Ca2+ puffs in generating repetitive Ca2+ oscillations. The EMBO Journal, 20(1&2): 65-76.

[19] MEANS S, SMITH AJ, SHEPHERD J, SHADID J, FOWLER J, WOJCIKIEWICZ RJH, MAZEL T, SMITH GD & WILSON BS. 2006. Reaction diffusion modeling of calcium dynamics with realistic ergeometry. Biophysical Journal, 91(2): 537-557.

[20] NAGAIAH Ch, RUDIGER S, WARNECKE G & FALCKE M. 2012. Adaptive space and time numerical simulation of reactiondiffusion models for intracellular calcium dynamics. Applied Mathematics and Computation, 218(20): 10194-10210.

[21] OKADA J-i, SUGIURA S, NISHIMURA S & HISADA T. 2005. Threedimensional simulation of calcium waves and contraction in car- diomyocytes using the finite element method. American Journal of Physiology - Cell Physiology, 288(3): C510-C522.

[22] PUTNEY JW & BIRD GSJ. 1993. The inositolphosphate-calcium signaling system in nonexcitable cells. Endocrine Reviews, 14(5): 610-631.

[23] RIDGWAY EB, GILKEY JC & JAFFE LF. 1977. Free calcium increases explosively in activating medaka eggs. Proc. Natl. Acad. Sci. USA, 74: 623-627.

[24] ROONEY TA, SASS EJ & THOMAS AP. 1989. Characterization of cytosolic calcium oscillations induced by phenylephrine and vasopressin in single fura-2-loaded hepatocytes. J. Biol. Chem.,264: 17131-17141.

[25] RUDIGER S, JUNG P & SHUAI JW. 2012. Termination of Ca2+ release for clustered IP3R channels. PLoS Comput. Biol., 8(5): e1002485.

[26] RUDIGER S, NAGAIAH Ch, WARNECKE G & SHUAI JW. 2010. Calcium domains around single and clustered IP3 receptors and their modulation by buffers. Biophysical Journal, 99(1): 3-12.

[27] RUDIGER S, SHUAI JW, HUISINGA W, NAGAIAH Ch, WARNECKE G, PARKER I & FALCKE M. 2007. Hybrid stochastic and determin- istic simulations of calcium blips. Bio-Phys. J., 93: 1847-1857.

[28] VAN DER VORST HA. 1994. Bi-CGSTAB: A fast and smoothly converging variant of bi-cg for the solution of nonsymmetric linear systems. SIAM J. Sci. Stat. Comput., 13: 631-644.