Research Article

Spatially compartmentalized phase regulation of a Ca²⁺-cAMP-PKA oscillatory circuit

Department of Biophysics and Biophysical Chemistry, The Johns Hopkins University School of Medicine, United States
Department of Pharmacology, University of California, San Diego, United States
Chemical Engineering Graduate Program, University of California, San Diego, United States
Department of Mechanical and Aerospace Engineering, University of California, San Diego, United States
Department of Chemistry and Biochemistry, University of California, San Diego, United States

Nov 17, 2020

Open access
Copyright information

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Signaling networks are spatiotemporally organized to sense diverse inputs, process information, and carry out specific cellular tasks. In β cells, Ca²⁺, cyclic adenosine monophosphate (cAMP), and Protein Kinase A (PKA) exist in an oscillatory circuit characterized by a high degree of feedback. Here, we describe a mode of regulation within this circuit involving a spatial dependence of the relative phase between cAMP, PKA, and Ca²⁺. We show that in mouse MIN6 β cells, nanodomain clustering of Ca²⁺-sensitive adenylyl cyclases (ACs) drives oscillations of local cAMP levels to be precisely in-phase with Ca²⁺ oscillations, whereas Ca²⁺-sensitive phosphodiesterases maintain out-of-phase oscillations outside of the nanodomain. Disruption of this precise phase relationship perturbs Ca²⁺ oscillations, suggesting the relative phase within an oscillatory circuit can encode specific functional information. This work unveils a novel mechanism of cAMP compartmentation utilized for localized tuning of an oscillatory circuit and has broad implications for the spatiotemporal regulation of signaling networks.

Introduction

Cyclic adenosine monophosphate (cAMP) and Ca²⁺ act as essential second messengers in almost every cell type and regulate many functional pathways within a cell, such as hormonal signal transduction, metabolism, and secretion (Clapham, 2007; Sassone-Corsi, 2012). In some cell types, including neurons, cardiomyocytes, and pancreatic β cells, these messengers’ concentrations oscillate intracellularly (Dupont et al., 2011; Dyachok et al., 2006), and the oscillations encode critical signaling information (e.g. signal strength, duration, and target diversity) into parameters such as frequency and amplitude (Berridge et al., 1998; De Pittà et al., 2008; Parekh, 2011). This phenomenon is perhaps best exemplified in β cells, where oscillations of Ca²⁺ drive pulsatile insulin secretion (Rorsman and Ashcroft, 2018) and oscillating cAMP levels (Nesher et al., 2002; Tengholm, 2012). Furthermore, Ca²⁺, cAMP, and the downstream cAMP-dependent kinase protein kinase (PKA) constitute a highly coordinated oscillatory circuit responsible for integrating metabolic and signaling information (Ni et al., 2011). In addition to temporal control, biochemical pathways are also spatially organized within the cell (Smith and Scott, 2002; White and Anderson, 2005). Both Ca²⁺ and cAMP are highly spatially compartmentalized and form signaling microdomains or nanodomains (Calebiro and Maiellaro, 2014; Petersen, 2002). While Ca²⁺ levels are locally controlled by channels, pumps, and intracellular buffering systems (Clapham, 2007; Stern, 1992), cAMP is thought to be regulated via controlled synthesis by adenylyl cyclases (ACs) and degradation by phosphodiesterases (PDEs) (Bender and Beavo, 2006; Defer et al., 2000). Despite extensive studies on cAMP compartmentation, the mechanisms that spatially constrain this mobile second messenger await to be fully elucidated (Saucerman et al., 2014; Lohse et al., 2017; Musheshe et al., 2018; Bock et al., 2020; Zhang et al., 2020). Furthermore, it is not clear how the spatial regulation of a second messenger influences its dynamic behaviors in the context of coordinated oscillations.

In this study, we investigated the spatiotemporal organization of the Ca²⁺-cAMP-PKA oscillatory circuit in MIN6 β cells and discovered that the relative oscillatory phase between cAMP/PKA and Ca²⁺ is spatially regulated within signaling nanodomains. By combining dynamic live-cell imaging, super-resolution microscopy, and computational modeling, we further found that fine-scale, compartment-specific perturbations of this precise phase regulation impact Ca²⁺ oscillations in β cells. These findings suggest that the relative phase in oscillatory signaling circuits, like the amplitude and frequency, represents yet another mode of informational encoding and processing, which is subjected to spatiotemporal regulation within the cell.

Results

The relative phase of β cell cAMP and Ca²⁺ oscillations is compartmentalized

In order to study the spatiotemporal relationship between key players of the Ca²⁺-cAMP-PKA circuit, we chose to focus our attention on an important class of molecular scaffolds, A-kinase anchoring proteins (AKAPs), which are responsible for recruiting PKA to specific substrates at distinct subcellular locations. In several excitable cell types, the plasma membrane (PM)-localized scaffold protein AKAP79 (rodent ortholog AKAP150) has been shown to organize a macromolecular complex with binding partners that include PKA, the voltage-gated Ca²⁺ channel Ca_V1.2, Protein Kinase C (PKC), the Ca²⁺/calmodulin-dependent protein phosphatase calcineurin, Ca²⁺-sensitive ACs, AMPA receptors, and many others (Gold et al., 2011). Due to the extensive and multivalent nature of AKAP79/150 (as the scaffold is commonly referred) and a report describing the functional impairment of glucose-stimulated insulin secretion (GSIS) in pancreatic β cells upon its knock-out (Hinke et al., 2012), we hypothesized that the AKAP79/150 scaffold might play an important role in the spatiotemporal regulation of the Ca²⁺-cAMP-PKA oscillatory circuit. Specifically, we were interested in testing if AKAP79/150 is able to create a spatially-distinct compartment in which recruitment of signaling effectors can locally fine-tune and reshape signaling dynamics within the circuit (Beene and Scott, 2007; Greenwald and Saucerman, 2011).

To test this hypothesis, we monitored intracellular cAMP and Ca²⁺ using the FRET-based cAMP biosensor (Ci/Ce)Epac2-camps (Everett and Cooper, 2013) and the red Ca²⁺ indicator RCaMP (Akerboom et al., 2013). We measured cAMP concentration changes in the immediate vicinity of AKAP150 in mouse MIN6 β cells by using (Ci/Ce)Epac2-camps fused to the full-length AKAP79 scaffold (gene AKAP5) (Figure 1a). Although human AKAP79 and the rodent AKAP150 share only 53% sequence identity overall, the interaction motifs and association between the AKAP scaffold and other key signaling players such as PKA, the voltage-gated calcium channel, adenylyl cyclases, and calcineurin are highly conserved (Willoughby et al., 2010; Zhang et al., 2016). The primary difference between these two closely related scaffolds is the presence of an internal repetitive amino acid sequence of unknown function in AKAP150 (Robertson et al., 2009). The functional equivalence of AKAP79 and AKAP150 has been demonstrated in several reciprocal knock-out and recovery experiments in which AKAP79 or AKAP150 was knocked out and expression of the other was shown to rescue a measured phenotype (Hoshi et al., 2005; Zhang et al., 2011). As a control, we also targeted the cAMP probe to the general plasma membrane by adding a lipid modification domain (Wachten et al., 2010). These targeted biosensors allowed us to compare the dynamics within the AKAP79/150-specific compartment versus the general plasma membrane compartment (Figure 1a).

Figure 1 with 2 supplements see all

Download asset Open asset

The phase of oscillating cAMP is shifted between the AKAP79/150 compartment and the general plasma membrane compartment, relative to Ca²⁺.

(a) Depiction of the AKAP79/150 and plasma membrane compartments, including the targeted cAMP biosensor (Ci/Ce)Epac2-camps to measure the compartment-specific cAMP signaling. Schematics of the lyn-(Ci/Ce)Epac2-camps and AKAP79-(Ci/Ce)Epac2-camps sensors are shown above. (b) Network diagram describing the key players in the β cell Ca²⁺-cAMP-PKA oscillatory circuit. (c) Representative single-cell trace of in-phase oscillating responses of AKAP79-(Ci/Ce)Epac2-camps and RCaMP, whole-cell fluorescence measured. Red trace is cAMP (cyan direct channel divided by CY-FRET channel) and black trace is Ca²⁺ (RFP). (d) Representative single-cell trace of out-of-phase oscillating responses of lyn-(Ci/Ce)Epac2-camps and RCaMP, whole-cell fluorescence measured. Blue trace is cAMP (cyan direct channel divided by CY-FRET channel) and black trace is Ca²⁺ (RFP). (e) Cross-correlation between the oscillatory Ca²⁺ and cAMP signals from the representative in-phase AKAP79 (red) and out-of-phase plasma membrane (PM, blue) responses from c, d. Time lag (sec) between the cAMP and Ca²⁺ signals for the two compartments (AKAP79/150, red, is 13 ± 3 sec n=60 and PM, blue, is 47 ± 4 s n=24). ****p<0.0001; unpaired two-tailed Student’s t-test.

Although both targeted sensors were trafficked to and distributed along the plasma membrane (Figure 1—figure supplement 1a,b), we observed notable differences in their respective cAMP signals relative to Ca²⁺ oscillations after triggering the circuit (Figure 1b) with tetraethylammonium chloride (TEA, 20 mM), a potent K⁺ channel blocker. cAMP oscillations measured within the AKAP79/150 compartment were in-phase with oscillating Ca²⁺, such that each transient spike in intracellular Ca²⁺ was closely associated with a transient increase in cAMP (Figure 1c) (n = 60 cells). This was in sharp contrast to cAMP oscillations measured within the general plasma membrane compartment, where each local Ca²⁺ peak corresponded to a local trough in cAMP (n = 24), followed by a slow reversal of both signals to a pre-stimulated baseline (Figure 1d). While these out-of-phase cAMP-Ca²⁺ oscillations were consistent with those observed in the cytoplasm of β cells (Landa et al., 2005; Ni et al., 2011), in-phase cAMP-Ca²⁺ oscillations had not previously been observed under these conditions. To quantify the cAMP-Ca²⁺ phase relationship, we measured the lag time by calculating the cross-correlation between the two normalized, oscillatory signals and finding the shortest delay yielding the maximum correlation (see Appendix 1 for details) (Figure 1e). In-phase cAMP oscillations corresponded to short lag times (typically <20 s), while out-of-phase oscillations mostly possessed longer lag times. Within the AKAP79/150 compartment, cAMP lagged behind Ca²⁺ by an average of only 13 ± 3 s (mean ± SEM, n = 60); however, cAMP within the general plasma membrane compartment oscillated with a lag time of 47 ± 4 s (n = 24) behind Ca²⁺ (Figure 1e). In order to measure these different cAMP dynamics within the same cells, we transfected MIN6 cells with two cAMP sensors of different targeting sequences and colors, AKAP79-(Ci/Ce)Epac2-camps and the red cytosolic cAMP sensor R-FlincA (Ohta et al., 2018). After stimulating the signaling circuit, we observed anti-correlated cAMP oscillations, indicating that cAMP oscillates with different phases between the cytosol/plasma membrane and the immediate vicinity of AKAP79/150 within the same cell (n = 25) (Figure 1—figure supplement 2). This stark difference in the cAMP-Ca²⁺ phase relationship suggests that the relative phase of this oscillatory circuit is compartmentalized and hints at differential regulation of the circuit between the AKAP79/150 compartment and the general plasma membrane compartment.

Oscillatory phase is regulated by balanced activities of Ca²⁺-sensitive ACs and PDEs

Given that in-phase cAMP oscillations were only observed within the AKAP79/150 compartment (Figure 1c) and that out-of-phase cAMP oscillations were observed in the general plasma membrane compartment (Figure 1d) and the cytoplasm (Ni et al., 2011), we hypothesized that Ca²⁺ oscillations are coupled to cAMP oscillations by a ubiquitous mechanism throughout the cell, while additional mechanisms specifically regulate the phase relationship within the AKAP79/150 compartment. We first sought to identify the component that is responsible for coupling cAMP dynamics to Ca²⁺ dynamics globally. Since TEA induces continuous Ca²⁺ oscillations, we determined the temporal relationship between Ca²⁺ and cAMP at the general plasma membrane more precisely by measuring the impulse response of the circuit following a transient membrane depolarization. After the addition of KCl (15 mM) followed by a subsequent washout to elicit a transient influx of Ca²⁺ (Dou et al., 2015), we observed a synchronous cAMP decrease (n = 20) followed by a return to baseline (Figure 2a). These data suggest that increasing cytosolic Ca²⁺ was coupled to a decrease in cAMP at the plasma membrane through Ca²⁺-sensitive AC or PDE activities. Given that Ca²⁺-inhibited ACs (AC5, AC6) have low specific activity in both the presence and absence of physiological Ca²⁺, as well as relatively low expression in the pancreas (Defer et al., 2000), we instead focused on probing the role of PDEs. The Ca²⁺-dependent PDE1 family in β cells, specifically PDE1C, has been implicated in modulating GSIS (Han et al., 1999). Indeed, acute addition of 8MM-IBMX (100 μM), a relatively selective PDE1 inhibitor, effectively uncoupled cAMP dynamics from Ca²⁺ oscillations (Figure 2b, Figure 2—figure supplement 1a) (n = 18), indicating that Ca²⁺-triggered activation of PDE1 mediates the transient cAMP decreases. We also observed that the overall increase in cAMP led to an increase in the Ca²⁺ oscillation frequency, consistent with the previously identified role of cAMP/PKA in regulating the Ca²⁺ oscillations (Ni et al., 2011). We tested the roles of two additional families of abundant PDEs in pancreatic β cells, PDE3 and PDE4, by acute pharmacologic inhibition. While treating cells with either milrinone (PDE3 inhibitor, 10 μM, n = 12) or rolipram (PDE4 inhibitor, 1 μM, n = 15) slightly increased cAMP levels, neither inhibitor had an effect on cAMP-Ca²⁺ coupling or relative phase (Figure 2—figure supplement 1b,c). These data suggest that PDE1 is the key component that couples Ca²⁺ and cAMP oscillations within this signaling circuit.

Figure 2 with 3 supplements see all

Download asset Open asset

The oscillation phase is regulated by a balance between Ca²⁺-sensitive AC and PDE activity.

(a) Impulse response of plasma membrane cAMP (blue) to a spike in Ca²⁺ entry (black), triggered by KCl-mediated membrane depolarization (wash in/out). The transient decrease in PM-cAMP is coupled to the transient increase in intracellular Ca²⁺. (b) Acute inhibition of Ca²⁺-sensitive PDE1 decouples the out-of-phase PM-cAMP oscillations from Ca²⁺ oscillations, as observed in this representative cell trace (Ca²⁺ – black, PM-cAMP – blue). (c) The oscillatory phase of cAMP can be manipulated by tuning the relative activity of Ca²⁺-sensitive PDE and AC, as demonstrated by a simplified mathematical model. The schematic shows the network architecture. The relative activities of AC and PDE, denoted as v(AC) and v(PDE), control the delay between Ca²⁺ and cAMP. Point (i), with high AC activity and low PDE, shows in phase oscillations of Ca²⁺ and cAMP, whereas point (ii), with high PDE and low AC activity shows out of phase oscillations. These dynamics are shown in the line graphs. (d) Knocking down AC8 is correlated with an increase in the time lag for oscillatory cAMP at the AKAP79/150 microdomain (37 ± 9 s, n = 11), indicating more cells exhibiting out-of-phase cAMP oscillations (representative cell trace, Ca²⁺ – black, AKAP79/150-cAMP – red). (e) Over-expressing AC8 is sufficient to reverse the phase at the PM to in-phase (23 ± 2 s, n = 56) (representative cell trace, Ca²⁺ – black, PM cAMP – blue). **p=0.0014, ****p<0.0001; unpaired two-tailed Student’s t-test.

How is the phase relationship between Ca^2+- and cAMP-regulated within distinct signaling compartments? In order to gain a more quantitative understanding of the regulation of the cAMP-Ca²⁺ phase relationship, we created a simplified well-mixed mathematical model involving Ca²⁺, cAMP, and Ca²⁺-driven PDE and AC activity components (Cooper et al., 1995; Figure 2c, see Appendix 1 for details). This simple circuit represents the key aspects of the oscillatory cAMP-Ca²⁺ circuit and is applicable to different signaling compartments. Opposite to the Ca²⁺-stimulated PDE1 (Ang and Antoni, 2002) is the Ca²⁺-stimulated AC8 (gene Adcy8) (Masada et al., 2009; Masada et al., 2012), an abundant Ca²⁺-sensitive transmembrane AC isoform in β cells that has been shown to mediate sustained insulin secretion and associate with the AKAP79/150 scaffold (Dou et al., 2015; Willoughby and Cooper, 2006; Willoughby et al., 2010). By computationally manipulating the activity of each arm, we found that cAMP can oscillate either out-of-phase or in-phase when a Ca²⁺ pulse train is used as an input (Figure 2c). In particular, when the relative activity of PDE1 is greater than the activity of AC8, Ca²⁺-driven cAMP degradation dominates, resulting in an out-of-phase cAMP-Ca²⁺ relationship. On the other hand, if the relative activity of AC8 is greater than that of PDE1, Ca²⁺-stimulated cAMP production is favored, and an in-phase relationship is observed, consistent with previous modeling studies (Fridlyand et al., 2007; Peercy et al., 2015).

Thus, our simplified model indicates that the phase relationship can be tuned by altering the relative strength between Ca²⁺-sensitive ACs and PDEs (Figure 2c). This model provided a blueprint for understanding the interplay between the Ca²⁺-stimulated AC/PDE balance and the cAMP-Ca²⁺ phase relationship within the AKAP79/150 compartment. Based on the findings from our model, we predicted that decreasing the relative contribution of AC8 will shift the cAMP-Ca²⁺ phase relationship from in-phase to out-of-phase. To test this prediction, we knocked down endogenous AC8 in MIN6 cells as previously done (Raoux et al., 2015) and observed that most cells exhibited out-of-phase cAMP oscillations within the AKAP79/150 compartment (average lag time 37 ± 9 s, n = 11) (Figure 2d), indicating an AC8-specific role in mediating the cAMP-Ca²⁺ phase signature.

Conversely, increasing the relative contribution of AC8, for example, by increasing the concentration of AC8 throughout the PM, should shift the cAMP-Ca²⁺ phase relationship from out-of-phase to in-phase. To test this prediction, we overexpressed full-length AC8 and examined the effect in the general plasma membrane compartment. Interestingly, we found that AC8 overexpression reversed the out-of-phase cAMP-Ca²⁺ phase relationship in a titratable manner, where the percentage of in-phase oscillating cells correlated with increasing amounts of co-transfected AC8 (average lag time 23 ± 2 s, n = 56) (Figure 2e, Figure 2—figure supplement 2a–d). In order to examine the effects of the targeted cAMP biosensors on AC8 expression, we also measured endogenous AC8 levels in cells transfected with either lyn-(Ci/Ce)Epac2-camps or AKAP79-(Ci/Ce)Epac2-camps and observed no significant difference (Figure 2—figure supplement 3a,b). These data demonstrate that higher levels of AC8 are sufficient to reverse the cAMP phase at the plasma membrane. In summary, these phase-manipulation experiments suggest that the cAMP-Ca²⁺ phase relationship is representative of a sensitive, compartmentalized balance between the Ca²⁺-stimulated activities of PDE1 and AC8.

Membrane-localized AKAP150:AC8 nanoclusters regulate cAMP-Ca²⁺ oscillatory phase

The close spatial juxtaposition between the AKAP79/150 and general plasma membrane compartments presents a significant challenge for cAMP compartmentation, in that cAMP oscillations must be distinctly regulated within these adjacent signaling domains. Indeed, how cAMP, a rapidly diffusing small molecule, is spatially compartmentalized in cells is not yet completely understood, especially given the relatively low catalytic efficiency of a single cAMP-producing AC and degrading PDE (Conti et al., 2014; Lohse et al., 2017; Bock et al., 2020; Zhang et al., 2020). Given that AKAP79/150 exists in nanoclusters at the plasma membrane in multiple cell types (Mo et al., 2017; Zhang et al., 2016) and associates with AC8 in β cells (Willoughby et al., 2010), we hypothesized that AC8 could form nanoclusters on the plasma membrane of MIN6 cells and compartmentalize cAMP dynamics. To test this hypothesis, we examined the spatial organization of AC8 and AKAP150 at the membrane using Stochastic Optical Reconstruction Microscopy (STORM). We found that AKAP150 molecules were organized in clusters with a mean radius of 127 ± 9 nm and an average nearest-neighbor spacing of 313 ± 20 nm between cluster centers (n = 20) (Figure 3a and Figure 3—figure supplement 1), consistent with several recent reports demonstrating AKAP79/150’s tendency to form nanoclusters in other cell types (Mo et al., 2017; Purkey et al., 2018; Tajada et al., 2017; Zhang et al., 2016). Thus, the AKAP79/150 compartment-specific cAMP phase is likely representative of the balanced cAMP generation and degradation within these AKAP clusters.

Figure 3 with 3 supplements see all

Download asset Open asset

AKAP150 and AC8 both form nanoclusters at the surface of MIN6 β cells.

(a) Representative super-resolution STORM image of the AKAP150 scaffold (scale 5 μm, inset 500 nm). Ripley-K analysis measures the average radii of the nanoclusters and indicates that AKAP150 forms clusters of 127 ± 9 nm, n = 20 (uniform random distribution control in inset). The nearest-neighbor distance distribution describes the distance between nanoclusters (average distance for AKAP150 is 313 ± 20 nm). (b) Representative super-resolution STORM image of Ca²⁺-sensitive AC8 (scale 5 μm, inset 500 nm) depicts AC8 nanoclusters of average radius 88 ± 8 nm and average nearest-neighbor distance 292 ± 16 nm, n = 16.

Due to the known interaction between AKAP79/150 and AC8 (Willoughby et al., 2010) and diffusion of membrane-localized signaling complexes, we next probed the spatial organization and mobility of AC8. We found AC8 also distributes non-uniformly at the plasma membrane and forms clusters with a mean radius of 88 ± 8 nm and an average nearest-neighbor spacing of 292 ± 16 nm between cluster centers (n = 16) (Figure 3b; Figure 3—figure supplement 1). In order to measure the degree of colocalization between AKAP150 and AC8 nanoclusters, we also performed 2-color STORM and observed 72% of AKAP150 localizations were co-clustered with AC8 localizations, suggesting the presence of AKAP150:AC8 co-clustered nanodomains (Figure 3—figure supplement 2a–e). Next, to determine the mobility of AC8 molecules within the timescale of the oscillatory circuit’s period, we performed Fluorescence Recovery After Photobleaching (FRAP) with EGFP-tagged AC8 and found AC8 diffuses slowly at the membrane, and that there exists a significant immobile fraction (D_AC8 = 0.019 ± 0.002 μm²/s, avg. % immobile = 42.3%, n = 16) (Figure 3—figure supplement 3).

With the evidence of the nanoscale organization of AKAP150 and AC8 on the plasma membrane, we further hypothesized that the increased spatial density of Ca²⁺-driven cAMP sources within the AKAP150 clusters, in conjunction with dispersed PDE1 in the cytosol (Bender and Beavo, 2006; Goraya et al., 2008), is important in compartmentalizing cAMP production and mediating the in-phase cAMP signal. To test this idea, we sought to build a mathematical framework to describe the spatial compartmentalization of the in- and out-of-phase cAMP-Ca²⁺ oscillations. Briefly, we expanded our network motif model (Figure 2c) by including a 3D spatial component with cAMP diffusion (D_cAMP = 60 μm²/s; Agarwal et al., 2016) and incorporating our previously published well-mixed β cell model (Ni et al., 2011). We used the AKAP79/150:AC8 cluster pattern measurements from the STORM imaging to set model parameters in a hexagonal prism domain (200 nm edge, 600 nm depth), with one immobile AKAP79/150:AC8 cluster centered in the domain for simulation (Figure 4a, see Appendix 1 for model development details). By localizing AC8 within the AKAP79/150:AC8 cluster on the plasma membrane face and leaving PDE1 well-mixed throughout the volume, we could simulate Ca²⁺-driven cAMP oscillations that were in-phase within the immediate vicinity of a cluster, but sharply transitioned out-of-phase outside the cluster. Specifically, during a Ca²⁺ influx event, Ca²⁺-triggered cAMP production dominated at the center of the AKAP79/150 cluster while Ca²⁺-triggered cAMP degradation was favored outside the cluster at the PM and in the center of the unit volume (Figure 4a). Not surprisingly, the regime that recapitulates this phase relationship is sensitive to the spatially-restricted AC8/PDE1 balance and the diffusivity of cAMP. For example, by increasing the AC8 concentration within the cluster, the Ca²⁺-driven cAMP generative flux could exceed the rate of cAMP degradation by PDE1 and thus cAMP and Ca²⁺ oscillations could exist in-phase within both compartments (Figure 4—figure supplement 1). Alternatively, assuming that AC8 clustering is driven by AKAP150:AC8 interactions, weakening this interaction would then reduce the AC8 cluster stabilization and lead to a redistribution of AC8 away from the nanoclusters and a decrease in the local concentration of AC8 within the clusters (Figure 4b). Without the high local concentration of AC8 driving a net positive cAMP production within an AKAP79/150 cluster, the spatial domain where cAMP oscillates in-phase with Ca²⁺ is predicted to shrink while the out-of-phase regime expands and can reverse the phase at the cluster center (Figure 4c).

Figure 4 with 1 supplement see all

Download asset Open asset

cAMP-Ca²⁺ phase relationship can be described by a 3D reaction-diffusion model involving clusters of AKAP79/150 and AC8.

(a) 3D reaction-diffusion model with a single AKAP79/150:AC8 co-cluster positioned at the PM in the β cell in a hexagonal prism volume. cAMP oscillates in-phase immediately within the AKAP79/150:AC8 nanocluster due to the high effective concentration of AC8, but out-of-phase at the PM or cytosol due to the presence of PDE1 (cAMP – red, Ca²⁺– blue). (b) Disruption of the AKAP79/150:AC8 interaction can redistribute AC8 from within the cluster to the PM, shown by the half-Gaussian cross-sections and representative AC8 concentration heatmaps at the PM. (c) Heatmap depicting the time lag (s) for AC8 distribution (% Gaussian) and spatial distance (nm) from cluster center along PM.

To test this prediction, we overexpressed the amino terminus of AC8 (AC8^1-106) required for interaction with AKAP79/150 (Willoughby et al., 2010) to compete with the binding of endogenous AC8 to the endogenous AKAP150 scaffold. The disruption of the AKAP150:AC8 interaction was validated using a proximity ligation assay (PLA) to visualize the AKAP150:AC8 interaction in situ. Compared to non-transfected cells, cells expressing the AC8^1-106 peptide had a 39 ± 4% reduction in the PLA signals, indicating a decrease in the AKAP150:AC8 interaction (Figure 5—figure supplement 1a,b). Furthermore, STORM imaging showed that overexpression of the AC8^1-106 peptide led to a decrease in the percentage of AC8 single-molecule localizations within AC8 nanoclusters (n = 9) (Figure 5a), consistent with the predicted redistribution of AC8 molecules (Figure 4b). To test the impact of loss of AC8 molecules from the nanoclusters on the oscillation phase, we measured AKAP79/150-localized cAMP in the presence of AC8^1-106 and observed a significant increase in the average lag time (43 ± 6 s, n = 33) (Figure 5b). This is due to a higher proportion of cells exhibiting out-of-phase cAMP oscillations, indicating that the AKAP79/150:AC8 competitor peptide was sufficient in reversing the phase relationship in the AKAP79/150 compartment. Interestingly, we also observed many cells displaying irregular Ca²⁺ oscillations, as indicated by a disruption in the periodic timing of individual cells’ Ca²⁺ peaks (Figure 5b, left). This nanoscale perturbation establishes the regulatory role of the AKAP79/150:AC8 interaction in mediating the compartmentalized cAMP-Ca²⁺ phase relationship.

Figure 5 with 1 supplement see all

Download asset Open asset

Disruption of the AKAP79/150:AC8 interaction is associated with a redistribution of AC8 at the PM and a phase shift of cAMP at the AKAP79/150 nanodomain.

(a) Over-expression of the N-terminus of AC8, which is necessary and sufficient for mediating the AKAP79/150:AC8 interaction, redistributes AC8 from within nanoclusters to the general PM, as seen in the STORM image (scale 5 μm, inset 500 nm) and measured by the percent of localizations that fall into nanoclusters. (b) Disruption of the AKAP79/150:AC8 interaction lengthens the time lag between the cAMP (red) and Ca²⁺ (black) signals at the AKAP79/150 compartment (avg. time lag in absence of disruptor is 13 ± 3 s, n = 60, and in presence of disruptor 43 ± 6 s, n = 33) due to more cells displaying out-of-phase cAMP oscillations. *p=0.0341, ****p<0.0001; unpaired two-tailed Students t-test.

AKAP79/150-mediated phase relationship is critical for regulating oscillatory Ca²⁺

Next we systematically examined the impact of perturbing the precisely regulated phase relationship within the AKPA79/150 compartment. Due to the modulatory role of PKA in the Ca²⁺-cAMP-PKA oscillatory circuit and the interaction between PKA and AKAP79/150, we wondered how the in-phase cAMP oscillations with respect to Ca²⁺ are translated into PKA activities and if spatial compartmentalization of the phase relationship is also maintained at the PKA activity level. Therefore, we extended our 3D model to include AKAP79/150-associated PKA (see Appendix 1 for model details). In this extended model, PKA activity oscillations exhibit distinct phase relationships with respect to Ca²⁺ within and outside of the AKAP79/150 compartment (Figure 6—figure supplement 1a). To test this prediction, we fused our FRET-based biosensor for PKA activity (AKAR4) (Depry et al., 2011) to either full-length AKAP79 or the PM-targeting motif and expressed the sensors in MIN6 cells. Upon TEA stimulation, PKA activity was observed to oscillate with a lag time of 25 ± 6 s (n = 15) within the AKAP79/150 compartment but with a lag time of 55 ± 8 s (n = 12) (Figure 6—figure supplement 1b–e) at the general plasma membrane, indicating that the compartmentalized phase relationship is preserved from cAMP to PKA.

Spatiotemporal organization of PKA signaling and its phosphorylation targets via AKAPs have been implicated in regulating several important pathways. For example, PKA has been shown to phosphorylate Ca_V1.2 in an AKAP79/150-dependent manner, and this modification can influence the open probability of the channel (Murphy et al., 2014), suggesting a mechanistic link between local cAMP/PKA activity and global oscillatory Ca²⁺. Thus, we sought to study the functional role of the spatially compartmentalized cAMP-Ca²⁺ phase relationship in regulating intracellular Ca²⁺ dynamics. We measured Ca²⁺ oscillations by RCaMP in the presence of either the EGFP-tagged AKAP79/150:AC8 disruptor peptide, AC8^1-106, or EGFP alone as a control. Population-wide differences in Ca²⁺ dynamics, such as strength and timing, were observed in AC8^1-106-transfected cells and visualized in heat maps depicting the normalized Ca²⁺ signal per cell versus time (Figure 6a). Interestingly, we found that the expression of the disruptor peptide was correlated with a significant decrease in the peak ratio between the second Ca²⁺ peak and the first Ca²⁺ peak (control average −1.6%, n = 270; AC8^1-108 average −10.8%, n = 562) post TEA addition, indicative of less sustained oscillations (Figure 6b,c). In addition to intracellular Ca²⁺ concentration, the precise timing of internal oscillatory events is critical for modulating β cell functions such as glucose homeostasis and pulsatile insulin secretion (Fridlyand et al., 2010). In the presence of the disruptor peptide, cells also exhibited a longer elapsed time between oscillatory Ca²⁺ peaks (control average 3.9 ± 0.1 min, n = 270; AC8^1-108 average 4.6 ± 0.1 min, n = 562), suggesting that the timing of the signaling circuit was disturbed (Figure 6b,c). It is well established that besides the precise timing, the regularity of cytoplasmic Ca²⁺ in β cells is crucial in mediating pulsatile insulin secretion from the pancreas (Gilon et al., 2002; Schmitz et al., 2002). By stratifying the disruptor peptide-expressing cell population into ‘low, ‘medium,’ and ‘high’ expressers and performing a blinded classification of responding cells based on the regularity of the Ca²⁺ oscillation (see Appendix 1 for details), we found a positive correlation between the percentage of cells exhibiting irregular oscillations and the expression level of the disruptor peptide (42% for low-expressing vs. 68% for high-expressing AC8^1-106 disruptor) (Figure 6d). Taken together, these data signify that the compartmentalized cAMP-Ca²⁺ phase relationship regulates the oscillatory Ca²⁺ signal and plays an important role in determining the pace, regularity, and sustainability of the Ca²⁺ oscillations.

Figure 6 with 2 supplements see all

Download asset Open asset

Ca²⁺ oscillatory dynamics are affected by expression of the disruptor peptide in β cells.

(a) Heatmaps depicting Ca²⁺ oscillations for 220 randomly selected cells co-expressing EGFP alone (control, left) or EGFP-tagged AC8^1-106 (AKAP79/150:AC8 disruptor, right), ordered by a mixed parameter describing the time lag between the first two Ca²⁺ peaks and the avg. timelag between all Ca²⁺ peaks. (b) Schematic describing two Ca²⁺ oscillatory parameters: the ratio between the first two Ca²⁺ peaks (R = P₂ / P₁) and the interpeak timing (t). (c) The peak ratio R is decreased in the presence of the AKAP79/150:AC8 disruptor (left), indicating less of a sustained Ca²⁺ oscillatory response. Over-expression of the disruptor also lengthens the timing between peaks (right). (d) The level of disruptor peptide expression is correlated with an increase in the percentage of cells exhibiting irregular oscillations (total n = 562). p=0.0054, ordinary one-way ANOVA.

Discussion

Biological oscillations represent a rich way of encoding information using the amplitude and frequency. Here, we show the phase in an oscillatory signaling circuit represents a novel mode of informational encoding for which the phase itself can be spatiotemporally regulated. In the case of the Ca²⁺-cAMP-PKA circuit, the oscillatory cAMP/PKA phase relative to a widespread Ca²⁺ signal is distinctly regulated within two adjacent membrane compartments through intracellular organization of scaffolds and signaling effectors. Localized perturbation of this spatial phase signature disrupts global Ca²⁺ oscillations and thus has far-reaching consequences on the functional landscape of the β cell.

Compartmentalization of cAMP/PKA signaling is instrumental in processing a diverse set of inputs and mediating specific cellular functions; however, the mechanistic details of compartmentalization await to be fully elucidated (Musheshe et al., 2018). Given the measured kinetic rates of most ACs and PDEs, the generation of local cAMP gradients around single enzymes appears unfeasible (Conti et al., 2014). Context-dependent discrepancies in some of the kinetics (i.e. differences between in vitro and in vivo measurements) and slower cAMP diffusion due to buffering have been shown to be important for cAMP compartmentalization (Bock et al., 2020; Zhang et al., 2020). Here, we propose the nanoscale organization of key cAMP effectors and regulators as a novel mechanism for cAMP compartmentation. Despite the slow rates measured for individual ACs, we computationally and experimentally describe conditions in which the generation of compartmentalized cAMP can emerge from the clustering of many relatively immobile AC8 enzymes at the membrane and bulk distribution of PDE1 in the cytoplasm. Alterations to this nanoscale organization lead to dysregulated Ca²⁺ oscillations, demonstrating the functional importance of maintaining this organization. This system also serves as a nanoscale demonstration of how a cell can translate a global signal (Ca²⁺) into a compartmentalized signal (cAMP/PKA activity) by local activation and global inhibition, a strategy that is likely utilized in many other cellular contexts (Levchenko et al., 2000; Purvis and Lahav, 2013).

In order to understand the functional impact of the nanoscale architecture, it is important to also consider the discrete scale of such signaling organization. According to our rough estimates, approximately 15–150 individual AC8 molecules comprise AKAP150:AC8 nanoclusters (see Materials and methods for STORM-based estimation). To investigate the role of small molecule numbers, we recast our previous network motif model as a stochastic system and simulated the impact of changing initial conditions of PDE1 and AC8. We found that even in these discrete simulations, the previously noted phase relationships are preserved. As the balance between PDE and AC is perturbed, the system moves from an in-phase to an out-of-phase cAMP response with respect to Ca²⁺. Interestingly, we also noted that increasing the number of AC8 molecules leads to a reduction of the noise in the cAMP response. Thus, we predict that tuning AC8 cluster size might be a biological design strategy of nanoclusters to maintain accurate responses to perturbations (Figure 6—figure supplement 2a–d).

The interplay between Ca²⁺ and cAMP is very intricate, and multiple mechanisms could contribute to the observed dysregulated signaling effects upon perturbation of the AKAP79/150-compartmentalized cAMP-Ca²⁺ phase (Figure 6a–d). AKAPs can recruit PKA to regulate channel activities (Dell'Acqua et al., 2006; Mo et al., 2017; Torres-Quesada et al., 2017), such as in the regulation of voltage-mediated Ca²⁺ entry via PKA-dependent phosphorylation of Ca_V1.2 (Murphy et al., 2014) or the modulation of store-operated Ca²⁺ entry by both PKA-dependent STIM1 and Orai1 phosphorylation (Thompson and Shuttleworth, 2015; Zhang et al., 2019). Additional levels of regulatory feedback within the Ca²⁺-cAMP-PKA oscillatory circuit have also been identified, such as a negative feedback loop involving PKA phosphorylation of AC8, thereby fine-tuning the circuit dynamics (Willoughby et al., 2012).

The precisely regulated Ca²⁺ and cAMP oscillations are likely functionally important as insulin secretion requires coordination between Ca²⁺, cAMP, and PKA. Localized cAMP/PKA signaling at the AKAP79/150 scaffold and the in-phase oscillations of cAMP and Ca²⁺ within this compartment might play a critical role in regulating insulin secretion due to close interactions between AKAP79/150 and insulin secretory granules via Ca_V1.2 (Barg et al., 2001). Several important processes and components of the secretory machinery have been identified as targets of PKA signaling here, such as PKA-dependent mobilization of granules (Renström et al., 1997) and modulation of the synaptosomal protein SNAP25 (Gao et al., 2016). These PKA phosphorylation events may be highly coordinated with oscillations of Ca²⁺, which drive the exocytosis process for pulsatile insulin secretion. In addition to PKA-dependent secretory control, cAMP has recently been implicated to play a direct role in fusion pore formation via the cAMP-regulated guanine exchange factor Epac (Guček et al., 2019). Compartmentalized phase regulation of the Ca²⁺-cAMP-PKA oscillatory circuit at the AKAP79/150 macromolecular complex is likely involved in the regulation of many β cell processes, in addition to insulin secretion, and more work will be needed to further establish the mechanisms involved in decoding the information embedded in the local phase relationship.

The Ca²⁺-cAMP-PKA oscillatory circuit in pancreatic β cells integrates many important regulators of cellular function, and the precise coordination of each is required for proper signaling control. Here, we have uncovered a spatiotemporal organization of the circuit where the oscillatory phase between cAMP/PKA and Ca²⁺ depends on the spatial proximity of the AKAP79/150 scaffold protein and AC8. The construction principles of this signaling nanodomain, including the spatial distributions of sinks and sources, likely represent a generalized strategy for the generation of other compartmentalized signals and provide a unique modality in which cells embed, process, and produce signaling information.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Cell line (Mus musculus)	MIN6	Dr. Jun-Ichi Miyazaki, Osaka University	RRID:CVCL_0431
Chemical compound, drug	Tetraethylammonium chloride (TEA)	Sigma	T2265	20 mM
Chemical compound, drug	8-Methoxymethyl-3-isobutyl- 1-methylzanthine (8MM-IBMX)	Sigma-Aldrich	M2547	100 μM
Chemical compound, drug	Milrinone	Alexis	Cat# ALX-270–083 M005	10 μM
Chemical compound, drug	Rolipram	Alexis	Cat# ALX-270–119	1 μM
Chemical compound, Drug	KCl	Sigma-Aldrich	P9541	15 mM
Commercial assay or kit	Lipofectamine-2000	Invitrogen	11668019
Recombinant DNA reagent	pCDNA3 AKAR4	PMID:20838685
Recombinant DNA reagent	(Ci/Ce)Epac2-camps	Dr. Dermot Cooper, University of Cambridge
Recombinant DNA reagent	AKAP79 (AKAP5, Homo sapiens)	Dr. John D. Scott, University of Washington
Recombinant DNA reagent	AC8 (Adcy8, Rattus norvegicus)	Dr. Dermot Cooper, University of Cambridge
Recombinant DNA reagent	ShAdcy8 Plasmid	Dr. Jochen Lang, University of Bordeaux	ShAdcy8 #2
Recombinant DNA reagent	RCaMP	Dr. Loren Looger, Janelia Farms
Recombinant DNA reagent	R-FlincA	Dr. Kazuki Horikawa, Tokushima University Graduate School
Commercial assay or kit	Duolink in situ red starter kit	Sigma	DUO92101
Commercial assay or kit	Duolink in situ Probemaker MINUS	Sigma	DUO92010
Commercial assay or kit	Duolink in situ Probemaker PLUS	Sigma	DUO92009
Antibody	Anti-AC8, rat polyclonal	Abcam	ab196686	1:2000
Antibody	Anti-AKAP150, rat polyclonal	Millipore Sigma	07–210	1:500
Antibody	Anti-AKAP79, mouse monoclonal	BD	610314	1:500
Antibody	Goat anti-rabbit AlexaFluor647	ThermoFisher Scientific	A21245	1:1000
Antibody	Goat anti-mouse AlexaFluor568	ThermoFisher Scientific	A11031	1:1000
Software, algorithm	FIJI	https://imagej.net/Fiji	RRID:SCR_014294
Software, algorithm	MetaFluor	https://www.moleculardevices.com/	RRID:SCR_002285
Software, algorithm	MATLAB	https://www.mathworks.com/	RRID:SCR_001622
Software, algorithm	GraphPad Prism	https://www.graphpad.com/ scientific-software/prism/	RRID:SCR_002798
Software, algorithm	COPASI	http://copasi.org/ Hoops et al., 2006	RRID:SCR_014260
Software, algorithm	Virtual Cell	https://vcell.org/ Cowan et al., 2012	RRID:SCR_007421

Gene construction

Request a detailed protocol

For AKAP79-(Ci/Ce)Epac2-camps, AKAP79 (gene AKAP5, from Dr. John D. Scott) was PCR amplified to have HindIII/BamHI digestion sites, and (Ci/Ce)Epac2-camps (from Dr. Dermot Cooper) was PCR amplified to have BamHI/EcoRI digestions sites. Both fragments were inserted into a pcDNA3 (Invitrogen) backbone for mammalian expression (cAMP sensor is C-terminal to AKAP79). For AKAP79-AKAR4, a similar approach was taken where AKAR4 was dropped between BamHI/EcoRI. For AC8 (gene Adcy8, from Dr. D. Cooper), AC8^1-108, Gibson Assembly was used to insert the genes into the pcDNA3 mammalian expression vector. The ShAdcy8 construct for Adcy8 knockdown was previously verified and a gift from Dr. Jochen Lang. RCaMP was a gift from Dr. Loren Looger.

Cell culture

Request a detailed protocol

MIN6 cells (a mouse insulinoma β cell line, gift from Dr. Jun-Ichi Miyazaki) were plated onto sterilized glass coverslips in 35 mm dishes and grown to 50–90% confluency in DMEM (10% FBS, 4.5 g/L glucose) at 37°C with 5% CO₂. Identity was verified by phenotype (Ca²⁺ oscillations, morphology, and insulin secretion). Recurrent mycoplasma testing was performed and results were negative. Cells were transfected using Lipofectamine 2000 (Invitrogen) for 20–48 hr before imaging.

Imaging

Request a detailed protocol

Cells were washed twice with Hank’s balanced salt solution buffer and maintained in the dark at room temperature. Cells were imaged on a Zeiss Axiovert 200M microscope with a cooled charge-coupled device camera (MicroMAX BFT512, Roper Scientific, Trenton, NJ) controlled by METAFLUOR 6.2 software (Universal Imaging, Downingtown, PA). Dual cyan/yellow emission from FRET used a 420DF20 excitation filter, a 450DRLP dichroic mirror, and two emission filters [475DF40 for CFP and 535DF25 for YFP]. RFP fluorescence was imaged using a 568D55 excitation filter, a 600DRLP dichroic mirror, and a 653DF95 emission filter. GFP fluorescence was imaged using a 480DF30 excitation filter, a 505DRLP dichroic mirror, and a 535DF45 emission filter. YFP fluorescence was imaged using a 495DF40 excitation filter, a 515DRLP dichroic mirror, and a 535DF25 emission filter. These filters were alternated by a filter-changer Lambda 10–2 (Sutter Instruments, Novato, CA). Exposure time was 50–500 ms, and images were taken every 10–30 s.

For confocal imaging, images were collected with a C2 plus on a Nikon Ti2 inverted microscope equipped with a Plan Apo lambda 60x oil immersion objective NA 1.4. YFP fluorescence was excited with the 488 nm line from a LU-N4 laser. Images were acquired with a DUVB detector collecting emission from 495 nm to 600 nm with a virtual spectral GaAsP detector controlled by NIS Elements software. The pinhole was set at 30 μm. Frame size was 1024 × 1024 pix.

For FRAP acquisition, cells were imaged using a C2 plus mounted on a Nikon Ti2 with a 488 nm laser (LU-N4 laser engine, Nikon Instruments) and sensitivity PMT detector (C2-DU3, Nikon) with 525/50 nm bandpass filter. The excitation light is reflected by a dichroic mirror 405/488/561/640 and focused through an Apo 100 × 1.49 NA objective. Image stacks were acquired in Galvano mode with unidirectional scanning with a 488 nm laser at 1.5% laser power with frame size 512 × 512 at scan zoom 2, one frame per second and 97.1 µm pinhole size. Small stimulation ROIs were drawn at the plasma membrane. The total FRAP series contained three images before bleaching (obtained with 2 s intervals), two cycles of ROI bleaching with the 488 nm laser at 100% laser power (5 frames at one frame per second), and two minutes of continuous acquisition to monitor fluorescence recovery.

Image analysis

Request a detailed protocol

Fluorescence images were background-corrected by subtracting the fluorescence intensity of background with no cells from the emission intensities of cells expressing fluorescent reporters. Cells for experiments using the targeted biosensors were manually segmented from microscopy images based on emission intensity (YFP) and RCaMP (RFP) using custom-written Java and MATLAB scripts and whole-cell fluorescence emission was measured. The ratio of yellow/cyan emission (AKAR4) or cyan/yellow emission (Epac2-camps) and RFP intensity were then calculated at different time points using MATLAB scripts. The values of all-time courses were normalized by dividing each by the average basal value before drug addition. FRAP measurements were background-subtracted, bleach-corrected, and normalized, and fit to a previously described diffusion-dominant model (Ellenberg et al., 1997; Phair et al., 2004). For the AKAP150:AC8 disruption experiments, cells were segmented with custom CellProfiler (Broad Institute) pipelines based on RCaMP (RFP) fluorescence emission, and both RFP and GFP fluorescence intensities (EGFP-AC8^1-106) per cell were measured every 30 s. Graphpad Prism eight was used for visualization and statistical analysis. Custom-written MATLAB scripts were used for correlation analysis and regularity classification.

Super-resolution imaging (STORM)

Request a detailed protocol

For fixed-cell stochastic optical reconstruction microscopy (STORM) imaging, cells were fixed with 4% paraformaldehyde (PFA) and 0.2% glutaraldehyde (GA) for 20 min and then washed with 100 mM glycine in HBSS to quench the free PFA. Cells were permeabilized and blocked in a permeabilization solution with 0.1% Triton X-100, 0.2% bovine serum albumin, 5% goat serum, and 0.01% sodium azide in HBSS. The cells were then incubated overnight at 4°C with an anti-AC8 antibody (Abcam, ab196686) at a 1:2000 dilution or an anti-AKAP150 (Millipore Sigma, 07–210) antibody at a 1:500 dilution, followed by 1 to 2 hr with goat anti-rabbit Alexa 647–conjugated antibody (ThermoFisher Scientific, A21245) at 1:1000 dilution. For the 2-color STORM experiments, cells were incubated with both the anti-AC8 antibody (1:2000 dilution) and an anti-AKAP79 antibody that recognizes a conserved epitope in AKAP150 (BD, 610314; raised against the C-terminal 248 residues of AKAP79 which shares 70% similarity with a region in AKAP150; 1:500 dilution), followed by a 2 hr incubation with a goat anti-rabbit Alexa 647-conjugated antibody (1:1000 dilution) and a goat anti-mouse Alexa 568-conjugated antibody (ThermoFisher Scientific, A11031; 1:1000 dilution). The cells were then post-fixed again in 4% PFA and 0.2% GA, quenched with 100 mM glycine in HBSS, and washed with HBSS to prepare for imaging. Immediately before imaging, the medium was changed to STORM-compatible buffer (50 mM Tris-HCl, pH 8.0, 10 mM NaCl, and 10% glucose) with glucose oxidase (560 μg/ml), catalase (170 μg/ml), and mercapto-ethylamide (7.7 mg/ml). STORM images were obtained using a Nikon Ti total internal reflection fluorescence (TIRF) microscope with N-STORM, an Andor IXON3 Ultra DU897 EMCCD, and a 100 × oil immersion TIRF objective. Photoactivation was driven by a Coherent 405 nm laser, while excitation was driven with a Coherent 647 nm laser or Coherent 561 nm laser.

All image analysis and image reconstruction were performed using both Nikon Elements analysis software and custom-written MATLAB scripts. Blinking correction was performed by implementing the pairwise Distance Distribution Correction (DDC) algorithm (Bohrer et al., 2019). Cluster property measurements were performed using Ripley-K analysis and custom mean-shift MATLAB code for segmentation, as described before (Mo et al., 2017). Co-clustering measurements were made using custom MATLAB scripts implementing Getis-Franklin pattern analysis, as described previously (Getis and Franklin, 1987; Mo et al., 2017). Briefly, we first calculated the degree of local clustering of AKAP150 molecules using the L(r) spatial statistic at r = 200 nm which counts the number of AKAP150 localizations within 200 nm (length scale of nanoclusters) of each AKAP150 localization, normalized appropriately. We then calculated Lcross(r) at r = 200 nm which counts the number of AC8 localizations within 200 nm of each AKAP150 localization, normalized appropriately. Plots of Lcross(200) vs. L(200) were generated and an L threshold of 150 was chosen based on the proportion of localizations found within clusters from the mean-shift segmentation analysis, and used to divide the plot into quadrants to calculate the proportion of AKAP150 localizations co-clustered with AC8 as previously described (Mo et al., 2017).

We obtained a very rough estimate of AC8 molecules within nanoclusters from our blink-corrected STORM datasets by assuming background localizations outside the cell are due to single secondary antibodies with conjugated dyes. We measured an average of 2 localizations per secondary antibody outside the cell and approximately 50 localizations per AC8 nanocluster within the cell. Estimating a labeling stoichiometry of 0.4–1.5 secondary antibodies per primary antibody and 0.4–1.5 primary antibodies per AC8 molecule (Zanacchi et al., 2017; Nieuwenhuizen et al., 2015), we can roughly approximate the number of AC8 molecules per nanocluster to be between 15 and 150 molecules/cluster (N = 3 cells).

Proximity ligation assay

Request a detailed protocol

Antibodies for AC8 and AKAP150, mentioned in STORM section, were buffer exchanged into DPBS and conjugated with MINUS or PLUS oligos, following the Sigma DuoLink in situ Probemaker kits. PLA experiments were performed using the Duolink in situ red kit for PLAs according to the provided protocol. The only protocol modification was to extend the amplification time by 50 min. Briefly, cells were fixed and permeabilized as in the STORM experiments before incubation with PLUS and MINUS oligo-conjugated primary antibodies for 30 min at 37°C each, with washes after each step. Ligation of the nucleotides and amplification of the strand occurred sequentially by incubating cells with first ligase then polymerase and detection solution. PLA experiments with AKAP95 antibodies from different species were used as positive controls in HEK293T cells, and experiments with just one oligo-labeled primary antibody or the other were our negative control. Images were acquired on a Nikon Ti Eclipse epifluorescence scope with z-control and maximum intensity projections were created. A cross section of the DAPI-stained nucleus (3.6–5 μm) was also acquired and the number of dots per cell was counted using the nucleus as reference.

Statistical analysis

Request a detailed protocol

Statistics, such as means, SEMs, and comparisons were calculated using Graphpad Prism 8 software. Experiments comparing multiple conditions were analyzed using one-way ANOVA or unpaired t-tests. Statistical significance was assessed using p<0.05 as a cutoff. Experiments were repeated >3 times. Individual cells with high expression of the AKAP79/150-fused biosensors were removed from analysis (see Appendix 1 for details).

Appendix 1

cAMP analysis and quantification

AKAP79-(Ci/Ce)Epac2-camps transfected cells displayed cAMP oscillations that were either in-phase or out-of-phase with their respective Ca²⁺ signal. This was in sharp contrast to responsive lyn-(Ci/Ce)Epac2-camps transfected cells where all cells yielded only out-of-phase oscillations. We found a strong correlation between the AKAP79-fused sensor expression level and the observed cAMP-Ca²⁺ phase relationship, with cells having lower levels of sensor present displaying predominantly in-phase cAMP oscillations and cells with higher levels of the AKAP79/150-fused biosensor exhibiting out-of-phase oscillations (Appendix 1—figure 1). Overexpression of the AKAP79 scaffold likely changed the stoichiometry of the signaling complexes and resulted in unsuccessful targeting of the biosensor to functional AKAP79/150 domains, and so in this manuscript we considered only TEA-responsive cells below an AKAP79 expression threshold determined by the YFP fluorescence (Appendix 1—figure 1).

Time lag calculation

Due to the heterogeneity of cellular Ca²⁺ and cAMP/PKA activity oscillatory responses in each cell (i.e. variations in frequency, amplitude, and regularity), we sought an applicable metric to describe the phase relationship. Here we measure the lag time (sec) between the Ca²⁺ signal trace and the cAMP/PKA activity signal trace. Specifically, we high-pass filtered both the Ca²⁺ and cAMP/PKA activity traces (approximately 20 min) to subtract out slowly varying baseline changes, normalized the traces so that the maximum intensity/FRET ratio was set to 1, and then computed the cross-correlation to measure the signal overlap for different lag times. To calculate the lag time, we identified peaks in the cross-correlation passing a peak prominence cutoff and found the absolute value of the shortest lag time corresponding to a peak maximum. For in-phase oscillations, the lag time was typically small (τ ≤20 s) due to the two signal traces oscillating in synchrony. On the other hand, out-of-phase oscillations typically corresponded to longer lag times (τ >20 s) due to the anti-phasic relationship seen in the peak timing and peak shape. Analysis was performed with custom scripts in MATLAB and Java, and pipelines in CellProfiler.

Quantification of nanodomain perturbation effects on global Ca²⁺

In order to measure the effects of AKAP79/150:AC8 disruption on Ca²⁺ dynamics, we transiently transfected and expressed AC8^1-106 in MIN6 cells and measured Ca²⁺ with RCaMP. For quantification of the interpeak timing and peak ratio, we first selected cells that responded to the TEA treatment, identified Ca²⁺ peaks passing a peak prominence cutoff, and finally calculated the avg. time between peak maxima and RFP intensity ratio between the second and first Ca²⁺ peak maxima. To find the percentage of cells with regular vs. irregular Ca²⁺ oscillations, we randomized all Ca²⁺ traces from the experimental and control samples and performed a blinded classification to sort the single cell traces as regular, irregular, or nonresponsive. Analysis was performed with custom scripts in MATLAB and Java, and pipelines in CellProfiler.

Computational modeling

Well-mixed system

Assumptions

Signaling components are present in large enough quantities that concentration changes are smooth and move in a deterministic fashion.
Well-mixed kinetic rate constants contain conversion factors between compartments, that is membrane to cytosol.
Binding interactions occur rapidly enough such that any kinetic parameter, k, remains constant on surfaces (Berry, 2002).
AKAP does not alter the activity of the catalytic subunits of Protein kinase A (PKA) instead it only affects localization.
Ca²⁺ independent activity of Adenylyl cyclase (AC) is of the same strength as inactive Ca²⁺ dependent Adenylyl cyclase (AC8) and stays at a constant value.

Modeling well-mixed chemical reactions

Mass-action kinetics

We generated an ordinary differential equation (ODE) for every species using mass-action kinetics for each binding interaction. The law of mass-action states that the rate of a chemical reaction is proportional to the product of the concentration of the reactants raised to the power of their stoichiometric coefficient (Guldberg and Waage, 1879). Mass-action kinetics rely on the assumption that the rate constant, k, is constant over time. For example, consider the one-reaction system:

A + 2 B ⇌ 3 C,

where the forward and backward rates are $k_{1}$ and $k_{2}$ . The differential equations describing the dynamics of species $A, B$ , and $C$ under mass-action kinetics are:

\frac{d [A]}{d t} = k_{2} [C]^{3} - k_{1} [A] [B]^{2} \frac{d [B]}{d t} = 2 k_{2} [C]^{3} - 2 k_{1} [A] [B]^{2} \frac{d [C]}{d t} = 3 k_{1} [A] [B]^{2} - 3 k_{2} [C]^{3} .

Michaelis-Menten kinetics

We used Michaelis-Menten kinetics to model the kinetics of enzyme-catalyzed reactions. When a reaction is catalyzed by an enzyme with kinetic properties $k_{c a t}$ and $K_{M}$ ,

S E \to P

then the reaction rate is given by

\frac{d [S]}{d t} = - \frac{k_{c a t} [E] [S]}{K_{M} + [S]} = \frac{d [P]}{d t} .

For Michaelis-Menten kinetics, concentrations of reactants and products must be in large enough quantities, and one of the following conditions must apply: the concentration of the substrate is much larger than the concentration of products, $[S] ≫ [p]$ , and/or the energy released in the reaction is very large, $Δ G ≪ 0$ .

Model development

We constructed a biochemical network to represent interactions between Ca²⁺ and cAMP in $β$ -cells, Figure 1b, (Fridlyand and Philipson, 2016; Ni et al., 2011) after a depolarization event. The computational model considered the dynamics of calcium, potassium, leaky, and calcium-sensitive potassium channels (Appendix 1—table 1). Importantly, we included feedback of PKA with KATP channels and the inclusion of Ca²⁺-sensitive ACs and PDEs. The model contains a total of 92 parameters with 11 free parameters. The values of the parameters are constrained through both previously peer-reviewed publication results (Ang and Antoni, 2002; Boras et al., 2014; Lai et al., 2015; Masada et al., 2009; Ni et al., 2011) and with new experimental results using obtained FRET measurements. To constrain source and sink activation rates, previously published literature of AC and PDE stimulus-response curves (of related isoforms) was utilized (Ang and Antoni, 2002; Masada et al., 2009). COPASI was used to calculate initial guesses for kinetic activation parameters. FRET measurements took precedence over binding curves, especially for the spatial model where further fitting routines were performed to refine the model. Model variations were performed to attain both semi-physiological concentrations (within the ranges of the sensor) and phase (period and relation) information. Predictions are made on qualitative behavior as opposed to quantitative as proper parameter fitting would require much more data than available for this system.

Values and reaction sets used in the well-mixed model can be found in Appendix 1—tables 1–5. The network of interactions was constructed using COPASI (version 4.23, build 184)

(http://www.nrcam.uchc.edu, http://copasi.org/) (Hoops et al., 2006). The model was built in COPASI to leverage the inbuilt fitting techniques for initial parameter guesses pre-FRET.

Well-mixed computational results

The network shown in Figure 1b has been shown to exhibit oscillations through cAMP variation due to the action of PKA on IP3 receptors and KATP plasma membrane channels. This network has been studied in many labs previous work (Fridlyand et al., 2007; Fridlyand et al., 2003; Han et al., 1999; Ni et al., 2011) and has been explored to show:

The system can be moved in and out-of-phase through tuning of AC (Fridlyand and Philipson, 2016).
The oscillation rate can be tuned through PKA feedback to KATP channels (Ni et al., 2011).
Both Ca²⁺ sensitive PDE and AC is necessary for oscillations to occur (Ni et al., 2011).

Our model results agree with the above findings.

Variations in the connection strength of sources and sinks also introduced another finding— changing the component’s (source or sink) variability in the regime of Ca²⁺ spiking (0.1-1.2 µM) a switch in-phase can also be observed, Appendix 1—figure 2. Compared to well-mixed results, Appendix 1—figure 2a, by decreasing the activation rate of CaMAC from 59.5 to 0.6 s⁻¹, a switch to out-of-phase can occur, Appendix 1—figure 2c. Yet, we notice how increasing basal PDE activity does not switch the phase Appendix 1—figure 2b, only after subsequently decreasing PDE and CaM association does the phase change Appendix 1—figure 2d. This is due to both changes being required for the activity variation of PDE to outperform that of AC. These findings lead to more questions about how source and sink activities relate in a variable Ca²⁺ regime.

Well-mixed reaction tables

Appendix 1—table 1

Voltage gated channel reactions.

#	Species	Expression	Parameters	Ref.
1	Membrane Voltage	$\frac{d V}{d t} = \frac{- I_{C a} - I_{K} - I_{L} - I_{K C a}}{C_{m}}$	$C_{m}$ = 5.3 pF	Ni et al., 2011
2	Ca²⁺ current	$I_{C a} = g_{C a} m_{\infty} (V - E_{C a})$	$g_{C a}$ = 600 pS, $E_{C a}$ = 100 mV	Ni et al., 2011
3	Fraction of open VGCC (at steady state)	$m_{\infty} = \frac{1}{2} (1 + t a n h (\frac{V - v_{1}}{v_{2}}))$	$v_{1}$ = -20 mV, $v_{2}$ = 24 mV	Ni et al., 2011
4	K⁺ current	$I_{K} = g_{K} w (V - E_{K})$	$g_{K}$ = 240 pS, $E_{C a}$ = -75 mV	Ni et al., 2011
5	Fraction of open K⁺ channels (at steady state)	$\frac{d w}{d t} = \frac{Φ (w_{\infty - w})}{τ}$	ϕ = 35 $\frac{1}{s}$	Ni et al., 2011
6	Time constant for K⁺ channel open probability	$τ = \frac{1}{c o s h (\frac{V - v_{3}}{2 v_{4}})}$	$v_{3}$ = -16 mV $v_{4}$ = 11.2 mV	Ni et al., 2011
7	leak current	$I_{L} = g_{L} (V - E_{L})$	$g_{L}$ = 150 pS, $E_{L}$ = -75 mV	Zhang et al., 2005
8	Ca²⁺ gated K⁺ current	$I_{K C a} = g_{K C a} \frac{C a}{C a + K_{K C a}} (V - E_{K})$	$g_{K C a}$ = 2000 pS, $E_{C a}$ = -75 mV, $K_{K C a}$ = 5 $μ M$	Ni et al., 2011

Appendix 1—table 2

Ca²⁺flux and reactions.

#	Reaction	Reaction flux	Kinetic parameters	Ref.
9	→ Ca	$j_{C a_{v}} (1 + k_{P K A_{V}} [P K A]) + j_{C a_{l}}$	$k_{P K A V}$ = 3000 s⁻¹ · µM⁻¹	Ni et al., 2011
10	$j_{C a_{V}}$	$f_{i} (- α I_{C a} - v_{L P M} [C a])$	$f_{i}$ = 1 x 10⁻⁵ $α$ = 0.0045 $μ M$ · f^A−1 · s⁻¹ $v_{L P M}$ = 75 s⁻¹	Ni et al., 2011
11	$j_{C a_{I}}$	$\frac{C k_{I P 3 R} [P K A] A [C a]}{1 + A [C a] + [B] [C a]^{2}} ([C a_{s t o r e s}] - [C a]) - \frac{V_{s} [C a]^{2}}{K_{s}^{2} + [C a]^{2}}$	K_s = 10 µM, Ca_stores = 1.56 µM, V_s = 0.1 µM · s⁻¹, k_IP3R = 0.05, A = 0.2869µM⁻¹, B = 2.869µM⁻², C = 0.2133	Ni et al., 2011
12	2Ca + CaM ↔ Ca²CaM	$k_{f} [C a] [C a M] - k_{r} [C a^{2} C a M]$	$k_{f}$ = 3.6 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 8 $s^{- 1}$	Lai et al., 2015
13	Ca + Ca²CaM ↔ Ca³CaM	$k_{f} [C a] [C a^{2} C a M] - k_{r} [C a^{3} C a M]$	$k_{f}$ =11 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 195 $s^{- 1}$	Lai et al., 2015
14	Ca + Ca³CaM ↔ Ca⁴CaM	$k_{f} [C a] [C a^{3} C a M] - k_{r} [C a^{4} C a M]$	$k_{f}$ = 59 $s^{- 1}$ ⋅µM⁻¹, $k_{r}$ = 500 $s^{- 1}$	Lai et al., 2015
15	AC + Ca²CaMCaM ↔ CaM · AC	$k_{f} [A C] [C a^{2} C a M] - k_{r} [C a M \cdot A C]$	$k_{f}$ = 1.7 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 10 $s^{- 1}$	Masada et al., 2009; Masada et al., 2012
16	CaM · AC + 2Ca ↔ AC*	$K_{c a t} \frac{[C a] [C a M \cdot A C]}{C a + K_{m}} - k_{r} [A C *]$	$K_{c a t}$ = 59.5 $s^{- 1}$ , $K_{m}$ = 0.1 $μ M$ , $k_{r}$ = 10 $s^{- 1}$	Masada et al., 2009; Masada et al., 2012
17	PDE + Ca²CaM ↔ CaM · PDE	$k_{f} [P D E] [C a^{2} C a M] - k_{r} [C a M \cdot P D E]$	$k_{f}$ = 435 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 1 $s^{- 1}$	Ang and Antoni, 2002
18	CaM · PDE + 2Ca ↔ PDE*	$K_{c a t} \frac{[C a] [C a M \cdot P D E]}{C a + K_{m}} - k_{r} [P D E *]$	$K_{c a t}$ = 1.81 $s^{- 1}$ , $K_{m}$ = 0.18 $μ M$ , $k_{r}$ = 1 $s^{- 1}$	Ang and Antoni, 2002
19	PDE + Ca⁴CaM ↔ PDE*	$k_{f} [P D E] [C a^{4} C a M] - k_{r} [P D E *]$	$k_{f}$ = 435 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 1 $s^{- 1}$	Ang and Antoni, 2002

*denotes the activated form

Appendix 1—table 3

cAMP reactions.

#	Reaction	Reaction flux	Kinetic parameters	Ref.
20	→ cAMP	$k_{b a s e} ([C a M \cdot A C] + [A C] + [A C_{i n d}]) + k_{a c t} [A C *]$	$k_{b a s e}$ = 0.1 $s^{- 1}$ , $k_{a c t}$ = 0.785 $s^{- 1}$	Leonid E. Fridlyand & Philipson, 2016; Ni et al., 2011
21	cAMP→	$k_{b a s e} [c A M P] \frac{[C a M \cdot P D E] + [P D E]}{[c A M P] + K_{m}} + k_{a c t} \frac{[c A M P] [P D E *]}{[c A M P] + K_{m}}$	$k_{b a s e}$ = 0.2 $s^{- 1}$ , $K_{m}$ = 0.6 $μ M$ $k_{a c t}$ = 2.5 $s^{- 1}$	Ni et al., 2011; Fridlyand and Philipson, 2016
22	cAMP→	$V_{i n d} \frac{[c A M P]}{[c A M P] + K_{m}}$	$V_{i n d}$ = 2.5 $μ M$ ⋅ s⁻¹, $K_{m}$ = 1.4 $μ M$	Ni et al., 2011; Fridlyand and Philipson, 2016
23	cAMP + R2 → R2_b	$k_{f} [c A M P] [R 2] - k_{r} [R 2_{b}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.00033 $s^{- 1}$	Boras et al., 2014
24	cAMP + R2_b → R2_ba	$k_{f} [c A M P] [R 2_{b}] - k_{r} [R 2_{b a}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.00105 $s^{- 1}$	Boras et al., 2014
25	cAMP + R2_b → R2_bb	$k_{f} [c A M P] [R 2_{b}] - k_{r} [R 2_{b b}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.00132 $s^{- 1}$	Boras et al., 2014
26	cAMP + R2_ba → R2_bba	$k_{f} [c A M P] [R 2_{b a}] - k_{r} [R 2_{b b a}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.0013 $s^{- 1}$	Boras et al., 2014
27	cAMP + R2_bb → R2_bba	$k_{f} [c A M P] [R 2_{b b}] - k_{r} [R 2_{b b a}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.00103 $s^{- 1}$	Boras et al., 2014
28	cAMP + R2_bba → R2_bbaa	$k_{f} [c A M P] [R 2_{b b a}] - k_{r} [R 2_{b b a a}]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 0.0114 $s^{- 1}$	Boras et al., 2014
29	PKA + R2 → R2C	$k_{f} [P K A] [R 2] - k_{r} [R 2 C]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 1.26E-7 $s^{- 1}$	Boras et al., 2014
30	PKA + R2_b → R2_bcC	$k_{f} [P K A] [R 2_{b}] - k_{r} [R 2_{b} C]$	$k_{f}$ = 1 $s^{- 1}$ ⋅ µM⁻¹, $k_{r}$ = 2.52E-7 $s^{- 1}$	Boras et al., 2014

Appendix 1—table 4

cAMP reactions (cont.).

#	Reaction	Reaction flux	Kinetic parameters	Ref.
31	PKA + R2_ba → R2_ba C	$k_{f} [P K A] [R 2_{b a}] - k_{r} [R 2_{b a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ $k_{r}$ = 3.4E-6 $s^{- 1}$	Boras et al., 2014
32	PKA + R2_bba → R2_bba C	$k_{f} [P K A] [R 2_{b b a}] - k_{r} [R 2_{b b a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.000936 $s^{- 1}$	Boras et al., 2014
33	PKA + R2_bbaa → R2_bbaa C	$k_{f} [P K A] [R 2_{b b a a}] - k_{r} [R 2_{b b a a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.645 $s^{- 1}$	Boras et al., 2014
34	cAMP + R2C → R2_b C	$k_{f} [c A M P] [R 2 C] - k_{r} [R 2_{b} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ , $k_{r}$ = 0.000659 s^-1	Boras et al., 2014
35	cAMP + R2_b C → R2_ba C	$k_{f} [c A M P] [R 2_{b} C] - k_{r} [R 2_{b a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ , $k_{r}$ = 0.0142 s^-1	Boras et al., 2014
36	cAMP + R2_ba C → R2_bbc C	$k_{f} [c A M P] [R 2_{b a} C] - k_{r} [R 2_{b b a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.0142 s^-1	Boras et al., 2014
37	cAMP + R2_bba C → R2_bbaa C	$k_{f} [c A M P] [R 2_{b b a} C] - k_{r} [R 2_{b b a a} C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 7.84 s^-1	Boras et al., 2014
38	PKA + R2_b C → R2_b C₂	$k_{f} [P K A] [R 2_{b} C] - k_{r} [R 2_{b} C_{2}]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.00324 s^-1	Boras et al., 2014
39	PKA + R2C → R2C₂	$k_{f} [P K A] [R 2 C] - k_{r} [R 2 C_{2}]$	$k_{f}$ = 1 $s^{- 1}$ · µM⁻¹, $k_{r}$ = 2.81E-6 s^-1	Boras et al., 2014
40	PKA + R2_ba C → R2_ba C₂	$k_{f} [P K A] [R 2_{b a} C] - k_{r} [R 2_{b a} C_{2}]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.666 s^-1	Boras et al., 2014
41	cAMP + R2C₂→ R2_b C₂	$k_{f} [c A M P] [R 2 C_{2}] - k_{r} [R 2_{b} C_{2}]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 0.762 s^-1	Boras et al., 2014
42	cAMP + R2_b C₂→ R2_ba C₂	$k_{f} [c A M P] [R 2_{b} C_{2}] - k_{r} [R 2_{b a} C_{2}]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹, $k_{r}$ = 2.91 s^-1	Boras et al., 2014

Appendix 1—table 5

Initial conditions.

#	Species	Initial value
IC1	AC8	1 µM
IC2	AC_ind	1 µM
IC3	Ca²⁺	1 µM
IC4	CaM	10 µM
IC5	cAMP	0.1 µM
IC6	PDE1	1 µM
IC7	PDE4	0.4 µM
IC8	R2C2	0.4 µM
IC9	V	−60 mV

Minimal model to explore Ca²⁺-cAMP phase behavior

Based on the well-mixed model, we propose a minimal circuit to understand the phase behavior of Ca²⁺-cAMP. Consider the system in Figure 2c, where Ca²⁺ is the stimulus which uses a pulse train to set influx times, AC is an activator, PDE1C is an inhibitor, and cAMP is the response element.

Assumptions

The system is in a state such that the change in cAMP allows for further Ca²⁺ influx in a semi-predictable manner, as represented by the pulse train. Therefore, this system is deemed to be stable if there exists a Ca²⁺ value that gives a stable solution for cAMP. All constants must be positive to remain physically relevant. We assume that there exists a constant independent source and sink within the system. The Ca²⁺-dependent and independent sources are localized homogeneously on the membrane, with both sinks located uniformly in the cytosol. For simplicity, we assume that the activation function for both $A C$ and $P D E 1 C$ are linear functions of Ca²⁺ of the form $a S + b$ .

Governing Equations

We define a stimulus (S, Ca²⁺) and a response element (R, cAMP). The well-mixed rate of change function for cAMP is then given by,

\frac{d R}{d t} = v_{1} (a_{1} S + b_{1}) - v_{2} (a_{2} S + b_{2}) R + v_{i p} - v_{i d} R

Here, $v_{1}$ denotes the velocity of cAMP production by AC, $v_{2}$ the rate of degradation by PDE, and $v_{i p}$ and $v_{i d}$ the rates of independent production and degradation respectively.

Analytical solutions for the minimal model

To analyze if the system lies in an in- or out-of-phase state we find the direction of the system change after initialization to $S_{0}$ [i.e. the basal stimulus (initial concentration of Ca²⁺)]. First, we must solve for $R_{0}$ (initial concentration of cAMP) by setting $\frac{d R}{d t} = 0$ , we find:

R_{0} = \frac{v_{1} (a_{1} S_{0} + b_{1}) + v_{i p}}{v_{2} (a_{2} S_{0} + b_{2}) + v_{i d}}

We start the system at equilibrium and prescribe a discontinuous pulse of $S$ from $S_{0}$ to $S_{h}$ , akin to VGCC opening allowing Ca²⁺ flux. Therefore, solving for the sign of R, we obtain

\begin{array}{ccc} \frac{d R}{d t} |_{t = 0} & = & v_{1} (a_{1} S_{h} + b_{1}) + v_{i p} - (v_{2} (a_{2} S_{h} + b_{2}) + v_{i d}) \frac{v_{1} (a_{1} S_{0} + b_{1}) + v_{i p}}{v_{2} (a_{2} S_{0} + b_{2}) + v_{i d}} \\ \frac{d R}{d t} |_{t = 0} & = & \frac{(v_{1} a_{1} (v_{2} b_{2} + v_{i d}) - v_{2} a_{2} (v_{1} b_{1} + v_{i p})) (S_{h} - S_{0})}{v_{2} (a_{2} S_{0} + b_{2}) + v_{i d}} \end{array}

Since we consider the sign we can then characterize the solution by

\begin{array}{ccc} \frac{v_{1} a_{1} (v_{2} b_{2} + v_{i d})}{v_{2} a_{2} (v_{1} b_{1} + v_{i p})} > 1 & in phase (low τ) \\ \frac{v_{1} a_{1} (v_{2} b_{2} + v_{i d})}{v_{2} a_{2} (v_{1} b_{1} + v_{i p})} = 1 & transition (flat line) \\ \frac{v_{1} a_{1} (v_{2} b_{2} + v_{i d})}{v_{2} a_{2} (v_{1} b_{1} + v_{i p})} < 1 & out of phase (high τ) \end{array}

This model was then computationally run to confirm the results with arbitrary parameters (Figure 2c). The system was simulated with all parameters set to 1 except $v_{1}$ and $v_{2}$ at values of 4 and 2, respectively and vice-versa. The pulse train is ON for 1 sec (S_h = 10) and OFF for 3 sec (S₀ = 1), and repeats until the simulation total time is reached.

Numerical implementation of the minimal model

We performed numerical simulations in MATLAB R2018b and checked against analytic solutions provided in the previous section. Minimal model solutions are shown in Figure 2c.

Stochastic implementation of the minimal model

To further explore the possible effects of noise and small numbers of molecules on the simplified four component model, we created a stochastic version of the model. Due to the nature of stochastic simulations, all reactions were rewritten to multistep interactions. In order to allow for oscillatory kinetics at the membrane, a Lotka-Volterra model (in a larger volume, for numerical stability) was utilized to simulate depolarization events shown in Appendix 1—table 8. The system was solved using the Gibson method in the Virtual Cell modeling platform. All Ca²⁺, AC, and PDE reactions are shown in Appendix 1—tables 6–9.

Appendix 1—table 6

Initial conditions for the stochastic model.

Species	Number of molecules in the system
Ca²⁺	10
PDE	45
PDE*	5
AC	10–200 (log scale)
s0	10000
s4	1000

Appendix 1—table 7

Compartment volumes for the stochastic model.

Compartment	Compartment size
Cytosol	0.083 μm³
Membrane	0.28 μm²
Oscillator	6563.0 μm²

Appendix 1—table 8

Lotka-Volterra model to initiate the membrane depolarization events.

#	Reaction	Parameters
1	s0 → 2s0	K_f = 1.7 s⁻¹
2	s0 + s4 → 2s4	K_f = 1 μm² ⋅ molecule^-1 ⋅ s^-1
3	s4 →	K_f = 1.7 s^-1

Appendix 1—table 9

simplified model reactions for the stochastic model.

#	Reaction	Parameters
4	s4 → > s4 + Ca²⁺	K_f = 150 s⁻¹
5	Ca²⁺→	K_f = 500 μm⁻² ⋅ molecule ⋅ s⁻¹ ⋅ μM⁻¹
6	AC + Ca²⁺ → AC*	K_f = 200 s⁻¹ ⋅ μM⁻¹ K_r = 300 s⁻¹
7	AC* → > cAMP + AC*	K_f = 300 s⁻¹
8	AC → cAMP + AC	K_f = 5 s⁻¹
9	PDE + Ca²⁺ → PDE*	K_f = 150 s⁻¹ ⋅ μM⁻¹ K_f = 100 s⁻¹
10	PDE + cAMP → PDE	K_f = 50 s⁻¹ ⋅ μM⁻¹
11	PDE* + cAMP → PDE*	K_f = 1000 s⁻¹ ⋅ μM⁻¹
12	→ cAMP	J = 9000 μm⁻² ⋅ molecule ⋅ s⁻¹
13	cAMP →	K_f = 50 s⁻¹

Solutions to the stochastic model are shown in Figure 6—figure supplement 2. The Virtual Cell Model, ‘Minimum Stochastic Model BetaCell’ by user ‘mgetz’, can be accessed within the VCell software (available at https://vcell.org). The Virtual Cell is supported by NIH Grant R24 GM134211 from the National Institute for General Medical Sciences (Cowan et al., 2012; Schaff et al., 1997).

Simulations of the full spatial systems

For computational simplification, simulations were performed with a Gaussian profile on the top boundary, the size of the domain and Gaussian profile were informed by STORM images, Figure 3. The system is a hexagonal prism with outer diameter of 0.4 µm and depth 0.6 µm with periodic boundary conditions in the x and y planes, see Figure 4a. The top plane is assumed to be the membrane and the bottom is a no-flux condition. For the membrane plane, a Gaussian profile was normalized such that the average value is 4 x 10⁻¹⁰ mol/m². The Ca²⁺ sensitive AC initial condition (Gaussian profile) is fixed for all simulation times by setting the diffusion constant to ≈0. To test the effect of localization, a percentage of AC mass was moved from the Gaussian profile into a uniform profile (a 25% Localization means that 25% of the total AC mass is within the Gaussian profile and the remaining 75% is within the uniform profile). Examples of surface profiles at multiple % localizations can be seen in Figure 4b.

Assumptions

Membrane patterns are pre-existing and not affected by a single signaling event such that no diffusion occurs.
Clustering events were approximated by a Gaussian profile in the center of the hexagonal prism.

Model development

Although the well-mixed system shows the ability to oscillate in and out-of-phase in a well-behaved manner, the spatial model will not produce this effect if a homogeneous boundary is present. Previous studies have looked at spatial gradients in the context of cAMP and PKA and found out that localization of some species is necessary for the system to demonstrate different behaviors (Yang et al., 2016). Therefore, when moving to a 3D spatial map, we must consider how the two solution regimes can be recovered. Experimental data suggest that AKAP dimerizes and may form oligomers (Gao et al., 2011; Gold et al., 2011), which is important for the function of these cells (Zhang et al., 2013). This could allow spatial instabilities like those seen in Haselwandter et al., 2015 used to describe post synaptic domains. The final Gaussian profile on the top boundary had the size of the domain and Gaussian profile informed by STORM images of AC clustering (Figure 3). Statistics of the images show that, on average, 90% of the AC sits within 54 nm of the cluster center. The system was determined to be 0.35 x 0.35 x 0.6 µm hexagonal prism with a Gaussian standard deviation of 25 nm.

The kinetic parameters used were the same as for the well-mixed model, except for a few cases in which tuning through surface/volume relationships was needed. Post fitting was performed to further refine the relationship between PDE and AC through use of obtained FRET measurements (see section 4). Due to the large computational expense of the model, the PDE interactions were reduced to two steps (Appendix 1—table 10).

Appendix 1—table 10

Ca²⁺ Flux and reactions modified from Appendix 1—table 2.

#	Reaction	Reaction flux	Kinetic parameters	Ref.
S1	→ Ca²⁺	$j_{C a_{V}} (1 + k_{P K A_{V}} [P K A]) + J_{C a_{I}}$	$k_{P K A V}$ = 100 $μ M^{- 1}$	FRET constraint
S2	$j_{C a_{V}}$	$f_{i} (- α I_{C a} - v_{L P M} [C a])$	$f_{i}$ = 1 x 10⁻⁶, $α$ = 4.15 x 10⁵mol · m⁻² · A⁻² · s⁻¹,, $v_{L P M}$ = 7.5 x 10⁻⁴m · s⁻¹	Ni et al., 2011
S3	$j_{C a_{I}}$	$\frac{C k_{I P A R} [P K A] A [C a]}{1 + A [C a] + [B] [C a]^{2}} ([C a_{s t o r e s}] - [C a]) - \frac{V_{s} [C a]^{2}}{K_{s}^{2} + [C a]^{2}}$	$K_{s}$ = 10 $μ M$ , $C a_{s t o r e s}$ = 1.56 $μ M$ , $V_{s}$ = 0.1 µM · s⁻¹, $k_{I P 3 R}$ = 0.05 $μ M^{- 1} s^{- 1}$ , A = 0.2869 $μ M^{- 1}$ , B = 2.869 $μ M^{- 2}$ , C = 0.2133	Ni et al., 2011
S4	AC + Ca²CaM ↔ CaM · AC	$k_{f} [A C] [C a^{2} C a M] - k_{r} [C a M \cdot A C]$	$k_{f}$ = 10.8 $s^{- 1}$ · µM⁻¹, $k_{r}$ = 10 $s^{- 1}$	FRET constraint
S5	CaM · AC + 2Ca2+ ↔ AC*	$K_{c a t} \frac{[C a] [C a M \cdot A C]}{C a + K_{m}} - k_{r} [A C *]$	$K_{c a t}$ = 90 $s^{- 1}$ , $K_{m}$ = 1 $μ M$ , $k_{r}$ = 10 $s^{- 1}$	FRET constraint
S6	PDE + Ca²CaM ↔ CaM · PDE	$k_{f} [P D E] [C a^{2} C a M] - k_{r} [C a M \cdot P D E]$	$k_{f}$ = 0.25 $s^{- 1}$ · µM⁻¹, $k_{r}$ = 1 $s^{- 1}$	FRET constraint
S7	CaM · PDE + 2Ca2+ ↔ PDE*	$K_{c a t} \frac{[C a] [C a M \cdot P D E]}{C a + K_{m}} - k_{r} [P D E *]$	$K_{c a t}$ = 60 $s^{- 1}$ , $K_{m}$ = 1 $μ M$ , $k_{r}$ = 1 $s^{- 1}$	FRET constraint
S8	PDE + Ca⁴CaM ↔ PDE*	$k_{f} [P D E] [C a^{4} C a M] - k_{r} [P D E *]$	$k_{f}$ = 0.25 $s^{- 1}$ · µM⁻¹, $k_{r}$ = 1 $s^{- 1}$	FRET constraint

Numerical simulation

The well-mixed network was imported into COMSOL Multiphysics5.4 (Build:295), to solve the spatial model with inhomogeneous boundary conditions at the plasma membrane.

Reaction-Diffusion Partial differential equations (PDEs)

Consider the same one-reaction system as the well-mixed system:

A + 2 B ⇌ 3 C,

where the forward and backward rates are $k_{1}$ and $k_{2}$ . But, now we spatially discretize the system, the partial differential equations describing the dynamics of species $A, B$ , and $C$ with unrestricted diffusion are then:

\begin{array}{cc} \frac{\partial [A]}{\partial t} & = D_{A} \nabla^{2} [A] + k_{2} {[C]}^{3} - k_{1} [A] {[B]}^{2} \\ \frac{\partial [B]}{\partial t} & = D_{B} \nabla^{2} [B] + 2 k_{2} {[C]}^{3} - 2 k_{1} [A] {[B]}^{2} \\ \frac{\partial [C]}{\partial t} & = D_{C} \nabla^{2} [C] + 3 k_{1} [A] {[B]}^{2} - 3 k_{2} {[C]}^{3} . \end{array}

Where $\nabla^{2}$ is the Laplacian and $D_{A}$ is a diagonal matrix of diffusion coefficients for component A. Yet within the cell, all reactions do not occur in the free volume. In fact, most interactions occur on a membrane surface, which requires a different boundary condition.

Reactions on the system boundary

Now let us assume the previous reaction occurs on the boundary and C is a membrane species on surface $Ω$

A + 2 B ⇌ 3 C,

This would mean within the volume there is only free diffusion,

\begin{array}{r} \frac{\partial [A]}{\partial t} & = D_{A} \nabla^{2} [A] \\ \frac{\partial [B]}{\partial t} & = D_{B} \nabla^{2} [B] \end{array}

At the surface we would have a Neumann boundary condition, that is, a defined flux occurring into the boundary normal. We define the reaction as a flux occurring at a membrane surface ( $Ω$ )

\begin{array}{cc} \nabla A (x, t) \cdot \hat{n} (x) |_{\partial Ω} & = k_{2} {[C]}^{3} - k_{1} [A] {[B]}^{2} \\ \nabla B (x, t) \cdot \hat{n} (x) |_{\partial Ω} & = 2 k_{2} {[C]}^{3} - 2 k_{1} [A] {[B]}^{2} \\ \nabla C (x, t) \cdot \hat{n} (x) |_{\partial Ω} & = 3 k_{1} [A] {[B]}^{2} - 3 k_{2} {[C]}^{3} \end{array}

Here, $\hat{n}$ is the unit normal to $Ω$ . We therefore have defined a flux between the volume and the membrane $Ω$ . Finally, we consider $Ω$ as a 2-dimensional surface existing at the system boundary with free diffusion of C;

\frac{\partial [C]}{\partial t} = D_{C} Δ [C]

where Δ is the surface Laplacian.

PKA spatial response mirrors cAMP

PKA with full diffusion values ( $\approx 10 \frac{μ m^{2}}{s}$ ) did not follow cAMP dynamics and only elicited a single global response. We asked if this was due to AKAP patterning at the surface, which should mirror the AC profile at the surface. Adding this interaction into the model did not allow any sizable spatial gradient of PKA to develop. Recent work has suggested the PKA catalytic subunit in the presence of non-excess cAMP is effectively activated but its diffusion is restricted (Smith et al., 2017). By varying the diffusion constant, we found PKA activity could follow cAMP dynamics in the nanocluster and PM compartments in our computational model for restricted diffusivities (Figure 6–figure supplement 1a). We experimentally tested this prediction using the AKAP79-fused and PM-targeted PKA activity sensors and found that indeed PKA activity did follow the cAMP-Ca²⁺ phase relationship within the two compartments (Figure 6–figure supplement 1b-d). This suggests that anchored PKA holoenzyme action is much more restricted than originally anticipated.

Comparisons with experimental data

Raw FRET data (Figure 1) was used for model refinement of the cAMP-Ca²⁺ phase relationship. The experimental and model data were qualitatively compared against each other to preserve oscillation time, amplitude, range, and waveform. Although a quantitative calibration was not performed to map FRET ratio to cAMP concentration, we utilized the fact that these sensors are known to measure cAMP in the range of ≈0.1-10 µM to constrain our model. Voltage-gated channel sensitivities were not tuned, and only connection strengths between CaM to ACs and PDEs, which are largely less constrained in comparison, were varied. Values modified from the well-mixed model values can be found in Appendix 1—tables 10–13.

Model validation and predictions

The model generated predictions relating concentration perturbations (AC, PDE, etc.), disruption of patterning (AC binding disruption), and their changes to the phase of the signal. The system, once moved to the spatial model, was allowed free parameters along the six-component axis for CaM connection and cAMP production strength of ACs and PDEs. This includes the flux differential between basal and activated ACs and PDEs.

Reaction tables of modified well-mixed parameters for the spatial model

Appendix 1—table 11

cAMP reactions modified from Appendix 1—table 3.

#	Reaction	Reaction flux	Kinetic parameters	Ref.
S9	→ cAMP	$k_{b a s e} ([C a M \cdot A C] + [A C] + [A C_{i n d}]) + k_{a c t} [A C *]$	$k_{b a s e}$ = 0.2 $s^{- 1}$ , $k_{a c t}$ = 23.55 $s^{- 1}$ , ACind = 3 x 10⁻⁸	FRET constraint
S10	cAMP →	$\begin{array}{ll} k_{b a s e} [c A M P] \frac{[C a M \cdot P D E] + [P D E]}{[c A M P] + K_{m}} + \\ k_{a c t} \frac{[c A M P] [P D E *]}{[c A M P] + K_{m}} \end{array}$	$k_{b a s e}$ = 0.6 $s^{- 1}$ , $K_{m}$ = 0.6 µM $k_{a c t}$ = 720 $s^{- 1}$	FRET constraint
S11	cAMP →	$V_{i n d} \frac{[c A M P]}{[c A M P] + K_{m}}$	$V_{i n d}$ = 0.25 µM · s⁻¹, $K_{m}$ = 1.4 µM	FRET constraint
S12	2cAMP + R2C2 → R2C + PKA	$k_{f} {[c A M P]}^{2} [R 2 C 2] - k_{r} [R 2 C] [P K A]$	$k_{f}$ = 20 min⁻¹ · µM⁻², $k_{r}$ = 12 min⁻¹ · µM⁻¹
S13	2cAMP + R2C → R2 + PKA	$k_{f} {[c A M P]}^{2} [R 2 C] - k_{r} [R 2] [P K A]$	$k_{f}$ = 20 min⁻¹ · µM⁻², $k_{r}$ = 12 min⁻¹ · µM⁻¹

Appendix 1—table 12

Additional AKAP interactions for the spatial model.

#	Reaction	Reaction flux	Kinetic parameters	Ref
S14	AKAP + R2 → AKAP-R2	$k_{f} [R 2] [A K A P - R 2 C 2] - k_{r} [A K A P - R 2]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ $k_{f}$ = 0.1 $s^{- 1}$	Est.
S15	AKAP + R2C → AKAP-R2C	$k_{f} [R 2 C] [A K A P - R 2 C 2] - k_{r} [A K A P - R 2 C]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ $k_{f}$ = 0.1 $s^{- 1}$	Est.
S16	AKAP + R2C2 → AKAP-R2C2	$k_{f} [R 2 C 2] [A K A P - R 2 C 2] - k_{r} [A K A P - R 2 C 2]$	$k_{f}$ = 1 s⁻¹ · µM⁻¹ $k_{f}$ = 0.1 $s^{- 1}$	Est.
S17	2cAMP + AKAP-R2C2 → AKAP-R2C + PKA	$k_{f} {[c A M P]}^{2} [A K A P - R 2 C 2] - k_{r} [A K A P - R 2 C] [P K A]$	$k_{f}$ = 20 min⁻¹ · µM⁻², $k_{r}$ = 12 min⁻¹ · µM⁻²
S18	2cAMP + AKAP-R2C → AKAP-R2 + PKA	$k_{f} {[c A M P]}^{2} [A K A P - R 2 C] - k_{r} [A K A P - R 2] [P K A]$	$k_{f}$ = 20 min⁻¹ · µM⁻², $k_{r}$ = 12 min⁻¹ · µM⁻²

Appendix 1—table 13

Initial conditions and diffusion coefficients used in the model.

#	Species	Initial value	Diffusion	Ref.
SIC1	Ca²⁺	0.001 µM	100 $\frac{μ m^{2}}{s}$	(Donahue and Abercrombie, 1987), Est. from steady state
SIC2	CaM	2.9 µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC3	Ca₂CaM	0.1 µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC4	Ca₃CaM	4 x 10⁻³ µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC5	Ca₄CaM	1 x 10⁻² µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC6	R₂	0.04 µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC7	R₂C₂	0.2 µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC8	PKA	0.05 µM	0.01 $\frac{μ m^{2}}{s}$	Est. from steady state, diffusion fitted
SIC9	PDE1	0.9 µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC10	PDE1_act	1 x 10⁻³ µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC11	CaMPDE1	1 x 10⁻³ µM	10 $\frac{μ m^{2}}{s}$	Est. from steady state*
SIC12	AC	4 x 10⁻¹⁰ $\frac{m o l}{m^{2}}$	0	Est.
SIC13	CaMAC	0 $\frac{m o l}{m^{2}}$	0	Est.
SIC14	AC_act	0 $\frac{m o l}{m^{2}}$	0	Est.
SIC15	AC_ind	4 x 10⁻¹⁰ $\frac{m o l}{m^{2}}$	1 $\frac{μ m^{2}}{s}$	Est.^†
SIC16	V	−60 mV	1 $\frac{μ m^{2}}{s}$	Fridlyand et al., 2003; Ni et al., 2011^†
SIC17	cAMP	4 x 10⁻⁶ µM	60 $\frac{μ m^{2}}{s}$	Agarwal et al., 2016, Est. from steady state

*For all cytosolic species without well-constrained diffusions, we used 10 $\frac{μ m^{2}}{s}$

^†For all membrane species without well-constrained diffusions we used 1 $\frac{μ m^{2}}{s}$

Appendix 1—figure 1

Download asset Open asset

Phase of cAMP correlates with the expression level of the AKAP79-(Ci/Ce)Epac2-camps.

Scatter plot of the time lag (sec) and the YFP donor channel intensity (normalized to non-saturating maximum) for each cell expressing AKAP79-(Ci/Ce)Epac2-camps. Cells with higher expression of the probe correlated with a longer time lag, therefore a YFP intensity threshold was designated for analysis purposes.

Appendix 1—figure 2

Download asset Open asset

Phase is driven by activity variability within the Ca²⁺ oscillatory regime.

The phase of the system can be switched by tuning the association of CaM to sources (ACs) and sinks (PDEs). (a) At base system conditions, the system acts in an in-phase manner. (b) Decreasing the rate of Ca²⁺ association to the AC-CaM complex causes the phase to switch to out-of-phase. (c) Increasing basal PDE activity does not allow a phase switch only after decreasing PDE and CaM association rates will the system allow a phase switch (d). A phase switch is controlled by the variability in the activity of source or sink. If the sink dominates, then the system is out-of-phase. If the source dominates, the system is in-phase.

Data availability

Data has been deposited at Dryad (https://doi.org/10.6075/J0NP22TK). Analysis scripts have been deposited on GitHub (https://github.com/btenner/calcium-cAMP-PKA; copy archived at https://archive.softwareheritage.org/swh:1:rev:5c4c971e275cdc3e97b891e4b507393b209f80a1/).

The following data sets were generated

1. Tenner B
(2020) Dryad Digital Repository
Spatially compartmentalized phase regulation of a Ca2+-cAMP-PKA oscillatory circuit.

https://doi.org/10.6075/J0NP22TK

References

(2016) Mechanisms restricting diffusion of intracellular cAMP
Scientific Reports 6:19577.

https://doi.org/10.1038/srep19577
- PubMed
- Google Scholar
1. Akerboom J
2. Carreras Calderón N
3. Tian L
4. Wabnig S
5. Prigge M
6. Tolö J
7. Gordus A
8. Orger MB
9. Severi KE
10. Macklin JJ
11. Patel R
12. Pulver SR
13. Wardill TJ
14. Fischer E
15. Schüler C
16. Chen TW
17. Sarkisyan KS
18. Marvin JS
19. Bargmann CI
20. Kim DS
21. Kügler S
22. Lagnado L
23. Hegemann P
24. Gottschalk A
25. Schreiter ER
26. Looger LL
(2013) Genetically encoded calcium indicators for multi-color neural activity imaging and combination with optogenetics
Frontiers in Molecular Neuroscience 6:2.

https://doi.org/10.3389/fnmol.2013.00002
- PubMed
- Google Scholar
1. Ang KL
2. Antoni FA
(2002) Reciprocal regulation of calcium dependent and calcium independent cyclic AMP hydrolysis by protein phosphorylation
Journal of Neurochemistry 81:422–433.

https://doi.org/10.1046/j.1471-4159.2002.00903.x
- PubMed
- Google Scholar
1. Barg S
2. Ma X
3. Eliasson L
4. Galvanovskis J
5. Göpel SO
6. Obermüller S
7. Platzer J
8. Renström E
9. Trus M
10. Atlas D
11. Striessnig J
12. Rorsman P
(2001) Fast exocytosis with few ca(2+) channels in insulin-secreting mouse pancreatic B cells
Biophysical Journal 81:3308–3323.

https://doi.org/10.1016/S0006-3495(01)75964-4
- PubMed
- Google Scholar
1. Beene DL
2. Scott JD
(2007) A-kinase anchoring proteins take shape
Current Opinion in Cell Biology 19:192–198.

https://doi.org/10.1016/j.ceb.2007.02.011
- PubMed
- Google Scholar
1. Bender AT
2. Beavo JA
(2006) Cyclic nucleotide phosphodiesterases: molecular regulation to clinical use
Pharmacological Reviews 58:488–520.

https://doi.org/10.1124/pr.58.3.5
- PubMed
- Google Scholar
(1998) Calcium--a life and death signal
Nature 395:645–648.

https://doi.org/10.1038/27094
- PubMed
- Google Scholar
1. Berry H
(2002) Monte carlo simulations of enzyme reactions in two dimensions: fractal kinetics and spatial segregation
Biophysical Journal 83:1891–1901.

https://doi.org/10.1016/S0006-3495(02)73953-2
- Google Scholar
1. Bock A
2. Annibale P
3. Konrad C
4. Hannawacker A
5. Anton SE
6. Maiellaro I
7. Zabel U
8. Sivaramakrishnan S
9. Falcke M
10. Lohse MJ
(2020) Optical mapping of cAMP signaling at the nanometer scale
Cell 182:1519–1530.

https://doi.org/10.1016/j.cell.2020.07.035
- PubMed
- Google Scholar
Preprint
1. Bohrer CH
2. Yang X
3. Weng X
4. Tenner B
5. Ross B
6. Mcquillen R
7. Zhang J
8. Roberts E
9. Xiao J
(2019) A pairwise distance distribution correction (DDC) algorithm for blinking-free super-resolution microscopy
bioRxiv.

https://doi.org/10.1101/768051
- Google Scholar
(2014) Using Markov state models to develop a mechanistic understanding of protein kinase A regulatory subunit riα activation in response to cAMP binding
Journal of Biological Chemistry 289:30040–30051.

https://doi.org/10.1074/jbc.M114.568907
- PubMed
- Google Scholar
1. Calebiro D
2. Maiellaro I
(2014) cAMP signaling microdomains and their observation by optical methods
Frontiers in Cellular Neuroscience 8:350.

https://doi.org/10.3389/fncel.2014.00350
- PubMed
- Google Scholar
1. Clapham DE
(2007) Calcium signaling
Cell 131:1047–1058.

https://doi.org/10.1016/j.cell.2007.11.028
- PubMed
- Google Scholar
1. Conti M
2. Mika D
3. Richter W
(2014) Cyclic AMP compartments and signaling specificity: role of cyclic nucleotide phosphodiesterases
Journal of General Physiology 143:29–38.

https://doi.org/10.1085/jgp.201311083
- PubMed
- Google Scholar
(1995) Adenylyl cyclases and the interaction between calcium and cAMP signalling
Nature 374:421–424.

https://doi.org/10.1038/374421a0
- PubMed
- Google Scholar
(2012) Spatial modeling of cell signaling networks
Methods in Cell Biology 110:195–221.

https://doi.org/10.1016/B978-0-12-388403-9.00008-4
- PubMed
- Google Scholar
(2008) Coexistence of amplitude and frequency modulations in intracellular calcium dynamics
Physical Review E 77:030903.

https://doi.org/10.1103/PhysRevE.77.030903
- PubMed
- Google Scholar
(2000) Tissue specificity and physiological relevance of various isoforms of adenylyl cyclase
American Journal of Physiology-Renal Physiology 279:F400–F416.

https://doi.org/10.1152/ajprenal.2000.279.3.F400
- PubMed
- Google Scholar
1. Dell'Acqua ML
2. Smith KE
3. Gorski JA
4. Horne EA
5. Gibson ES
6. Gomez LL
(2006) Regulation of neuronal PKA signaling through AKAP targeting dynamics
European Journal of Cell Biology 85:627–633.

https://doi.org/10.1016/j.ejcb.2006.01.010
- PubMed
- Google Scholar
1. Depry C
2. Allen MD
3. Zhang J
(2011) Visualization of PKA activity in plasma membrane microdomains
Mol. BioSyst. 7:52–58.

https://doi.org/10.1039/C0MB00079E
- PubMed
- Google Scholar
1. Donahue BS
2. Abercrombie RF
(1987) Free diffusion coefficient of ionic calcium in cytoplasm
Cell Calcium 8:437–448.

https://doi.org/10.1016/0143-4160(87)90027-3
- PubMed
- Google Scholar
1. Dou H
2. Wang C
3. Wu X
4. Yao L
5. Zhang X
6. Teng S
7. Xu H
8. Liu B
9. Wu Q
10. Zhang Q
11. Hu M
12. Wang Y
13. Wang L
14. Wu Y
15. Shang S
16. Kang X
17. Zheng L
18. Zhang J
19. Raoux M
20. Lang J
21. Li Q
22. Su J
23. Yu X
24. Chen L
25. Zhou Z
(2015) Calcium influx activates adenylyl cyclase 8 for sustained insulin secretion in rat pancreatic beta cells
Diabetologia 58:324–333.

https://doi.org/10.1007/s00125-014-3437-z
- PubMed
- Google Scholar
(2011) Calcium oscillations
Cold Spring Harbor Perspectives in Biology 3:a004226.

https://doi.org/10.1101/cshperspect.a004226
- PubMed
- Google Scholar
(2006) Oscillations of cyclic AMP in hormone-stimulated insulin-secreting beta-cells
Nature 439:349–352.

https://doi.org/10.1038/nature04410
- PubMed
- Google Scholar
(1997) Nuclear membrane dynamics and reassembly in living cells: targeting of an inner nuclear membrane protein in interphase and mitosis
Journal of Cell Biology 138:1193–1206.

https://doi.org/10.1083/jcb.138.6.1193
- PubMed
- Google Scholar
1. Everett KL
2. Cooper DM
(2013) An improved targeted cAMP sensor to study the regulation of adenylyl cyclase 8 by Ca2+ entry through voltage-gated channels
PLOS ONE 8:e75942.

https://doi.org/10.1371/journal.pone.0075942
- PubMed
- Google Scholar
(2003)
Modeling of Ca2+ flux in pancreatic beta-cells: role of the plasma membrane and intracellular stores

American Journal of Physiology. Endocrinology and Metabolism 285:E138–E154.
- Google Scholar
(2007) Regulation of cAMP dynamics by Ca²⁺ and G protein-coupled receptors in the pancreatic beta-cell: a computational approach
American Journal of Physiology-Cell Physiology 293:C1924–C1933.

https://doi.org/10.1152/ajpcell.00555.2006
- PubMed
- Google Scholar
(2010) Bursting and calcium oscillations in pancreatic beta-cells: specific pacemakers for specific mechanisms
American Journal of Physiology-Endocrinology and Metabolism 299:E517–E532.

https://doi.org/10.1152/ajpendo.00177.2010
- PubMed
- Google Scholar
1. Fridlyand LE
2. Philipson LH
(2016) Pancreatic beta cell G-Protein coupled receptors and second messenger interactions: a systems biology computational analysis
PLOS ONE 11:e0152869.

https://doi.org/10.1371/journal.pone.0152869
- PubMed
- Google Scholar
1. Gao S
2. Wang HY
3. Malbon CC
(2011) AKAP12 and AKAP5 form higher-order hetero-oligomers
Journal of Molecular Signaling 6:8.

https://doi.org/10.1186/1750-2187-6-8
- PubMed
- Google Scholar
1. Gao J
2. Hirata M
3. Mizokami A
4. Zhao J
5. Takahashi I
6. Takeuchi H
7. Hirata M
(2016) Differential role of SNAP-25 phosphorylation by protein kinases A and C in the regulation of SNARE complex formation and exocytosis in PC12 cells
Cellular Signalling 28:425–437.

https://doi.org/10.1016/j.cellsig.2015.12.014
- PubMed
- Google Scholar
1. Getis A
2. Franklin J
(1987) Second-Order neighborhood analysis of mapped point patterns
Ecology 68:473–477.

https://doi.org/10.2307/1938452
- Google Scholar
1. Gilon P
2. Ravier MA
3. Jonas JC
4. Henquin JC
(2002) Control mechanisms of the oscillations of insulin secretion in vitro and in vivo
Diabetes 51 Suppl 1:S144–S151.

https://doi.org/10.2337/diabetes.51.2007.S144
- PubMed
- Google Scholar
1. Gold MG
2. Stengel F
3. Nygren PJ
4. Weisbrod CR
5. Bruce JE
6. Robinson CV
7. Barford D
8. Scott JD
(2011) Architecture and dynamics of an A-kinase anchoring protein 79 (AKAP79) signaling complex
PNAS 108:6426–6431.

https://doi.org/10.1073/pnas.1014400108
- PubMed
- Google Scholar
1. Goraya TA
2. Masada N
3. Ciruela A
4. Willoughby D
5. Clynes MA
6. Cooper DM
(2008) Kinetic properties of Ca2+/calmodulin-dependent phosphodiesterase isoforms dictate intracellular cAMP dynamics in response to elevation of cytosolic Ca2+
Cellular Signalling 20:359–374.

https://doi.org/10.1016/j.cellsig.2007.10.024
- PubMed
- Google Scholar
1. Greenwald EC
2. Saucerman JJ
(2011) Bigger, better, faster: principles and models of AKAP anchoring protein signaling
Journal of Cardiovascular Pharmacology 58:462–469.

https://doi.org/10.1097/FJC.0b013e31822001e3
- PubMed
- Google Scholar
1. Guček A
2. Gandasi NR
3. Omar-Hmeadi M
4. Bakke M
5. Døskeland SO
6. Tengholm A
7. Barg S
(2019) Fusion pore regulation by cAMP/Epac2 controls cargo release during insulin exocytosis
eLife 8:e41711.

https://doi.org/10.7554/eLife.41711
- PubMed
- Google Scholar
1. Guldberg CM
2. Waage P
(1879) Ueber die chemische affinität. § 1. einleitung
Journal Für Praktische Chemie 19:69–114.

https://doi.org/10.1002/prac.18790190111
- Google Scholar
1. Han P
2. Werber J
3. Surana M
4. Fleischer N
5. Michaeli T
(1999) The calcium/calmodulin-dependent phosphodiesterase PDE1C down-regulates glucose-induced insulin secretion
Journal of Biological Chemistry 274:22337–22344.

https://doi.org/10.1074/jbc.274.32.22337
- PubMed
- Google Scholar
(2015) Self-assembly and plasticity of synaptic domains through a reaction-diffusion mechanism
Physical Review E 92:032705.

https://doi.org/10.1103/PhysRevE.92.032705
- Google Scholar
1. Hinke SA
2. Navedo MF
3. Ulman A
4. Whiting JL
5. Nygren PJ
6. Tian G
7. Jimenez-Caliani AJ
8. Langeberg LK
9. Cirulli V
10. Tengholm A
11. Dell'Acqua ML
12. Santana LF
13. Scott JD
(2012) Anchored phosphatases modulate glucose homeostasis
The EMBO Journal 31:3991–4004.

https://doi.org/10.1038/emboj.2012.244
- PubMed
- Google Scholar
1. Hoops S
2. Sahle S
3. Gauges R
4. Lee C
5. Pahle J
6. Simus N
7. Singhal M
8. Xu L
9. Mendes P
10. Kummer U
(2006) COPASI--a COmplex PAthway SImulator
Bioinformatics 22:3067–3074.

https://doi.org/10.1093/bioinformatics/btl485
- PubMed
- Google Scholar
(2005) Distinct enzyme combinations in AKAP signalling complexes permit functional diversity
Nature Cell Biology 7:1066–1073.

https://doi.org/10.1038/ncb1315
- PubMed
- Google Scholar
(2015) Modulation of calmodulin lobes by different targets: an allosteric model with hemiconcerted conformational transitions
PLOS Computational Biology 11:e1004063.

https://doi.org/10.1371/journal.pcbi.1004063
- PubMed
- Google Scholar
1. Landa LR
2. Harbeck M
3. Kaihara K
4. Chepurny O
5. Kitiphongspattana K
6. Graf O
7. Nikolaev VO
8. Lohse MJ
9. Holz GG
10. Roe MW
(2005) Interplay of Ca²⁺ and cAMP signaling in the insulin-secreting MIN6 beta-cell line
The Journal of Biological Chemistry 280:31294–31302.

https://doi.org/10.1074/jbc.M505657200
- PubMed
- Google Scholar
(2000) Scaffold proteins may biphasically affect the levels of mitogen-activated protein kinase signaling and reduce its threshold properties
PNAS 97:5818–5823.

https://doi.org/10.1073/pnas.97.11.5818
- PubMed
- Google Scholar
1. Lohse C
2. Bock A
3. Maiellaro I
4. Hannawacker A
5. Schad LR
6. Lohse MJ
7. Bauer WR
(2017) Experimental and mathematical analysis of cAMP nanodomains
PLOS ONE 12:e0174856.

https://doi.org/10.1371/journal.pone.0174856
- PubMed
- Google Scholar
(2009) Distinct mechanisms of regulation by ^Ca2+/calmodulin of type 1 and 8 adenylyl cyclases support their different physiological roles
Journal of Biological Chemistry 284:4451–4463.

https://doi.org/10.1074/jbc.M807359200
- PubMed
- Google Scholar
1. Masada N
2. Schaks S
3. Jackson SE
4. Sinz A
5. Cooper DM
(2012) Distinct mechanisms of calmodulin binding and regulation of adenylyl cyclases 1 and 8
Biochemistry 51:7917–7929.

https://doi.org/10.1021/bi300646y
- PubMed
- Google Scholar
1. Mo GC
2. Ross B
3. Hertel F
4. Manna P
5. Yang X
6. Greenwald E
7. Booth C
8. Plummer AM
9. Tenner B
10. Chen Z
11. Wang Y
12. Kennedy EJ
13. Cole PA
14. Fleming KG
15. Palmer A
16. Jimenez R
17. Xiao J
18. Dedecker P
19. Zhang J
(2017) Genetically encoded biosensors for visualizing live-cell biochemical activity at super-resolution
Nature Methods 14:427–434.

https://doi.org/10.1038/nmeth.4221
- PubMed
- Google Scholar
(2014) AKAP-anchored PKA maintains neuronal L-type calcium channel activity and NFAT transcriptional signaling
Cell Reports 7:1577–1588.

https://doi.org/10.1016/j.celrep.2014.04.027
- PubMed
- Google Scholar
(2018) cAMP: from Long-Range second messenger to nanodomain signalling
Trends in Pharmacological Sciences 39:209–222.

https://doi.org/10.1016/j.tips.2017.11.006
- PubMed
- Google Scholar
1. Nesher R
2. Anteby E
3. Yedovizky M
4. Warwar N
5. Kaiser N
6. Cerasi E
(2002) Beta-cell protein kinases and the dynamics of the insulin response to glucose
Diabetes 51 Suppl 1:S68–S73.

https://doi.org/10.2337/diabetes.51.2007.S68
- PubMed
- Google Scholar
1. Ni Q
2. Ganesan A
3. Aye-Han NN
4. Gao X
5. Allen MD
6. Levchenko A
7. Zhang J
(2011) Signaling diversity of PKA achieved via a Ca2+-cAMP-PKA oscillatory circuit
Nature Chemical Biology 7:34–40.

https://doi.org/10.1038/nchembio.478
- PubMed
- Google Scholar
(2015) Quantitative localization microscopy: effects of photophysics and labeling stoichiometry
PLOS ONE 10:e0127989.

https://doi.org/10.1371/journal.pone.0127989
- PubMed
- Google Scholar
1. Ohta Y
2. Furuta T
3. Nagai T
4. Horikawa K
(2018) Red fluorescent cAMP Indicator with increased affinity and expanded dynamic range
Scientific Reports 8:1866.

https://doi.org/10.1038/s41598-018-20251-1
- PubMed
- Google Scholar
1. Parekh AB
(2011) Decoding cytosolic Ca2+ oscillations
Trends in Biochemical Sciences 36:78–87.

https://doi.org/10.1016/j.tibs.2010.07.013
- PubMed
- Google Scholar
(2015) Modeling of glucose-induced cAMP oscillations in pancreatic β cells: camp rocks when metabolism rolls
Biophysical Journal 109:439–449.

https://doi.org/10.1016/j.bpj.2015.06.024
- PubMed
- Google Scholar
1. Petersen OH
(2002) Calcium signal compartmentalization
Biological Research 35:177–182.

https://doi.org/10.4067/S0716-97602002000200008
- PubMed
- Google Scholar
(2004) Measurement of dynamic protein binding to chromatin in vivo, using photobleaching microscopy
Methods in Enzymology 375:393–414.

https://doi.org/10.1016/s0076-6879(03)75025-3
- PubMed
- Google Scholar
1. Purkey AM
2. Woolfrey KM
3. Crosby KC
4. Stich DG
5. Chick WS
6. Aoto J
7. Dell'Acqua ML
(2018) AKAP150 palmitoylation regulates synaptic incorporation of Ca²⁺-Permeable AMPA Receptors to Control LTP
Cell Reports 25:974–987.

https://doi.org/10.1016/j.celrep.2018.09.085
- PubMed
- Google Scholar
1. Purvis JE
2. Lahav G
(2013) Encoding and decoding cellular information through signaling dynamics
Cell 152:945–956.

https://doi.org/10.1016/j.cell.2013.02.005
- PubMed
- Google Scholar
1. Raoux M
2. Vacher P
3. Papin J
4. Picard A
5. Kostrzewa E
6. Devin A
7. Gaitan J
8. Limon I
9. Kas MJ
10. Magnan C
11. Lang J
(2015) Multilevel control of glucose homeostasis by adenylyl cyclase 8
Diabetologia 58:749–757.

https://doi.org/10.1007/s00125-014-3445-z
- PubMed
- Google Scholar
(1997) Protein kinase A-dependent and -independent stimulation of exocytosis by cAMP in mouse pancreatic B-cells
The Journal of Physiology 502 ( Pt 1:105–118.

https://doi.org/10.1111/j.1469-7793.1997.105bl.x
- PubMed
- Google Scholar
(2009) Regulation of postsynaptic structure and function by an A-kinase anchoring protein-membrane-associated guanylate kinase scaffolding complex
Journal of Neuroscience 29:7929–7943.

https://doi.org/10.1523/JNEUROSCI.6093-08.2009
- PubMed
- Google Scholar
1. Rorsman P
2. Ashcroft FM
(2018) Pancreatic β-Cell electrical activity and insulin secretion: of mice and men
Physiological Reviews 98:117–214.

https://doi.org/10.1152/physrev.00008.2017
- PubMed
- Google Scholar
1. Sassone-Corsi P
(2012) The cyclic AMP pathway
Cold Spring Harbor Perspectives in Biology 4:a011148.

https://doi.org/10.1101/cshperspect.a011148
- PubMed
- Google Scholar
(2014) Mechanisms of cyclic AMP compartmentation revealed by computational models
Journal of General Physiology 143:39–48.

https://doi.org/10.1085/jgp.201311044
- PubMed
- Google Scholar
1. Schaff J
2. Fink CC
3. Slepchenko B
4. Carson JH
5. Loew LM
(1997) A general computational framework for modeling cellular structure and function
Biophysical Journal 73:1135–1146.

https://doi.org/10.1016/S0006-3495(97)78146-3
- PubMed
- Google Scholar
1. Schmitz O
2. Brock B
3. Hollingdal M
4. Juhl CB
5. Pørksen N
(2002)
High frequency insulin pulsatility and type 2 diabetes: from physiology and pathophysiology to clinical pharmacology

Diabetes & Metabolism 28:4S14–420.
- PubMed
- Google Scholar
1. Smith FD
2. Esseltine JL
3. Nygren PJ
4. Veesler D
5. Byrne DP
6. Vonderach M
7. Strashnov I
8. Eyers CE
9. Eyers PA
10. Langeberg LK
11. Scott JD
(2017) Local protein kinase A action proceeds through intact holoenzymes
Science 356:1288–1293.

https://doi.org/10.1126/science.aaj1669
- Google Scholar
1. Smith FD
2. Scott JD
(2002) Signaling complexes: junctions on the intracellular information super highway
Current Biology 12:R32–R40.

https://doi.org/10.1016/S0960-9822(01)00646-7
- PubMed
- Google Scholar
1. Stern MD
(1992) Buffering of calcium in the vicinity of a channel pore
Cell Calcium 13:183–192.

https://doi.org/10.1016/0143-4160(92)90046-U
- PubMed
- Google Scholar
1. Tajada S
2. Moreno CM
3. O'Dwyer S
4. Woods S
5. Sato D
6. Navedo MF
7. Santana LF
(2017) Distance constraints on activation of TRPV4 channels by AKAP150-bound pkcα in arterial myocytes
Journal of General Physiology 149:639–659.

https://doi.org/10.1085/jgp.201611709
- PubMed
- Google Scholar
1. Tengholm A
(2012) Cyclic AMP dynamics in the pancreatic β-cell
Upsala Journal of Medical Sciences 117:355–369.

https://doi.org/10.3109/03009734.2012.724732
- PubMed
- Google Scholar
1. Thompson JL
2. Shuttleworth TJ
(2015) Anchoring protein AKAP79-mediated PKA phosphorylation of STIM1 determines selective activation of the ARC channel, a store-independent orai channel
The Journal of Physiology 593:559–572.

https://doi.org/10.1113/jphysiol.2014.284182
- PubMed
- Google Scholar
(2017) The many faces of compartmentalized PKA signalosomes
Cellular Signalling 37:1–11.

https://doi.org/10.1016/j.cellsig.2017.05.012
- PubMed
- Google Scholar
1. Wachten S
2. Masada N
3. Ayling LJ
4. Ciruela A
5. Nikolaev VO
6. Lohse MJ
7. Cooper DM
(2010) Distinct pools of cAMP centre on different isoforms of adenylyl cyclase in pituitary-derived GH3B6 cells
Journal of Cell Science 123:95–106.

https://doi.org/10.1242/jcs.058594
- PubMed
- Google Scholar
1. White MA
2. Anderson RG
(2005) Signaling networks in living cells
Annual Review of Pharmacology and Toxicology 45:587–603.

https://doi.org/10.1146/annurev.pharmtox.45.120403.095807
- PubMed
- Google Scholar
1. Willoughby D
2. Masada N
3. Wachten S
4. Pagano M
5. Halls ML
6. Everett KL
7. Ciruela A
8. Cooper DM
(2010) AKAP79/150 interacts with AC8 and regulates ^Ca2+-dependent cAMP synthesis in pancreatic and neuronal systems
Journal of Biological Chemistry 285:20328–20342.

https://doi.org/10.1074/jbc.M110.120725
- PubMed
- Google Scholar
(2012) A key phosphorylation site in AC8 mediates regulation of ca(2+)-dependent cAMP dynamics by an AC8-AKAP79-PKA signalling complex
Journal of Cell Science 125:5850–5859.

https://doi.org/10.1242/jcs.111427
- PubMed
- Google Scholar
1. Willoughby D
2. Cooper DM
(2006) Ca2+ stimulation of adenylyl cyclase generates dynamic oscillations in cyclic AMP
Journal of Cell Science 119:828–836.

https://doi.org/10.1242/jcs.02812
- PubMed
- Google Scholar
1. Yang PC
2. Boras BW
3. Jeng MT
4. Docken SS
5. Lewis TJ
6. McCulloch AD
7. Harvey RD
8. Clancy CE
(2016) A computational modeling and simulation approach to investigate mechanisms of subcellular cAMP compartmentation
PLOS Computational Biology 12:e1005005.

https://doi.org/10.1371/journal.pcbi.1005005
- PubMed
- Google Scholar
(2017) A DNA origami platform for quantifying protein copy number in super-resolution
Nature Methods 14:789–792.

https://doi.org/10.1038/nmeth.4342
- PubMed
- Google Scholar
1. Zhang J
2. Hupfeld CJ
3. Taylor SS
4. Olefsky JM
5. Tsien RY
(2005) Insulin disrupts β-adrenergic signalling to protein kinase A in adipocytes
Nature 437:569–573.

https://doi.org/10.1038/nature04140
- Google Scholar
1. Zhang J
2. Bal M
3. Bierbower S
4. Zaika O
5. Shapiro MS
(2011) AKAP79/150 signal complexes in G-protein modulation of neuronal ion channels
Journal of Neuroscience 31:7199–7211.

https://doi.org/10.1523/JNEUROSCI.4446-10.2011
- PubMed
- Google Scholar
1. Zhang M
2. Patriarchi T
3. Stein IS
4. Qian H
5. Matt L
6. Nguyen M
7. Xiang YK
8. Hell JW
(2013) Adenylyl cyclase anchoring by a kinase anchor protein AKAP5 (AKAP79/150) Is important for postsynaptic β-Adrenergic signaling
Journal of Biological Chemistry 288:17918–17931.

https://doi.org/10.1074/jbc.M112.449462
- Google Scholar
(2016) Clustering and functional coupling of diverse ion channels and signaling proteins revealed by Super-resolution STORM microscopy in neurons
Neuron 92:461–478.

https://doi.org/10.1016/j.neuron.2016.09.014
- PubMed
- Google Scholar
1. Zhang X
2. Pathak T
3. Yoast R
4. Emrich S
5. Xin P
6. Nwokonko RM
7. Johnson M
8. Wu S
9. Delierneux C
10. Gueguinou M
11. Hempel N
12. Putney JW
13. Gill DL
14. Trebak M
(2019) A calcium/cAMP signaling loop at the ORAI1 mouth drives channel inactivation to shape NFAT induction
Nature Communications 10:1971.

https://doi.org/10.1038/s41467-019-09593-0
- PubMed
- Google Scholar
1. Zhang JZ
2. Lu TW
3. Stolerman LM
4. Tenner B
5. Yang JR
6. Zhang JF
7. Falcke M
8. Rangamani P
9. Taylor SS
10. Mehta S
11. Zhang J
(2020) Phase separation of a PKA regulatory subunit controls cAMP compartmentation and oncogenic signaling
Cell 182:1531–1544.

https://doi.org/10.1016/j.cell.2020.07.043
- PubMed
- Google Scholar

Article and author information

Author details

Brian Tenner
1. Department of Biophysics and Biophysical Chemistry, The Johns Hopkins University School of Medicine, Baltimore, United States
2. Department of Pharmacology, University of California, San Diego, La Jolla, United States
Contribution
Conceptualization, Resources, Data curation, Formal analysis, Supervision, Validation, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

Competing interests
No competing interests declared
Michael Getz

Chemical Engineering Graduate Program, University of California, San Diego, La Jolla, United States

Contribution
Software, Formal analysis, Validation, Visualization, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9886-8454
Brian Ross

Department of Pharmacology, University of California, San Diego, La Jolla, United States

Contribution
Data curation, Software, Investigation, Visualization, Methodology, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9020-627X
Donya Ohadi

Department of Mechanical and Aerospace Engineering, University of California, San Diego, La Jolla, United States

Contribution
Data curation, Software, Formal analysis, Validation, Writing - review and editing

Competing interests
No competing interests declared
Christopher H Bohrer

Department of Biophysics and Biophysical Chemistry, The Johns Hopkins University School of Medicine, Baltimore, United States

Contribution
Data curation, Software, Formal analysis, Methodology, Writing - review and editing

Competing interests
No competing interests declared
Eric Greenwald

Department of Pharmacology, University of California, San Diego, La Jolla, United States

Contribution
Conceptualization, Data curation, Software, Formal analysis, Writing - review and editing

Competing interests
No competing interests declared
Sohum Mehta

Department of Pharmacology, University of California, San Diego, La Jolla, United States

Contribution
Conceptualization, Resources, Supervision, Investigation, Visualization, Project administration, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4764-8579
Jie Xiao

Department of Biophysics and Biophysical Chemistry, The Johns Hopkins University School of Medicine, Baltimore, United States

Contribution
Project administration, Writing - review and editing

Competing interests
Reviewing editor, eLife

"This ORCID iD identifies the author of this article:" 0000-0003-1433-5437
Padmini Rangamani
1. Chemical Engineering Graduate Program, University of California, San Diego, La Jolla, United States
2. Department of Mechanical and Aerospace Engineering, University of California, San Diego, La Jolla, United States
Contribution
Conceptualization, Supervision, Investigation, Project administration, Writing - review and editing

For correspondence
padmini.rangamani@eng.ucsd.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5953-4347
Jin Zhang
1. Department of Pharmacology, University of California, San Diego, La Jolla, United States
2. Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, United States
Contribution
Conceptualization, Supervision, Funding acquisition, Project administration, Writing - review and editing

For correspondence
jzhang32@ucsd.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7145-7823

Funding

National Institutes of Health (R01 DK073368)

Jin Zhang

U.S. Department of Defense (DOD AFOSR FA9550-18-1-0051)

Padmini Rangamani
Jin Zhang

Office of Naval Research (ONR N00014-17-1-2628)

Padmini Rangamani

National Institutes of Health (R35 CA197622)

Jin Zhang

National Institutes of Health (5T32 GM007231)

Christopher H Bohrer

National Institutes of Health (R01 GM112008)

Jie Xiao

National Institutes of Health (R35 GM127075)

Jie Xiao

Johns Hopkins University (Johns Hopkins Discovery Award)

Jie Xiao
Jin Zhang

Johns Hopkins University (Hamilton Innovation Research Award)

Jie Xiao

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. John D Scott for providing AKAP79, Dr. Dermot Cooper for providing (Ci/Ce)Epac2-camps and AC8, Dr. Loren Looger for providing RCaMP, Dr. Kazuki Horikawa for providing R-FlincA and Dr. Jochen Lang for providing the shAC8 construct. We thank Dr. Susan Taylor for critical reading of the manuscript. We also thank the members of the Nikon Imaging Center at UCSD for help with fluorescence microscopy experiments, and the Neuroscience Imaging Center at the Moore’s Cancer Center for help with the STORM acquisition and analysis. This work was supported by NIH R01 DK073368 and R35 CA197622 (to JZ), DOD AFOSR FA9550-18-1-0051 (to PR and JZ), ONR N00014-17-1-2628 (to PR), NIH R01 GM112008 and R35 GM127075 (to JX), NIH T32GM007231 (to CHB), Johns Hopkins Discovery Award (to JX and JZ) and Hamilton Innovation Rsearch Award (to JX).

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.