Roles of GRK and PDE4 activities in the regulation of beta2 adrenergic signaling

Computational modelling of dynamic cAMP responses to GPCR agonists for exploration of GLP-1R ligand effects in pancreatic β-cells and neurons

The glucagon-like peptide-1 receptor (GLP-1R) is a class B G protein-coupled receptor (GPCR) which plays important physiological roles in insulin release and promoting fullness. GLP-1R agonists initiate cellular responses by cyclic AMP (cAMP) pathway signal transduction. Understanding of the potential of GLP-1R agonists in the treatment of type 2 diabetes may be advanced by considering the cAMP dynamics for agonists at GLP-1R in both pancreatic β-cells (important in insulin release) and neurons (important in appetite regulation). Receptor desensitisation in the cAMP pathway is known to be an important regulatory mechanism, with different ligands differentially promoting G protein activation and desensitisation. Here, we use mathematical modelling to quantify and understand experimentally obtained cAMP timecourses for two GLP-1R agonists, exendin-F1 (ExF1) and exendin-D3 (ExD3), which give markedly different signals in β-cells and neurons. We formulate an ordinary differential equation (ODE) model for the dynamics of cAMP signalling in response to G protein-coupled receptor (GPCR) ligands, encompassing ligand binding, receptor activation, G protein activation, desensitisation and second messenger generation. We validate our model initially by fitting to timecourse data for HEK293 cells, then proceed to parameterise the model for β-cells and neurons. Through numerical simulation and sensitivity studies, our analysis adds support to the hypothesis that ExF1 offers more potential glucose regulation benefit than ExD3 over long timescales via signalling in pancreatic β-cells, but that there is little difference between the two ligands in the potential appetite suppression effects offered via long-time signalling in neurons on the same timescales.

Computational investigation of the dynamic control of cAMP signaling by PDE4 isoform types

Cyclic adenosine monophosphate (cAMP) is a generic signaling molecule that, through precise control of its signaling dynamics, exerts distinct cellular effects. Consequently, aberrant cAMP signaling can have detrimental effects. Phosphodiesterase 4 (PDE4) enzymes profoundly control cAMP signaling and comprise different isoform types of which the enzymatic activity is modulated by differential feedback mechanisms. Because these feedback dynamics are non-linear and occur coincidentally, their effects are difficult to examine experimentally, but can be well simulated computationally. Through understanding the role of PDE4 isoform types in regulating cAMP signaling, PDE4-targeted therapeutic strategies can be better specified. Here, we established a computational model to study how feedback mechanisms on different PDE4 isoform types lead to dynamic, isoform-specific control of cAMP signaling. Ordinary differential equations describing cAMP dynamics were implemented in the VirtualCell (VCell) environment. Simulations indicated that long PDE4 isoforms exert the most profound control on oscillatory cAMP signaling, as opposed to the PDE4-mediated control of single cAMP input pulses. Moreover, elevating cAMP levels or decreasing PDE4 levels revealed different effects on downstream signaling. Together these results underline that cAMP signaling is distinctly regulated by different PDE4 isoform types and that this isoform-specificity should be considered in both computational and experimental follow-up studies to better define PDE4 enzymes as therapeutic targets in diseases in which cAMP signaling is aberrant.

PDE4DIP in health and diseases

Cyclic-AMP (cAMP), the first second messenger to be identified, is synthesized, and is universally utilized as a second messenger, and plays important roles in integrity, and function of organs, including heart. Through its coupling with other intracellular messengers, cAMP facilitates excitation-contraction coupling, increases heart rate and conduction velocity. It is degraded by a class of enzymes called cAMP-dependent phosphodiesterase (PDE), with PDE3 and PDE4 being the predominant isoforms in the heart. This highly diverse class of enzymes degrade cAMP and through anchoring proteins generates dynamic microdomains to target specific proteins and control specific cell functions in response to various stimuli. The impaired function of the anchoring protein either by inherited genetic mutations or acquired injuries results in altered intracellular targeting, and blunted responsiveness to stimulating pathways and contributes to pathological cardiac remodeling, cardiac arrhythmias and reduced cell survival. Recent genetic studies provide compelling evidence for an association between the variants in the anchoring protein PDE4DIP and atrial fibrillation, stroke, and heart failure.

Ion Channel-Based Reporters for cAMP Detection

In the last 20 years tremendous progress has been made in the development of single cell cAMP sensors. Sensors are based upon cAMP binding proteins that have been modified to transduce cAMP concentrations into electrical or fluorescent readouts that can be readily detected using patch clamp amplifiers, photomultiplier tubes, or cameras. Here, we describe two complementary approaches for the detection and measurement of cAMP signals near the plasma membrane of cells using cyclic nucleotide (CNG) channel-based probes. These probes take advantage of the ability of CNG channels to transduce small changes in cAMP concentration into ionic flux through channel pores that can be readily detected by measuring Ca²⁺ and/or Mn²⁺ influx or by measuring ionic currents.

Identification of phosphodiesterase-4 as the therapeutic target of arctigenin in alleviating psoriatic skin inflammation

Introduction

Arctigenin, derived from Arctium lappa L., has multiple pharmacological activities, including immunoregulatory, anti-diabetic, anti-tumor, and neuroprotective effects. Nevertheless, the potential therapeutic target of arctigenin in modulating inflammation remains undefined.

Objectives

In the present study, we identified that arctigenin was a phosphodiesterase-4 (PDE4) selective inhibitor for the first time. Further investigations were performed to fully uncover the effects and mechanism of arctigenin on experimental murine psoriasis model.

Methods

Crystal structure determination, PDEs enzyme assay, and isothermal titration calorimetry were included to illustrate the binding specialty, inhibitory effects, and selectivity of arctigenin on PDE4D. The anti-inflammatory effects were conducted in LPS-activated human peripheral blood mononuclear cells (PBMCs) and RAW264.7 cells. Imiquimod-induced murine psoriasis was performed to uncover the therapeutic effects and mechanism of arctigenin in vivo.

Results

Arctigenin could bind to the catalytic domain of PDE4D via formation of hydrogen bonds as well as π-π stacking interactions between the dibenzyl butyrolactone of arctigenin and several residues of PDE4D. Accordingly, arctigenin showed prominent anti-inflammation in human PBMCs and murine RAW264.7 cells. PDE4 inhibition by arctigenin resulted in elevation of intracellular cyclic adenosine monophosphate (cAMP) and phosphorylation of cAMP-response element binding protein (CREB), which were largely blocked through intervention of protein kinase A (PKA) activity by H89 treatment or reduction of protein expression by siRNA transfection. Moreover, we first identified that a topical application of arctigenin ameliorated experimental psoriatic manifestations in imiquimod-induced murine psoriasis model by decreasing adhesion and chemotaxis of several inflammatory cells. Further proteomics analysis revealed that arctigenin could rectify the immune dysfunction and hyperactivation of keratinocytes in the inflamed skin microenvironments, which might be largely related to the expression of Keratins.

Conclusion

The research provided credible clew that inhibition of PDE4 by arctigenin might function as the potential therapeutic approach for the treatment of psoriasis.

Phosphodiesterase 4 mediates interleukin-8-induced heterologous desensitization of the β2 -adrenergic receptor

Acute respiratory distress syndrome (ARDS) is a life-threatening illness characterized by decreased alveolar-capillary barrier function, pulmonary edema consisting of proteinaceous fluid, and inhibition of net alveolar fluid transport responsible for resolution of pulmonary edema. There is currently no pharmacotherapy that has proven useful to prevent or treat ARDS, and two trials using beta-agonist therapy to treat ARDS demonstrated no effect. Prior studies indicated that IL-8-induced heterologous desensitization of the beta2-adrenergic receptor (β₂ -AR) led to decreased beta-agonist-induced mobilization of cyclic adenosine monophosphate (cAMP). Interestingly, phosphodiesterase (PDE) 4 inhibitors have been used in human airway diseases characterized by low intracellular cAMP levels and increases in specific cAMP hydrolyzing activity. Therefore, we hypothesized that PDE4 would mediate IL-8-induced heterologous internalization of the β₂ -AR and that PDE4 inhibition would restore beta-agonist-induced functions. We determined that CINC-1 (a functional IL-8 analog in rats) induces internalization of β₂ -AR from the cell surface, and arrestin-2, PDE4, and β₂ -AR form a complex during this process. Furthermore, we determined that cAMP associated with the plasma membrane was adversely affected by β₂ -AR heterologous desensitization. Additionally, we determined that rolipram, a PDE4 inhibitor, reversed CINC-1-induced derangements of cAMP and also caused β₂ -AR to successfully recycle back to the cell surface. Finally, we demonstrated that rolipram could reverse CINC-1-mediated inhibition of beta-agonist-induced alveolar fluid clearance in a murine model of trauma-shock. These results indicate that PDE4 plays a role in CINC-1-induced heterologous internalization of the β₂ -AR; PDE4 inhibition reverses these effects and may be a useful adjunct in particular ARDS patients.

Epistatic interaction of PDE4DIP and DES mutations in familial atrial fibrillation with slow conduction

The genetic causes of atrial fibrillation (AF) with slow conduction are unknown. Eight kindreds with familial AF and slow conduction, including a family affected by early-onset AF, heart block, and incompletely penetrant nonischemic dilated cardiomyopathy (DCM) underwent whole exome sequencing. A known pathogenic mutation in the desmin (DES) gene resulting in p.S13F substitution (NM_001927.3:c.38C>T) at a PKC phosphorylation site was identified in all four members of the kindred with early-onset AF and heart block, while only two developed DCM. Higher penetrance for AF and heart block prompted a genetic screening for DES modifier(s). A deleterious mutation in the phosphodiesterase-4D-interacting-protein (PDE4DIP) gene resulting in p.A123T substitution (NM_001002811:c.367G>A) was identified that segregated with early-onset AF, heart block, and the DES mutation. Three additional novel deleterious PDE4DIP mutations were identified in four other unrelated kindreds. Characterization of PDE4DIP^A123T in vitro suggested impaired compartmentalization of PKA and PDE4D characterized by reduced colocalization with PDE4D, increased cAMP activation leading to higher PKA phosphorylation of the β2-adrenergic-receptor, and decreased PKA phosphorylation of desmin after isoproterenol stimulation. Our findings identify PDE4DIP as a novel gene for slow AF and unravel its epistatic interaction with DES mutations in development of conduction disease and arrhythmia.

Analyzing kinetic signaling data for G-protein-coupled receptors

In classical pharmacology, bioassay data are fit to general equations (e.g. the dose response equation) to determine empirical drug parameters (e.g. EC50 and Emax), which are then used to calculate chemical parameters such as affinity and efficacy. Here we used a similar approach for kinetic, time course signaling data, to allow empirical and chemical definition of signaling by G-protein-coupled receptors in kinetic terms. Experimental data are analyzed using general time course equations (model-free approach) and mechanistic model equations (mechanistic approach) in the commonly-used curve-fitting program, GraphPad Prism. A literature survey indicated signaling time course data usually conform to one of four curve shapes: the straight line, association exponential curve, rise-and-fall to zero curve, and rise-and-fall to steady-state curve. In the model-free approach, the initial rate of signaling is quantified and this is done by curve-fitting to the whole time course, avoiding the need to select the linear part of the curve. It is shown that the four shapes are consistent with a mechanistic model of signaling, based on enzyme kinetics, with the shape defined by the regulation of signaling mechanisms (e.g. receptor desensitization, signal degradation). Signaling efficacy is the initial rate of signaling by agonist-occupied receptor (kτ), simply the rate of signal generation before it becomes affected by regulation mechanisms, measurable using the model-free analysis. Regulation of signaling parameters such as the receptor desensitization rate constant can be estimated if the mechanism is known. This study extends the empirical and mechanistic approach used in classical pharmacology to kinetic signaling data, facilitating optimization of new therapeutics in kinetic terms.

Control of βAR- and N-methyl-D-aspartate (NMDA) Receptor-Dependent cAMP Dynamics in Hippocampal Neurons

Norepinephrine, a neuromodulator that activates β-adrenergic receptors (βARs), facilitates learning and memory as well as the induction of synaptic plasticity in the hippocampus. Several forms of long-term potentiation (LTP) at the Schaffer collateral CA1 synapse require stimulation of both βARs and N-methyl-D-aspartate receptors (NMDARs). To understand the mechanisms mediating the interactions between βAR and NMDAR signaling pathways, we combined FRET imaging of cAMP in hippocampal neuron cultures with spatial mechanistic modeling of signaling pathways in the CA1 pyramidal neuron. Previous work implied that cAMP is synergistically produced in the presence of the βAR agonist isoproterenol and intracellular calcium. In contrast, we show that when application of isoproterenol precedes application of NMDA by several minutes, as is typical of βAR-facilitated LTP experiments, the average amplitude of the cAMP response to NMDA is attenuated compared with the response to NMDA alone. Models simulations suggest that, although the negative feedback loop formed by cAMP, cAMP-dependent protein kinase (PKA), and type 4 phosphodiesterase may be involved in attenuating the cAMP response to NMDA, it is insufficient to explain the range of experimental observations. Instead, attenuation of the cAMP response requires mechanisms upstream of adenylyl cyclase. Our model demonstrates that Gs-to-Gi switching due to PKA phosphorylation of βARs as well as Gi inhibition of type 1 adenylyl cyclase may underlie the experimental observations. This suggests that signaling by β-adrenergic receptors depends on temporal pattern of stimulation, and that switching may represent a novel mechanism for recruiting kinases involved in synaptic plasticity and memory.

Noradrenaline is a stress related molecule that facilitates learning and memory when released in the hippocampus. The facilitation of memory is related to modulation of synaptic plasticity, but the mechanisms underlying this modulation are not well understood. We utilize a combination of live cell imaging and computational modeling to discover how noradrenergic receptor stimulation interacts with other molecules, such as calcium, required for synaptic plasticity and memory storage. Though prior work has shown that noradrenergic receptors and calcium interact synergistically to elevate intracellular second messengers when combined simultaneously, our results demonstrate that prior stimulation of noradrenergic receptors inhibits the elevation of intracellular second messengers. Our results further demonstrate that the inhibition may be caused by the noradrenergic receptor switching signaling pathways, thereby recruiting a different set of memory kinases. This switching represents a novel mechanism for recruiting molecules involved in synaptic plasticity and memory.

Channel-based reporters for cAMP detection

In the last 15 years, tremendous progress has been made in the development of single-cell cAMP sensors. Sensors are based upon cAMP-binding proteins that have been modified to transduce cAMP concentrations into electrical or fluorescent readouts that can be readily detected using patch clamp amplifiers, photomultiplier tubes, or cameras. Here we describe two complementary approaches for the detection and measurement of cAMP signals near the plasma membrane of cells. These probes take advantage of the ability of cyclic nucleotide-gated (CNG) channels to transduce small changes in cAMP concentrations into ionic flux through channel pores that can be readily detected by measuring Ca(2+) and/or Mn(2+) influx or by measuring ionic currents.

Heterogeneity of pulmonary endothelial cyclic nucleotide response to Pseudomonas aeruginosa ExoY infection

Here, we tested the hypothesis that a promiscuous bacterial cyclase synthesizes purine and pyrimidine cyclic nucleotides in the pulmonary endothelium. To test this hypothesis, pulmonary endothelial cells were infected with a strain of the Gram-negative bacterium Pseudomonas aeruginosa that introduces only exoenzyme Y (PA103 ΔexoUexoT::Tc pUCPexoY; ExoY(+)) via a type III secretion system. Purine and pyrimidine cyclic nucleotides were simultaneously detected using mass spectrometry. Pulmonary artery (PAECs) and pulmonary microvascular (PMVECs) endothelial cells both possess basal levels of four different cyclic nucleotides in the following rank order: cAMP > cUMP ≈ cGMP ≈ cCMP. Endothelial gap formation was induced in a time-dependent manner following ExoY(+) intoxication. In PAECs, intercellular gaps formed within 2 h and progressively increased in size up to 6 h, when the experiment was terminated. cGMP concentrations increased within 1 h postinfection, whereas cAMP and cUMP concentrations increased within 3 h, and cCMP concentrations increased within 4 h postinfection. In PMVECs, intercellular gaps did not form until 4 h postinfection. Only cGMP and cUMP concentrations increased at 3 and 6 h postinfection, respectively. PAECs generated higher cyclic nucleotide levels than PMVECs, and the cyclic nucleotide levels increased earlier in response to ExoY(+) intoxication. Heterogeneity of the cyclic nucleotide signature in response to P. aeruginosa infection exists between PAECs and PMVECs, suggesting the intracellular milieu in PAECs is more conducive to cNMP generation.

The genetically encoded tool set for investigating cAMP: more than the sum of its parts

Intracellular fluctuations of the second messenger cyclic AMP (cAMP) are regulated with spatial and temporal precision. This regulation is supported by the sophisticated arrangement of cyclases, phosphodiesterases, anchoring proteins, and receptors for cAMP. Discovery of these nuances to cAMP signaling has been facilitated by the development of genetically encodable tools for monitoring and manipulating cAMP and the proteins that support cAMP signaling. In this review, we discuss the state-of-the-art in development of different genetically encoded tools for sensing cAMP and the activity of its primary intracellular receptor protein kinase A (PKA). We introduce sequences for encoding adenylyl cyclases that enable cAMP levels to be artificially elevated within cells. We chart the evolution of sequences for selectively modifying protein-protein interactions that support cAMP signaling, and for driving cAMP sensors and manipulators to different subcellular locations. Importantly, these different genetically encoded tools can be applied synergistically, and we highlight notable instances that take advantage of this property. Finally, we consider prospects for extending the utility of the tool set to support further insights into the role of cAMP in health and disease.

Estimating the magnitude of near-membrane PDE4 activity in living cells

Recent studies have demonstrated that functionally discrete pools of phosphodiesterase (PDE) activity regulate distinct cellular functions. While the importance of localized pools of enzyme activity has become apparent, few studies have estimated enzyme activity within discrete subcellular compartments. Here we present an approach to estimate near-membrane PDE activity. First, total PDE activity is measured using traditional PDE activity assays. Second, known cAMP concentrations are dialyzed into single cells and the spatial spread of cAMP is monitored using cyclic nucleotide-gated channels. Third, mathematical models are used to estimate the spatial distribution of PDE activity within cells. Using this three-tiered approach, we observed two pharmacologically distinct pools of PDE activity, a rolipram-sensitive pool and an 8-methoxymethyl IBMX (8MM-IBMX)-sensitive pool. We observed that the rolipram-sensitive PDE (PDE4) was primarily responsible for cAMP hydrolysis near the plasma membrane. Finally, we observed that PDE4 was capable of blunting cAMP levels near the plasma membrane even when 100 μM cAMP were introduced into the cell via a patch pipette. Two compartment models predict that PDE activity near the plasma membrane, near cyclic nucleotide-gated channels, was significantly lower than total cellular PDE activity and that a slow spatial spread of cAMP allowed PDE activity to effectively hydrolyze near-membrane cAMP. These results imply that cAMP levels near the plasma membrane are distinct from those in other subcellular compartments; PDE activity is not uniform within cells; and localized pools of AC and PDE activities are responsible for controlling cAMP levels within distinct subcellular compartments.

Salmeterol Efficacy and Bias in the Activation and Kinase-Mediated Desensitization of β2-Adrenergic Receptors

Salmeterol is a long-acting β2-adrenergic receptor (β2AR) agonist that is widely used as a bronchodilator for the treatment of persistent asthma and chronic obstructive pulmonary disease in conjunction with steroids. Previous studies demonstrated that salmeterol showed weak efficacy for activation of adenylyl cyclase; however, its efficacy in the complex desensitization of the β2AR remains poorly understood. In this work, we provide insights into the roles played by the G protein-coupled receptor kinase/arrestin and protein kinase A in salmeterol-mediated desensitization through bioluminescence resonance energy transfer (BRET) studies of liganded-β2AR binding to arrestin and through kinetic studies of cAMP turnover. First, BRET demonstrated a much reduced efficacy for salmeterol recruitment of arrestin to β2AR relative to isoproterenol. The ratio of BRETISO/BRETSALM after 5-minute stimulation was 20 and decreased to 5 after 35 minutes, reflecting a progressive decline in BRETISO and a stable BRETSALM. Second, to assess salmeterol efficacy for functional desensitization, we examined the kinetics of salmeterol-induced cAMP accumulation (0-30 minutes) in human airway smooth muscle cells in the presence and absence of phosphodiesterase inhibition. Analysis of shaping of cAMP turnover for both agonists demonstrated significant salmeterol desensitization, although it was reduced relative to isoproterenol. Using an isoproterenol rescue protocol after either short-term (10 minutes) or long-term (2 and 14 hours) salmeterol pretreatments, we found that salmeterol progressively depressed isoproterenol stimulation but did not prevent subsequent rescue by isoproterenol and additional isoproterenol-mediated desensitization. Our findings reveal a complex efficacy for functional desensitization, demonstrating that although salmeterol shows weak efficacy for adenylyl cyclase activation and G protein-coupled receptor kinase/arrestin-mediated desensitization, it acts as a strong agonist in highly amplified protein kinase A-mediated events.

Mathematical modeling improves EC50 estimations from classical dose-response curves

The β-adrenergic response is impaired in failing hearts. When studying β-adrenergic function in vitro, the half-maximal effective concentration (EC50 ) is an important measure of ligand response. We previously measured the in vitro contraction force response of chicken heart tissue to increasing concentrations of adrenaline, and observed a decreasing response at high concentrations. The classical interpretation of such data is to assume a maximal response before the decrease, and to fit a sigmoid curve to the remaining data to determine EC50 . Instead, we have applied a mathematical modeling approach to interpret the full dose-response curve in a new way. The developed model predicts a non-steady-state caused by a short resting time between increased concentrations of agonist, which affect the dose-response characterization. Therefore, an improved estimate of EC50 may be calculated using steady-state simulations of the model. The model-based estimation of EC50 is further refined using additional time-resolved data to decrease the uncertainty of the prediction. The resulting model-based EC50 (180-525 nm) is higher than the classically interpreted EC50 (46-191 nm). Mathematical modeling thus makes it possible to re-interpret previously obtained datasets, and to make accurate estimates of EC50 even when steady-state measurements are not experimentally feasible.

DATABASE

The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database, and may be accessed at http://jjj.bio.vu.nl/database/nyman.

Modulation of the cAMP response by Gαi and Gβγ: a computational study of G protein signaling in immune cells

Cyclic AMP is important for the resolution of inflammation, as it promotes anti-inflammatory signaling in several immune cell lines. In this paper, we present an immune cell specific model of the cAMP signaling cascade, paying close attention to the specific isoforms of adenylyl cyclase (AC) and phosphodiesterase that control cAMP production and degradation, respectively, in these cells. The model describes the role that G protein subunits, including Gαs, Gαi, and Gβγ, have in regulating cAMP production. Previously, Gαi activation has been shown to increase the level of cAMP in certain immune cell types. This increase in cAMP is thought to be mediated by βγ subunits which are released upon Gα activation and can directly stimulate specific isoforms of AC. We conduct numerical experiments in order to explore the mechanisms through which Gαi activation can increase cAMP production. An important conclusion of our analysis is that the relative abundance of different G protein subunits is an essential determinant of the cAMP profile in immune cells. In particular, our model predicts that limited availability of βγ subunits may both (i) enable immune cells to link inflammatory Gαi signaling to anti-inflammatory cAMP production thereby creating a balanced immune response to stimulation with low concentrations of PGE2, and (ii) prohibit robust anti-inflammatory cAMP signaling in response to stimulation with high concentrations of PGE2.

cAMP-PKA inhibition of SK3 channel reduced both Ca2+ entry and cancer cell migration by regulation of SK3-Orai1 complex

SK3 channel mediates the migration of various cancer cells. When expressed in breast cancer cells, SK3 channel forms a complex with Orai1, a voltage-independent Ca(2+) channel. This SK3-Orai1 complex associates within lipid rafts where it controls a constitutive Ca(2+) entry leading to cancer cell migration and bone metastases development. Since cAMP was found to modulate breast cancer cell migration, we hypothesized that this could be explained by a modulation of SK3 channel activity. Herein, we study the regulation of SK3 channel by the cAMP-PKA pathway and the consequences for SK3-dependent Ca(2+) entry and cancer cell migration. We established that the beta-adrenergic receptor agonist, isoprenaline, or the direct adenylyl cyclase activator forskolin alone or in combination with the PDE4 inhibitor, CI-1044, decreased SK3 channel activity without modifying the expression of SK3 protein at the plasma membrane. Forskolin and CI-1044 reduced the SK3-dependent constitutive Ca(2+) entry and the SK3-dependent migration of MDA-MB-435s cells. PKA inhibition with KT 5720 reduced: (1) the effect of forskolin and CI-1044 by 50 % on Ca(2+) entry and (2) SK3 activity by inhibiting the serine phosphorylation of SK3. These cAMP-elevating agents displaced Orai1 protein outside lipid rafts in contrast to SK3, which remained in the lipid rafts fractions. All together, these results show that activation of the cAMP-PKA pathway decreases SK3 channel and SK3-Orai1 complex activities, leading to a decrease in both Ca(2+) entry and cancer cell migration. This work supports the potential use of cAMP-elevating agents to reduce cancer cell migration and may provide novel opportunities to address/prevent bone metastasis.

Biochemical and Cellular Specificity of Peptide Inhibitors of G Protein-Coupled Receptor Kinases

Identifying novel allosteric inhibitors of G protein-coupled receptor kinases (GRKs) would be of considerable use in limiting both the extent of desensitization of GPCRs as well as downstream positive regulation through GRKs. Several peptides have previously been identified as inhibitors of specific GRKs, but to date there have been few comparisons of the selectivities of these materials on the seven GRKs, modifications to allow cell penetration, or off-target activities. The goal of this study was to determine if a panel of peptides mimicking domains on either GPCRs or GRKs would exhibit selective inhibition of GRKs 2, 5, 6 and 7 phosphorylation of rhodopsin. Peptides included sequences from GRK5; helices 3, 9, and 10 (α3, α9, and α10) in the RH domain, and the N-terminal peptide (N-Ter), as well as the intracellular loop 1 (iL1) of the β2-adrenergic receptor (β2AR), and the G_α transducin C-tail (TCT). While some selectivity for individual GRKs was found, overall selectivity was limited and often not reflective of structural predictions. Off-target effects were probed by determining peptide inhibition of adenylyl cyclase (AC) and PKA, and while peptides had no effect on AC activity, N-Ter, iL1, and α10 were potent inhibitors of PKA. To probe inhibition of GRK activity in intact cells, we synthesized TAT-tagged peptides, and found that TAT-α9-R169A and TAT-TCT inhibited isoproterenol-stimulated GRK phosphorylation of the β2AR; however, the TAT peptides also inhibited isoproterenol and forskolin stimulation of AC activity. Our findings demonstrate potent peptide inhibition of GRK activities in vitro, highlight the differences in the environments of biochemical and cell-based assays, and illustrate the care that must be exercised in interpreting results of either assay alone.

Protective effect of high-dose montelukast on salbutamol-induced homologous desensitisation in airway smooth muscle

Montelukast (MK) is a potent cysteinyl-leukotriene receptor antagonist that causes dose-related improvements in chronic asthma. We sought to determine whether MK was able to prevent salbutamol-induced tolerance in airway smooth muscle. Homologous β2-adrenoceptor desensitisation models were established in guinea-pigs and in human bronchial smooth muscle cells (BSMC) by chronic salbutamol administration. Characterisation tools included measurement of the response of tracheal smooth muscle tissues to salbutamol, analysis of gene expression and receptor trafficking, evaluation of intracellular cAMP levels and phosphodiesterase (PDE) activity in human bronchial smooth muscle cells. Salbutamol-induced β2-adrenoceptor desensitisation was characterised by β2-agonist hyporesponsiveness (-30%, p < 0.001) in desensitised tracheal smooth muscle, as compared to controls. MK, given intraperitoneally at 5 mg/kg/day for 6 consecutive days, completely restored tissue responsiveness to salbutamol. Prolonged salbutamol treatment significantly decreased cAMP synthesis, induced a complete removal of the β2-adrenoceptor from plasma membrane with a parallel increase in the cytosol and increased PDE4D5 gene transcription and PDE activity in human bronchial smooth muscle cells. In homologously desensitised BSMC, MK 30 μM for 24 h was able to prevent salbutamol subsensitivity and such an effect was associated with inhibition of salbutamol-induced PDE4 activity and restoration of membrane β2-adrenoceptor expression and function. These findings suggest the presence of a favourable interaction between MK and β2-adrenoceptor agonists that might improve the therapeutic index of bronchodilators in patients with chronic respiratory diseases.

Selective inhibition of phosphodiesterase 1 relaxes urinary bladder smooth muscle: role for ryanodine receptor-mediated BK channel activation

The large conductance voltage- and Ca(2+)-activated K(+) (BK) channel is a major regulator of detrusor smooth muscle (DSM) excitability and contractility. Recently, we showed that nonselective phosphodiesterase (PDE) inhibition reduces guinea pig DSM excitability and contractility by increasing BK channel activity. Here, we investigated how DSM excitability and contractility changes upon selective inhibition of PDE type 1 (PDE1) and the underlying cellular mechanism involving ryanodine receptors (RyRs) and BK channels. PDE1 inhibition with 8-methoxymethyl-3-isobutyl-1-methylxanthine (8MM-IBMX; 10 μM) increased the cAMP levels in guinea pig DSM cells. Patch-clamp experiments on freshly isolated DSM cells showed that 8MM-IBMX increased transient BK currents and the spontaneous transient hyperpolarization (STH) frequency by ∼2.5- and ∼1.8-fold, respectively. 8MM-IBMX hyperpolarized guinea pig and human DSM cell membrane potential and significantly decreased the intracellular Ca(2+) levels in guinea pig DSM cells. Blocking BK channels with 1 μM paxilline or inhibiting RyRs with 30 μM ryanodine abolished the STHs and the 8MM-IBMX inhibitory effects on the DSM cell membrane potential. Isometric DSM tension recordings showed that 8MM-IBMX significantly reduced the spontaneous phasic contraction amplitude, muscle force integral, duration, frequency, and tone of DSM isolated strips. The electrical field stimulation-induced DSM contraction amplitude, muscle force integral, and duration were also attenuated by 10 μM 8MM-IBMX. Blocking BK channels with paxilline abolished the 8MM-IBMX effects on DSM contractions. Our data provide evidence that PDE1 inhibition relaxes DSM by raising cellular cAMP levels and subsequently stimulates RyRs, which leads to BK channel activation, membrane potential hyperpolarization, and decrease in intracellular Ca(2+) levels.

β-Adrenergic receptor stimulation increases surface NKCC2 expression in rat thick ascending limbs in a process inhibited by phosphodiesterase 4

The thick ascending limb of the loop of Henle (THAL) reabsorbs ∼30% of the filtered NaCl in a process mediated by the apical Na-K-2Cl cotransporter NKCC2. Stimulation of β-adrenergic receptors in the THAL enhances NaCl reabsorption and increases intracellular cAMP. We found that intracellular cAMP stimulates NKCC2 trafficking to the apical membrane via protein kinase A (PKA). Several cAMP-specific phosphodiesterases (PDE) have been identified in rat THALs, and PDE4 decreases cAMP generated by β-adrenergic stimulation in other cells. However, it is not known whether β-adrenergic receptors activation stimulates NKCC2 trafficking. Thus we hypothesized that β-adrenergic receptor stimulation enhances THAL apical membrane NKCC2 expression via the PKA pathway and PDE4 blunts this effect. THAL suspensions were obtained from Sprague-Dawley rats, and surface NKCC2 expression was measured by surface biotinylation and Western blot. Incubation of THALs with the β-adrenergic receptor agonist isoproterenol at 0.5 and 1.0 μM increased surface NKCC2 by 17 ± 1 and 29 ± 5% respectively (P < 0.05). Preventing cAMP degradation with 3-isobutyl-methylxanthine (IBMX; a nonselective phosphodiesterase inhibitor) enhanced isoproterenol-stimulated surface NKCC2 expression to 51 ± 7% (P < 0.05 vs. isoproterenol). The β-adrenergic receptor antagonist propranolol or the PKA inhibitor H-89 completely blocked isoproterenol + IBMX-induced increase on surface NKCC2, while propranolol or H-89 alone had no effect. Selective inhibition of PDE4 with rolipram (20 μM) potentiated the effect of isoproterenol on surface NKCC2 and increased cAMP levels. We concluded that β-adrenergic receptor stimulation enhances surface NKCC2 expression in the THALs via PKA and PDE4 blunts this effect.

Calcineurin regulates homologous desensitization of natriuretic peptide receptor-A and inhibits ANP-induced testosterone production in MA-10 cells

Receptor desensitization is a ubiquitous regulatory mechanism that defines the activatable pool of receptors, and thus, the ability of cells to respond to environmental stimuli. In recent years, the molecular mechanisms controlling the desensitization of a variety of receptors have been established. However, little is known about the molecular mechanisms that underlie desensitization of natriuretic peptide receptors, including natriuretic peptide receptor-A (NPR-A). Here we report that calcineurin (protein phosphatase 2B, PP2B, PPP3C) regulates homologous desensitization of NPR-A in murine Leydig tumor (MA-10) cells. We demonstrate that both pharmacological inhibition of calcineurin activity and siRNA-mediated suppression of calcineurin expression potentiate atrial natriuretic peptide (ANP)-induced cGMP synthesis. Treatment of MA-10 cells with inhibitors of other phosphoprotein phosphatases had little or no effect on ANP-induced cGMP accumulation. In addition, overexpression of calcineurin blunts ANP-induced cGMP synthesis. We also present data indicating that the inhibition of calcineurin potentiates ANP-induced testosterone production. To better understand the contribution of calcineurin in the regulation of NPR-A activity, we examined the kinetics of ANP-induced cGMP signals. We observed transient ANP-induced cGMP signals, even in the presence of phosphodiesterase inhibitors. Inhibition of both calcineurin and phosphodiesterase dramatically slowed the decay in the response. These observations are consistent with a model in which calcineurin mediated dephosphorylation and desensitization of NPR-A is associated with significant inhibition of cGMP synthesis. PDE activity hydrolyzes cGMP, thus lowering intracellular cGMP toward the basal level. Taken together, these data suggest that calcineurin plays a previously unrecognized role in the desensitization of NPR-A and, thereby, inhibits ANP-mediated increases in testosterone production.

A-kinase anchoring proteins regulate compartmentalized cAMP signaling in airway smooth muscle

A-kinase anchoring proteins (AKAPs) have emerged as important regulatory molecules that can compartmentalize cAMP signaling transduced by β2-adrenergic receptors (β(2)ARs); such compartmentalization ensures speed and fidelity of cAMP signaling and effects on cell function. This study aimed to assess the role of AKAPs in regulating global and compartmentalized β(2)AR signaling in human airway smooth muscle (ASM). Transcriptome and proteomic analyses were used to characterize AKAP expression in ASM. Stable expression or injection of peptides AKAP-IS or Ht31 was used to disrupt AKAP-PKA interactions, and global and compartmentalized cAMP accumulation stimulated by β-agonist was assessed by radioimmunoassay and membrane-delineated flow through cyclic nucleotide-gated channels, respectively. ASM expresses multiple AKAP family members, with gravin and ezrin among the most readily detected. AKAP-PKA disruption had minimal effects on whole-cell cAMP accumulation stimulated by β-agonist (EC(50) and B(max)) concentrations, but significantly increased the duration of plasma membrane-delineated cAMP (τ=251±51 s for scrambled peptide control vs. 399±79 s for Ht31). Direct PKA inhibition eliminated decay of membrane-delineated cAMP levels. AKAPs coordinate compartmentalized cAMP signaling in ASM cells by regulating multiple elements of β(2)AR-mediated cAMP accumulation, thereby representing a novel target for manipulating β(2)AR signaling and function in ASM.

Regulation of nuclear PKA revealed by spatiotemporal manipulation of cyclic AMP

Understanding how specific cAMP signals are organized and relayed to their effectors in different compartments of the cell to achieve functional specificity requires molecular tools that allow precise manipulation of cAMP in these compartments. Here we characterize a new method using bicarbonate-activatable and genetically targetable soluble adenylyl cyclase (sAC) to control the location, kinetics and magnitude of the cAMP signal. Using this live-cell cAMP manipulation in conjunction with fluorescence imaging and mechanistic modeling, we uncover the activation of a resident pool of PKA holoenzyme in the nuclei of HEK-293 cells, modifying the existing dogma of cAMP-PKA signaling in the nucleus. Furthermore, we show that phosphodiesterases (PDE) and A-Kinase Anchoring Proteins (AKAP) are critical in shaping nuclear PKA responses. Collectively, our data suggests a new model where AKAP-localized PDEs tune an activation threshold for nuclear PKA holoenzyme, thereby converting spatially distinct second messenger signals to temporally controlled nuclear kinase activity.

A biphasic and brain-region selective down-regulation of cyclic adenosine monophosphate concentrations supports object recognition in the rat

BACKGROUND

We aimed to further understand the relationship between cAMP concentration and mnesic performance.

METHODS AND FINDINGS

Rats were injected with milrinone (PDE3 inhibitor, 0.3 mg/kg, i.p.), rolipram (PDE4 inhibitor, 0.3 mg/kg, i.p.) and/or the selective 5-HT4R agonist RS 67333 (1 mg/kg, i.p.) before testing in the object recognition paradigm. Cyclic AMP concentrations were measured in brain structures linked to episodic-like memory (i.e. hippocampus, prefrontal and perirhinal cortices) before or after either the sample or the testing phase. Except in the hippocampus of rolipram treated-rats, all treatment increased cAMP levels in each brain sub-region studied before the sample phase. After the sample phase, cAMP levels were significantly increased in hippocampus (1.8 fold), prefrontal (1.3 fold) and perirhinal (1.3 fold) cortices from controls rat while decreased in prefrontal cortex (∼0.83 to 0.62 fold) from drug-treated rats (except for milrinone+RS 67333 treatment). After the testing phase, cAMP concentrations were still increased in both the hippocampus (2.76 fold) and the perirhinal cortex (2.1 fold) from controls animals. Minor increase were reported in hippocampus and perirhinal cortex from both rolipram (respectively, 1.44 fold and 1.70 fold) and milrinone (respectively 1.46 fold and 1.56 fold)-treated rat. Following the paradigm, cAMP levels were significantly lower in the hippocampus, prefrontal and perirhinal cortices from drug-treated rat when compared to controls animals, however, only drug-treated rats spent longer time exploring the novel object during the testing phase (inter-phase interval of 4 h).

CONCLUSIONS

Our results strongly suggest that a "pre-sample" early increase in cAMP levels followed by a specific lowering of cAMP concentrations in each brain sub-region linked to the object recognition paradigm support learning efficacy after a middle-term delay.

Assessment of cellular mechanisms contributing to cAMP compartmentalization in pulmonary microvascular endothelial cells

Cyclic AMP signals encode information required to differentially regulate a wide variety of cellular responses; yet it is not well understood how information is encrypted within these signals. An emerging concept is that compartmentalization underlies specificity within the cAMP signaling pathway. This concept is based on a series of observations indicating that cAMP levels are distinct in different regions of the cell. One such observation is that cAMP production at the plasma membrane increases pulmonary microvascular endothelial barrier integrity, whereas cAMP production in the cytosol disrupts barrier integrity. To better understand how cAMP signals might be compartmentalized, we have developed mathematical models in which cellular geometry as well as total adenylyl cyclase and phosphodiesterase activities were constrained to approximate values measured in pulmonary microvascular endothelial cells. These simulations suggest that the subcellular localizations of adenylyl cyclase and phosphodiesterase activities are by themselves insufficient to generate physiologically relevant cAMP gradients. Thus, the assembly of adenylyl cyclase, phosphodiesterase, and protein kinase A onto protein scaffolds is by itself unlikely to ensure signal specificity. Rather, our simulations suggest that reductions in the effective cAMP diffusion coefficient may facilitate the formation of substantial cAMP gradients. We conclude that reductions in the effective rate of cAMP diffusion due to buffers, structural impediments, and local changes in viscosity greatly facilitate the ability of signaling complexes to impart specificity within the cAMP signaling pathway.

β2-Agonists upregulate PDE4 mRNA but not protein or activity in human airway smooth muscle cells from asthmatic and nonasthmatic volunteers

β(2)-Adrenergic receptor (β2AR) agonists induce airway relaxation via cAMP. Phosphodiesterase (PDE)s degrade and regulate cAMP, and in airway smooth muscle (ASM) cells PDE4D degrades cAMP. Long-acting β(2)-agonists are now contraindicated as monotherapy for asthma, and increased PDE4D has been speculated to contribute to this phenomenon. In this study we investigated the expression of PDE4D in asthmatic and nonasthmatic ASM cells and its regulation by formoterol and budesonide. Primary ASM cells from people with or without asthma were stimulated with transforming growth factor (TGF)-β(1), formoterol, and/or budesonide. PDE4D mRNA was assessed by real-time PCR, or PCR to assess splice variant production. PDE4D protein was assessed by Western blotting, and we investigated the effect of formoterol on cAMP production and PDE activity. Interleukin (IL)-6 was assessed using ELISA. PDE4D mRNA was dose dependently upregulated by formoterol, with a single splice variant, PDE4D5, present. Formoterol did not induce PDE4D protein at time points between 3 to 72 h, whereas it did induce and increase IL-6 secretion. We pretreated cells with actinomycin D and a proteasome inhibitor, MG132, and found no evidence of alterations in mRNA, protein expression, or degradation of PDE4D. Finally PDE activity was not altered by formoterol. This study shows, for the first time, that PDE4D5 is predominantly expressed in human ASM cells from people with and without asthma and that formoterol does not upregulate PDE4D protein production. This leads us to speculate that continual therapy with β2AR agonists is unlikely to cause PDE4-mediated tachyphylaxis.

All preconditioning-related G protein-coupled receptors can be demonstrated in the rabbit cardiomyocyte

G protein-coupled receptors for adenosine (A(1), A(3), A(2A), and A(2B)), bradykinin (B(1)) and opioids (δ) are all involved in the mechanism of ischemic preconditioning. Although the heart is comprised of many tissue types, it has been assumed that preconditioning's protective signaling occurs in the cardiomyocyte. We critically tested that hypothesis by testing for the presence of each of these receptors in isolated adult rabbit ventricular myocytes that had been transfected with cyclic nucleotide-gated (CNG) ion channels. Because subsarcolemmal cyclic adenosine monophosphate (cAMP) opens the CNG channels, we could monitor cAMP levels within a single cardiomyocyte by measuring channel current with a patch pipette. The presence of a receptor would be confirmed if we could alter cAMP in the cell with a selective agonist to the receptor being studied. Superfusion with the β-adrenergic G(s)-coupled receptor agonist isoproterenol (50 nmol/L) transiently increased cAMP levels and, therefore, channel current. Pretreatment with selective agonists to A(1) or A(3) adenosine receptors (ARs) that are G(i)-coupled markedly attenuated the response to isoproterenol, indicating inhibition of adenylyl cyclase by increased G(i) activity. Agonists to bradykinin or δ-opioid receptors also attenuated isoproterenol's response. A(2A)AR and A(2B)AR are G(s)-coupled. The A(2A)AR-selective agonist CGS21680 increased current through CNG channels but only in the presence of phosphodiesterase (PDE) inhibitors, indicating low surface receptor activity and high intracellular PDE activity. As we previously reported, BAY 60-6583, an A(2B)AR-selective agonist which mimics preconditioning's protection in rabbit heart, neither increased nor decreased membrane current in transfected cardiomyocytes, suggesting the absence or a markedly limited number of A(2B)AR in the sarcolemma. However, reverse transcription polymerase chain reaction (RT-PCR) of purified cardiomyocytes yielded an A(2B)AR band, implying that rabbit cardiomyocytes do indeed express A(2B)AR. These data reveal that all receptors reported to be involved in ischemic preconditioning do exist on or within the cardiomyocyte.

A2B adenosine receptors inhibit superoxide production from mitochondrial complex I in rabbit cardiomyocytes via a mechanism sensitive to Pertussis toxin

BACKGROUND AND PURPOSE

A(2B) adenosine receptors protect against ischaemia/reperfusion injury by activating survival kinases including extracellular signal-regulated kinase (ERK) and phosphatidylinositol 3-kinase (PI3K). However, the underlying mechanism(s) and signalling pathway(s) remain undefined.

EXPERIMENTAL APPROACH

HEK 293 cells stably transfected with human A(2B) adenosine receptors (HEK-A(2B) ) and isolated adult rabbit cardiomyocytes were used to assay phosphorylation of ERK by Western blot and cation flux through cAMP-gated channels by patch clamp methods. Generation of reactive oxygen species (ROS) by mitochondria was measured with a fluorescent dye.

KEY RESULTS

In HEK-A(2B) cells, the selective A(2B) receptor agonist Bay 60-6583 (Bay 60) increased ERK phosphorylation and cAMP levels, detected by current through cAMP-gated ion channels. However, increased cAMP or its downstream target protein kinase A was not involved in ERK phosphorylation. Pertussis toxin (PTX) blocked ERK phosphorylation, suggesting receptor coupling to G(i) or G(o) proteins. Phosphorylation was also blocked by inhibition of PI3K (with wortmannin) or of ERK kinase (MEK1/2, with PD 98059) but not by inhibition of NO synthase (NOS). In cardiomyocytes, Bay 60 did not affect cAMP levels but did block the increased superoxide generation induced by rotenone, a mitochondrial complex I inhibitor. This effect of Bay 60 was inhibited by PD 98059, wortmannin or PTX. Inhibition of NOS blocked superoxide production because NOS is downstream of ERK.

CONCLUSION AND IMPLICATIONS

Activation of A(2B) adenosine receptors reduced superoxide generation from mitochondrial complex I through G(i/o) , ERK, PI3K, and NOS, all of which have been implicated in ischaemic preconditioning.

Local control of β-adrenergic stimulation: Effects on ventricular myocyte electrophysiology and Ca(2+)-transient

Local signaling domains and numerous interacting molecular pathways and substrates contribute to the whole-cell response of myocytes during β-adrenergic stimulation (βARS). We aimed to elucidate the quantitative contribution of substrates and their local signaling environments during βARS to the canine epicardial ventricular myocyte electrophysiology and calcium transient (CaT). We present a computational compartmental model of βARS and its electrophysiological effects. Novel aspects of the model include localized signaling domains, incorporation of β1 and β2 receptor isoforms, a detailed population-based approach to integrate the βAR and Ca(2+)/Calmodulin kinase (CaMKII) signaling pathways and their effects on a wide range of substrates that affect whole-cell electrophysiology and CaT. The model identifies major roles for phosphodiesterases, adenylyl cyclases, PKA and restricted diffusion in the control of local cAMP levels and shows that activation of specific cAMP domains by different receptor isoforms allows for specific control of action potential and CaT properties. In addition, the model predicts increased CaMKII activity during βARS due to rate-dependent accumulation and increased Ca(2+) cycling. CaMKII inhibition, reduced compartmentation, and selective blockade of β1AR is predicted to reduce the occurrence of delayed afterdepolarizations during βARS. Finally, the relative contribution of each PKA substrate to whole-cell electrophysiology is quantified by comparing simulations with and without phosphorylation of each target. In conclusion, this model enhances our understanding of localized βAR signaling and its whole-cell effects in ventricular myocytes by incorporating receptor isoforms, multiple pathways and a detailed representation of multiple-target phosphorylation; it provides a basis for further studies of βARS under pathological conditions.

Receptor-mediated enhancement of beta adrenergic drug activity by ascorbate in vitro and in vivo

Rationale

Previous in vitro research demonstrated that ascorbate enhances potency and duration of activity of agonists binding to alpha 1 adrenergic and histamine receptors.

Objectives

Extending this work to beta 2 adrenergic systems in vitro and in vivo.

Methods

Ultraviolet spectroscopy was used to study ascorbate binding to adrenergic receptor preparations and peptides. Force transduction studies on acetylcholine-contracted trachealis preparations from pigs and guinea pigs measured the effect of ascorbate on relaxation due to submaximal doses of beta adrenergic agonists. The effect of inhaled albuterol with and without ascorbate was tested on horses with heaves and sheep with carbachol-induced bronchoconstriction.

Measurements

Binding constants for ascorbate binding to beta adrenergic receptor were derived from concentration-dependent spectral shifts. Dose- dependence curves were obtained for the relaxation of pre-contracted trachealis preparations due to beta agonists in the presence and absence of varied ascorbate. Tachyphylaxis and fade were also measured. Dose response curves were determined for the effect of albuterol plus-and-minus ascorbate on airway resistance in horses and sheep.

Main Results

Ascorbate binds to the beta 2 adrenergic receptor at physiological concentrations. The receptor recycles dehydroascorbate. Physiological and supra-physiological concentrations of ascorbate enhance submaximal epinephrine and isoproterenol relaxation of trachealis, producing a 3–10-fold increase in sensitivity, preventing tachyphylaxis, and reversing fade. In vivo, ascorbate improves albuterol's effect on heaves and produces a 10-fold enhancement of albuterol activity in “asthmatic” sheep.

Conclusions

Ascorbate enhances beta-adrenergic activity via a novel receptor-mediated mechanism; increases potency and duration of beta adrenergic agonists effective in asthma and COPD; prevents tachyphylaxis; and reverses fade. These novel effects are probably caused by a novel mechanism involving phosphorylation of aminergic receptors and have clinical and drug-development applications.

Modeling cardiac β-adrenergic signaling with normalized-Hill differential equations: comparison with a biochemical model

Background

New approaches are needed for large-scale predictive modeling of cellular signaling networks. While mass action and enzyme kinetic approaches require extensive biochemical data, current logic-based approaches are used primarily for qualitative predictions and have lacked direct quantitative comparison with biochemical models.

Results

We developed a logic-based differential equation modeling approach for cell signaling networks based on normalized Hill activation/inhibition functions controlled by logical AND and OR operators to characterize signaling crosstalk. Using this approach, we modeled the cardiac β₁-adrenergic signaling network, including 36 reactions and 25 species. Direct comparison of this model to an extensively characterized and validated biochemical model of the same network revealed that the new model gave reasonably accurate predictions of key network properties, even with default parameters. Normalized Hill functions improved quantitative predictions of global functional relationships compared with prior logic-based approaches. Comprehensive sensitivity analysis revealed the significant role of PKA negative feedback on upstream signaling and the importance of phosphodiesterases as key negative regulators of the network. The model was then extended to incorporate recently identified protein interaction data involving integrin-mediated mechanotransduction.

Conclusions

The normalized-Hill differential equation modeling approach allows quantitative prediction of network functional relationships and dynamics, even in systems with limited biochemical data.

The role of type 4 phosphodiesterases in generating microdomains of cAMP: large scale stochastic simulations

Cyclic AMP (cAMP) and its main effector Protein Kinase A (PKA) are critical for several aspects of neuronal function including synaptic plasticity. Specificity of synaptic plasticity requires that cAMP activates PKA in a highly localized manner despite the speed with which cAMP diffuses. Two mechanisms have been proposed to produce localized elevations in cAMP, known as microdomains: impeded diffusion, and high phosphodiesterase (PDE) activity. This paper investigates the mechanism of localized cAMP signaling using a computational model of the biochemical network in the HEK293 cell, which is a subset of pathways involved in PKA-dependent synaptic plasticity. This biochemical network includes cAMP production, PKA activation, and cAMP degradation by PDE activity. The model is implemented in NeuroRD: novel, computationally efficient, stochastic reaction-diffusion software, and is constrained by intracellular cAMP dynamics that were determined experimentally by real-time imaging using an Epac-based FRET sensor (H30). The model reproduces the high concentration cAMP microdomain in the submembrane region, distinct from the lower concentration of cAMP in the cytosol. Simulations further demonstrate that generation of the cAMP microdomain requires a pool of PDE4D anchored in the cytosol and also requires PKA-mediated phosphorylation of PDE4D which increases its activity. The microdomain does not require impeded diffusion of cAMP, confirming that barriers are not required for microdomains. The simulations reported here further demonstrate the utility of the new stochastic reaction-diffusion algorithm for exploring signaling pathways in spatially complex structures such as neurons.

Quantitative modeling of GRK-mediated beta2AR regulation

We developed a unified model of the GRK-mediated β2 adrenergic receptor (β2AR) regulation that simultaneously accounts for six different biochemical measurements of the system obtained over a wide range of agonist concentrations. Using a single deterministic model we accounted for (1) GRK phosphorylation in response to various full and partial agonists; (2) dephosphorylation of the GRK site on the β2AR; (3) β2AR internalization; (4) recycling of the β2AR post isoproterenol treatment; (5) β2AR desensitization; and (6) β2AR resensitization. Simulations of our model show that plasma membrane dephosphorylation and recycling of the phosphorylated receptor are necessary to adequately account for the measured dephosphorylation kinetics. We further used the model to predict the consequences of (1) modifying rates such as GRK phosphorylation of the receptor, arrestin binding and dissociation from the receptor, and receptor dephosphorylation that should reflect effects of knockdowns and overexpressions of these components; and (2) varying concentration and frequency of agonist stimulation “seen” by the β2AR to better mimic hormonal, neurophysiological and pharmacological stimulations of the β2AR. Exploring the consequences of rapid pulsatile agonist stimulation, we found that although resensitization was rapid, the β2AR system retained the memory of the previous stimuli and desensitized faster and much more strongly in response to subsequent stimuli. The latent memory that we predict is due to slower membrane dephosphorylation, which allows for progressive accumulation of phosphorylated receptor on the surface. This primes the receptor for faster arrestin binding on subsequent agonist activation leading to a greater extent of desensitization. In summary, the model is unique in accounting for the behavior of the β2AR system across multiple types of biochemical measurements using a single set of experimentally constrained parameters. It also provides insight into how the signaling machinery can retain memory of prior stimulation long after near complete resensitization has been achieved.

The β2 adrenergic receptor (β2AR) is involved in regulating many cellular processes such as smooth muscle relaxation in the airways and the vasculature. Drugs that activate the β2AR are used in treating asthma and chronic obstructive pulmonary disorder (COPD), and prolonged use of these drugs leads to the loss of their effects. Thus, a dynamic model of how the β2AR responds to different drugs is fundamental to their rational use. In this study a consensus model of G protein coupled receptor kinase (GRK)-mediated receptor regulation was formulated based on quantitative measures of six processes involved in β2AR regulation. This model was then used to simulate the consequences of manipulating key rates associated with the GRK-mediated β2AR regulation, leading to predictions which will provide a useful framework for further tests and elaborations of the model in basic and clinical research.

A mathematical analysis of second messenger compartmentalization

Intracellular compartmentalization of second messengers can lead to microdomains of elevated concentration that are thought to be involved in ensuring signaling specificity. Most experimental evidence for this compartmentalization involves the second messenger adenosine monophosphate (cAMP), which is degraded by phosphodiesterases (PDEs). One possible way of creating these compartments, supported by recent experiments, is to spatially separate the source of cAMP from regions of elevated PDE concentration. To quantify this possibility, we study here a simplified geometry in two dimensions (2D) and in three dimensions (3D), containing a cAMP point source and regions with different degradation constants. Using the symmetry of our geometry, we are able to derive steady state solutions for the cAMP concentration as a function of the system parameters. Furthermore, we show, using analytics as well as direct numerical simulations, that for physiologically relevant time scales the steady state solution has been reached. Our results indicate that elevating the degradation constant throughout the cell, except for a small microdomain surrounding the source, requires an unphysiologically high cellular PDE concentration. On the other hand, a tight spatial relationship of localized PDEs with the cAMP source can result in functional microdomains while maintaining a physiologically plausible cellular PDE concentration.

Endogenous Gs-coupled receptors in smooth muscle exhibit differential susceptibility to GRK2/3-mediated desensitization

Although G protein-coupled receptor (GPCR) kinases (GRKs) have been shown to mediate desensitization of numerous GPCRs in studies using cellular expression systems, their function under physiological conditions is less well understood. In the current study, we employed various strategies to assess the effect of inhibiting endogenous GRK2/3 on signaling and function of endogenously expressed G s-coupled receptors in human airway smooth muscle (ASM) cells. GRK2/3 inhibition by expression of a Gbetagamma sequestrant, a GRK2/3 dominant-negative mutant, or siRNA-mediated knockdown increased intracellular cAMP accumulation mediated via beta-agonist stimulation of the beta-2-adrenergic receptor (beta 2AR). Conversely, neither 5'-( N-ethylcarboxamido)-adenosine (NECA; activating the A2b adenosine receptor) nor prostaglandin E2 (PGE 2; activating EP2 or EP4 receptors)-stimulated cAMP was significantly increased by GRK2/3 inhibition. Selective knockdown using siRNA suggested the majority of PGE 2-stimulated cAMP in ASM was mediated by the EP2 receptor. Although a minor role for EP3 receptors in influencing PGE 2-mediated cAMP was determined, the GRK2/3-resistant nature of EP2 receptor signaling in ASM was confirmed using the EP2-selective agonist butaprost. Somewhat surprisingly, GRK2/3 inhibition did not augment the inhibitory effect of the beta-agonist on mitogen-stimulated increases in ASM growth. These findings demonstrate that with respect to G s-coupled receptors in ASM, GRK2/3 selectively attenuates beta 2AR signaling, yet relief of GRK2/3-dependent beta 2AR desensitization does not influence at least one important physiological function of the receptor.

Critical role of PDE4D in beta2-adrenoceptor-dependent cAMP signaling in mouse embryonic fibroblasts

One of the defining properties of beta2-adrenergic receptor (beta(2)AR) signaling is the transient and rapidly reversed accumulation of cAMP. Here we have investigated the contribution of different PDE4 proteins to the generation of this transient response. To this aim, mouse embryonic fibroblasts deficient in PDE4A, PDE4B, or PDE4D were generated, and the regulation of PDE activity, the accumulation of cAMP, and CREB phosphorylation in response to isoproterenol were monitored. Ablation of PDE4D, but not PDE4A or PDE4B, had a major effect on the beta-agonist-induced PDE activation, with only a minimal increase in PDE activity being retained in PDE4D knock-out (KO) cells. Accumulation of cAMP was markedly enhanced, and the kinetics of cAMP accumulation were altered in their properties in PDE4DKO but not PDE4BKO cells. Modest effects were observed in PDE4AKO mouse embryonic fibroblasts. The return to basal levels of both cAMP accumulation and CREB phosphorylation was greatly delayed in the PDE4DKO cells, suggesting that PDE4D is critical for dissipation of the beta2AR stimulus. This effect of PDE4D ablation was in large part due to inactivation of a negative feedback mechanism consisting of the PKA-mediated activation of PDE4D in response to elevated cAMP levels, as indicated by experiments using the cAMP-dependent protein kinase inhibitors H89 and PKI. Finally, PDE4D ablation affected the kinetics of beta2AR desensitization as well as the interaction of the receptor with Galphai. These findings demonstrate that PDE4D plays a major role in shaping the beta2AR signal.