Abstract | ||
---|---|---|
Simulation of electrophysiology and intracellular Ca2+ dynamics in cardiomyocytes comprises fast stochastic dynamics in tiny subcompartments, partial differential equations (PDEs) with stochastic source terms for concentration fields, and the globally coupling membrane potential. We use highly unstructured meshes appropriate for the spatial heterogeneity of intracellular Ca2+ release, and adaptive time-stepping algorithms appropriate for the simulation of stochastic channel opening and closing in subcompartments like the Ca2+ release units (CRUs). A set of reaction-diffusion equations describes the behavior of the intracellular concentration fields on length scales from tens of nanometers to cell size (tens of micrometers) and milliseconds to tens of seconds. Detailed highly stochastic CRU models drive source functions in the PDE model. These CRU models cover dynamics with time scales below 1 ms and length scales from a few to a few hundred nm. The spatially detailed Ca2+ dynamics and the cardiomyocyte membrane potential interact. Membrane potential introduces a global spatial coupling across the whole cell due to its large coupling length. Its dynamics are consequently represented by a set of ordinary differential equations (ODEs). We developed an efficient adaptive finite element simulator interface for the numerical simulation of this multiphysics and multiscale problem. The use of stationary Green functions within the CRU models and highly unstructured meshes for the PDE integration allows for bridging of many orders of magnitude of spatial scale, to represent accurately the Ca2+ concentration dynamics from within a single CRU up to the level of the whole cell. The time scale separation between fast stochastic CRU dynamics and slower PDE dynamics limits efficiency with traditional approaches. We present new methods to circumvent that problem. We demonstrate large-scale numerical results for a 436 CRU cellular subdomain using many-core parallel machines. |
Year | DOI | Venue |
---|---|---|
2018 | 10.1137/17M1121639 | MULTISCALE MODELING & SIMULATION |
Keywords | Field | DocType |
calcium cycling,multiscale modeling,calcium release units,membrane potential model,reaction-diffusion equations,hybrid algorithm,finite element method,adaptive Runge-Kutta methods,high-performance computing | Myocyte,Calcium,Computer simulation,Stochastic dynamics,Mathematical analysis,Finite element method,Multiscale modeling,Reaction–diffusion system,Partial differential equation,Mathematics | Journal |
Volume | Issue | ISSN |
16 | 3 | 1540-3459 |
Citations | PageRank | References |
0 | 0.34 | 2 |
Authors | ||
6 |
Name | Order | Citations | PageRank |
---|---|---|---|
Chamakuri Nagaiah | 1 | 16 | 3.50 |
Wilhelm Neubert | 2 | 0 | 0.34 |
Stephen H. Gilbert | 3 | 15 | 4.55 |
Janine Vierheller | 4 | 0 | 0.34 |
G. Warnecke | 5 | 85 | 17.57 |
M Falcke | 6 | 13 | 2.75 |