Applications of FLIKA, a Python-based image processing and analysis platform, for studying local events of cellular calcium signaling

Functional investigation of a putative calcium-binding site involved in the inhibition of inositol 1,4,5-trisphosphate receptor activity

The regulation of inositol 1,4,5-trisphosphate (IP3) receptor (IP3R) activity is thought to define the spatiotemporal patterns of Ca2+ signals necessary for the appropriate activation of downstream effectors. The binding of both IP3 and Ca2+ is obligatory for IP3R channel opening. Ca2+ however regulates IP3R activity in a biphasic manner. Ca2+ binding to a high-affinity pocket formed by the third armadillo repeat domain and linker domain promotes IP3R channel opening without altering the Ca2+ dependency for channel inactivation. These data suggest that a distinct low-affinity Ca2+-binding site is responsible for the reduction in IP3R activity at higher [Ca2+]. We mutated a cluster of acidic residues in the second armadillo repeat domain and central linker domain of IP3R type 1, reported to coordinate Ca2+ in the cryo-EM structures of the IP3R type 3. This "CD Ca2+-binding site" is well conserved in all IP3R subtypes. CD site Ca2+-binding mutants where the negatively charged glutamic acid residues were mutated to alanine exhibited enhanced sensitivity to IP3-generating agonists. Ca2+-binding mutants displayed spontaneous elemental Ca2+ puffs, and the number of IP3-induced Ca2+ puffs was augmented in cells stably expressing Ca2+-binding site mutants. The inhibitory effect of high [Ca2+] on single-channel open probability (Po) was reduced in mutant channels, and this effect was dependent on [ATP]. This indicates that Ca2+ binding to the putative CD Ca2+ inhibitory site facilitates the reduction in IP3R channel activation at subsaturating, likely physiological cytosolic [ATP], and suggest that at higher [ATP], additional Ca2+-binding motifs may contribute to the biphasic regulation of IP3-induced Ca2+ release.

Visualizing PIEZO1 Localization and Activity in hiPSC-Derived Single Cells and Organoids with HaloTag Technology

PIEZO1 is critical to numerous physiological processes, transducing diverse mechanical stimuli into electrical and chemical signals. Recent studies underscore the importance of visualizing endogenous PIEZO1 activity and localization to understand its functional roles. To enable physiologically and clinically relevant studies on human PIEZO1, we genetically engineered human induced pluripotent stem cells (hiPSCs) to express a HaloTag fused to endogenous PIEZO1. Combined with advanced imaging, our chemogenetic platform allows precise visualization of PIEZO1 localization dynamics in various cell types. Furthermore, the PIEZO1-HaloTag hiPSC technology facilitates the non-invasive monitoring of channel activity across diverse cell types using Ca2+-sensitive HaloTag ligands, achieving temporal resolution approaching that of patch clamp electrophysiology. Finally, we used lightsheet imaging of hiPSC-derived neural organoids to achieve molecular scale imaging of PIEZO1 in three-dimensional tissue organoids. Our advances offer a novel platform for studying PIEZO1 mechanotransduction in human cells and tissues, with potential for elucidating disease mechanisms and targeted therapeutic development.

Dual regulation of IP3 receptors by IP3 and PIP2 controls the transition from local to global Ca2+ signals

The spatial organization of inositol 1,4,5-trisphosphate (IP3)-evoked Ca2+ signals underlies their versatility. Low stimulus intensities evoke Ca2+ puffs, localized Ca2+ signals arising from a few IP3 receptors (IP3Rs) within a cluster tethered beneath the plasma membrane. More intense stimulation evokes global Ca2+ signals. Ca2+ signals propagate regeneratively as the Ca2+ released stimulates more IP3Rs. How is this potentially explosive mechanism constrained to allow local Ca2+ signaling? We developed methods that allow IP3 produced after G-protein coupled receptor (GPCR) activation to be intercepted and replaced by flash photolysis of a caged analog of IP3. We find that phosphatidylinositol 4,5-bisphosphate (PIP2) primes IP3Rs to respond by partially occupying their IP3-binding sites. As GPCRs stimulate IP3 formation, they also deplete PIP2, relieving the priming stimulus. Loss of PIP2 resets IP3R sensitivity and delays the transition from local to global Ca2+ signals. Dual regulation of IP3Rs by PIP2 and IP3 through GPCRs controls the transition from local to global Ca2+ signals.

Functional investigation of a putative calcium-binding site involved in the inhibition of inositol 1,4,5-trisphosphate receptor activity

A wide variety of factors influence inositol 1,4,5-trisphosphate (IP 3 ) receptor (IP 3 R) activity resulting in modulation of intracellular Ca 2+ release. This regulation is thought to define the spatio-temporal patterns of Ca 2+ signals necessary for the appropriate activation of downstream effectors. The binding of both IP 3 and Ca 2+ are obligatory for IP 3 R channel opening, however, Ca 2+ regulates IP 3 R activity in a biphasic manner. Mutational studies have revealed that Ca 2+ binding to a high-affinity pocket formed by the ARM3 domain and linker domain promotes IP 3 R channel opening without altering the Ca 2+ dependency for channel inactivation. These data suggest a distinct low-affinity Ca 2+ binding site is responsible for the reduction in IP 3 R activity at higher [Ca 2+ ]. We determined the consequences of mutating a cluster of acidic residues in the ARM2 and central linker domain reported to coordinate Ca 2+ in cryo-EM structures of the IP 3 R type 3. This site is termed the "CD Ca 2+ binding site" and is well-conserved in all IP 3 R sub-types. We show that the CD site Ca 2+ binding mutants where the negatively charged glutamic acid residues are mutated to alanine exhibited enhanced sensitivity to IP 3 -generating agonists. Ca 2+ binding mutants displayed spontaneous elemental Ca 2+ events (Ca 2+ puffs) and the number of IP 3 -induced Ca 2+ puffs was significantly augmented in cells stably expressing Ca 2+ binding site mutants. When measured with "on-nucleus" patch clamp, the inhibitory effect of high [Ca 2+ ] on single channel-open probability (P o ) was reduced in mutant channels and this effect was dependent on [ATP]. These results indicate that Ca 2+ binding to the putative CD Ca 2+ inhibitory site facilitates the reduction in IP 3 R channel activation when cytosolic [ATP] is reduced and suggest that at higher [ATP], additional Ca 2+ binding motifs may contribute to the biphasic regulation of IP 3 -induced Ca 2+ release.

Dissociation of inositol 1,4,5-trisphosphate from IP3 receptors contributes to termination of Ca2+ puffs

Ca²⁺ puffs are brief, localized Ca²⁺ signals evoked by physiological stimuli that arise from the coordinated opening of a few clustered inositol 1,4,5-trisphosphate receptors (IP₃Rs). However, the mechanisms that control the amplitude and termination of Ca²⁺ puffs are unresolved. To address these issues, we expressed SNAP-tagged IP₃R3 in HEK cells without endogenous IP₃Rs and used total internal reflection fluorescence microscopy to visualize the subcellular distribution of IP₃Rs and the Ca²⁺ puffs that they evoke. We first confirmed that SNAP-IP₃R3 were reliably identified and that they evoked normal Ca²⁺ puffs after photolysis of a caged analog of IP₃. We show that increased IP₃R expression caused cells to assemble more IP₃R clusters, each of which contained more IP₃Rs, but the mean amplitude of Ca²⁺ puffs (indicative of the number of open IP₃Rs) was unaltered. We thus suggest that functional interactions between IP₃Rs constrain the number of active IP₃Rs within a cluster. Furthermore, Ca²⁺ puffs evoked by IP₃R with reduced affinity for IP₃ had undiminished amplitude, but the puffs decayed more quickly. The selective effect of reducing IP₃ affinity on the decay times of Ca²⁺ puffs was not mimicked by exposing normal IP₃R to a lower concentration of IP₃. We conclude that distinct mechanisms constrain recruitment of IP₃Rs during the rising phase of a Ca²⁺ puff and closure of IP₃Rs during the falling phase, and that only the latter is affected by the rate of IP₃ dissociation.

Algorithm for biological second messenger analysis with dynamic regions of interest

Physiological function is regulated through cellular communication that is facilitated by multiple signaling molecules such as second messengers. Analysis of signal dynamics obtained from cell and tissue imaging is difficult because of intricate spatially and temporally distinct signals. Signal analysis tools based on static region of interest analysis may under- or overestimate signals in relation to region of interest size and location. Therefore, we developed an algorithm for biological signal detection and analysis based on dynamic regions of interest, where time-dependent polygonal regions of interest are automatically assigned to the changing perimeter of detected and segmented signals. This approach allows signal profiles to be rigorously and precisely tracked over time, eliminating the signal distortion observed with static methods. Integration of our approach with state-of-the-art image processing and particle tracking pipelines enabled the isolation of dynamic cellular signaling events and characterization of biological signaling patterns with distinct combinations of parameters including amplitude, duration, and spatial spread. Our algorithm was validated using synthetically generated datasets and compared with other available methods. Application of the algorithm to volumetric time-lapse hyperspectral images of cyclic adenosine monophosphate measurements in rat microvascular endothelial cells revealed distinct signal heterogeneity with respect to cell depth, confirming the utility of our approach for analysis of 5-dimensional data. In human tibial arteries, our approach allowed the identification of distinct calcium signal patterns associated with atherosclerosis. Our algorithm for automated detection and analysis of second messenger signals enables the decoding of signaling patterns in diverse tissues and identification of pathologic cellular responses.

Satellite Image Processing by Python and R Using Landsat 9 OLI/TIRS and SRTM DEM Data on Côte d'Ivoire, West Africa

In this paper, we propose an advanced scripting approach using Python and R for satellite image processing and modelling terrain in Côte d'Ivoire, West Africa. Data include Landsat 9 OLI/TIRS C2 L1 and the SRTM digital elevation model (DEM). The EarthPy library of Python and 'raster' and 'terra' packages of R are used as tools for data processing. The methodology includes computing vegetation indices to derive information on vegetation coverage and terrain modelling. Four vegetation indices were computed and visualised using R: the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index 2 (EVI2), Soil-Adjusted Vegetation Index (SAVI) and Atmospherically Resistant Vegetation Index 2 (ARVI2). The SAVI index is demonstrated to be more suitable and better adjusted to the vegetation analysis, which is beneficial for agricultural monitoring in Côte d'Ivoire. The terrain analysis is performed using Python and includes slope, aspect, hillshade and relief modelling with changed parameters for the sun azimuth and angle. The vegetation pattern in Côte d'Ivoire is heterogeneous, which reflects the complexity of the terrain structure. Therefore, the terrain and vegetation data modelling is aimed at the analysis of the relationship between the regional topography and environmental setting in the study area. The upscaled mapping is performed as regional environmental analysis of the Yamoussoukro surroundings and local topographic modelling of the Kossou Lake. The algorithms of the data processing include image resampling, band composition, statistical analysis and map algebra used for calculation of the vegetation indices in Côte d'Ivoire. This study demonstrates the effective application of the advanced programming algorithms in Python and R for satellite image processing.

KRAP is required for diffuse and punctate IP3-mediated Ca2+ liberation and determines the number of functional IP3R channels within clusters

KRas-induced actin-interacting protein (KRAP) has been identified as crucial for the appropriate localization and functioning of the inositol trisphosphate receptors (IP₃Rs) that mediate Ca²⁺ release from the endoplasmic reticulum. Here, we used siRNA knockdown of KRAP expression in HeLa and HEK293 cells to examine the roles of KRAP in the generation of IP₃-mediated local Ca²⁺ puffs and global, cell-wide Ca²⁺ signals. High resolution Ca²⁺ imaging revealed that the mean amplitude of puffs was strongly reduced by KRAP knockdown, whereas the Ca²⁺ flux during openings of individual IP₃R channels was little affected. In both control and KRAP knockdown cells the numbers of functional channels in the clusters underlying puff sites were stochastically distributed following a Poisson relationship, but the mean number of functional channels per site was reduced by about two thirds by KRAP knockdown. We conclude that KRAP is required for activity of IP₃R channels at puff sites and stochastically 'licenses' the function of individual channels on a one-to-one basis, rather than determining the functioning of the puff site as a whole. In addition to puff activity ('punctate' Ca²⁺ release), global, cell-wide Ca²⁺ signals evoked by higher levels of IP₃ are further composed from a discrete 'diffuse' mode of Ca²⁺ release. By applying fluctuation analysis to isolate the punctate component during global Ca²⁺ signals, we find that KRAP knockdown suppresses to similar extents punctate and diffuse Ca²⁺ release in wild-type cells and in HEK293 cells exclusively expressing type 1 and type 3 IP3Rs. Thus, KRAP appears essential for the functioning of the IP₃Rs involved in diffuse Ca²⁺ release as well as the clustered IP₃Rs that generate local Ca²⁺ puffs.

Termination of Ca2+ puffs during IP3-evoked global Ca2+ signals

We previously described that cell-wide cytosolic Ca²⁺ transients evoked by inositol trisphosphate (IP₃) are generated by two modes of Ca²⁺ liberation from the ER; 'punctate' release via an initial flurry of transient Ca²⁺ puffs from local clusters of IP₃ receptors, succeeded by a spatially and temporally 'diffuse' Ca²⁺ liberation. Those findings were derived using statistical fluctuation analysis to monitor puff activity which is otherwise masked as global Ca²⁺ levels rise. Here, we devised imaging approaches to resolve individual puffs during global Ca²⁺ elevations to better investigate the mechanisms terminating the puff flurry. We find that puffs contribute about 40% (∼90 attomoles) of the total Ca²⁺ liberation, largely while the global Ca²⁺ signal rises halfway to its peak. The major factor terminating punctate Ca²⁺ release is an abrupt decline in puff frequency. Although the amplitudes of large puffs fall during the flurry, the amplitudes of more numerous small puffs remain steady, so overall puff amplitudes decline only modestly (∼30%). The Ca²⁺ flux through individual IP₃ receptor/channels does not measurably decline during the flurry, or when puff activity is depressed by pharmacological lowering of Ca²⁺ levels in the ER lumen, indicating that the termination of punctate release is not a simple consequence of reduced driving force for Ca²⁺ liberation. We propose instead that the gating of IP₃ receptors at puff sites is modulated such that their openings become suppressed as the bulk [Ca²⁺] in the ER lumen falls during global Ca²⁺ signals.

A protocol for detecting elemental calcium signals (Ca2+ puffs) in mammalian cells using total internal reflection fluorescence microscopy

This protocol outlines steps to visualize and detect Ca2+ puffs following photo-liberation of caged inositol-1,4,5-trisphosphate (IP3) from HEK-293 cells expressing only the native IP3R type 1 receptor using total internal reflection fluorescence (TIRF) microscopy. TIRF microscopy offers high axial resolution and allows imaging at high speed, with a higher signal-to-background ratio. Additionally, we shed light on commonly encountered pitfalls, which should be considered while recording Ca2+ puffs using TIRF microscopy. For complete details on the use and execution of this protocol, please refer to Emrich et al. (2021) and Lock et al. (2015a).

Mechanisms and functions of calcium microdomains produced by ORAI channels, d-myo-inositol 1,4,5-trisphosphate receptors, or ryanodine receptors

With the discovery of local Ca²⁺ signals in the 1990s the concept of 'elementary Ca²⁺ signals' and 'fundamental Ca²⁺ signals' was developed. While 'elementary Ca²⁺signals' relate to optical signals gained by activity of small clusters of Ca²⁺channels, 'fundamental signals' describe such optical signals that arise from opening of single Ca²⁺channels. In this review, we discuss (i) concepts of local Ca²⁺ signals and Ca²⁺ microdomains, (ii) molecular mechanisms underlying Ca²⁺ microdomains, (iii) functions of Ca²⁺ microdomains, and (iv) mathematical modelling of Ca²⁺ microdomains. We focus on Ca²⁺ microdomains produced by ORAI channels, D-myo-inositol 1,4,5-trisphosphate receptors, or ryanodine receptors. In summary, research on local Ca²⁺ signals in different cell models aims to better understand how cells use the Ca²⁺ toolkit to produce Ca²⁺ microdomains as relevant signals for specific cellular responses, but also how local Ca²⁺ signals as building blocks merge into global Ca²⁺ signaling.

PunctaSpecks: A tool for automated detection, tracking, and analysis of multiple types of fluorescently labeled biomolecules

Recent advances in imaging technology and fluorescent probes have made it possible to gain information about the dynamics of subcellular processes at unprecedented spatiotemporal scales. Unfortunately, a lack of automated tools to efficiently process the resulting imaging data encoding fine details of the biological processes remains a major bottleneck in utilizing the full potential of these powerful experimental techniques. Here we present a computational tool, called PunctaSpecks, that can characterize fluorescence signals arising from a wide range of biological molecules under normal and pathological conditions. Among other things, the program can calculate the number, areas, life-times, and amplitudes of fluorescence signals arising from multiple sources, track diffusing fluorescence sources like moving mitochondria, and determine the overlap probability of two processes or organelles imaged using indicator dyes of different colors. We have tested PunctaSpecks on synthetic time-lapse movies containing mobile fluorescence objects of various sizes, mimicking the activity of biomolecules. The robustness of the software is tested by varying the level of noise along with random but known pattern of appearing, disappearing, and movement of these objects. Next, we use PunctaSpecks to characterize protein-protein interaction involved in store-operated Ca2+ entry through the formation and activation of plasma membrane-bound ORAI1 channel and endoplasmic reticulum membrane-bound stromal interaction molecule (STIM), the evolution of inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ signals from sub-micrometer size local events into global waves in human cortical neurons, and the activity of Alzheimer's disease-associated β amyloid pores in the plasma membrane. The tool can also be used to study other dynamical processes imaged through fluorescence molecules. The open source algorithm allows for extending the program to analyze more than two types of biomolecules visualized using markers of different colors.

ER-luminal [Ca2+] regulation of InsP3 receptor gating mediated by an ER-luminal peripheral Ca2+-binding protein

Modulating cytoplasmic Ca2+ concentration ([Ca2+]i) by endoplasmic reticulum (ER)-localized inositol 1,4,5-trisphosphate receptor (InsP3R) Ca2+-release channels is a universal signaling pathway that regulates numerous cell-physiological processes. Whereas much is known regarding regulation of InsP3R activity by cytoplasmic ligands and processes, its regulation by ER-luminal Ca2+ concentration ([Ca2+]ER) is poorly understood and controversial. We discovered that the InsP3R is regulated by a peripheral membrane-associated ER-luminal protein that strongly inhibits the channel in the presence of high, physiological [Ca2+]ER. The widely-expressed Ca2+-binding protein annexin A1 (ANXA1) is present in the nuclear envelope lumen and, through interaction with a luminal region of the channel, can modify high-[Ca2+]ER inhibition of InsP3R activity. Genetic knockdown of ANXA1 expression enhanced global and local elementary InsP3-mediated Ca2+ signaling events. Thus, [Ca2+]ER is a major regulator of InsP3R channel activity and InsP3R-mediated [Ca2+]i signaling in cells by controlling an interaction of the channel with a peripheral membrane-associated Ca2+-binding protein, likely ANXA1.

IP3 mediated global Ca2+ signals arise through two temporally and spatially distinct modes of Ca2+ release

The 'building-block' model of inositol trisphosphate (IP3)-mediated Ca2+ liberation posits that cell-wide cytosolic Ca2+ signals arise through coordinated activation of localized Ca2+ puffs generated by stationary clusters of IP3 receptors (IP3Rs). Here, we revise this hypothesis, applying fluctuation analysis to resolve Ca2+ signals otherwise obscured during large Ca2+ elevations. We find the rising phase of global Ca2+ signals is punctuated by a flurry of puffs, which terminate before the peak by a mechanism involving partial ER Ca2+ depletion. The continuing rise in Ca2+, and persistence of global signals even when puffs are absent, reveal a second mode of spatiotemporally diffuse Ca2+ signaling. Puffs make only small, transient contributions to global Ca2+ signals, which are sustained by diffuse release of Ca2+ through a functionally distinct process. These two modes of IP3-mediated Ca2+ liberation have important implications for downstream signaling, imparting spatial and kinetic specificity to Ca2+-dependent effector functions and Ca2+ transport.

A Primary Cell and Organoid Platform for Evaluating Pharmacological Responses in Mammary Epithelial Cells

An essential process in predicting the in vivo pharmacological activity of a candidate molecule involves the evaluation of target responses using established model systems. While these models largely comprise immortalized cells, which are often serially passaged as monolayers on uniformly stiff substrates and are modified to overexpress one or more components of the pathway-of-interest, the importance of cell identity, heterogeneity, and three-dimensional (3D) context to target response is gaining increasing attention. Here, we assess intracellular calcium responses in mouse mammary epithelial cells in three distinct model systems: 3D primary organoids, 2D primary epithelial cells, and 2D immortalized cells. Specifically, we assess intracellular calcium responses to a number of extracellular signals implicated in the regulation of basal (or myoepithelial) cell function. These findings provide further insights into cell type and context-specific pharmacological responses in mammary epithelial cells and highlight the opportunities and challenges in the adoption of architecturally complex and heterogeneous in vitro assays in pharmacological research.

Spatial-temporal patterning of Ca2+ signals by the subcellular distribution of IP3 and IP3 receptors

The patterning of cytosolic Ca²⁺ signals in space and time underlies their ubiquitous ability to specifically regulate numerous cellular processes. Signals mediated by liberation of Ca²⁺ sequestered in the endoplasmic reticulum (ER) through inositol trisphosphate receptor (IP₃R) channels constitute a hierarchy of events; ranging from openings of individual IP₃ channels, through the concerted openings of several clustered IP₃Rs to generate local Ca²⁺ puffs, to global Ca²⁺ waves and oscillations that engulf the entire cell. Here, we review recent progress in elucidating how this hierarchy is shaped by an interplay between the functional gating properties of IP₃Rs and their spatial distribution within the cell. We focus in particular on the subset of IP₃Rs that are organized in stationary clusters and are endowed with the ability to preferentially liberate Ca²⁺.