1
|
Mathematical modeling of intracellular calcium in presence of receptor: a homeostatic model for endothelial cell. Biomech Model Mechanobiol 2023; 22:217-232. [PMID: 36219362 DOI: 10.1007/s10237-022-01643-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 09/20/2022] [Indexed: 11/02/2022]
Abstract
Calcium is a ubiquitous molecule and second messenger that regulates many cellular functions ranging from exocytosis to cell proliferation at different time scales. In the vasculature, a constant adenosine triphosphate (ATP) concentration is maintained because of ATP released by red blood cells (RBCs). These ATP molecules continuously react with purinergic receptors on the surface of endothelial cells (ECs). Consequently, a cascade of chemical reactions are triggered that result in a transient cytoplasmic calcium (Ca[Formula: see text]), followed by return to its basal concentration. The mathematical models proposed in the literature are able to reproduce the transient peak. However, the trailing concentration is always higher than the basal cytoplasmic Ca[Formula: see text] concentrations, and the Ca[Formula: see text] concentration in endoplasmic reticulum (ER) remains lower than its initial concentration. This means that the intracellular homeostasis is not recovered. We propose, herein, a minimal model of calcium kinetics. We find that the desensitization of EC surface receptors due to phosphorylation and recycling plays a vital role in maintaining calcium homeostasis in the presence of a constant stimulus (ATP). The model is able to capture several experimental observations such as refilling of Ca[Formula: see text] in the ER, variation of cytoplasmic Ca[Formula: see text] transient peak in ECs, the resting cytoplasmic Ca[Formula: see text] concentration, the effect of removing ATP from the plasma on Ca[Formula: see text] homeostasis, and the saturation of cytoplasmic Ca[Formula: see text] transient peak with increase in ATP concentration. Direct confrontation with several experimental results is conducted. This work paves the way for systematic studies on coupling between blood flow and chemical signaling, and should contribute to a better understanding of the relation between (patho)physiological conditions and Ca[Formula: see text] kinetics.
Collapse
|
2
|
Modeling of Endothelial Calcium Responses within a Microfluidic Generator of Spatio-Temporal ATP and Shear Stress Signals. MICROMACHINES 2021; 12:mi12020161. [PMID: 33562260 PMCID: PMC7914997 DOI: 10.3390/mi12020161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Revised: 01/27/2021] [Accepted: 02/04/2021] [Indexed: 11/18/2022]
Abstract
Intracellular calcium dynamics play essential roles in the proper functioning of cellular activities. It is a well known important chemosensing and mechanosensing process regulated by the spatio-temporal microenvironment. Nevertheless, how spatio-temporal biochemical and biomechanical stimuli affect calcium dynamics is not fully understood and the underlying regulation mechanism remains missing. Herein, based on a developed microfluidic generator of biochemical and biomechanical signals, we theoretically analyzed the generation of spatio-temporal ATP and shear stress signals within the microfluidic platform and investigated the effect of spatial combination of ATP and shear stress stimuli on the intracellular calcium dynamics. The simulation results demonstrate the capacity and flexibility of the microfluidic system in generating spatio-temporal ATP and shear stress. Along the transverse direction of the microchannel, dynamic ATP signals of distinct amplitudes coupled with identical shear stress are created, which induce the spatio-temporal diversity in calcium responses. Interestingly, to the multiple combinations of stimuli, the intracellular calcium dynamics reveal two main modes: unimodal and oscillatory modes, showing significant dependence on the features of the spatio-temporal ATP and shear stress stimuli. The present study provides essential information for controlling calcium dynamics by regulating spatio-temporal biochemical and biomechanical stimuli, which shows the potential in directing cellular activities and understanding the occurrence and development of disease.
Collapse
|
3
|
Meng G, Li C, Sun H, Lee I. Multiple calcium patterns of rat osteoblasts under fluidic shear stress. J Orthop Res 2018; 36:2039-2051. [PMID: 29266507 DOI: 10.1002/jor.23843] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/27/2017] [Accepted: 12/14/2017] [Indexed: 02/04/2023]
Abstract
The intracellular calcium ([Ca2+ ]i ) response induced by external forces can be diverse and complex. Using primary osteoblasts from Wistar rats, we found multiple patterns of [Ca2+ ]i responses induced by fluidic shear stress (Fss), including homogeneous non-oscillation and heterogeneous oscillations. These multiple-patterned [Ca2+ ]i responses could be influenced by Fss intensity, cell density, and cell differentiation. Our real-time measurements with free calcium, ATP, ATP without calcium, suramin, apyrase, and thapsigargin confirmed homogeneous [Ca2+ ]i patterns and/or heterogeneous [Ca2+ ]i oscillations with respect to the combined degree of external Ca2+ influx, and intracellular Ca2+ release. Our theoretical model supported diverse Fss-induced calcium activities as well. We concluded that a singular factor of Ca2+ influx or release dominated to produce smooth homogeneous patterns, but combined factors produced oscillatory heterogeneous patterns. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:2039-2051, 2018.
Collapse
Affiliation(s)
- Guixian Meng
- The Key Laboratory of Weak-Light Nonlinear Photonics, Ministry of Education, TEDA Institute of Applied Physics and School of Physics, Nankai University, Tianjin, China.,Academy of Laboratory, Jilin Medical University, Jilin, China
| | - Cunbo Li
- The Key Laboratory of Weak-Light Nonlinear Photonics, Ministry of Education, TEDA Institute of Applied Physics and School of Physics, Nankai University, Tianjin, China
| | - Haiying Sun
- The Key Laboratory of Weak-Light Nonlinear Photonics, Ministry of Education, TEDA Institute of Applied Physics and School of Physics, Nankai University, Tianjin, China
| | - Imshik Lee
- The Key Laboratory of Weak-Light Nonlinear Photonics, Ministry of Education, TEDA Institute of Applied Physics and School of Physics, Nankai University, Tianjin, China
| |
Collapse
|
4
|
Simulation of intracellular $$\hbox {Ca}^{2+}$$ Ca 2 + transients in osteoblasts induced by fluid shear stress and its application. Biomech Model Mechanobiol 2016; 16:509-520. [DOI: 10.1007/s10237-016-0833-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Accepted: 09/09/2016] [Indexed: 10/21/2022]
|
5
|
Vainio I, Abu Khamidakh A, Paci M, Skottman H, Juuti-Uusitalo K, Hyttinen J, Nymark S. Computational Model of Ca2+ Wave Propagation in Human Retinal Pigment Epithelial ARPE-19 Cells. PLoS One 2015; 10:e0128434. [PMID: 26070134 PMCID: PMC4466493 DOI: 10.1371/journal.pone.0128434] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2014] [Accepted: 04/27/2015] [Indexed: 12/21/2022] Open
Abstract
Objective Computational models of calcium (Ca2+) signaling have been constructed for several cell types. There are, however, no such models for retinal pigment epithelium (RPE). Our aim was to construct a Ca2+ signaling model for RPE based on our experimental data of mechanically induced Ca2+ wave in the in vitro model of RPE, the ARPE-19 monolayer. Methods We combined six essential Ca2+ signaling components into a model: stretch-sensitive Ca2+ channels (SSCCs), P2Y2 receptors, IP3 receptors, ryanodine receptors, Ca2+ pumps, and gap junctions. The cells in our epithelial model are connected to each other to enable transport of signaling molecules. Parameterization was done by tuning the above model components so that the simulated Ca2+ waves reproduced our control experimental data and data where gap junctions were blocked. Results Our model was able to explain Ca2+ signaling in ARPE-19 cells, and the basic mechanism was found to be as follows: 1) Cells near the stimulus site are likely to conduct Ca2+ through plasma membrane SSCCs and gap junctions conduct the Ca2+ and IP3 between cells further away. 2) Most likely the stimulated cell secretes ligand to the extracellular space where the ligand diffusion mediates the Ca2+ signal so that the ligand concentration decreases with distance. 3) The phosphorylation of the IP3 receptor defines the cell’s sensitivity to the extracellular ligand attenuating the Ca2+ signal in the distance. Conclusions The developed model was able to simulate an array of experimental data including drug effects. Furthermore, our simulations predict that suramin may interfere ligand binding on P2Y2 receptors or accelerate P2Y2 receptor phosphorylation, which may partially be the reason for Ca2+ wave attenuation by suramin. Being the first RPE Ca2+ signaling model created based on experimental data on ARPE-19 cell line, the model offers a platform for further modeling of native RPE functions.
Collapse
Affiliation(s)
- Iina Vainio
- Department of Electronics and Communications Engineering, Tampere University of Technology, Tampere, Finland
- Institute of Biosciences and Medical Technology, Tampere University of Technology, Tampere, Finland
- * E-mail:
| | - Amna Abu Khamidakh
- Department of Electronics and Communications Engineering, Tampere University of Technology, Tampere, Finland
- Institute of Biosciences and Medical Technology, Tampere University of Technology, Tampere, Finland
| | - Michelangelo Paci
- Department of Electronics and Communications Engineering, Tampere University of Technology, Tampere, Finland
- Institute of Biosciences and Medical Technology, Tampere University of Technology, Tampere, Finland
| | - Heli Skottman
- Institute of Biosciences and Medical Technology, University of Tampere, Tampere, Finland
| | - Kati Juuti-Uusitalo
- Institute of Biosciences and Medical Technology, University of Tampere, Tampere, Finland
| | - Jari Hyttinen
- Department of Electronics and Communications Engineering, Tampere University of Technology, Tampere, Finland
- Institute of Biosciences and Medical Technology, Tampere University of Technology, Tampere, Finland
| | - Soile Nymark
- Department of Electronics and Communications Engineering, Tampere University of Technology, Tampere, Finland
- Institute of Biosciences and Medical Technology, Tampere University of Technology, Tampere, Finland
| |
Collapse
|
6
|
Modelling the coupling between intracellular calcium release and the cell cycle during cortical brain development. J Theor Biol 2014; 347:17-32. [DOI: 10.1016/j.jtbi.2014.01.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2013] [Revised: 11/28/2013] [Accepted: 01/03/2014] [Indexed: 01/28/2023]
|
7
|
Ji L, Cheng M, Zhang H. Stochastic contribution for the coding of agonist induced calcium oscillation in hepatocytes. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2013; 42:607-12. [PMID: 23653093 DOI: 10.1007/s00249-013-0908-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2013] [Revised: 04/13/2013] [Accepted: 04/23/2013] [Indexed: 10/26/2022]
Abstract
The influence of stochastic inositol-1,4,5-trisphosphate receptor (IP₃R) dynamics and their clustering have been extensively investigated to explore the mechanism through which the stochastic molecular event finally shape the intracellular calcium signaling. Most of the previous works employed simplified models which take the concentration of IP₃ instead of that of the agonist as the stimulation intensity. However, the IP₃ level is not linearly dependent on the agonist concentration in stimulus induced signaling systems because there are feedback links in the transduction network. In this work, we include both the IP₃R dynamics and the typical agonist induced signaling transduction cascade in the model to investigate the essential influence of stochastic IP₃R dynamics on the coding of the stimulus induced calcium signal. Simulation results reveal two distinct oscillation areas under different stimulation levels. The signal is optimally modulate by the IP₃R cluster number in the weak stimulated area while affected by the stimulus intensity in the strong stimulated area. Different dependences of coefficient of variance (CV) on the number of clusters are obtained in these two areas, which explains the disagreement in the previous reported results. Besides, the transition between these areas explains the significant CV reduction observed in experiments.
Collapse
Affiliation(s)
- Lin Ji
- Department of Chemistry, Capital Normal University, Beijing 100048, China.
| | | | | |
Collapse
|
8
|
Sun GX, Wang LJ, Xiang C, Qin KR. A dynamic model for intracellular calcium response in fibroblasts induced by electrical stimulation. Math Biosci 2013; 244:47-57. [PMID: 23624256 DOI: 10.1016/j.mbs.2013.04.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2012] [Revised: 04/13/2013] [Accepted: 04/15/2013] [Indexed: 11/27/2022]
Abstract
Regulation of intracellular calcium ion concentration ([Ca(2+)]in) in fibroblasts induced by exogenous electrical stimulation could be applied to control gene expressions selectively which in turn modulate the function of the fibroblasts. Regarding the mechanism for electric-field-induced Ca(2+) influx via voltage-gated Ca(2+) channels and/or stretch-activated cation channels in the fibroblasts, a dynamic mathematical model is proposed to quantify the [Ca(2+)]in dynamics in response to direct current or alternating current electric fields. Simulation results demonstrate that the changes in [Ca(2+)]in predicted by our dynamic model are consistent with the experimental data in the literature. The proposed dynamic model could provide not only more insights into the electric-field-induced intracellular Ca(2+) response but also a quantitative way to regulate the [Ca(2+)]in dynamics by controlling the external electrical stimulation.
Collapse
Affiliation(s)
- Guo-Xin Sun
- Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, No. 2, Linggong Rd., Dalian 116024, China
| | | | | | | |
Collapse
|
9
|
Su JH, Xu F, Lu XL, Lu TJ. Fluid flow induced calcium response in osteoblasts: mathematical modeling. J Biomech 2011; 44:2040-6. [PMID: 21665208 DOI: 10.1016/j.jbiomech.2011.05.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2010] [Revised: 05/03/2011] [Accepted: 05/03/2011] [Indexed: 11/28/2022]
Abstract
Fluid flow in the bone lacuno-canalicular network can induce dynamic fluctuation of intracellular calcium concentration ([Ca(2+)](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(2+) channels and extracellular messengers. In this study we developed a mathematical model to describe [Ca(2+)](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(2+)](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(2+)](i) baseline (73% increase in [Ca(2+)](i) with extracellular ATP concentration varying between 0 and 10 μM), as compared to that induced by SS (25% variation in [Ca(2+)](i) with SS varying from 0 to 3.5 Pa); (ii) Pretreatment with ATP-medium results in different [Ca(2+)](i) response as compared to the control group (ATP-free medium) under SS; (iii) Relative [Ca(2+)](i) fluctuation over baseline is more reliable to show the [Ca(2+)](i) response process than the absolute [Ca(2+)](i) response peak. The developed model may improve the experimental design and facilitate our understanding of the mechanotransduction process in osteoblasts.
Collapse
Affiliation(s)
- J H Su
- Biomedical Engineering and Biomechanics Center, School of Aerospace, Xi'an Jiaotong University, 710049 Xi'an, PR China
| | | | | | | |
Collapse
|
10
|
Jovic A, Howell B, Cote M, Wade SM, Mehta K, Miyawaki A, Neubig RR, Linderman JJ, Takayama S. Phase-locked signals elucidate circuit architecture of an oscillatory pathway. PLoS Comput Biol 2010; 6:e1001040. [PMID: 21203481 PMCID: PMC3009597 DOI: 10.1371/journal.pcbi.1001040] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2010] [Accepted: 11/27/2010] [Indexed: 12/26/2022] Open
Abstract
This paper introduces the concept of phase-locking analysis of oscillatory cellular signaling systems to elucidate biochemical circuit architecture. Phase-locking is a physical phenomenon that refers to a response mode in which system output is synchronized to a periodic stimulus; in some instances, the number of responses can be fewer than the number of inputs, indicative of skipped beats. While the observation of phase-locking alone is largely independent of detailed mechanism, we find that the properties of phase-locking are useful for discriminating circuit architectures because they reflect not only the activation but also the recovery characteristics of biochemical circuits. Here, this principle is demonstrated for analysis of a G-protein coupled receptor system, the M3 muscarinic receptor-calcium signaling pathway, using microfluidic-mediated periodic chemical stimulation of the M3 receptor with carbachol and real-time imaging of resulting calcium transients. Using this approach we uncovered the potential importance of basal IP3 production, a finding that has important implications on calcium response fidelity to periodic stimulation. Based upon our analysis, we also negated the notion that the Gq-PLC interaction is switch-like, which has a strong influence upon how extracellular signals are filtered and interpreted downstream. Phase-locking analysis is a new and useful tool for model revision and mechanism elucidation; the method complements conventional genetic and chemical tools for analysis of cellular signaling circuitry and should be broadly applicable to other oscillatory pathways. Key to robust discernment of cell circuit architecture is to have as many distinct response features as possible for comparison and evaluation. One under-appreciated characteristic of oscillatory circuits is that under periodic stimulation, these systems will exhibit responses synchronized to this stimulatory input, a phenomenon termed phase-locking. We demonstrate that phase-locked response characteristics vary noticeably depending on circuit activation and recovery properties; these response characteristics thereby provide a unique set of criteria for oscillatory circuit architecture analysis. The concept is validated through experiments on an oscillatory calcium pathway in mammalian cells; the experimental setup allowed us to explore, for the first time, the properties of chemically induced phase-locking of intracellular signals. Observations of this phenomenon were then used to test the predictions of several existing mathematical models of calcium signaling. Most of the models we evaluated were unable to match all our experimental observations, suggesting that current models are missing mechanistic elements in the context of calcium signaling for the cell type and receptor/stimulant tested. The observations of phase-locking further led us to identify one simple mechanistic modification that would account for all the experimental observations. The techniques and methodology presented should be broadly applicable to a variety of biological oscillators.
Collapse
Affiliation(s)
- Andreja Jovic
- Biomedical Engineering Department, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Bryan Howell
- Biomedical Engineering Department, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Michelle Cote
- Biomedical Engineering Department, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Susan M. Wade
- Pharmacology Department, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Khamir Mehta
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Atsushi Miyawaki
- Laboratory for Cell Function and Dynamics, Advanced Technology Development Center, Brain Science Institute, Wako City, Saitama, Japan
| | - Richard R. Neubig
- Pharmacology Department, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Jennifer J. Linderman
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan, United States of America
- * E-mail: (ST); (JJL)
| | - Shuichi Takayama
- Biomedical Engineering Department, University of Michigan, Ann Arbor, Michigan, United States of America
- Macromolecular Science and Engineering Department, University of Michigan, Ann Arbor, Michigan, United States of America
- * E-mail: (ST); (JJL)
| |
Collapse
|