Inositol (1,4,5)-trisphosphate receptor microarchitecture shapes Ca2+ puff kinetics

Three-Dimensional Model of Sub-Plasmalemmal Ca2+ Microdomains Evoked by the Interplay Between ORAI1 and InsP3 Receptors

Ca²⁺ signaling plays an essential role in T cell activation, which is a key step to start an adaptive immune response. During the transition from a quiescent to a fully activated state, Ca²⁺ microdomains characterized by reduced spatial and temporal extents are observed in the junctions between the plasma membrane (PM) and the endoplasmic reticulum (ER). Such Ca²⁺ responses can also occur in response to T cell adhesion to other cells or extracellular matrix proteins in otherwise unstimulated T cells. These non-TCR/CD3-dependent Ca²⁺ microdomains rely on d-myo-inositol 1,4,5-trisphosphate (IP₃) signaling and subsequent store operated Ca²⁺ entry (SOCE) via the ORAI/STIM system. The detailed molecular mechanism of adhesion-dependent Ca²⁺ microdomain formation remains to be fully elucidated. We used mathematical modeling to investigate the spatiotemporal characteristics of T cell Ca²⁺ microdomains and their molecular regulators. We developed a reaction-diffusion model using COMSOL Multiphysics to describe the evolution of cytosolic and ER Ca²⁺ concentrations in a three-dimensional ER-PM junction. Equations are based on a previously proposed realistic description of the junction, which is extended to take into account IP₃ receptors (IP₃R) that are located next to the junction. The first model only considered the ORAI channels and the SERCA pumps. Taking into account the existence of preformed clusters of ORAI1 and STIM2, ORAI1 slightly opens in conditions of a full ER. These simulated Ca²⁺ microdomains are too small as compared to those observed in unstimulated T cells. When considering the opening of the IP₃Rs located near the junction, the local depletion of ER Ca²⁺ allows for larger Ca²⁺ fluxes through the ORAI1 channels and hence larger local Ca²⁺ concentrations. Computational results moreover show that Ca²⁺ diffusion in the ER has a major impact on the Ca²⁺ changes in the junction, by affecting the local Ca²⁺ gradients in the sub-PM ER. Besides pointing out the likely involvement of the spontaneous openings of IP₃Rs in the activation of SOCE in conditions of T cell adhesion prior to full activation, the model provides a tool to investigate how Ca²⁺ microdomains extent and interact in response to T cell receptor activation.

Using two dyes to observe the competition of Ca2+ trapping mechanisms and their effect on intracellular Ca2+ signals

The specificity and universality of intracellular [Formula: see text] signals rely on the variety of spatio-temporal patterns that the [Formula: see text] concentration can display. [Formula: see text] liberation through inositol 1,4,5-trisphosphate receptors ([Formula: see text]) is key for this variety. In this paper, we study how the competition between buffers of different kinetics affects [Formula: see text] signals that involve [Formula: see text] release through [Formula: see text]. The study also provides insight into the underlying spatial distribution of the channels that participate in the signals. Previous works on the effects of [Formula: see text] buffers have drawn conclusions 'indirectly' by observing the [Formula: see text]-bound dye distributions in the presence of varying concentrations of exogenous buffers and using simulations to interpret the results. In this paper, we make visible the invisible by observing the signals simultaneously with two dyes, [Formula: see text] and [Formula: see text], each of which plays the role of a slow or fast [Formula: see text] buffer, respectively. Our observations obtained for different concentrations of [Formula: see text] highlight the dual role that fast buffers exert on the dynamics, either reducing the intracluster channel coupling or preventing channel inhibition and allowing the occurrence of relatively long cycles of [Formula: see text] release. Our experiments also show that signals with relatively high [Formula: see text] release rates remain localized in the presence of large [Formula: see text] concentrations, while the mean speed of the elicited waves increases. We interpret this as a consequence of the more effective uncoupling between [Formula: see text] clusters as the slow dye concentration increases. Combining the analysis of the experiments with numerical simulations, we also conclude that [Formula: see text] release not only occurs within the close vicinity of the centers of the clearly identifiable release sites ([Formula: see text] clusters) but there are also functional [Formula: see text] in between them.

IP3 receptor signaling and endothelial barrier function

The endothelium, a monolayer of endothelial cells lining vessel walls, maintains tissue-fluid homeostasis by restricting the passage of the plasma proteins and blood cells into the interstitium. The ion Ca²⁺, a ubiquitous secondary messenger, initiates signal transduction events in endothelial cells that is critical to control of vascular tone and endothelial permeability. The ion Ca²⁺ is stored inside the intracellular organelles and released into the cytosol in response to environmental cues. The inositol 1,4,5-trisphosphate (IP₃) messenger facilitates Ca²⁺ release through IP₃ receptors which are Ca²⁺-selective intracellular channels located within the membrane of the endoplasmic reticulum. Binding of IP₃ to the IP₃Rs initiates assembly of IP₃R clusters, a key event responsible for amplification of Ca²⁺ signals in endothelial cells. This review discusses emerging concepts related to architecture and dynamics of IP₃R clusters, and their specific role in propagation of Ca²⁺ signals in endothelial cells.

Calcium waves in a grid of clustered channels with synchronous IP3 binding and unbinding

Calcium signals in cells occur at multiple spatial scales and variable temporal duration. However, a physical explanation for transitions between long-lasting global oscillations and localized short-term elevations (puffs) of cytoplasmic Ca²⁺ is still lacking. Here we introduce a phenomenological, coarse-grained model for the calcium variable, which is represented by ordinary differential equations. Due to its small number of parameters, and its simplicity, this model allows us to numerically study the interplay of multi-scale calcium concentrations with stochastic ion channel gating dynamics even in larger systems. We apply this model to a single cluster of inositol trisphosphate (IP ₃) receptor channels and find further evidence for the results presented in earlier work: a single cluster may be capable of producing different calcium release types, where long-lasting events are accompanied by unbinding of IP ₃ from the receptor (Rückl et al., PLoS Comput. Biol. 11, e1003965 (2015)). Finally, we show the practicability of the model in a grid of 64 clusters which is computationally intractable with previous high-resolution models. Here long-lasting events can lead to synchronized oscillations and waves, while short events stay localized. The frequency of calcium releases as well as their coherence can thereby be regulated by the amplitude of IP ₃ stimulation. Finally the model allows for a new explanation of oscillating [IP ₃], which is not based on metabolic production and degradation of IP ₃.

Luminal Ca(2+) dynamics during IP3R mediated signals

The role of cytosolic Ca(2+) on the kinetics of Inositol 1,4,5-triphosphate receptors (IP3Rs) and on the dynamics of IP3R-mediated Ca(2+) signals has been studied at large both experimentally and by modeling. The role of luminal Ca(2+) has not been investigated with that much detail although it has been found that it is relevant for signal termination in the case of Ca(2+) release through ryanodine receptors. In this work we present the results of observing the dynamics of luminal and cytosolic Ca(2+) simultaneously in Xenopus laevis oocytes. Combining observations and modeling we conclude that there is a rapid mechanism that guarantees the availability of free Ca(2+) in the lumen even when a relatively large Ca(2+) release is evoked. Comparing the dynamics of cytosolic and luminal Ca(2+) during a release, we estimate that they are consistent with a 80% of luminal Ca(2+) being buffered. The rapid availability of free luminal Ca(2+) correlates with the observation that the lumen occupies a considerable volume in several regions across the images.

On the Adjacency Matrix of RyR2 Cluster Structures

In the heart, electrical stimulation of cardiac myocytes increases the open probability of sarcolemmal voltage-sensitive Ca2+ channels and flux of Ca2+ into the cells. This increases Ca2+ binding to ligand-gated channels known as ryanodine receptors (RyR2). Their openings cause cell-wide release of Ca2+, which in turn causes muscle contraction and the generation of the mechanical force required to pump blood. In resting myocytes, RyR2s can also open spontaneously giving rise to spatially-confined Ca2+ release events known as "sparks." RyR2s are organized in a lattice to form clusters in the junctional sarcoplasmic reticulum membrane. Our recent work has shown that the spatial arrangement of RyR2s within clusters strongly influences the frequency of Ca2+ sparks. We showed that the probability of a Ca2+ spark occurring when a single RyR2 in the cluster opens spontaneously can be predicted from the precise spatial arrangements of the RyR2s. Thus, "function" follows from "structure." This probability is related to the maximum eigenvalue (λ1) of the adjacency matrix of the RyR2 cluster lattice. In this work, we develop a theoretical framework for understanding this relationship. We present a stochastic contact network model of the Ca2+ spark initiation process. We show that λ1 determines a stability threshold for the formation of Ca2+ sparks in terms of the RyR2 gating transition rates. We recapitulate these results by applying the model to realistic RyR2 cluster structures informed by super-resolution stimulated emission depletion (STED) microscopy. Eigendecomposition of the linearized mean-field contact network model reveals functional subdomains within RyR2 clusters with distinct sensitivities to Ca2+. This work provides novel perspectives on the cardiac Ca2+ release process and a general method for inferring the functional properties of transmembrane receptor clusters from their structure.

Abortive and propagating intracellular calcium waves: analysis from a hybrid model

The functional properties of inositol(1,4,5)-triphosphate (IP3) receptors allow a variety of intracellular Ca(2+) phenomena. In this way, global phenomena, such as propagating and abortive Ca(2+) waves, as well as local events such as puffs, have been observed. Several experimental studies suggest that many features of global phenomena (e.g., frequency, amplitude, speed wave) depend on the interplay of biophysical processes such as diffusion, buffering, efflux and influx rates, which in turn depend on parameters such as buffer concentration, Ca(2+) pump density, cytosolic IP3 level, and intercluster distance. Besides, it is known that cells are able to modify some of these parameters in order to regulate the Ca(2+) signaling. By using a hybrid model, we analyzed different features of the hierarchy of calcium events as a function of two relevant parameters for the calcium signaling, the intercluster distance and the pump strength or intensity. In the space spanned by these two parameters, we found two modes of calcium dynamics, one dominated by abortive calcium waves and the other by propagating waves. Smaller distances between the release sites promote propagating calcium waves, while the increase of the efflux rate makes the transition from propagating to abortive waves occur at lower values of intercluster distance. We determined the frontier between these two modes, in the parameter space defined by the intercluster distance and the pump strength. Furthermore, we found that the velocity of simulated calcium waves accomplishes Luther's law, and that an effective rate constant for autocatalytic calcium production decays linearly with both the intercluster distance and the pump strength.

Ca²⁺ images obtained in different experimental conditions shed light on the spatial distribution of IP₃ receptors that underlie Ca²⁺ puffs

Many intracellular Ca(2+) signals involve Ca(2+) release from the endoplasmic reticulum through inositol 1,4,5-trisphosphate receptors (IP3Rs). The open probability of IP3Rs depends on cytosolic Ca(2+) so that these signals involve Ca(2+)-induced Ca(2+)-release (CICR). IP3Rs are organized in clusters. The signals they mediate are observed using single-wavelength dyes and, often, a slow Ca(2+) buffer (EGTA) is added to disrupt CICR between clusters and keep the signals spatially restricted. It is assumed that the presence of the dye or of EGTA does not alter the intra-cluster Ca(2+) dynamics. In this paper we analyze this issue combining experiments and numerical simulations. We compare the properties of local signals known as puffs observed with different dyes and EGTA concentrations. We determine that although the dye or EGTA does not alter the intra-cluster dynamics, the set of observable events is different depending on the degree of inter-cluster uncoupling of the experiment. An analysis of the observations shows that the events that are missed for insufficient inter-cluster uncoupling are those of fastest amplitude growth rate. This agrees with a spatial organization in which the largest amplitude events correspond to clusters with densely packed active IP3Rs.

Fluorescence fluctuations and equivalence classes of Ca²⁺ imaging experiments

release into the cytosol through inositol 1,4,5-trisphosphate receptors (IP₃Rs) plays a relevant role in numerous physiological processes. IP₃R-mediated signals involve -induced -release (CICR) whereby release through one open IP₃R induces the opening of other channels. IP₃Rs are apparently organized in clusters. The signals can remain localized (i.e., puffs) if CICR is limited to one cluster or become waves that propagate between clusters. puffs are the building blocks of waves. Thus, there is great interest in determining puff properties, especially in view of the current controversy on the spatial distribution of activatable IP₃Rs. puffs have been observed in intact cells with optical techniques proving that they are intrinsically stochastic. Obtaining a correct picture of their dynamics then entails being able to detect the whole range of puff sizes. puffs are observed using visible single-wavelength dyes, slow exogenous buffers (e.g., EGTA) to disrupt inter-cluster CICR and UV-photolyzable caged IP₃. Single-wavelength dyes increase their fluorescence upon calcium binding producing images that are strongly dependent on their kinetic, transport and photophysical properties. Determining the artifacts that the imaging setting introduces is particularly relevant when trying to analyze the smallest signals. In this paper we introduce a method to estimate the expected signal-to-noise ratio of imaging experiments that use single-wavelength dyes. The method is based on the Number and Brightness technique. It involves the performance of a series of experiments and their subsequent analysis in terms of a fluorescence fluctuation model with which the model parameters are quantified. Using the model, the expected signal-to-noise ratio is then computed. Equivalence classes between different experimental conditions that produce images with similar signal-to-noise ratios can then be established. The method may also be used to estimate the smallest signals that can reliably be observed with each setting.

Buffer regulation of calcium puff sequences

Puffs are localized Ca(2 +) signals that arise in oocytes in response to inositol 1,4,5-trisphosphate (IP3). They are the result of the liberation of Ca(2 +) from the endoplasmic reticulum through the coordinated opening of IP3 receptor/channels clustered at a functional release site. The presence of buffers that trap Ca(2 +) provides a mechanism that enriches the spatio-temporal dynamics of cytosolic calcium. The expression of different types of buffers along the cell's life provides a tool with which Ca(2 +) signals and their responses can be modulated. In this paper we extend the stochastic model of a cluster of IP3R-Ca(2 +) channels introduced previously to elucidate the effect of buffers on sequences of puffs at the same release site. We obtain analytically the probability laws of the interpuff time and of the number of channels that participate of the puffs. Furthermore, we show that under typical experimental conditions the effect of buffers can be accounted for in terms of a simple inhibiting function. Hence, by exploring different inhibiting functions we are able to study the effect of a variety of buffers on the puff size and interpuff time distributions. We find the somewhat counter-intuitive result that the addition of a fast Ca(2 +) buffer can increase the average number of channels that participate of a puff.

'Eventless' InsP3-dependent SR-Ca2+ release affecting atrial Ca2+ sparks

Augmented inositol 1,4,5-trisphosphate receptor (InsP3R) function has been linked to a variety of cardiac pathologies, including cardiac arrhythmia. The contribution of inositol 1,4,5-trisphosphate-induced Ca(2+) release (IP3ICR) in excitation-contraction coupling (ECC) under physiological conditions, as well as under cellular remodelling, remains controversial. Here we test the hypothesis that local IP3ICR directly affects ryanodine receptor (RyR) function and subsequent Ca(2+)-induced Ca(2+) release in atrial myocytes. IP3ICR was evoked by UV-flash photolysis of caged InsP3 under whole-cell configuration of the voltage-clamp technique in atrial myocytes isolated from C57/BL6 mice. Photolytic release of InsP3 was accompanied by a significant increase in the Ca(2+) release event frequency (4.14 ± 0.72 vs. 6.20 ± 0.76 events (100 μm)(-1) s(-1)). These individual photolytically triggered Ca(2+) release events were identified as Ca(2+) sparks, which originated from RyR openings. This was verified by Ca(2+) spark analysis and pharmacological separation between RyR and InsP3R-dependent sarcoplasmic reticulum (SR)-Ca(2+) release (2-aminoethoxydiphenyl borate, xestospongin C, tetracaine). Significant SR-Ca(2+) flux but eventless SR-Ca(2+) release through InsP3R were characterized using SR-Ca(2+) leak/SR-Ca(2+) load measurements. These results strongly support the idea that IP3ICR can effectively modulate RyR openings and Ca(2+) spark probability. We conclude that eventless and highly efficient InsP3-dependent SR-Ca(2+) flux is the main mechanism of functional cross-talk between InsP3Rs and RyRs, which may be an important factor in the modulation of ECC sensitivity.

IP3, a small molecule with a powerful message

Research conducted over the past two decades has provided convincing evidence that cell death, and more specifically apoptosis, can exceed single cell boundaries and can be strongly influenced by intercellular communication networks. We recently reported that gap junctions (i.e. channels directly connecting the cytoplasm of neighboring cells) composed of connexin43 or connexin26 provide a direct pathway to promote and expand cell death, and that inositol 1,4,5-trisphosphate (IP3) diffusion via these channels is crucial to provoke apoptosis in adjacent healthy cells. However, IP3 itself is not sufficient to induce cell death and additional factors appear to be necessary to create conditions in which IP3 will exert proapoptotic effects. Although IP3-evoked Ca(2+) signaling is known to be required for normal cell survival, it is also actively involved in apoptosis induction and progression. As such, it is evident that an accurate fine-tuning of this signaling mechanism is crucial for normal cell physiology, while a malfunction can lead to cell death. Here, we review the role of IP3 as an intracellular and intercellular cell death messenger, focusing on the endoplasmic reticulum-mitochondrial synapse, followed by a discussion of plausible elements that can convert IP3 from a physiological molecule to a killer substance. Finally, we highlight several pathological conditions in which anomalous intercellular IP3/Ca(2+) signaling might play a role. This article is part of a Special Issue entitled:12th European Symposium on Calcium.

Intracellular calcium signals display an avalanche-like behavior over multiple lengthscales

Many natural phenomena display “self-organized criticality” (SOC), (Bak et al., 1987). This refers to spatially extended systems for which patterns of activity characterized by different lengthscales can occur with a probability density that follows a power law with pattern size. Differently from power laws at phase transitions, systems displaying SOC do not need the tuning of an external parameter. Here we analyze intracellular calcium (Ca²⁺) signals, a key component of the signaling toolkit of almost any cell type. Ca²⁺ signals can either be spatially restricted (local) or propagate throughout the cell (global). Different models have suggested that the transition from local to global signals is similar to that of directed percolation. Directed percolation has been associated, in turn, to the appearance of SOC. In this paper we discuss these issues within the framework of simple models of Ca²⁺ signal propagation. We also analyze the size distribution of local signals (“puffs”) observed in immature Xenopus Laevis oocytes. The puff amplitude distribution obtained from observed local signals is not Gaussian with a noticeable fraction of large size events. The experimental distribution of puff areas in the spatio-temporal record of the image has a long tail that is approximately log-normal. The distribution can also be fitted with a power law relationship albeit with a smaller goodness of fit. The power law behavior is encountered within a simple model that includes some coupling among individual signals for a wide range of parameter values. An analysis of the model shows that a global elevation of the Ca²⁺ concentration plays a major role in determining whether the puff size distribution is long-tailed or not. This suggests that Ca²⁺-clearing from the cytosol is key to determine whether IP₃-mediated Ca²⁺ signals can display a SOC-like behavior or not.

Inositol 1,4,5-trisphosphate and its receptors

Activation of cells by many extracellular agonists leads to the production of inositol 1,4,5-trisphosphate (IP₃). IP₃ is a global messenger that easily diffuses in the cytosol. Its receptor (IP₃R) is a Ca(2+)-release channel located on intracellular membranes, especially the endoplasmic reticulum (ER). The IP₃R has an affinity for IP(3) in the low nanomolar range. A prime regulator of the IP₃R is the Ca(2+) ion itself. Cytosolic Ca(2+) is considered as a co-agonist of the IP₃R, as it strongly increases IP(3)R activity at concentrations up to about 300 nM. In contrast, at higher concentrations, cytosolic Ca(2+) inhibits the IP₃R. Also the luminal Ca(2+) sensitizes the IP₃R. In higher organisms three genes encode for an IP₃R and additional diversity exists as a result of alternative splicing mechanisms and the formation of homo- and heterotetramers. The various IP₃R isoforms have a similar structure and a similar function, but due to differences in their affinity for IP₃, their variable sensitivity to regulatory parameters, their differential interaction with associated proteins, and the variation in their subcellular localization, they participate differently in the formation of intracellular Ca(2+) signals and this affects therefore the physiological consequences of these signals.

Mean field strategies induce unrealistic non-linearities in calcium puffs

Mean field models are often useful approximations to biological systems, but sometimes, they can yield misleading results. In this work, we compare mean field approaches with stochastic models of intracellular calcium release. In particular, we concentrate on calcium signals generated by the concerted opening of several clustered channels (calcium puffs). To this end we simulate calcium puffs numerically and then try to reproduce features of the resulting calcium distribution using mean field models were all the channels open and close simultaneously. We show that an unrealistic non-linear relationship between the current and the number of open channels is needed to reproduce the simulated puffs. Furthermore, a single channel current which is five times smaller than the one of the stochastic simulations is also needed. Our study sheds light on the importance of the stochastic kinetics of the calcium release channel activity to estimate the release fluxes.