fluid flow induced calcium response in osteoblasts: mathematical modeling

fluid flow in the bone lacuno–canalicular network can induce dynamic fluctuation of intracellular calcium concentration ([ca²⁺]i) in osteoblasts, which plays an important role in bone remodeling. there has been limited progress in the mathematical modeling of this process probably due to its complexity, which is controlled by various factors such as ca²⁺ channels and extracellular messengers. in this study we developed a mathematical model to describe [ca²⁺]i response induced by fluid shear stress (ss) by integrating the major factors involved and analyzed the effects of different experimental setups (e.g. [ca²⁺]i baseline, pretreatment with atp). in this model we considered the atp release process and the activities of multiple ion channels and purinergic receptors. the model was further verified quantitatively by comparing the simulation results with experimental data reported in literature. the results showed that: (i) extracellular atp concentration has more significant effect on [ca²⁺]i baseline (73% increase in [ca²⁺]i with extracellular atp concentration varying between 0 and 10 μm), as compared to that induced by ss (25% variation in [ca²⁺]i with ss varying from 0 to 3.5 pa); (ii) pretreatment with atp-medium results in different [ca²⁺]i response as compared to the control group (atp-free medium) under ss; (iii) relative [ca2+]i fluctuation over baseline is more reliable to show the [ca²⁺]i response process than the absolute [ca²⁺]i response peak. the developed model may improve the experimental design and facilitate our understanding of the mechanotransduction process in osteoblasts.