Kinetic diversity in G-protein-coupled receptor signalling

Systems modeling of oncogenic G-protein and GPCR signaling reveals unexpected differences in downstream pathway activation

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gαq/11 and CysLT2R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT2R was impaired at activating the FAK/YAP/TAZ pathway relative to Gαq/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

Neomorphic Gαo mutations gain interaction with Ric8 proteins in GNAO1 encephalopathies

GNAO1 mutated in pediatric encephalopathies encodes the major neuronal G protein Gαo. Of the more than 80 pathogenic mutations, most are single amino acid substitutions spreading across the Gαo sequence. We performed extensive characterization of Gαo mutants, showing abnormal GTP uptake and hydrolysis and deficiencies in binding Gβγ and RGS19. Plasma membrane localization of Gαo was decreased for a subset of mutations that leads to epilepsy; dominant interactions with GPCRs also emerged for the more severe mutants. Pathogenic mutants massively gained interaction with Ric8A and, surprisingly, Ric8B proteins, relocalizing them from cytoplasm to Golgi. Of these 2 mandatory Gα-subunit chaperones, Ric8A is normally responsible for the Gαi/Gαo, Gαq, and Gα12/Gα13 subfamilies, and Ric8B solely responsible for Gαs/Gαolf. Ric8 mediates the disease dominance when engaging in neomorphic interactions with pathogenic Gαo through imbalance of the neuronal G protein signaling networks. As the strength of Gαo-Ric8B interactions correlates with disease severity, our study further identifies an efficient biomarker and predictor for clinical manifestations in GNAO1 encephalopathies. Our work uncovers the neomorphic molecular mechanism of mutations underlying pediatric encephalopathies and offers insights into other maladies caused by G protein malfunctioning and further genetic diseases.

Kinetic Model for the Desensitization of G Protein-Coupled Receptor

Desensitization of G-protein-coupled receptors (GPCR) is a general regulatory mechanism adopted by biological organisms against overstimulation of G protein signaling. Although the details of the mechanism are extensively studied, it is not easy to gain an overarching understanding of the process constituted by a multitude of molecular events with vastly differing time scales. To offer a semiquantitative yet predictive understanding of the mechanism, we formulate a kinetic model for the G protein signaling and desensitization by considering essential biochemical steps from ligand binding to receptor internalization. The internalization, followed by receptor depletion from the plasma membrane, attenuates the downstream signal. Together with the kinetic model and its full numerics of the expression derived for the dose-response relation, an approximated form of the expression clarifies the role played by the individual biochemical processes and allows us to identify four distinct regimes for the downregulation that emerge from the balance between phosphorylation, dephosphorylation, and the cellular level of β-arrestin.

Clinical Cases and the Molecular Profiling of a Novel Childhood Encephalopathy-Causing GNAO1 Mutation P170R

De novo mutations in GNAO1, the gene encoding the major neuronal G protein Gαo, cause a spectrum of pediatric encephalopathies with seizures, motor dysfunction, and developmental delay. Of the >80 distinct missense pathogenic variants, many appear to uniformly destabilize the guanine nucleotide handling of the mutant protein, speeding up GTP uptake and deactivating GTP hydrolysis. Zinc supplementation emerges as a promising treatment option for this disease, as Zn2+ ions reactivate the GTP hydrolysis on the mutant Gαo and restore cellular interactions for some of the mutants studied earlier. The molecular etiology of GNAO1 encephalopathies needs further elucidation as a prerequisite for the development of efficient therapeutic approaches. In this work, we combine clinical and medical genetics analysis of a novel GNAO1 mutation with an in-depth molecular dissection of the resultant protein variant. We identify two unrelated patients from Norway and France with a previously unknown mutation in GNAO1, c.509C>G that results in the production of the Pro170Arg mutant Gαo, leading to severe developmental and epileptic encephalopathy. Molecular investigations of Pro170Arg identify this mutant as a unique representative of the pathogenic variants. Its 100-fold-accelerated GTP uptake is not accompanied by a loss in GTP hydrolysis; Zn2+ ions induce a previously unseen effect on the mutant, forcing it to lose the bound GTP. Our work combining clinical and molecular analyses discovers a novel, biochemically distinct pathogenic missense variant of GNAO1 laying the ground for personalized treatment development.

Restoration of the GTPase activity and cellular interactions of Gαo mutants by Zn2+ in GNAO1 encephalopathy models

De novo point mutations in GNAO1, gene encoding the major neuronal G protein Gα_o, have recently emerged in patients with pediatric encephalopathy having motor, developmental, and epileptic dysfunctions. Half of clinical cases affect codons Gly²⁰³, Arg²⁰⁹, or Glu²⁴⁶; we show that these mutations accelerate GTP uptake and inactivate GTP hydrolysis through displacement Gln²⁰⁵ critical for GTP hydrolysis, resulting in constitutive GTP binding by Gα_o. However, the mutants fail to adopt the activated conformation and display aberrant interactions with signaling partners. Through high-throughput screening of approved drugs, we identify zinc pyrithione and Zn²⁺ as agents restoring active conformation, GTPase activity, and cellular interactions of the encephalopathy mutants, with negligible effects on wild-type Gα_o. We describe a Drosophila model of GNAO1 encephalopathy where dietary zinc restores the motor function and longevity of the mutant flies. Zinc supplements are approved for diverse human neurological conditions. Our work provides insights into the molecular etiology of GNAO1 encephalopathy and defines a potential therapy for the patients.

The critical balance between dopamine D2 receptor and RGS for the sensitive detection of a transient decay in dopamine signal

In behavioral learning, reward-related events are encoded into phasic dopamine (DA) signals in the brain. In particular, unexpected reward omission leads to a phasic decrease in DA (DA dip) in the striatum, which triggers long-term potentiation (LTP) in DA D2 receptor (D2R)-expressing spiny-projection neurons (D2 SPNs). While this LTP is required for reward discrimination, it is unclear how such a short DA-dip signal (0.5-2 s) is transferred through intracellular signaling to the coincidence detector, adenylate cyclase (AC). In the present study, we built a computational model of D2 signaling to determine conditions for the DA-dip detection. The DA dip can be detected only if the basal DA signal sufficiently inhibits AC, and the DA-dip signal sufficiently disinhibits AC. We found that those two requirements were simultaneously satisfied only if two key molecules, D2R and regulators of G protein signaling (RGS) were balanced within a certain range; this balance has indeed been observed in experimental studies. We also found that high level of RGS was required for the detection of a 0.5-s short DA dip, and the analytical solutions for these requirements confirmed their universality. The imbalance between D2R and RGS is associated with schizophrenia and DYT1 dystonia, both of which are accompanied by abnormal striatal LTP. Our simulations suggest that D2 SPNs in patients with schizophrenia and DYT1 dystonia cannot detect short DA dips. We finally discussed that such psychiatric and movement disorders can be understood in terms of the imbalance between D2R and RGS.

High expression of PSMC2 promotes gallbladder cancer through regulation of GNG4 and predicts poor prognosis

Gallbladder cancer (GBC) is a common malignant tumor of the biliary tract, which accounts for 80-95% of biliary tumors worldwide, and is the leading cause of biliary malignant tumor-related death. This study identified PSMC2 as a potential regulator in the development of GBC. We showed that PSMC2 expression in GBC tissues is significantly higher than that in normal tissues, while high PSMC2 expression was correlated with more advanced tumor grade and poorer prognosis. The knockdown of PSMC2 in GBC cells induced significant inhibition of cell proliferation, colony formation and cell motility, while the promotion of cell apoptosis. The construction and observation of the mice xenograft model also confirmed the inhibitory effects of PSMC2 knockdown on GBC development. Moreover, our mechanistic study recognized GNG4 as a potential downstream target of PSMC2, knockdown of which could aggravate the tumor suppression induced by PSMC2 knockdown in vitro and in vivo. In conclusion, for the first time, PSMC2 was revealed as a tumor promotor in the development of GBC, which could regulate cell phenotypes of GBC cells through the interaction with GNG4, and maybe a promising therapeutic target in GBC treatment.

A novel molecular mechanism responsible for phosphorylation-independent desensitization of G protein-coupled receptors exemplified by the dopamine D3 receptor

GRK-mediated receptor phosphorylation followed by association with β-arrestins has been proposed to be the molecular mechanism involved in the desensitization of G protein-coupled receptors (GPCRs). However, this mechanism does not explain the desensitization of some GPCRs, such as dopamine D₃ receptor (D₃R), which does not undergo GRK-mediated phosphorylation. Loss-of-function approaches and mutants of dopamine D₂ receptor and D₃R, which exhibit different desensitization properties, were used to identify the cellular components and processes responsible for desensitization. D₃R mediated the recruitment of Mdm2 to the cytosol, which resulted in the constitutive ubiquitination of β-arrestin2 in the resting state. Under desensitization conditions, cytosolic Mdm2 returned to the nucleus, resulting in the deubiquitination of cytosolic β-arrestins. Deubiquitinated β-arrestins formed a tight complex with Gβγ, thereby sequestering it, causing interference in D₃R signaling. In conclusion, this study shows that β-arrestins, depending on their ubiquitination status, control the G protein cycling by regulating their interactions with Gβγ. This is a novel mechanism proposed to explain how certain GPCRs can undergo desensitization without receptor phosphorylation.

High capacity in G protein-coupled receptor signaling

G protein-coupled receptors (GPCRs) constitute a large family of receptors that activate intracellular signaling pathways upon detecting specific extracellular ligands. While many aspects of GPCR signaling have been uncovered through decades of studies, some fundamental properties, like its channel capacity—a measure of how much information a given transmission system can reliably transduce—are still debated. Previous studies concluded that GPCRs in individual cells could transmit around one bit of information about the concentration of the ligands, allowing only for a reliable on or off response. Using muscarinic receptor-induced calcium response measured in individual cells upon repeated stimulation, we show that GPCR signaling systems possess a significantly higher capacity. We estimate the channel capacity of this system to be above two, implying that at least four concentration levels of the agonist can be distinguished reliably. These findings shed light on the basic principles of GPCR signaling.

G protein-coupled receptors (GPCRs) activate intracellular signalling pathways upon extracellular stimulation. Here authors record single cell responses of GPCR signalling which allows the direct estimation of its channel capacity for each cell along with the reproducibility of its response.

A Model for Direction Sensing in Dictyostelium discoideum: Ras Activity and Symmetry Breaking Driven by a Gβγ-Mediated, Gα2-Ric8 -- Dependent Signal Transduction Network

Chemotaxis is a dynamic cellular process, comprised of direction sensing, polarization and locomotion, that leads to the directed movement of eukaryotic cells along extracellular gradients. As a primary step in the response of an individual cell to a spatial stimulus, direction sensing has attracted numerous theoretical treatments aimed at explaining experimental observations in a variety of cell types. Here we propose a new model of direction sensing based on experiments using Dictyostelium discoideum (Dicty). The model is built around a reaction-diffusion-translocation system that involves three main component processes: a signal detection step based on G-protein-coupled receptors (GPCR) for cyclic AMP (cAMP), a transduction step based on a heterotrimetic G protein Gα2βγ, and an activation step of a monomeric G-protein Ras. The model can predict the experimentally-observed response of cells treated with latrunculin A, which removes feedback from downstream processes, under a variety of stimulus protocols. We show that [Formula: see text] cycling modulated by Ric8, a nonreceptor guanine exchange factor for [Formula: see text] in Dicty, drives multiple phases of Ras activation and leads to direction sensing and signal amplification in cAMP gradients. The model predicts that both [Formula: see text] and Gβγ are essential for direction sensing, in that membrane-localized [Formula: see text], the activated GTP-bearing form of [Formula: see text], leads to asymmetrical recruitment of RasGEF and Ric8, while globally-diffusing Gβγ mediates their activation. We show that the predicted response at the level of Ras activation encodes sufficient 'memory' to eliminate the 'back-of-the wave' problem, and the effects of diffusion and cell shape on direction sensing are also investigated. In contrast with existing LEGI models of chemotaxis, the results do not require a disparity between the diffusion coefficients of the Ras activator GEF and the Ras inhibitor GAP. Since the signal pathways we study are highly conserved between Dicty and mammalian leukocytes, the model can serve as a generic one for direction sensing.

A Quantitative Model of the GIRK1/2 Channel Reveals That Its Basal and Evoked Activities Are Controlled by Unequal Stoichiometry of Gα and Gβγ

G protein-gated K⁺ channels (GIRK; Kir3), activated by Gβγ subunits derived from G_i/o proteins, regulate heartbeat and neuronal excitability and plasticity. Both neurotransmitter-evoked (I_evoked) and neurotransmitter-independent basal (I_basal) GIRK activities are physiologically important, but mechanisms of I_basal and its relation to I_evoked are unclear. We have previously shown for heterologously expressed neuronal GIRK1/2, and now show for native GIRK in hippocampal neurons, that I_basal and I_evoked are interrelated: the extent of activation by neurotransmitter (activation index, R_a) is inversely related to I_basal. To unveil the underlying mechanisms, we have developed a quantitative model of GIRK1/2 function. We characterized single-channel and macroscopic GIRK1/2 currents, and surface densities of GIRK1/2 and Gβγ expressed in Xenopus oocytes. Based on experimental results, we constructed a mathematical model of GIRK1/2 activity under steady-state conditions before and after activation by neurotransmitter. Our model accurately recapitulates I_basal and I_evoked in Xenopus oocytes, HEK293 cells and hippocampal neurons; correctly predicts the dose-dependent activation of GIRK1/2 by coexpressed Gβγ and fully accounts for the inverse I_basal-R_a correlation. Modeling indicates that, under all conditions and at different channel expression levels, between 3 and 4 Gβγ dimers are available for each GIRK1/2 channel. In contrast, available Gα_i/o decreases from ~2 to less than one Gα per channel as GIRK1/2's density increases. The persistent Gβγ/channel (but not Gα/channel) ratio support a strong association of GIRK1/2 with Gβγ, consistent with recruitment to the cell surface of Gβγ, but not Gα, by GIRK1/2. Our analysis suggests a maximal stoichiometry of 4 Gβγ but only 2 Gα_i/o per one GIRK1/2 channel. The unique, unequal association of GIRK1/2 with G protein subunits, and the cooperative nature of GIRK gating by Gβγ, underlie the complex pattern of basal and agonist-evoked activities and allow GIRK1/2 to act as a sensitive bidirectional detector of both Gβγ and Gα.

Many neurotransmitters and hormones inhibit the electric activity of excitable cells (such as cardiac cells and neurons) by activating a K⁺ channel, GIRK (G protein-gated Inwardly Rectifying K⁺ channel). GIRK channels also possess constitutive “basal” activity which contributes to regulation of neuronal and cardiac excitability and certain disorders, but the mechanism of this activity and its interrelation with the neurotransmitter-evoked activity are poorly understood. In this work we show that key features of basal and neurotransmitter-evoked activities are similar in cultured hippocampal neurons and in two model systems (mammalian HEK293 cells and Xenopus oocytes). Using experimental data of the neuronal GIRK1/2 channel function upon changes in GIRK and G protein concentrations, we constructed a mathematical model that quantitatively accounts for basal and evoked activity, and for the inverse correlation between the two. Our analysis suggests a novel and unexpected mechanism of interaction of GIRK1/2 with the G protein subunits, where the tetrameric GIRK channel can assemble with 4 molecules of the Gβγ subunits but only 2 molecules of Gα. GIRK is a prototypical effector of Gβγ, and the unequal stoichiometry of interaction with G protein subunits may have general implications for G protein signaling.

The structure of dynamic GPCR signaling networks

G-protein-coupled receptors (GPCRs) stimulate signaling networks that control a variety of critical physiological processes. Static information on the map of interacting signaling molecules at the basis of many cellular processes exists, but little is known about the dynamic operation of these networks. Here we focus on two questions. First, Is the network architecture underlying GPCR-activated cellular processes unique in comparison with others such as transcriptional networks? We discuss how spatially localized GPCR signaling requires uniquely organized networks to execute polarized cell responses. Second, What approaches overcome challenges in deciphering spatiotemporally dynamic networks that govern cell behavior? We focus on recently developed microfluidic and optical approaches that allow GPCR signaling pathways to be triggered and perturbed with spatially and temporally variant input while simultaneously visualizing molecular and cellular responses. When integrated with mathematical modeling, these approaches can help identify design principles that govern cell responses to extracellular signals. We outline why optical approaches that allow the behavior of a selected cell to be orchestrated continually are particularly well suited for probing network organization in single cells.

Heterotrimeric Go protein links Wnt-Frizzled signaling with ankyrins to regulate the neuronal microtubule cytoskeleton

Drosophila neuromuscular junctions (NMJs) represent a powerful model system with which to study glutamatergic synapse formation and remodeling. Several proteins have been implicated in these processes, including components of canonical Wingless (Drosophila Wnt1) signaling and the giant isoforms of the membrane-cytoskeleton linker Ankyrin 2, but possible interconnections and cooperation between these proteins were unknown. Here, we demonstrate that the heterotrimeric G protein Go functions as a transducer of Wingless-Frizzled 2 signaling in the synapse. We identify Ankyrin 2 as a target of Go signaling required for NMJ formation. Moreover, the Go-ankyrin interaction is conserved in the mammalian neurite outgrowth pathway. Without ankyrins, a major switch in the Go-induced neuronal cytoskeleton program is observed, from microtubule-dependent neurite outgrowth to actin-dependent lamellopodial induction. These findings describe a novel mechanism regulating the microtubule cytoskeleton in the nervous system. Our work in Drosophila and mammalian cells suggests that this mechanism might be generally applicable in nervous system development and function.

A model of cell biological signaling predicts a phase transition of signaling and provides mathematical formulae

A biological signal is transmitted by interactions between signaling molecules in the cell. To date, there have been extensive studies regarding signaling pathways using numerical simulation of kinetic equations that are based on equations of continuity and Fick’s law. To obtain a mathematical formulation of cell signaling, we propose a stability kinetic model of cell biological signaling of a simple two-parameter model based on the kinetics of the diffusion-limiting step. In the present model, the signaling is regulated by the binding of a cofactor, such as ATP. Non-linearity of the kinetics is given by the diffusion fluctuation in the interaction between signaling molecules, which is different from previous works that hypothesized autocatalytic reactions. Numerical simulations showed the presence of a critical concentration of the cofactor beyond which the cell signaling molecule concentration is altered in a chaos-like oscillation with frequency, which is similar to a discontinuous phase transition in physics. Notably, we found that the frequency is given by the logarithm function of the difference of the outside cofactor concentration from the critical concentration. This implies that the outside alteration of the cofactor concentration is transformed into the oscillatory alteration of cell inner signaling. Further, mathematical stability kinetic analysis predicted a discontinuous dynamic phase transition in the critical state at which the cofactor concentration is equivalent to the critical concentration. In conclusion, the present model illustrates a unique feature of cell signaling, and the stability analysis may provide an analytical framework of the cell signaling system and a novel formulation of biological signaling.

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.

Double suppression of the Gα protein activity by RGS proteins

Regulator of G protein signaling (RGS) proteins accelerate GTP hydrolysis on G protein α subunits, restricting their activity downstream from G protein-coupled receptors. Here we identify Drosophila Double hit (Dhit) as a dual RGS regulator of Gαo. In addition to the conventional GTPase-activating action, Dhit possesses the guanine nucleotide dissociation inhibitor (GDI) activity, slowing the rate of GTP uptake by Gαo; both activities are mediated by the same RGS domain. These findings are recapitulated using homologous mammalian Gαo/i proteins and RGS19. Crystal structure and mutagenesis studies provide clues into the molecular mechanism for this unprecedented GDI activity. Physiologically, we confirm this activity in Drosophila asymmetric cell divisions and HEK293T cells. We show that the oncogenic Gαo mutant found in breast cancer escapes this GDI regulation. Our studies identify Dhit and its homologs as double-action regulators, inhibiting Gαo/i proteins both through suppression of their activation and acceleration of their inactivation through the single RGS domain.

G-protein signaling leverages subunit-dependent membrane affinity to differentially control βγ translocation to intracellular membranes

Activation of G-protein heterotrimers by receptors at the plasma membrane stimulates βγ-complex dissociation from the α-subunit and translocation to internal membranes. This intermembrane movement of lipid-modified proteins is a fundamental but poorly understood feature of cell signaling. The differential translocation of G-protein βγ-subunit types provides a valuable experimental model to examine the movement of signaling proteins between membranes in a living cell. We used live cell imaging, mathematical modeling, and in vitro measurements of lipidated fluorescent peptide dissociation from vesicles to determine the mechanistic basis of the intermembrane movement and identify the interactions responsible for differential translocation kinetics in this family of evolutionarily conserved proteins. We found that the reversible translocation is mediated by the limited affinity of the βγ-subunits for membranes. The differential kinetics of the βγ-subunit types are determined by variations among a set of basic and hydrophobic residues in the γ-subunit types. G-protein signaling thus leverages the wide variation in membrane dissociation rates among different γ-subunit types to differentially control βγ-translocation kinetics in response to receptor activation. The conservation of primary structures of γ-subunits across mammalian species suggests that there can be evolutionary selection for primary structures that confer specific membrane-binding affinities and consequent rates of intermembrane movement.

Yellow submarine of the Wnt/Frizzled signaling: submerging from the G protein harbor to the targets

The Wnt/Frizzled signaling pathway plays multiple functions in animal development and, when deregulated, in human disease. The G-protein coupled receptor (GPCR) Frizzled and its cognate heterotrimeric Gi/o proteins initiate the intracellular signaling cascades resulting in cell fate determination and polarization. In this review, we summarize the knowledge on the ligand recognition, biochemistry, modifications and interacting partners of the Frizzled proteins viewed as GPCRs. We also discuss the effectors of the heterotrimeric Go protein in Frizzled signaling. One group of these effectors is represented by small GTPases of the Rab family, which amplify the initial Wnt/Frizzled signal. Another effector is the negative regulator of Wnt signaling Axin, which becomes deactivated in response to Go action. The discovery of the GPCR properties of Frizzled receptors not only provides mechanistic understanding to their signaling pathways, but also paves new avenues for the drug discovery efforts.

The Wnt/Frizzled GPCR signaling pathway

G protein-coupled receptors (GPCRs) represent the biggest transmembrane receptor family. The Frizzled group of GPCRs is evolutionarily conserved and serves to transduce signals from the Wnt-type lipoglycoprotein growth factors. The Wnt/Frizzled signaling cascades are repeatedly used during animal development and are mostly silent in the adult. Improper activation of these cascades, e.g. through somatic mutation, underlies cancer development in various tissues. Our research over the past years has identified the trimeric G proteins as crucial transducers of the Wnt/Frizzled cascades in insect and mammalian cells. The current mini-review summarizes our findings on the role of G proteins in Wnt/Frizzled signaling, as well as on identification of other signaling intermediates in this physiologically and pathologically important type of intracellular signal transduction.

CB2 cannabinoid receptor targets mitogenic Gi protein-cyclin D1 axis in osteoblasts

CB2 is a Gi protein-coupled receptor activated by endo- and phytocannabinoids, thus inhibiting stimulated adenylyl cyclase activity. CB2 is expressed in bone cells and Cb2 null mice show a marked age-related bone loss. CB2-specific agonists both attenuate and rescue ovariectomy-induced bone loss. Activation of CB2 stimulates osteoblast proliferation and bone marrow derived colony-forming units osteoblastic. Here we show that selective and nonselective CB2 agonists are mitogenic in MC3T3 E1 and newborn mouse calvarial osteoblastic cultures. The CB2 mitogenic signaling depends critically on the stimulation of Erk1/2 phosphorylation and de novo synthesis of MAP kinase-activated protein kinase 2 (Mapkapk2) mRNA and protein. Further downstream, CB2 activation enhances CREB transcriptional activity and cyclin D1 mRNA expression. The CB2-induced stimulation of CREB and cyclin D1 is inhibitable by pertussis toxin, the MEK-Erk1/2 inhibitors PD098059 and U0126, and Mapkapk2 siRNA. These data demonstrate that in osteoblasts CB2 targets a Gi protein-cyclin D1 mitogenic axis. Erk1/2 phosphorylation and Mapkapk2 protein synthesis are critical intermediates in this axis.

Competing activities of heterotrimeric G proteins in Drosophila wing maturation

Drosophila genome encodes six alpha-subunits of heterotrimeric G proteins. The Galphas alpha-subunit is involved in the post-eclosion wing maturation, which consists of the epithelial-mesenchymal transition and cell death, accompanied by unfolding of the pupal wing into the firm adult flight organ. Here we show that another alpha-subunit Galphao can specifically antagonize the Galphas activities by competing for the Gbeta13F/Ggamma1 subunits of the heterotrimeric Gs protein complex. Loss of Gbeta13F, Ggamma1, or Galphas, but not any other G protein subunit, results in prevention of post-eclosion cell death and failure of the wing expansion. However, cell death prevention alone is not sufficient to induce the expansion defect, suggesting that the failure of epithelial-mesenchymal transition is key to the folded wing phenotypes. Overactivation of Galphas with cholera toxin mimics expression of constitutively activated Galphas and promotes wing blistering due to precocious cell death. In contrast, co-overexpression of Gbeta13F and Ggamma1 does not produce wing blistering, revealing the passive role of the Gbetagamma in the Galphas-mediated activation of apoptosis, but hinting at the possible function of Gbetagamma in the epithelial-mesenchymal transition. Our results provide a comprehensive functional analysis of the heterotrimeric G protein proteome in the late stages of Drosophila wing development.

Cellular modelling: experiments and simulation to develop a physiological model of G-protein control of muscarinic K+ channels in mammalian atrial cells

The first model of G-protein-K(ACh) channel interaction was developed 14 years ago and then expanded to include both the receptor-G-protein cycle and G-protein-K(ACh) channel interaction. The G-protein-K(ACh) channel interaction used the Monod-Wyman-Changeux allosteric model with the idea that one K(ACh) channel is composed of four subunits, each of which binds one active G-protein subunit (G(betagamma)). The receptor-G-protein cycle used a previous model to account for the steady-state relationship between K(ACh) current and intracellular guanosine-5-triphosphate at various extracellular concentrations of acetylcholine (ACh). However, simulations of the activation and deactivation of K(ACh) current upon ACh application or removal were much slower than experimental results. This clearly indicates some essential elements were absent from the model. We recently found that regulators of G-protein signalling are involved in the control of K(ACh) channel activity. They are responsible for the voltage-dependent relaxation behaviour of K(ACh) channels. Based on this finding, we have improved the receptor-G-protein cycle model to reproduce the relaxation behaviour. With this modification, the activation and deactivation of K(ACh) current in the model are much faster and now fall within physiological ranges.

Snazer: the simulations and networks analyzer

BACKGROUND

Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.From a pure dynamical perspective, simulation represents a relevant way-out. Many simulator tools capitalized on the "noisy" behavior of some systems and used formal models to represent cellular activities as temporal trajectories. Statistical methods have been applied to a fairly large number of replicated trajectories in order to infer knowledge.A tool which both graphically manipulates reactive models and deals with sets of simulation time-course data by aggregation, interpretation and statistical analysis is missing and could add value to simulators.

RESULTS

We designed and implemented Snazer, the simulations and networks analyzer. Its goal is to aid the processes of visualizing and manipulating reactive models, as well as to share and interpret time-course data produced by stochastic simulators or by any other means.

CONCLUSIONS

Snazer is a solid prototype that integrates biological network and simulation time-course data analysis techniques.

Design of versatile biochemical switches that respond to amplitude, duration, and spatial cues

Cells often mount ultrasensitive (switch-like) responses to stimuli. The design principles underlying many switches are not known. We computationally studied the switching behavior of GTPases, and found that this first-order kinetic system can show ultrasensitivity. Analytical solutions indicate that ultrasensitive first-order reactions can yield switches that respond to signal amplitude or duration. The three-component GTPase system is analogous to the physical fermion gas. This analogy allows for an analytical understanding of the functional capabilities of first-order ultrasensitive systems. Experiments show amplitude- and time-dependent Rap GTPase switching in response to Cannabinoid-1 receptor signal. This first-order switch arises from relative reaction rates and the concentrations ratios of the activator and deactivator of Rap. First-order ultrasensitivity is applicable to many systems where threshold for transition between states is dependent on the duration, amplitude, or location of a distal signal. We conclude that the emergence of ultrasensitivity from coupled first-order reactions provides a versatile mechanism for the design of biochemical switches.

Coordinate regulation of G protein signaling via dynamic interactions of receptor and GAP

Signal output from receptor–G-protein–effector modules is a dynamic function of the nucleotide exchange activity of the receptor, the GTPase-accelerating activity of GTPase-activating proteins (GAPs), and their interactions. GAPs may inhibit steady-state signaling but may also accelerate deactivation upon removal of stimulus without significantly inhibiting output when the receptor is active. Further, some effectors (e.g., phospholipase C-β) are themselves GAPs, and it is unclear how such effectors can be stimulated by G proteins at the same time as they accelerate G protein deactivation. The multiple combinations of protein–protein associations and interacting regulatory effects that allow such complex behaviors in this system do not permit the usual simplifying assumptions of traditional enzyme kinetics and are uniquely subject to systems-level analysis. We developed a kinetic model for G protein signaling that permits analysis of both interactive and independent G protein binding and regulation by receptor and GAP. We evaluated parameters of the model (all forward and reverse rate constants) by global least-squares fitting to a diverse set of steady-state GTPase measurements in an m1 muscarinic receptor–G_q–phospholipase C-β1 module in which GTPase activities were varied by ∼10⁴-fold. We provide multiple tests to validate the fitted parameter set, which is consistent with results from the few previous pre-steady-state kinetic measurements. Results indicate that (1) GAP potentiates the GDP/GTP exchange activity of the receptor, an activity never before reported; (2) exchange activity of the receptor is biased toward replacement of GDP by GTP; (3) receptor and GAP bind G protein with negative cooperativity when G protein is bound to either GTP or GDP, promoting rapid GAP binding and dissociation; (4) GAP indirectly stabilizes the continuous binding of receptor to G protein during steady-state GTPase hydrolysis, thus further enhancing receptor activity; and (5) receptor accelerates GDP/GTP exchange primarily by opening an otherwise closed nucleotide binding site on the G protein but has minimal effect on affinity (K_assoc = k_assoc/k_dissoc) of G protein for nucleotide. Model-based simulation explains how GAP activity can accelerate deactivation >10-fold upon removal of agonist but still allow high signal output while the receptor is active. Analysis of GTPase flux through distinct reaction pathways and consequent accumulation of specific GTPase cycle intermediates indicate that, in the presence of a GAP, the receptor remains bound to G protein throughout the GTPase cycle and that GAP binds primarily during the GTP-bound phase. The analysis explains these behaviors and relates them to the specific regulatory phenomena described above. The work also demonstrates the applicability of appropriately data-constrained system-level analysis to signaling networks of this scale.

Throughout the eukaryotes, G proteins convey information from receptors for diverse stimuli—neurotransmitters, hormones, light, odors, and pheromones—to intracellular regulatory proteins collectively known as effectors. G proteins function by transiting a dynamic cycle of activation and deactivation. Receptors accelerate activation, which allows G proteins to regulate effectors, and receptors thus increase signal output. GTPase-activating proteins, GAPs, accelerate deactivation. GAPs can thus attenuate signaling, but GAPs can also accelerate signal termination when stimulus is removed without inhibiting signal output while stimulus is present. Surprisingly, some effectors are also GAPs for the G proteins that activate them, essentially turning off their activator. We developed a mathematical model that describes control of G protein signaling by receptor and GAP and used experimental data to determine its important parameters. We show that GAPs actually potentiate G protein activation by receptor, a previously unsuspected effect. Further, GAPs indirectly stabilize receptor–G protein binding during stimulation, which we had previously predicted based on inconsistencies among other experimental results. The present results elucidate how GAPs can independently control signaling kinetics and amplitude and thus clarify how effectors can both respond to G proteins and act as G protein GAPs.

Bistability and oscillations in the Huang-Ferrell model of MAPK signaling

Physicochemical models of signaling pathways are characterized by high levels of structural and parametric uncertainty, reflecting both incomplete knowledge about signal transduction and the intrinsic variability of cellular processes. As a result, these models try to predict the dynamics of systems with tens or even hundreds of free parameters. At this level of uncertainty, model analysis should emphasize statistics of systems-level properties, rather than the detailed structure of solutions or boundaries separating different dynamic regimes. Based on the combination of random parameter search and continuation algorithms, we developed a methodology for the statistical analysis of mechanistic signaling models. In applying it to the well-studied MAPK cascade model, we discovered a large region of oscillations and explained their emergence from single-stage bistability. The surprising abundance of strongly nonlinear (oscillatory and bistable) input/output maps revealed by our analysis may be one of the reasons why the MAPK cascade in vivo is embedded in more complex regulatory structures. We argue that this type of analysis should accompany nonlinear multiparameter studies of stationary as well as transient features in network dynamics.

Molecular studies of cell communication systems lead to models with multiple free parameters. Analysis of dynamical behavior of these models presents considerable challenge. We have developed a computational approach for the efficient exploration of dynamic behavior in such models and applied this method to the model of the Mitogen Activated Protein Kinase cascade, a signaling network conserved in all eukaryotes. Previous analysis of this model suggested that it works as a reversible switch. We have shown that it can also function as an irreversible switch and as a clock.