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

Nikolai P Kaimachnikov; Polina D Kotova; Ekaterina N Kochkina; Olga A Rogachevskaja; Alexander A Khokhlov; Marina F Bystrova; Stanislav S Kolesnikov

doi:10.1016/j.bbadva.2021.100012

Modeling of Ca²⁺ transients initiated by GPCR agonists in mesenchymal stromal cells

BBA Adv. 2021 Apr 21:1:100012. doi: 10.1016/j.bbadva.2021.100012. eCollection 2021.

Authors

Nikolai P Kaimachnikov¹, Polina D Kotova¹, Ekaterina N Kochkina¹, Olga A Rogachevskaja¹, Alexander A Khokhlov¹, Marina F Bystrova¹, Stanislav S Kolesnikov¹

Affiliation

¹ Institute of Cell Biophysics, Pushchino Scientific Center for Biological Research of the Russian Academy of Sciences, 3 Institutskaya Street, Pushchino, Moscow Region, 142290, Russia.

Abstract

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.

Keywords: Agonist-induced Ca2+ signaling; Ca2+-induced Ca2+ release; IP3 receptors; Mathematical modeling; Mesenchymal stromal cells.