Dual mechanisms of Ca2+ oscillations in hepatocytes

Numerical Analysis of Compound Biochemical Calcium Oscillations Process in Hepatocyte Cells

The hepatocyte cells regulate the wide range of liver function by moderating cellular activities such as lipid, protein metabolism, carbohydrate, and interact with other cells for proliferation and maintenance. In hepatocyte cells, the concentration of calcium uptake is quite extensive from various agonists such as activeG α ${G_\alpha}$ subunit, active phospholipase C, free calcium in the cytosol, and endoplasmic reticulum. The overproduction and degradation of calcium signals can cause homeostasis, liver inflammation, and liver diseases. The spatiotemporal behavior of calcium oscillation reveals the physiological role of these cellular entities in understanding the process of production and degradation. No computational attempt has been registered to date on the compound calcium regulation of these cellular entities including the memory of cells. Hence, the authors proposed a fractional order compartmental model that systematically simulates the exchange of calcium intake in cellular entities. The nonlinear equations of the rate of changes in the activeG α ${G_\alpha}$ subunit, active phospholipase C, free calcium in the cytosol, and endoplasmic reticulum are coupled to form a nonlinear fractional order initial value problem. The existence and uniqueness, stability analysis of the model is performed that validate the theoretical results and explore the dynamic behaviour of calcium oscillation in each compartment.

A Tale of two receptors

Calcium (Ca²⁺) oscillations in hepatocytes have a wide dynamic range. In particular, recent experimental evidence shows that agonist stimulation of the P2Y family of receptors leads to qualitatively diverse Ca²⁺ oscillations. We present a new model of Ca²⁺ oscillations in hepatocytes based on these experiments to investigate the mechanisms controlling P2Y-activated Ca²⁺ oscillations. The model accounts for Ca²⁺ regulation of the IP₃ receptor (IP₃R), the positive feedback from Ca²⁺ on phospholipase C (PLC) and the P2Y receptor phosphorylation by protein kinase C (PKC). Furthermore, PKC is shown to control multiple cellular substrates. Utilising the model, we suggest the activity and intensity of PLC and PKC necessary to explain the qualitatively diverse Ca²⁺ oscillations in response to P2Y receptor activation.

Modeling of Ca2+ transients initiated by GPCR agonists in mesenchymal stromal cells

The integrative study that included experimentation and mathematical modeling was carried out to analyze dynamic aspects of transient Ca²⁺ signaling induced by brief pulses of GPCR agonists in mesenchymal stromal cells from the human adipose tissue (AD-MSCs). The experimental findings argued for IP₃/Ca²⁺-regulated Ca²⁺ release via IP₃ receptors (IP₃Rs) as a key mechanism mediating agonist-dependent Ca²⁺ transients. The consistent signaling circuit was proposed to formalize coupling of agonist binding to Ca²⁺ mobilization for mathematical modeling. The model properly simulated the basic phenomenology of agonist transduction in AD-MSCs, which mostly produced single Ca²⁺ spikes upon brief stimulation. The spike-like responses were almost invariantly shaped at different agonist doses above a threshold, while response lag markedly decreased with stimulus strength. In AD-MSCs, agonists and IP₃ uncaging elicited similar Ca²⁺ transients but IP₃ pulses released Ca²⁺ without pronounced delay. This suggested that IP₃ production was rate-limiting in agonist transduction. In a subpopulation of AD-MSCs, brief agonist pulses elicited Ca²⁺ bursts crowned by damped oscillations. With properly adjusted parameters of IP₃R inhibition by cytosolic Ca²⁺, the model reproduced such oscillatory Ca²⁺ responses as well. GEM-GECO1 and R-CEPIA1er, the genetically encoded sensors of cytosolic and reticular Ca²⁺, respectively, were co-expressed in HEK-293 cells that also responded to agonists in an "all-or-nothing" manner. The experimentally observed Ca²⁺ signals triggered by ACh in both compartments were properly simulated with the suggested signaling circuit. Thus, the performed modeling of the transduction process provides sufficient theoretical basis for deeper interpretation of experimental findings on agonist-induced Ca²⁺ signaling in AD-MSCs.