T cell Ca²⁺ microdomains through the lens of computational modeling

Abstract

Cellular Ca²⁺ signaling is highly organized in time and space. Locally restricted and short-lived regions of Ca²⁺ increase, called Ca²⁺ microdomains, constitute building blocks that are differentially arranged to create cellular Ca²⁺ signatures controlling physiological responses. Here, we focus on Ca²⁺ microdomains occurring in restricted cytosolic spaces between the plasma membrane and the endoplasmic reticulum, called endoplasmic reticulum-plasma membrane junctions. In T cells, these microdomains have been finely characterized. Enough quantitative data are thus available to develop detailed computational models of junctional Ca²⁺ dynamics. Simulations are able to predict the characteristics of Ca²⁺ increases at the level of single channels and in junctions of different spatial configurations, in response to various signaling molecules. Thanks to the synergy between experimental observations and computational modeling, a unified description of the molecular mechanisms that create Ca²⁺ microdomains in the first seconds of T cell stimulation is emerging.

Ca²⁺ signaling research: from global to local

Ca²⁺ signaling is the major intracellular signaling system found in almost all cell types in the animal kingdom, but also extending to plant cells and fungi (1, 2). Ca²⁺ signaling requires Ca²⁺ gradients across organellar membranes and across the plasma membrane (PM). ATP-driven Ca²⁺ pumps, such as sarco-/endoplasmic reticular Ca²⁺ ATPase (SERCA) and plasma membrane Ca²⁺ ATPase (PMCA), build up such Ca²⁺ gradients allowing for rapid Ca²⁺ release from organellar Ca²⁺ stores and/or Ca²⁺ entry across the plasma membrane upon cell stimulation.

Global Ca²⁺ signaling in the form of an increase in the free cytosolic Ca²⁺ concentration ([Ca²⁺]_i) throughout the cytosol, most often extending into the nucleus, was first described in the 1980s (3). Such experiments were first performed in cell suspensions analyzed in a fluorimeter; however, with the advent of microscopic imaging technology, global Ca²⁺ signaling was also observed in single cells (4). With the increasing spatio-temporal resolution of microscopes and cameras, analyses of the dynamics of local Ca²⁺ signaling has become possible. In particular, the advent of confocal microscopy and the development of super-resolution imaging methods were major steps to allow for very high resolution in both space and time. Further, novel deconvolution approaches allow to increase signal-to-noise ratio even in dim raw images, adding another important step to disentangle fast local Ca²⁺ signaling events (5).

Recently, we reviewed functions and mechanisms of formation of Ca²⁺ microdomains (6). In summary, Ca²⁺ microdomains in which the rise in concentration reaches a few hundred nanomolar and extending on 1-3 microns were found in different cell types. Pioneering work by Peter Lipp, Martin Bootman and colleagues revealed elementary Ca²⁺ signal events, likely representing Ca²⁺ release from clusters of either ryanodine receptors (RYR), termed ‘sparks’ or d-myo-inositol 1,4,5-trisphosphate receptors (IP₃R), termed ‘puffs’. Smaller fundamental Ca²⁺ signal events that relate to Ca²⁺ release from single RYR, termed ‘quarks’ or IP₃R, termed ‘blips’ were also described (7–10). In addition to these local Ca²⁺ release events, Ca²⁺ sparklets were defined as local signals resulting from brief openings of plasma membrane Ca²⁺ channels (11–13).

More recently, Ca²⁺ microdomains occurring in confined spaces between closely apposed membranes have been described. Locally elevated Ca²⁺ concentrations arise around mitochondria-associated ER membranes (MAM) that provide privileged platforms for Ca²⁺ transfer from the ER into mitochondria (14). Similarly, sharp and short-lived Ca²⁺ increases can sometimes be observed at the ER-PM contact sites or junctions. These signals are dependent on Ca²⁺ entry from the extracellular medium through store-operated Ca²⁺ entry (SOCE). The latter mechanism relies on the stromal interaction molecules (STIM1 and/or STIM2) that are ER Ca²⁺ sensors regulating the activity of the Orai1 PM Ca²⁺ channels. Thus, upon local depletion of the ER Ca²⁺ store, Ca²⁺ dissociates from STIM, which allows recruitment of Orai1 to the junction and their gating (15–18). These channels are also referred to as “Ca²⁺ release activated channels (CRAC)”. The ER-PM junctions are areas in the cytosol where the ER gets really close to the PM, with a depth of approximately 15 nm and an extension of ~200 nm (18).

As proposed by 19, Ca²⁺ microdomains can be defined as “subcellular, rapid and localized high Ca²⁺ concentration regions that develop near open Ca²⁺ channels at the plasma membrane or internal stores”. Key factors that determine their specificity in triggering different signaling pathways in different cell types are their amplitude, frequency and duration. Variability in these factors leads to specific forms of global Ca²⁺ increases. The detailed characterization of Ca²⁺ microdomains is limited by the spatio-temporal resolution of available imaging techniques. Computational modelling has therefore much been used to further investigate their molecular origin. These studies rely on the theoretical tools that have been developed to describe global Ca²⁺ signaling (20), although they face specific problems related to the confined geometry and the steep gradients that need to be simulated (18). Moreover, stochastic approaches are sometimes required to take into account the local fluctuations in concentrations due to the small numbers of Ca²⁺ channels present in the microdomains (21–25).

Local Ca²⁺ entry via SOCE activates different cellular responses in different cell types depending on their spatio-temporal profiles. These include transcription factor stimulation in B cells by high amplitude Ca²⁺ spikes, NFAT activation in T cells by a sustained low amplitude signal and exocytosis stimulation in non-excitable cells by spontaneous repetitive short Ca²⁺ signals. It is known that the modulation of such localized Ca²⁺ signals involves clusters of ER (IP₃R or RYR) or PM (Orai1) Ca²⁺ channels, but the detailed interplay between these channels, including their molecular effectors, remain to be fully characterized.

Ca²⁺ microdomains in T cells

Although most of the Ca²⁺ toolbox molecules, e.g. second messenger forming and degrading enzymes, Ca²⁺ mobilizing second messengers and their receptors, Ca²⁺ channels, Ca²⁺ storing organelles, and Ca²⁺ pumps are often ubiquitously expressed, the specific configuration of a particular cell type determines the exact mechanism of both local and global Ca²⁺ signaling. Thus, in this review, we will concentrate on a particular cell type that has been extensively characterized regarding mechanisms involved in the formation of Ca²⁺ microdomains, mammalian T-lymphocytes. The main molecules involved are listed in Table 1, for ease of reading. Over the past 10 years, Ca²⁺ microdomains at the ER-PM junctions have been extensively studied in these cells. Despite disadvantages of this cell model, e.g. small diameter of about 6-7 μm and a large nucleus as compared to the cytosol, as well as spherical shape with many stacked z-layers that potentially result in image blurring, development of a high-resolution Ca²⁺ imaging method combined with image refining by off-line deconvolution resulted in spatio-temporal resolution of 368 nm and 25 ms (26). Further, activation of the T cell receptor (TCR)/CD3 complex and co-stimulation of CD28 by small beads coated with monoclonal antibodies against CD3 and CD28 allowed for point-shaped activation of few TCR/CD3 complexes, thereby mimicking activation during formation of an immune synapse (26).

So far, two different signaling mechanisms involved in formation of Ca²⁺ microdomains in T cells were identified: (i) adhesion-dependent Ca²⁺ microdomains (ADCM) evoked by integrin-mediated IP₃ signaling and subsequent SOCE (Figure 1A), or (ii) Ca²⁺ microdomains evoked by TCR/CD3/CD28 stimulation (TDCM), proceeding via nicotinic acid adenine dinucleotide phosphate (NAADP) formation by dual NADPH oxidase 2 (DUOX2), and activation of RYR1 via the NAADP receptor hematological and neurological expressed 1-like protein (HN1L)/jupiter microtubule associated homolog 2 (JPT2). Hereafter, ADCM and TDCM will refer to the spatially restricted and short-lived Ca²⁺ increases occurring in the so-called microdomains. TDMC are amplified by SOCE (27) (Figure 1B) and purinergic signaling via P2X4 and P2X7 (28).

Figure 1

Long before the observation of ADCM, Weissmann and colleagues described global Ca²⁺ responses in Jurkat T cells that were triggered by integrin-mediated adhesion to extracellular matrix (ECM) proteins (29). These Ca²⁺ signals were dependent on the phosphorylation of focal adhesion kinase (FAK). FAK is thought to be a key component of the integrin-mediated signaling cascade (30) and evokes Ca²⁺ signals selectively at membrane rafts during the cell–matrix adhesion process (31). Stimulation of FAK subsequently activates phospholipase C-γ (PLC-γ) (32), which in turn generates IP₃ leading to a Ca²⁺ release from the ER by targeting IP₃ receptors (33). Local decrease in ER Ca²⁺ concentration then activates SOCE. Indeed, locally restricted Ca²⁺ microdomains in the absence of TCR stimulation were shown to be dependent on the expression of Orai1 and STIM1/2 (27). Interestingly, the NAADP/RYR1 axis was not involved in this process in primary T cells (27). However, another Ca²⁺ mobilizing second messenger, termed cyclic ADP-ribose (cADPR), which evokes Ca²⁺ signals through ryanodine receptor type 3 (34), was involved in global Ca²⁺ responses mediated by integrins in Jurkat T cells (35). Purinergic signaling via P2X4, is also involved in the formation of TCR independent Ca²⁺ microdomains in primary T cells (28). However, whether integrin mediated Ca²⁺ signaling evokes ATP release followed by purinergic signaling via P2X4 in an autocrine manner is not yet clear. Concerning the physiological role of adhesion-dependent Ca²⁺ microdomains, it has been suggested that they may facilitate full TCR-mediated activation (36) given that T cells form adhesive interactions with endothelial cells or ECM proteins during their migration from blood vessels into inflamed tissues (37).

TDCM require functional NAADP signaling in T cells (Figure 1B). A crucial role of NAADP in T cells was discovered by showing global and local Ca²⁺ signaling upon NAADP microinjection (38, 39; Dammermann & Guse; 2005; 26, 40). RYR was identified as a major Ca²⁺ channel in NAADP signaling in T cells (39; 41; 26, 27, 40). In addition, NAADP evoked Ca²⁺ entry secondary to Ca²⁺ release was also shown for the first time in T cells (39). The mechanism of NAADP evoked release of Ca²⁺ from intracellular Ca²⁺ stores started to being disentangled by the discovery of small cytosolic NAADP receptor/binding proteins, firstly as 22/23 kDa band labeled by a NAADP-photoaffinity probe (42) and finally by identifying the NAADP receptor/binding protein as HN1L/JPT2 in both T cells (43) and erythrocytes (44). Of note, Lsm12 was identified as another NAADP binding protein (45). HN1L/JPT2 co-localized with RYR in T cells as shown by STED super-resolution microscopy and co-immunoprecipitation (43), suggesting that NAADP bound to HN1L/JPT2 directly activates RYR, as hypothesized in the unifying hypothesis for NAADP action (46). The latter also includes, as further potential link between NAADP and activation of ion channels, the possibility of activation of two-pore channels (46). In addition to the Ca²⁺ releasing mechanism of NAADP, formation and degradation of NAADP are rather important for this signaling pathway. Along these lines, it was reported for T cells that NAADP formation within seconds upon TCR/CD3 stimulation (47) proceeds by oxidation of NAADPH to NAADP, catalyzed by DUOX2 (48). Rapid reduction of NAADP to NAADPH may be catalyzed in intact cells by glucose 6-phosphate dehydrogenase, that, together with DUOX2, constitutes a redox cycle for fast formation of NAADP and its ‘storage’ as inactive reduced NAADPH (48). Full degradation of NAADP is catalyzed by CD38 (to 2-phospho-ADPR) or by alkaline phosphatase (to NAAD) (49, 50). Crucial importance of NAADP signaling for Ca²⁺ microdomain formation in T cells was demonstrated by gene deletion or silencing of Duox2, Hn1l/Jpt2, or Ryr1 (26, 27, 43, 48).

Gene deletion of the major proteins of SOCE in T cells, STIM1 and 2, and Orai1, revealed that in addition to NAADP evoked release of Ca²⁺, also SOCE is essentially required for the formation of the local Ca²⁺ microdomains (27). Very recently, another amplification system for Ca²⁺ increases inside microdomains in T cells was discovered as autocrine ATP evoked activation of purinergic Ca²⁺ channels P2X4 and P2X7 (28).

TDCM are observed as soon as a few hundreds of millisecond post TCR/CD3/CD28 stimulation and last for approx. 15 to 25s (Figure 1B). Among the other initial signaling pathway active in T cells in this very early period of time, tyrosine phosphorylation of target proteins is certainly outstanding, e.g. by the tyrosine kinases ZAP70 and/or p56^lck (51). Experimental evidence for their involvement in TDCM was obtained: NAADP formation evoked by TCR/CD3 stimulation was completely blocked by tyrosine kinase inhibitor genestein (47). In contrast, during the first 15s, levels of IP₃ and diacylglycerol, both products formed from phosphatidylinositol 4,5-bisphosphate by phospholipase C-γ (52) are still quite low (53, 54) and do not significantly impact on Ca²⁺ microdomains, as demonstrated, at least for IP₃ evoked Ca²⁺ release, by IP₃ antagonism (26). Later on, between approx. 20 and 30s post TCR/CD3/CD28 stimulation, TDCM merge into global Ca²⁺ signaling. Major mechanisms involved in global Ca²⁺ signaling are (i) IP₃ evoked Ca²⁺ release between approx. 30s and 10 to 15 min (53), (ii) cyclic ADP-ribose mediated Ca²⁺ release started between 5 and 10 min and lasting for at least 60 min (55). SOCE via STIM1 and Orai1 is the major Ca²⁺ entry mechanism secondary to ER Ca²⁺ store depletion by IP₃ and cyclic ADP-ribose [reviewed in (56).

ER-PM Ca²⁺ microdomains in T cells: modeling assumptions

As detailed above, current techniques in microscopy allow to analyze both global Ca²⁺ dynamics and localized Ca²⁺ events. These observations however raise numerous questions about the molecular mechanisms involved in the foundation of Ca²⁺ microdomains. Some of these questions are best addressed by computational modelling, which is made possible thanks to an impressive quantitative characterization of the properties of ER-PM Ca²⁺ microdomains. In particular, models allow to predict the number and spatial arrangement of the channels inside the junction, which remains inaccessible with current experimental techniques.

A few studies have focused on modelling Ca²⁺ microdomains due to SOCE. In a series of pioneering studies (57–59), Ca²⁺ changes at the immunological synapse have been modelled in 1D and later in 3D. These early models focused on the impact of mitochondria on the Ca²⁺ influx through Orai, but most of them did not specifically address the detailed conditions related to the small-scale microdomains in the ER-PM junctions. To assess the impact of mitochondria relocalization to the immunological synapse on the spatial distribution of Ca²⁺ inside the cell, (60) developed a whole cell model considering specific geometries for mitochondria. The model predicted that the experimentally observed locations of this Ca²⁺-exchanging organelle, moving close to the PM, can account for the increase in Ca²⁺ around a 100 nm radius cluster of ORAI channels forming a Ca²⁺ microdomain. In the simulations, such local Ca²⁺ increases occur in conditions of massive ER Ca²⁺ depletion that correspond to situations encountered a few minutes after T cell stimulation, when mitochondria start to move along the cytoskeleton. In this study, the channel locations in the PM, ER or mitochondrial membranes are not specified and Ca²⁺ fluxes are homogeneously distributed across these membranes. The detailed impact of Ca²⁺ channel locations on the spatio-temporal profiles of Ca²⁺ changes inside the ER-PM junctions were analyzed in 1D by (18) and modelled in 3D by (61, 62) and (63, 64). In this section, we describe these models and their main relying assumptions, while their physiological outcomes are presented in the next section.

Model building can be divided in three steps. First, the biologically relevant variables that the model aims at simulating must be selected. Second, the spatial geometry in which these variables are evolving have to be defined, including the locations of organelles, channels, pumps, transporters, etc. Finally, well-chosen rate equations establish how the selected variables evolve in time in the elected sub-cellular configuration. As compared to global cell model, simulations of microdomains call for specificities especially when located around the mouth of a channel (20).

Concerning the choice of the variables, the formation of Ca²⁺ microdomains involve Ca²⁺ movements in the junctions, in the ER portion just beneath the junction (triggering SOCE) and in the extracellular medium. Given that the latter can be seen as an infinite source of Ca²⁺, its concentration can be considered as constant. Thus, Ca²⁺ concentrations in the junction and in the sub-PM ER compartment represent key variables. Because they dynamically regulate the junctional Ca²⁺ concentration, buffers should in principle be taken into account. Most models however ignore the effect of buffers because of the short time and length scales involved in microdomain formation. Because of the high Ca²⁺ concentrations and the small volume of the junction, local Ca²⁺ buffers are expected to be saturated (62, 65, 66). Similarly, Ca²⁺ diffusion is assumed to occur in an unbuffered cytosol (D_c = 220 μm²s^-1). When considering the relation between cell adhesion/stimulation and Ca²⁺ signaling, IP₃ and NAADP should also be involved because they modify ER and junctional Ca²⁺ by affecting the opening states of the IP₃Rs and RYRs, respectively. Although NAADP mobilizes lysosomal Ca²⁺ (67), the involvement of this organelle in the formation of T cells Ca²⁺ microdomains remains to be established (68). Instead of considering [IP₃] and [NAADP] explicitly, their effects can be modelled by changing the opening states of their receptors directly (63, 64). This simplification is however not valid when considering the transition from local to global Ca²⁺ signaling which involves considering diffusion and metabolism of the messengers. These variables and their interplays are schematically shown in Figure 1.

To predict the exact spatio-temporal dynamics underlying Ca²⁺ increases in microdomains, Ca²⁺ fluxes and diffusion must be simulated in a realistic computational description of the ER-PM junction configuration. A typical configuration, used with slight variations in (61, 62) and 63); 64), is schematized in Figure 2A. ER-PM junctions are regions of the cell where both membranes are closely apposed and correspond to the locations where STIM and Orai co-localize upon store depletion. In resting T cells, STIM2 pre-formed complexes with Orai localize at the ER-PM junctions (27). The distance between the two membranes has been estimated to be 10-20 nm in depth and ~200 nm in length (15, 18, 69, 70). This flat cylinder defines the junctional space, surrounded by the cytoplasmic space where Ca²⁺ concentration is assumed to remain constant given the limited amount of Ca²⁺ ions involved in the formation of microdomain. Similarly, the luminal space is arbitrarily divided in two regions: one close to the membrane facing the ER in which [Ca²⁺] changes dictate SOCE dynamics and one further away where [Ca²⁺] is assumed to remain constant due to fast replenishment from the rest of the ER.

Figure 2

62) (A) Three-dimensional schematic view of the model geometry. (B) Ca²⁺ profiles at the ER membrane after 1 ms opening of 5 Orai channels under an Orai clustered (left) configuration with an inter-channel distance of 39 nm and a non-clustered (right) one with a 65 nm distance. (C) Ca²⁺ spatial signature along the dotted line in B for the clustered (blue line) and non-clustered (green line) configuration. Small decreases in junctional Ca²⁺ concentration can be seen near SERCAs. (D) Ca²⁺ spatial signature along the dotted line at the PM for the clustered and non-clustered configuration. All numerical simulations, initially run by McIvor et al. (62) using Green’s functions, are reproduced with COMSOL Multiphysics.

Five Orai channels are placed in the PM, as estimated for Jurkat T cells (18). They are arranged on a ring with an inter-channel distance of 47 nm (61). It should be noticed that the number of Orai channels in the junction is cell-dependent, as one channel per junction is reported for HeLa cells (71). The number of SERCA pumps present in the junction remains unknown and was arbitrarily set to 10 in in (61, 62 and (63, 64). In the whole cell model of (60), 10⁵ SERCA’s spread over the entire surface area of the ER, which would also correspond to ~10 pumps per ER-PM junction assuming a homogenous distribution. However, model outcomes are practically independent of the numbers of SERCA’s that are found to exert little influence on the spatio-temporal dynamics of microdomain formation. In line with this finding, SERCAs were not included in the first studies of Ca²⁺ changes at the immunological synapse (57–59). IP₃Rs and RYRs with a well-defined spatial distribution were explicitly considered by Gil et al. (63, 64), because of the focus of these studies on SOCE possibly induced by moderate, local increases of IP₃ and NAADP. The spatial arrangement of the IP₃Rs (see Figure 3A) mimics the observations of (72) in HeLa cells, where clusters of IP₃Rs reside alongside the junction and facing the PM. In the absence of specific information, the same arrangement was also considered for RYRs. Indeed, when RYRs are placed inside the junction, which would correspond to the situation encountered in the cardiac dyadic clefts (73), simulated Ca²⁺ microdomains do not correspond to observations [see next section, (64).

Figure 3

63). (B) Evolution of the average Ca²⁺ concentration reached in the ER-PM junction with the number of open IP₃Rs (green line) or RYRs (red line). SOCE relies on STIM2/2 homotetramers when local depletion in ER Ca²⁺ results from IP₃Rs opening and on STIM1/2 heterotetramers when local depletion in ER Ca²⁺ results from RYRs opening. Three IP₃Rs are needed to reproduce experimental observed adhesion-dependent Ca²⁺ microdomains (purple diamond) and seven RYRs are needed to reproduce experimental observed TCR/CD3-dependent Ca²⁺ microdomains (blue circle) showed in Diercks et al. (27). Simulations are performed as in Gil et al. (63, 64) and (2022) (C) Ca²⁺ profiles at the PM in the proposed T cell transition mechanism between quiescent adherent cells (t < 0, unstimulated T cell) to early activation (0 < t < 15s, first 15 seconds after stimulation). See text for detail. Simulations are performed as in Gil et al. (63, 64).

How Ca²⁺ ions evolve in space and time in the ER-PM junction is fixed by partial differential equations (PDE), or by a stochastic reaction-diffusion scheme such as, for example, the Gillespie’s algorithm (74). Given a basal [Ca²⁺] in the 10-100 nM range and a volume of ~10^-17 - 10^-15 liter, the number of Ca²⁺ ions in the ER-PM junction is between 0.06 and 60. Because the magnitude of the concentration fluctuations is inversely proportional to the number of ions, it is expected that stochastic fluctuations arising from molecular interactions and diffusion do not average out, thus calling for the requirement of stochastic simulations. All the studies mentioned here above however used PDE’s. This huge simplification is based on previous studies that have compared stochastic and deterministic modelling of Ca²⁺ dynamics in other confined cellular subspaces, such as a dendritic spine or a cardiac dyadic cleft. These works have indeed shown that for an appropriate range of physiological relevant parameters, stochastic Ca²⁺ signaling is well described using a deterministic model (75). This conclusion holds as long as the gating of the channels is not influenced by the ions that they transport, as it would be the case in the presence of Ca²⁺-induced Ca²⁺ release (76–78). Fast and mobile buffers also increase the magnitude of the fluctuations in the microdomain (79). When modeling Ca²⁺ microdomains occurring at the ER-PM junction on time periods of a few seconds following T cell stimulation, these conditions are met, and a deterministic description is appropriate. The resolution of the PDE’s remains challenging for several reasons, including the steep concentration gradients or the existence of widely different time scales. The numerous computational studies of Ca²⁺ “puffs” and “sparks” (21–24, 80, 81) have also paved the way for the modelling of ER-PM Ca²⁺ microdomains. However, for puffs and sparks, there is no barrier for the diffusion of Ca²⁺ flowing through the ER Ca²⁺ channels. In contrast, Ca²⁺ microdomains in dyadic clefts share with ER-PM junctions the characteristic of being embedded between two closely apposed membranes and many conceptual problems related to the deterministic Ca²⁺ dynamics in such restricted spaces have been formally addressed in this case (75, 82). It is also possible to resort to compartmental model. Wieder et al. (78) have summarized the evolution of numerical simulation techniques developed until recently, which accurately approach the different spatial and temporal scales found in the analysis of Ca²⁺ microdomain formation. For T cell Ca²⁺ microdomains, McIvor et al. (62) resorted to Green functions while Gil et al. (63, 64) took profit of the high accuracy of the COMSOL Multiphysics software (83) [COMSOL Multiphysics^®]. Thanks to its powerful mesh dissection, solution methods and abundant post process operation, this software offers a trustable and relatively easily tractable way of modelling multiscale cellular processes.

Modelling of the formation of T cell Ca²⁺ microdomains: interplay between Orai1, IP₃R and RYR1

Together with the increasingly accurate observations of Ca²⁺ dynamics in the ER-PM junctions, modelling studies have allowed to dissect Ca²⁺ signaling in these localized domains, and in particular the nature, number and spatial arrangement of the Ca²⁺ transporting channels that are at play at various stages following T cell stimulation. Samanta et al. [2015] were the first to provide a quantitative analysis of the Ca²⁺ profiles upon SOCE activation in a realistic 3D ER-PM junction. They showed that Orai clustering enables the Ca²⁺ sensors located on the ER surface to be exposed to Ca²⁺ levels considerably higher in the junction than in the nearby cytoplasm. In the same line, (62) focused on the Ca²⁺ spatial signatures in the ER-PM junction created by two different spatial arrangements of Orai channels and their impact in the ER refilling process: firstly, with an inter-channel distance of 39 nm (clustered) and secondly, of 65 nm (non-clustered), as shown in Figure 2B. They concluded that the spatial arrangement of the Orai channels significantly impacts Ca²⁺ changes in the junction. When clustered, junctional Ca²⁺ concentration rises from an initial value of 100 nM to 7 μM in the center of the junction, while this increase is limited to 4 μM in the non-clustered configuration (Figure 2C). In contrast, these changes do not impact on the rates of ER refilling because of the low K_D of SERCA pumps (~300 nM). Orai channels are assumed to remain open with a constant current of 2.1 fA (18) across each channel. At the mouth of the Orai channel Ca²⁺ concentrations around 60 μM are predicted (Figure 2D). The Ca²⁺ junctional concentration reaches its steady state already after 0.25ms, which highlights the temporal reactivity created by this spatial arrangement of channels and membranes.

Ca²⁺ microdomains in T cells appear under conditions of a full ER opposite to the conditions explored by (62). Yet they appear to be opposed to the PM suggesting that they originate directly at ER-PM junctions. Thus, their computational investigation has been performed in the same geometry as that proposed in (62) and schematized in Figure 2A. Diercks et al. (27) quantitatively characterized the two different types of Ca²⁺ increases in microdomains that precede global Ca²⁺ signaling in T cells. The first ones, termed ADCM, evoked by an interplay between IP₃R and/or SOCE and triggered by adhesion, are characterized by a duration of 44 ± 4 ms and a junctional average amplitude of 290 ± 12 nM. The second ones, termed TDCM, relying on RYR and SOCE and triggered by early TCR/CD3 stimulation, are much more frequent and have a duration of 64 ± 3 ms and a junctional average amplitude of 325 ± 11 nM [see 6, for review].

Simulation of the formation of adhesion-dependent Ca²⁺ microdomains

In weakly adhesive, but otherwise unstimulated T cells, STIM2 largely co-localizes with Orai1, which suggests that STIM2 might be pre-clustered and ready to activate Orai1. Considering the high sensitivity of this STIM isoform to small changes in luminal Ca²⁺ concentration (84), it is highly plausible that SOCE contributes to the generation of ADCM. Still, our simulations predicted that Orai1/STIM2 activation in conditions of full ER cannot account on its own for the observed increase in junctional Ca²⁺ (63) indicating that probably an ER Ca²⁺ channel is also involved. Given that RYR1s are not involved in ADCM (27), an interplay between ORAI1 and IP₃R, as suggested in (72), was proposed. This was found to be in agreement with the observation that adhesion triggers the phosphorylation of FAK and hence PLC activation and IP₃ synthesis (27). To computationally evaluate the possible involvement of IP₃Rs in the onset of ADCM, the eight IP₃Rs located around the junction are opened gradually one at a time in a deterministic way at given time points (Figure 3B, green curve). Depending on the luminal Ca²⁺ concentration around its closest IP₃R, each Orai1 opens in one of its four conductance states. This dependency of the open state on ER Ca²⁺ concentration follows a Hill function determined by STIM2 Ca²⁺ affinity. In order to reproduce ADCM observed experimentally (27), three open IP₃Rs are needed together with a cluster of five Orai1 channels bound to STIM2 (Figure 3B, purple diamond). The cluster of Orai1 channels or of IP₃R alone cannot account for the observed magnitude of the Ca²⁺ rise in the microdomain. Another interesting prediction of the model relates to the effect of changing the Ca²⁺ diffusion coefficient inside the ER. By default, a low value for this coefficient (D_ER = 10 μm²s^-1) was considered, to account for molecular crowding in the lumen and for the tortuosity of the tubular ER network (81). When D_ER = 100 μm²s^-1 is taken in the simulations, the amplitude of the Ca²⁺ increase and the spatial extent of the microdomains are much increased: opening of one IP₃R already leads to a Ca²⁺ increase in the microdomain that exceeds experimental observations. Because of the faster replenishment at the channel mouth, the flux through the IP₃R is indeed larger, which has a significant effect on junctional Ca²⁺ concentration. In other words, with the low value for D_ER, Ca²⁺ release via IP₃R is limited by local replenishment. SOCE activation does not compensate for release via the IP₃R because of the ~30 times larger current of the latter.

Simulation of the formation of TCR-dependent Ca²⁺ microdomains (TDCM)

In a similar way, Gil et al. (64) investigated how the interplay between Orai1 and RYR1 accounts for the formation of Ca²⁺ microdomains observed during the first 15 seconds after TCR/CD3 stimulation, termed TCR-dependent Ca²⁺ microdomains (TDCM). Following stimulation, NAADP signaling activates RYR1 and releases Ca²⁺ from the ER to the cytosol (40). Very recently, we identified hematological and neurological expressed 1-like protein (HN1L)/jupiter microtubule associated homolog 2 (JPT2) as NAADP receptor/binding protein that mediates activation of RYR1 upon NAADP binding (43); however, this step has not yet been included into the mathematical simulations. In parallel to increases of endogenous NAADP in the first seconds of TCR/CD3 stimulation, a higher co-localization of STIM1 with ORAI1 is observed in comparison to ADCM during which STIM2 co-localization predominates (27). This suggests that shortly after stimulation the main trigger of SOCE is Orai1 binding to STIM1 or most probably to heterotetramers of STIM1/2. Gil et al. (64) considered eight RYR1s around the junction, in the same way as described previously for IP₃Rs. Indeed, if considered inside the junction as in cardiac myocytes (73), one receptor would already saturate the signal, contrary to experimental observations. To consider the effect of STIM1/2 binding, the dependency of Orai1 opening to local luminal Ca²⁺ concentration at the mouth of RYR1 is modified, with a lower Ca²⁺ sensitivity to depletion than that of STIM2 but higher than that of STIM1. In order to reproduce the observed TDCM, the opening of seven RYR1s together with a cluster of five Orai1 channels bounded to STIM1/2 is needed (Figure 3B, blue circle). Moreover, simulations foresee that Orai1 bound to STIM2 homotetramers or to STIM1/2 heterotetramers do not activate further than 21% of their full conductance under early local depletion conditions of the sub-PM ER, either across IP₃R or RYR1. Thus, despite the involvement of different ER Ca²⁺ channels, the Ca²⁺ microdomains induced by adhesion of T cells before TCR/CD3 stimulation (ADCM) or by early NAADP signaling (TDCM) display similar characteristics in terms of amplitude and spatial extent because of the change in STIM isoform involved. Moreover, the model predicts that the flux at the mouth of the RYR1 or IP₃R, which determines the extent of local Ca²⁺ depletion, rapidly becomes similar despite the largely different conductance of these channels, due to the slow Ca²⁺ diffusion in the ER.

Towards a comprehensive model of Ca²⁺ signaling in the early phases of T cell stimulation

Based on these findings, T cell transition from quiescent to early activation can be described as schematized in Figure 3C. The resting state of a T cell (Figure 3C, left panel) is defined as an absence of any stimulatory contacts whatsoever, e.g. no contacts to cells, adhesive surfaces, or soluble activators; in this period, e.g. a T cell floating in gentle bloodstream, only basal (stochastic) openings of Orai1 channels are observed (Figure 3C, left panel). Upon weak (and likely transient) adhesion of the T cell, e.g. to proteins of the extracellular matrix during crossing of the basement membrane on the way to inflamed tissue areas, integrins are activated and cause stimulation of IP₃R. STIM2 recruits STIM1 close to the junction under such low stimulation conditions (ADCM) (Figure 3C, middle panel). Following TCR/CD3 stimulation (TDCM), NAADP signaling activates RYR1 and newly formed STIM1/2 heterotetramers need greater local luminal depletion to trigger a quantitatively similar Ca²⁺ influx through Orai1 Ca²⁺ channels (Figure 3C, right panel). To best describe T cell transition and Ca²⁺ microdomains evolution in the ER-PM junction, a computational model should include both the stochastic character of channels opening and closing, and the description of Ca²⁺ diffusion and buffering. This would require to further include a molecular description of STIM/Orai1 interactions considering the formation of molecular complexes of different sizes and compositions (85–87). Moreover, the opening of IP₃Rs or RYR1s in the cluster located near the junction should be described stochastically, as performed to describe Ca²⁺ puffs (88–90) or sparks (91). Mitochondria and PMCA’s should also be included, as in (57) and (60) The main challenges of such a model involve the formulation of the temporal shift from STIM1 to STIM2 pre-formed complexes with Orai1, and the interaction of ER Ca²⁺ channels IP₃R and RYR1 when active simultaneously during the transition time. Importantly, current T cell Ca²⁺ imaging techniques (27, 43, 48) are not able to resolve single junction contributions to the formation of Ca²⁺ microdomains. The observed Ca²⁺ increases at the PM upon adhesion or TCR/CD3 stimulation may in fact result from the input of several junctions depending on how close they are positioned to each other. Mathematical modeling may cast some light into the possible conformations of the different channels participating to the formation of Ca²⁺ microdomains in a cluster of ER-PM junctions. This would also require a dynamical description of Ca²⁺ buffering that affects inter-cluster communication. Finally, continuous modelling of the evolution of a Ca²⁺ increase in the form of a microdomain to global Ca²⁺ signaling will require a more complex multiscale whole-cell model. Such a study will build on the existing whole cell model taking into account specific spatial geometries for mitochondria and the ER (60). Simulated geometries will also include well-defined numbers and locations of the different types of Ca²⁺ channels and pumps, as they have been inferred from the small-scale modeling of the ER-PM junctions. Alternatively, it might be as valid to develop a compartmental model where the microdomain is modelled as a separate compartment in which the concentrations change must faster than in the rest of the cytoplasm (92, 93). The global model should not only be hybrid (i.e. stochastic and deterministic), but its numerical simulation should be adaptive in time since it must consider different time and spatial scales (82). It needs to be able to resolve both the initial single-channel Ca²⁺ gating and the later shift to a modeling of the global Ca²⁺ increase due to clusters of Ca²⁺ channels. Despite all these modelling challenges, such global simulations are expected to enhance our understanding of the localized regulation of Ca²⁺ dynamics, and hence of Ca²⁺ signaling in which small local effects can be either damped or amplified at the cellular level.

Concluding remarks

Computational modeling is increasingly used in cell biology to optimize the information that can be gained from the ever-improving techniques in imaging microscopy. In the field of ER-PM Ca²⁺ microdomains, this synergy has allowed to predict the number and spatial arrangement of Orai channels in the junction (18, 61), and to highlight the relatively minor role played by ER refilling by SERCA pumps in controlling the amplitude of the Ca²⁺ increase and the spatial spread of Ca²⁺ microdomains (62). When considering the ER Ca²⁺ channels surrounding the junction, a comparison between simulation results and high-resolution Ca²⁺ imaging uncovered the main role played by ER local depletion induced by weak IP₃ or NAADP signaling in triggering the formation of adhesion-dependent and TCR/CD3-dependent Ca²⁺ microdomains in T cells (63, 64). In the future, it will be necessary to resort to observations-based multiscale modeling to dissect the relationship between local and global Ca²⁺ signaling, which could emphasize the key role played by the PMCA’s and the mitochondria in the T cell Ca²⁺ dynamics beyond the first seconds after stimulation (60). Despite the corresponding computational challenge, this effort is necessary to reveal how the versatile Ca²⁺ signaling pathway fine-tunes T cell immune responses.

References

• 1

  BerridgeMJBootmanMDLippP. Calcium–a life and death signal. Nature (1998) 395:645–8. doi: 10.1038/27094

  • CrossRef
  • Google Scholar

• 2

  YangYXiePLiYBiYPruskyDB. Updating insights into the regulatory mechanisms of calcineurin-activated transcription factor crz1 in pathogenic fungi. J Fungi (Basel). (2022) 8:1082. doi: 10.3390/jof8101082

  • CrossRef
  • Google Scholar

• 3

  TsienRYRinkTJPoenieM. Measurement of cytosolic free Ca²⁺ in individual small cells using fluorescence microscopy with dual excitation wavelengths. Cell Calcium (1985) 6:145–57. doi: 10.1016/0143-4160(85)90041-7

  • CrossRef
  • Google Scholar

• 4

  WilliamsDAFogartyKETsienRYFayFS. Calcium gradients in single smooth muscle cells revealed by the digital imaging microscope using Fura-2. Nature (1985) 318:558–61. doi: 10.1038/318558a0

  • CrossRef
  • Google Scholar

• 5

  WoelkLMKannabiranSABrockVJGeeCELohrCGuseAHet al. Time-dependent image restoration of low-SNR live-cell ca² fluorescence microscopy data. Int J Mol Sci (2021) 22:11792. doi: 10.3390/ijms222111792

  • CrossRef
  • Google Scholar

• 6

  GuseAHGil MontoyaDCDiercksBP. Mechanisms and functions of calcium microdomains produced by ORAI channels, d-myo-inositol 1,4,5-trisphosphate receptors, or ryanodine receptors. Pharmacol Ther (2021) 223:107804. doi: 10.1016/j.pharmthera.2021.107804

  • CrossRef
  • Google Scholar

• 7

  ChengHLedererWJCannellMB. Calcium sparks: elementary events underlying excitation-contraction coupling in heart muscle. Science (1993) 262:740–4. doi: 10.1126/science.8235594

  • CrossRef
  • Google Scholar

• 8

  ParkerIChoiJYaoY. Elementary events of InsP₃-induced Ca²⁺ liberation in Xenopus oocytes: hot spots, puffs and blips. Cell Calcium (1996) 20:105–21. doi: 10.1016/S0143-4160(96)90100-1

  • CrossRef
  • Google Scholar

• 9

  BootmanMNiggliEBerridgeMLippP. Imaging the hierarchical Ca²⁺ signalling system in HeLa cells. J Physiol (1997) 499:307–14. doi: 10.1113/jphysiol.1997.sp021928

  • CrossRef
  • Google Scholar

• 10

  LippPNiggliE. Fundamental calcium release events revealed by two-photon excitation photolysis of caged calcium in Guinea-pig cardiac myocytes. J Physiol (1998) 508:801–9. doi: 10.1111/j.1469-7793.1998.801bp.x

  • CrossRef
  • Google Scholar

• 11

  ZouHLifshitzLMTuftRAFogartyKESingerJJ. Imaging Ca²⁺ entering the cytoplasm through a single opening of a plasma membrane cation channel. J Gen Physiol (1999) 114:575–88. doi: 10.1085/jgp.114.4.575

  • CrossRef
  • Google Scholar

• 12

  WangSQSongLSLakattaEGChengH. Ca²⁺ signalling between single L-type Ca²⁺ channels and ryanodine receptors in heart cells. Nature (2001) 410:592–6. doi: 10.1038/35069083

  • CrossRef
  • Google Scholar

• 13

  SonkusareSKBonevADLedouxJLiedtkeWKotlikoffMIHeppnerTJet al. Elementary Ca²⁺ signals through endothelial TRPV4 channels regulate vascular function. Science. (2012) 336:597–601. doi: 10.1126/science.1216283

  • CrossRef
  • Google Scholar

• 14

  MarchiSBittremieuxMMissiroliSMorgantiCPatergnaniSSbanoLet al. Endoplasmic reticulum-mitochondria communication through ca²⁺ Signaling: the importance of mitochondria-associated membranes (MAMs). Adv Exp Med Biol (2017) 997:49–67. doi: 10.1007/978-981-10-4567-7_4

  • CrossRef
  • Google Scholar

• 15

  LiouJKimMLHeoWDJonesJMyersJFerrellJJr.et al. STIM is a ca²⁺ Sensor essential for ca²⁺-store-depletion-triggered ca²⁺ Influx’. Curr Biol (2005) 15:1235–41. doi: 10.1016/j.cub.2005.05.055

  • CrossRef
  • Google Scholar

• 16

  FeskeSGwackYPrakriyaMSrikanthSPuppelS-HTanasaBet al. A mutation in orai1 causes immune deficiency by abrogating CRAC channel function. Nature (2006) 441:179–285. doi: 10.1038/nature04702

  • CrossRef
  • Google Scholar

• 17

  PrakriyaMFeskeSGwackYSrikanthSRaoAHoganP. Orai1 is an essential pore subunit of the CRAC channel. Nature (2006) 443:230–3. doi: 10.1038/nature05122

  • CrossRef
  • Google Scholar

• 18

  HoganPG. The STIM1-ORAI1 microdomain. Cell Calcium (2015) 58:357–67. doi: 10.1016/j.ceca.2015.07.001

  • CrossRef
  • Google Scholar

• 19

  BarakPParekhAB. Signaling through ca²⁺ Microdomains from store-operated CRAC channels. Cold Spring Harb. Perspect Biol (2020) 12:a035097. doi: 10.1101/cshperspect.a035097

  • CrossRef
  • Google Scholar

• 20

  DupontGFalckeMKirkVSneydJ. Models of calcium signalling. (Switzerland: Springer International Publishing) (2016).

  • Google Scholar

• 21

  SwillensSDupontGCombettesLChampeilP. From calcium blips to calcium puffs: theoretical analysis of the requirements for interchannel communication. Proc Natl Acad Sci U S A. (1999) 96:13750–5. doi: 10.1073/pnas.96.24.13750

  • CrossRef
  • Google Scholar

• 22

  GroffJRSmithGD. Ryanodine receptor allosteric coupling and the dynamics of calcium sparks. Biophys J (2008) 95:135–54. doi: 10.1529/biophysj.107.119982

  • CrossRef
  • Google Scholar

• 23

  ThulRThurleyKFalckeM. Toward a predictive model of Ca²⁺ puffs. Chaos (2009) 19:037108. doi: 10.1063/1.3183809

  • CrossRef
  • Google Scholar

• 24

  QiHHuangYRüdigerSShuaiJ. Frequency and relative prevalence of calcium blips and puffs in a model of small IP3R clusters. Biophys J (2014) 106:2353–63. doi: 10.1016/j.bpj.2014.04.027

  • CrossRef
  • Google Scholar

• 25

  QiHLiLShuaiJ. Optimal microdomain crosstalk between endoplasmic reticulum and mitochondria for Ca²⁺ oscillations. Sci Rep (2015) 5:7984. doi: 10.1038/srep07984

  • CrossRef
  • Google Scholar

• 26

  WolfIMDiercksBPGattkowskiECzarniakFKempskiJWernerRet al. Frontrunners of T cell activation: Initial, localized Ca²⁺ signals mediated by NAADP and the type 1 ryanodine receptor. Sci Signal (2015) 8:ra102. doi: 10.1126/scisignal.aab0863

  • CrossRef
  • Google Scholar

• 27

  WeißMHernandezLCGil MontoyaDCLöhndorfAKrügerAKopdagMet al. Adhesion to laminin-1 and collagen IV induces the formation of Ca²⁺ microdomains that sensitize mouse T cells for activation. Sci Signal (2023) 16(790):eabn9405. doi: 10.1126/scisignal.abn9405

  • CrossRef
  • Google Scholar

• 28

  BrockVJWolfIMAEr-LukowiakMLoryNStählerTWoelkLMet al. P2X4 and P2X7 are essential players in basal T cell activity and Ca²⁺ signaling milliseconds after T cell activation. Sci Adv (2022) 8:eabl9770. doi: 10.1126/sciadv.abl9770

  • CrossRef
  • Google Scholar

• 29

  WeismannMGuseAHSorokinLBrökerBFrieserMHallmannRet al. Integrin-mediated intracellular ca²⁺ Signaling in jurkat T lymphocytes. J Immunol (1997) 158:1618–27. doi: 10.4049/jimmunol.158.4.1618

  • CrossRef
  • Google Scholar

• 30

  GuanJLRole of focal adhesion kinase in integrin signaling. Intern J Biochem Cell Biol (1997) 29:1085–96. doi: 10.1016/s1357-2725(97)00051-4

  • CrossRef
  • Google Scholar

• 31

  KimT-JLeiLSeongJSuhJ-SJangY-KJungSHet al. Matrix rigidity-dependent regulation of ca²⁺ at plasma membrane microdomains by FAK visualized by fluorescence resonance energy transfer. Adv Sci (2019) 6:1801290. doi: 10.1002/advs.201801290

  • CrossRef
  • Google Scholar

• 32

  ZhangXChattopadhyayAJiQSOwenJDRuestPJCarpenterGet al. Focal adhesion kinase promotes phospholipase C-gamma1 activity. Proc Natl Acad Sci U S A. (1999) 96:9021–6. doi: 10.1073/pnas.96.16.9021

  • CrossRef
  • Google Scholar

• 33

  BosanacIAlattiaJ-RMalTChanJTalaricoSTongFet al. Structure of the inositol 1,4,5-trisphosphate receptor binding core in complex with its ligand’. Nature (2002) 420:696–700. doi: 10.1038/nature01268

  • CrossRef
  • Google Scholar

• 34

  KunerthSLanghorstMFSchwarzmannNGuXHuangLYangZet al. Amplification and propagation of pacemaker Ca²⁺ signals by cyclic ADP-ribose and the type 3 ryanodine receptor in T cells. J Cell Sci (2004) 117:2141–9. doi: 10.1242/jcs.01063

  • CrossRef
  • Google Scholar

• 35

  SchöttelndreierHPotterBVMayrGWGuseAH. Mechanisms involved in alpha6beta1-integrin-mediated ca²⁺ Signalling. Cell Signal (2001) 13:895–99. doi: 10.1016/S0898-6568(01)00225-X

  • CrossRef
  • Google Scholar

• 36

  RandriamampitaCBoullaGRevyPLemaitreFTrautmannA. T cell adhesion lowers the threshold for antigen detection. Eur J Immunol (2003) 33:1215–23. doi: 10.1002/eji.200323844

  • CrossRef
  • Google Scholar

• 37

  PribilaJJonathanTQualeAMuellerKShimizuY. Integrins and T cell–mediated immunity, annu. Rev Immunol (2004) 22:157–80. doi: 10.1146/annurev.immunol.22.012703.104649

  • CrossRef
  • Google Scholar

• 38

  BergIPotterBVMayrGWGuseAH. Nicotinic acid adenine dinucleotide phosphate (NAADP⁺) is an essential regulator of T-lymphocyte Ca²⁺ signaling. J Cell Biol (2000) 150:581–8. doi: 10.1083/jcb.150.3.581

  • CrossRef
  • Google Scholar

• 39

  LanghorstMFSchwarzmannNGuseAH. Ca²⁺ release via ryanodine receptors and Ca²⁺ entry: major mechanisms in NAADP-mediated Ca²⁺ signaling in T-lymphocytes. Cell Signal (2004) 16:1283–9. doi: 10.1016/j.cellsig.2004.03.013

  • CrossRef
  • Google Scholar

• 40

  DammermannWZhangBNebelMCordiglieriCOdoardiFKirchbergerTet al. NAADP-mediated Ca²⁺ signaling via type 1 ryanodine receptor in T cells revealed by a synthetic NAADP antagonist. Proc Natl Acad Sci U S A. (2009) 106:10678–83. doi: 10.1073/pnas.0809997106

  • CrossRef
  • Google Scholar

• 41

  DammermannWGuseAH. Functional ryanodine receptor expression is required for NAADP-mediated local Ca²⁺ signaling in T-lymphocytes. J Biol Chem (2005) 280:21394–9. doi: 10.1074/jbc.M413085200

  • CrossRef
  • Google Scholar

• 42

  WalsethTFLin-MoshierYWeberKMarchantJSSlamaJTGuseAH. Nicotinic acid adenine dinucleotide 2'-phosphate (NAADP) binding proteins in T-lymphocytes. Messenger (Los Angel). (2012) 1:86–94. doi: 10.1166/msr.2012.1008

  • CrossRef
  • Google Scholar

• 43

  RoggenkampHGKhansahibIHernandez CLCZhangYLodyginDKrügerAet al. HN1L/JPT2: A signaling protein that connects NAADP generation to Ca²⁺ microdomain formation. Sci Signal (2021) 14:eabd5647. doi: 10.1126/scisignal.abd5647

  • CrossRef
  • Google Scholar

• 44

  GunaratneGSBrailoiuEHeSUnterwaldEMPatelSSlamaJTet al. Essential requirement for JPT2 in NAADP-evoked Ca²⁺ signaling. Sci Signal (2021) 14:eabd5605. doi: 10.1126/scisignal.abd5605

  • CrossRef
  • Google Scholar

• 45

  ZhangJGuanXShahKYanJ. Lsm12 is an NAADP receptor and a two-pore channel regulatory protein required for calcium mobilization from acidic organelles. Nat Commun (2021) 12:4739. doi: 10.1038/s41467-021-24735-z

  • CrossRef
  • Google Scholar

• 46

  GuseAH. Linking NAADP to ion channel activity: a unifying hypothesis. Sci Signal (2012) 5:pe18. doi: 10.1126/scisignal.2002890

  • CrossRef
  • Google Scholar

• 47

  GasserABruhnSGuseAH. Second messenger function of nicotinic acid adenine dinucleotide phosphate revealed by an improved enzymatic cycling assay. J Biol Chem (2006) 281:16906–13. doi: 10.1074/jbc.M601347200

  • CrossRef
  • Google Scholar

• 48

  GuFKrügerARoggenkampHGAlpersRLodyginDJaquetVet al. Dual NADPH oxidases DUOX1 and DUOX2 synthesize NAADP and are necessary for Ca²⁺ signaling during T cell activation. Sci Signal (2021) 14:eabe3800. doi: 10.1126/scisignal.abe3800

  • CrossRef
  • Google Scholar

• 49

  SchmidFBruhnSWeberKMittrückerHWGuseAH. CD38: a NAADP degrading enzyme. FEBS Lett (2011) 585:3544–8. doi: 10.1016/j.febslet.2011.10.017

  • CrossRef
  • Google Scholar

• 50

  SchmidFFliegertRWestphalTBaucheAGuseAH. Nicotinic acid adenine dinucleotide phosphate (NAADP) degradation by alkaline phosphatase. J Biol Chem (2012) 287:32525–34. doi: 10.1074/jbc.M112.362715

  • CrossRef
  • Google Scholar

• 51

  WilliamsBLSchreiberKLZhangWWangeRLSamelsonLELeibsonPJet al. Genetic evidence for differential coupling of Syk family kinases to the T-cell receptor: reconstitution studies in a ZAP-70-deficient Jurkat T-cell line. Mol Cell Biol (1998) 18:1388–99. doi: 10.1128/MCB.18.3.1388

  • CrossRef
  • Google Scholar

• 52

  IrvinBJWilliamsBLNilsonAEMaynorHOAbrahamRT. Pleiotropic contributions of phospholipase C-gamma1 (PLC-gamma1) to T-cell antigen receptor-mediated signaling: reconstitution studies of a PLC-gamma1-deficient Jurkat T-cell line. Mol Cell Biol (2000) 20:9149–691611. doi: 10.1128/MCB.20.24.9149-9161.2000

  • CrossRef
  • Google Scholar

• 53

  GuseAHRothEEmmrichF. Intracellular ca²⁺ pools in jurkat T-lymphocytes. Biochem J (1993) 291:447–51. doi: 10.1042/bj2910447

  • CrossRef
  • Google Scholar

• 54

  QuannEJMerinoEFurutaTHuseM. Localized diacylglycerol drives the polarization of the microtubule-organizing center in T cells. Nat Immunol (2009) 10:627–35. doi: 10.1038/ni.1734

  • CrossRef
  • Google Scholar

• 55

  GuseAHda SilvaCPBergISkapenkoALWeberKHeyerPet al. Regulation of calcium signalling in T lymphocytes by the second messenger cyclic ADP-ribose. Nature (1999) 398:70–3. doi: 10.1038/18024

  • CrossRef
  • Google Scholar

• 56

  TrebakMKinetJP. Calcium signalling in T cells. Nat Rev Immunol (2019) 19:154–69. doi: 10.1038/s41577-018-0110-7

  • CrossRef
  • Google Scholar

• 57

  KummerowCJunkerCKruseKRiegerHQuintanaAHothM. The immunological synapse controls local and global calcium signals in T lymphocytes. Immunol Rev (2009) 231:132–47. doi: 10.1111/j.1600-065X.2009.00811.x

  • CrossRef
  • Google Scholar

• 58

  QuintanaAPascheMJunkerCAl-AnsaryDRiegerHKummerowCet al. Calcium microdomains at the immunological synapse: how ORAI channels, mitochondria and calcium pumps generate local calcium signals for efficient T-cell activation. EMBO J (2011) 30:3895–912. doi: 10.1038/emboj.2011.289

  • CrossRef
  • Google Scholar

• 59

  PeglowMNiemeyerBHothMRiegerH. Interplay of channels, pumps and organelle location in calcium microdomain formation. New J Phys (2013) 15:055022. doi: 10.1088/1367-2630/15/5/055022

  • CrossRef
  • Google Scholar

• 60

  MaccariIZhaoRPeglowMSchwarzKHornakIPascheMet al. Cytoskeleton rotation relocates mitochondria to the immunological synapse and increases calcium signals. Cell Calcium. (2016) 60:309–21. doi: 10.1016/j.ceca.2016.06.007

  • CrossRef
  • Google Scholar

• 61

  SamantaKKarPMiramsGRParekhAB. Ca²⁺ Channel re-localization to plasma-membrane microdomains strengthens activation of ca²⁺-dependent nuclear gene expression. Cell Rep (2015) 12:203–16. doi: 10.1016/j.celrep.2015.06.018

  • CrossRef
  • Google Scholar

• 62

  McIvorECoombesSThulR. Three-dimensional spatio-temporal modelling of store operated Ca²⁺ entry: Insights into ER refilling and the spatial signature of Ca²⁺ signals. Cell Calcium. (2018) 73:11–24. doi: 10.1016/j.ceca.2018.03.006

  • CrossRef
  • Google Scholar

• 63

  GilDGuseAHDupontG. Three-dimensional model of sub-plasmalemmal ca²⁺ Microdomains evoked by the interplay between ORAI1 and insP₃ receptors. Front Immunol (2021) 12:659790. doi: 10.3389/fimmu.2021.659790

  • CrossRef
  • Google Scholar

• 64

  GilDDiercksBPGuseAHDupontG. Three-dimensional model of sub-plasmalemmal ca²⁺ Microdomains evoked by T cell receptor/CD3 complex stimulation. Front Mol Biosci (2022) 9:811145. doi: 10.3389/fmolb.2022.811145

  • CrossRef
  • Google Scholar

• 65

  DarganSSchwallerBParkerI. Spatiotemporal patterning of IP₃-mediated Ca²⁺ signals in Xenopus oocytes by Ca²⁺-binding proteins. J Physiol (2004) 556:447–61. doi: 10.1113/jphysiol.2003.059204

  • CrossRef
  • Google Scholar

• 66

  PetersenOCourjaretRMachacaK. Ca²⁺ tunnelling through the ER lumen as a mechanism for delivering Ca²⁺ entering via store-operated Ca²⁺ channels to specific target sites. J Physiol (2017) 595:2999–3014. doi: 10.1113/JP272772

  • CrossRef
  • Google Scholar

• 67

  MarchantJGunaratneGCaiXSlamaJPatelS. NAADP-binding proteins and their identity. Trends Biochem Sci (2022) 47:235–49. doi: 10.1016/j.tibs.2021.10.008

  • CrossRef
  • Google Scholar

• 68

  WalsethTGuseAH. NAADP: from discovery to mechanism. Front Immunol (2021) 12:703326. doi: 10.3389/fimmu.2021.703326

  • CrossRef
  • Google Scholar

• 69

  WuMBuchananJLuikRLewisR. Ca²⁺ store depletion causes STIM1 to accumulate in ER regions closely associated with the plasma membrane. J Cell Biol (2006) 174:803–13. doi: 10.1083/jcb.200604014

  • CrossRef
  • Google Scholar

• 70

  VarnaiPTothBTothDHunyadyLBallaT. Visualization and manipulation of plasma membrane-endoplasmic reticulum contact sites indicates the presence of additional molecular components within the STIM1-Orai1 complex. J Biol Chem (2007) 282:29678–90. doi: 10.1074/jbc.M704339200

  • CrossRef
  • Google Scholar

• 71

  ShenYThillaiappanNTaylorC. The store-operated Ca²⁺ entry complex comprises a small cluster of STIM1 associated with one Orai1 channel. Proc Natl Acad Sci USA (2021) 118:e2010789118. doi: 10.1073/pnas.2010789118

  • CrossRef
  • Google Scholar

• 72

  ThillaiappanNBChavdaAPToveySCProleDLTaylorCW. Ca²⁺ signals initiate at immobile IP₃ receptors adjacent to ER-plasma membrane junctions. Nat Commun (2017) 8:1505. doi: 10.1038/s41467-017-01644-8

  • CrossRef
  • Google Scholar

• 73

  JonesPPMacQuaideNLouchWE. Dyadic plasticity in cardiomyocytes. Front Physiol (2018) 9:1773. doi: 10.3389/fphys.2018.01773

  • CrossRef
  • Google Scholar

• 74

  GillespieD. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J Comput Phys (1976) 22:403–34. doi: 10.1016/0021-9991(76)90041-3

  • CrossRef
  • Google Scholar

• 75

  HakeJLinesG. Stochastic binding of Ca²⁺ ions in the dyadic cleft: Continuous versus random walk description of diffusion. Biophys J (2008) 94:4184–201. doi: 10.1529/biophysj.106.103523

  • CrossRef
  • Google Scholar

• 76

  WinslowRLTanskanenAChenMGreensteinJL. Multiscale modeling of calcium signaling in the cardiac dyad. Ann NY. Acad Sci (2006) 1080:362–75. doi: 10.1196/annals.1380.027

  • CrossRef
  • Google Scholar

• 77

  WiederNFinkRvon WegnerF. Exact stochastic simulation of a calcium microdomain reveals the impact of Ca²⁺ fluctuations on IP3R gating. Biophys J (2015) 108:557–67. doi: 10.1016/j.bpj.2014.11.3458

  • CrossRef
  • Google Scholar

• 78

  WiederNFinkRHAvon WegnerF. Simulation strategies for calcium microdomains and calcium noise. Adv Exp Med Biol (2020) 1131:771–97. doi: 10.1007/978-3-030-12457-1_31

  • CrossRef
  • Google Scholar

• 79

  WeinbergSHSmithGD. The influence of Ca²⁺ buffers on free [Ca²⁺] fluctuations and the effective volume of Ca²⁺ microdomains. Biophys J (2014) 106:2693–709. doi: 10.1016/j.bpj.2014.04.045

  • CrossRef
  • Google Scholar

• 80

  ThulRFalckeM. Release currents of IP₃ receptor channel clusters and concentration profiles. Biophys J (2004) 86:2660–73. doi: 10.1016/S0006-3495(04)74322-2

  • CrossRef
  • Google Scholar

• 81

  MeansSSmithAJShepherdJShadidJFowlerJWojcikiewiczRJet al. Reaction diffusion modeling of calcium dynamics with realistic ER geometry. Biophys J (2006) 91:537–57. doi: 10.1529/biophysj.105.075036

  • CrossRef
  • Google Scholar

• 82

  VierhellerJNeubertWFalckeMGilbertSChamakuriN. A multiscale computational model of spatially resolved calcium cycling in cardiac myocytes: from detailed cleft dynamics to the whole cell concentration profiles. Front Physiol (2015) 6:255. doi: 10.3389/fphys.2015.00255

  • CrossRef
  • Google Scholar

• 83

  COMSOL multiphysics® v. 5.4. COMSOL AB, Stockholm, Sweden. Available at: www.comsol.com.

  • Google Scholar

• 84

  SubediKPOngHLSonGYLiuXAmbudkarIS. STIM2 induces activated conformation of STIM1 to control orai1 function in ER-PM junctions. Cell Rep (2018) 23:522–34. doi: 10.1016/j.celrep.2018.03.065

  • CrossRef
  • Google Scholar

• 85

  KilchTAlansaryDPeglowMDörrKRychkovGRiegerHet al. Mutations of the Ca^2+-sensing stromal interaction molecule STIM1 regulate Ca²⁺ influx by altered oligomerization of STIM1 and by destabilization of the Ca²⁺ channel Orai1. J Biol Chem (2013) 288:1653–64. doi: 10.1074/jbc.M112.417246

  • CrossRef
  • Google Scholar

• 86

  MelunisJHershbergU. A spatially heterogeneous Gillespie algorithm modeling framework that enables individual molecule history and tracking. Engineer. Appl Artif Intell (2017) 62:304–11. doi: 10.1016/j.engappai.2016.09.010

  • CrossRef
  • Google Scholar

• 87

  SchmidtBAlansaryDBogeskiINiemeyerBARiegerH. Reaction-diffusion model for STIM-ORAI interaction: The role of ROS and mutations. J Theor Biol (2019) 470:64–75. doi: 10.1016/j.jtbi.2019.02.010

  • CrossRef
  • Google Scholar

• 88

  FalckeMMoeinMTilūnaitėAThulRSkupinA. On the phase space structure of IP₃ induced Ca²⁺ signalling and concepts for predictive modeling. Chaos (2018) 28:045115. doi: 10.1063/1.5021073

  • CrossRef
  • Google Scholar

• 89

  VoorsluijsVDawsonSPDe DeckerYDupontG. Deterministic limit of intracellular calcium spikes. Phys Rev Lett (2019) 122:088101. doi: 10.1103/PhysRevLett.122.088101

  • CrossRef
  • Google Scholar

• 90

  PowellJFalckeMSkupinABellamyTCKypraiosTThulR. A statistical view on calcium oscillations. Adv Exp Med Biol (2020) 1131:799–826. doi: 10.1007/978-3-030-12457-1_32

  • CrossRef
  • Google Scholar

• 91

  IzuLTMeansSAShadidJNChen-IzuYBalkeCW. Interplay of ryanodine receptor distribution and calcium dynamics. Biophys J (2006) 91:95–112. doi: 10.1529/biophysj.105.077214

  • CrossRef
  • Google Scholar

• 92

  GreensteinJLHinchRWinslowRL. Mechanisms of excitation-contraction coupling in an integrative model of the cardiac ventricular myocyte. Biophys J (2006) 90:77–91. doi: 10.1529/biophysj.105.065169

  • CrossRef
  • Google Scholar

• 93

  PennyCJKilpatrickBSHanJMSneydJPatelS. A computational model of lysosome-ER Ca²⁺ microdomains. J Cell Sci (2014) 127:2934–43. doi: 10.1242/jcs.149047

  • CrossRef
  • Google Scholar