The gastrin and cholecystokinin receptors mediated signaling network: a scaffold for data analysis and new hypotheses on regulatory mechanisms
BMC Systems Biology volume 9, Article number: 40 (2015)
The gastrointestinal peptide hormones cholecystokinin and gastrin exert their biological functions via cholecystokinin receptors CCK1R and CCK2R respectively. Gastrin, a central regulator of gastric acid secretion, is involved in growth and differentiation of gastric and colonic mucosa, and there is evidence that it is pro-carcinogenic. Cholecystokinin is implicated in digestion, appetite control and body weight regulation, and may play a role in several digestive disorders.
We performed a detailed analysis of the literature reporting experimental evidence on signaling pathways triggered by CCK1R and CCK2R, in order to create a comprehensive map of gastrin and cholecystokinin-mediated intracellular signaling cascades. The resulting signaling map captures 413 reactions involving 530 molecular species, and incorporates the currently available knowledge into one integrated signaling network. The decomposition of the signaling map into sub-networks revealed 18 modules that represent higher-level structures of the signaling map. These modules allow a more compact mapping of intracellular signaling reactions to known cell behavioral outcomes such as proliferation, migration and apoptosis. The integration of large-scale protein-protein interaction data to this literature-based signaling map in combination with topological analyses allowed us to identify 70 proteins able to increase the compactness of the map. These proteins represent experimentally testable hypotheses for gaining new knowledge on gastrin- and cholecystokinin receptor signaling. The CCKR map is freely available both in a downloadable, machine-readable SBML-compatible format and as a web resource through PAYAO (http://sblab.celldesigner.org:18080/Payao11/bin/).
We have demonstrated how a literature-based CCKR signaling map together with its protein interaction extensions can be analyzed to generate new hypotheses on molecular mechanisms involved in gastrin- and cholecystokinin-mediated regulation of cellular processes.
Gastrin and cholecystokinin (CCK) are gastrointestinal peptide hormones that share a common C-terminal pentapeptide amide . Gastrin, produced in G-cells of the gastric antrum, is the central regulator of gastric acid secretion but also regulates growth and differentiation of gastric and colonic mucosa . CCK, produced primarily in I-cells of the small intestine, is involved in physiological processes such as digestion, appetite control and body weight regulation . The scientific interest in these hormones is further strengthened by their roles in several diseases. Indeed, CCK has been implicated in acute pancreatitis [4–6], obesity [7, 8], irritable bowel syndrome  and gall bladder disease [10, 11]. Gastrin is known to be pro-carcinogenic, affecting proliferation, angiogenesis and apoptosis , and is a risk co-factor for gastric carcinogenesis and atrophy upon Helicobacter pylori infection [12, 13]. In order to efficiently study and understand the molecular mechanisms triggered by gastrin and cholecystokinin, detailed knowledge concerning the signaling pathways they regulate is paramount.
Information concerning intracellular signaling is commonly retrieved from databases such as Reactome  and KEGG . However, none of these resources currently specify which specific molecular events are known to take place in response to gastrin or CCK. Researchers in need of such knowledge therefore must spend significant time reviewing current literature in order to gain an exhaustive and up-to-date understanding of the signaling network. A comprehensive map of gastrin and CCK intracellular signaling pathways would significantly assist in the study of normal or aberrant cholecystokinin receptor (CCKR) signaling.
In the past decade, several manually constructed maps of signaling events have been published [16–25], each providing solid foundations for a systems understanding of the signaling mechanisms. The present work extends this approach to the domain of CCKR signaling by providing a comprehensive literature-based CCKR signaling map that comprises 530 molecular species and 431 reactions, considerably extending previously compiled knowledge on CCK2R signaling [10, 26] including CCK1R downstream events.
Partitioning the total CCKR signaling map into sub-networks using the BiNoM tool  resulted in 18 modules that coordinate with each other to elicit the diverse intracellular signaling responses to gastrin and/or CCK. Finally, we used the CCKR map as a scaffold for protein-protein interaction (PPI) data integration assisted by PathExpand , in order to predict novel components of the signaling network. This resulted in the identification of 70 new proteins tightly connected to the CCKR signaling map, making them prime candidates for future experimental work aimed at further extending knowledge on regulation of CCKR mediated signaling mechanisms.
Construction of the CCKR map from literature
The CCKR map was constructed using CellDesigner 4.2, a structured diagram editor for drawing gene-regulatory and biochemical networks, following the Systems Biology Graphical Notation (SBGN) standard for process diagrams  and the Systems Biology Mark-up Language (SBML) for model representation . The MIRIAM (Minimum Information Requested In the Annotation of Models) guidelines were followed to characterize each species in the map .
Knowledge encoded in the CCKR map was obtained from scientific publications that were identified by searching for different combinations of cholecystokinin (CCK)/CCK1R and gastrin (G-17)/CCK2R in PubMed or through various literature mining tools, e.g. LitInspector (http://www.genomatix.de/solutions/genomatix-software-suite.html)  and iHOP (http://www.ihop-net.org/UniPub/iHOP/) . We adopted two main criteria as guidelines for including a scientific paper for the information extraction, namely that it must contain:
evidence that the reported signaling event is mediated by the specific interaction of CCK or gastrin (G-17) with their receptors CCK1R and CCK2R.
signaling information to allow for linkage of a new CCKR map component to its upstream and/or downstream regulators in CCKR signaling.
CellDesigner species and reaction “note” features were used to record PubMed IDs (PMID), cell-type specific information for each reaction, and the interacting components in the CCKR map.
Final curation and quality control was done in a collaborative effort involving five different research group members who collectively used the community curation platform PAYAO (http://www.payaologue.org) , enabling efficient exchange of comments and tags. Consensus and critical comments from each annotator about the precise representation of reactions, components, and their cellular localization were discussed and implemented (Fig. 1). Finally, the CCKR map was published with an open source license for the whole scientific community through PAYAO  (map available at: http://sblab.celldesigner.org:18080/Payao11/bin/). The CCKR map is also available in the SBML data exchange format (Additional file 1), and we are also in the process of submitting this map to the PANTHER database .
We welcome comments and tags from a community of curators to continue increasing the quality of this CCKR signaling map and will keep it up to date with our growing biological understanding.
Analysis of the CCKR map
Network topology studies were performed by importing the CellDesigner-generated SBML file of the CCKR map to Cytoscape version 2.8  using the BiNoM plugin , after removing connections downstream of the transcription factor layer. BiNoM considers both ‘reaction’ and ‘species’ of a CellDesigner map as a node. The Cytoscape version of the CCKR map consisted of 807 nodes (475 species and 332 reactions) with 963 edges. Node degree and closeness centrality were calculated using the ‘Network Analysis’ plugin  in Cytoscape, assuming the network to be undirected. The number of nodes connected directly to a particular node defines its node degree k, and the node degrees of all nodes in the network represent the ‘degree distribution’ of the network. Nodes with degree k > 5 were termed ‘hubs’.
BiNoM decomposition of CCKR map into modules
We used the ‘prune the graph’ function of BiNoM to automatically separate the strongly connected component (SCC), i.e. the central cyclic motif of the map, from more loosely connected upstream and downstream species. The SCC was further decomposed into smallest sub-networks with the function ‘extract material components’. Next, some manual curation work was performed: sub-networks with 50 % or more overlapping nodes were clustered together while the large sub-networks were decomposed further until each sub-network contained a unique central cyclic motif. Any unconnected nodes were removed from the sub-networks, yielding a total of 18 sub-networks or modules. Next, these 18 modules were merged together and the resulting network was compared with the initial map to check for completeness in terms of any missing interactions or nodes. The complete map and its 18 decomposed modules are available as a Cytoscape session file (Additional file 2).
Protein-protein interaction based expansion of the CCKR map
PPI data were downloaded using PSICQUIC (all databases, version June 2012), and filtered for binary physical interactions based on PSI-MI controlled vocabulary experimental method descriptions, following the procedure in  (Charles E. Chapple, personal communication). We then tested each of the 4119 proteins found to interact with CCKR signaling map proteins with the PathExpand method  to seek for tightly linked protein interactors that enhance the network compactness of the CCKR map by leading to higher node degree, betweenness and average local clustering as well as decreasing the shortest path lengths (for details see ) .
CCK and gastrin impinge on cellular functioning by binding to two different G protein-coupled receptors, CCK1R and CCK2R, respectively, located at the surface of multiple cell types in peripheral organs such as the gastrointestinal tract, the pancreas, and the gall bladder . Today, no comprehensive resource exists that compiles current knowledge on CCKR activated signaling pathways. Gastrin has a strong preference for CCK2R, while CCK can activate both receptors with similar affinities . Most cell types responsive to one or both peptide hormones express only one CCK-receptor variant. However, a range of normal and cancer cells in whole organisms as well as model cell lines (for instance the rat pancreatic acinar cell derived cell-line AR42J ) express both CCK1R and CCK2R.
The CCKR signaling map
We present a CCKR signaling map built with biological inferences from more than 250 scientific publications (including original articles and reviews), and based on experiments performed in 37 different cell lines representing a wide array of cell types (Fig. 2, Methods).
The CCKR map encompasses 199 unique proteins and their relationships to complexes, genes, and small molecules (including some inhibitors, e.g. rapamycin, wortmannin), describing a total of 530 molecular species and 431 reactions experimentally documented to play a role in CCKR signaling. The reactions include molecular state transitions (e.g. activation, phosphorylation), transport, heterodimer associations/dissociations as well as transcriptional regulation of target genes by transcription factors, which are reported to be associated with the regulation of gastrointestinal hormone responsive biological processes (Table 1).
Figure 2a displays a detailed view of the map depicting the two ligands gastrin and CCK, their biogenesis and processing, their binding to the two G protein coupled receptors CCK1R and CCK2R and the ensuing signal transduction pathways including activation of transcription factors and downstream target genes reported to be triggered in response to gastrin and/or CCK. Moreover, a number of cellular processes influenced by specific CCKR signaling events are indicated, including proliferation, migration, differentiation, anti-apoptosis and inflammation. Figure 2b is an aid for the reader to navigate in Fig. 2a, as it indicates the position of major pathways involved, including protein kinase C (PKC)-dependent activation of MAPK cascades, PI3K-mTOR signaling, protein kinase A (PKA)-dependent pathway, β-catenin, and Rho-GTPase-Bcl cascades.
Network topology analysis of the CCKR map indicated scale-free characteristics, with the majority of the proteins connected to only few other network proteins (Fig. 3). The 6 most highly connected proteins include four protein kinases AKT1, SRC, PKC and PAK1, and the small GTPase HRAS (Fig. 3, inserted table).
Signaling pathways shared by CCK1R and CCK2R
Downstream signaling mechanisms shared by CCK1R and CCK2R include the trimeric guanine nucleotide binding alpha q protein (G-protein αq), PKC dependent phosphorylation of adaptor protein SHC - transforming protein C, and its association with GRB2/SOS1 leading to activation of the HRAS/RAF/MAPK1/3 cascade [10, 41, 42]. Other shared pathways involve MAP3K11-mediated regulation of MAPK8, −9 and −10 and p38MAPK (MAPK14), PRKD1 (in PKC signaling), PI3K, AKT1, FAK1, JAK2-STAT3 and Rho GTPase [10, 41–44]. Both receptors activate PKC isoforms PKCα, −δ, −ε, −θ and −ζ [10, 45–47]. Transcription factors reported downstream of both CCK1R and CCK2R receptors include NFκB, CREB1, ELK1 and AP1.
CCK1R-specific signaling pathways
Two trimeric G-proteins appear to be regulated only by CCK1R. One is G alpha S (Gαs) , which leads to Protein kinase A activation via adenylate cyclase catalyzed cAMP production, and the other is G alpha 13 (Gα13) , involved in downstream activation of RHOA . The NOS1 signal transduction pathway downstream of CCK1R [50, 51] regulates Ca2+ signaling pathways by opening ryanodine receptors and two-pore channels that release calcium from the endoplasmic reticulum and endolysosomes, respectively [52–54]. Other signaling cascades, so far only reported for CCK1R, include the Ca2+/calmodulin (CaM) dependent activation of CaMKIV and calcineurin (CaN) – NFAT1 signaling pathways [55–57].
CCK2R-specific signaling pathways
CCK2R activates EGFR via PKC activated MMP3, which cleaves membrane-attached pro-HBEGF into mature HBEGF [58, 59]. PKC isoforms PKC-β and PKC-η have been reported only downstream of CCK2R signaling [60, 61]. CCKR2 specific activation of β-catenin and E-cadherin is mediated by PAK1 [62, 63], and CCK2R specific modulation of BCL-protein family signaling regulates mitochondrial cytochrome C release [64, 65]. CCK2R, but not CCK1R, is reported to activate MAPK7 , an upstream regulator of transcription factors MEF-B,-C and D, and the PKC-η target PRKD2 , which enhances nuclear export of HDAC7 thereby relieving transcriptional repression of target genes such as NR4A1 .
Segmentation of the CCKR signaling map into modules
The complete CCKR signaling map (Fig. 2) is very large. We therefore set out to identify structural and functional subdomains of this map. We used the BiNoM software to import the CCKR map in Cytoscape and built a modular view of the CCKR map similar to what has been previously reported for the E2F-Retinoblastoma signaling pathway [27, 68]. The resulting 18 modules (Methods, Table 2, Fig. 4a, Additional file 2) are fully based on the underlying detailed map and help to comprehend higher order map structure, navigate through the map and work on functional modules. Each of the modules represents a structural and functional signaling subunit, combining a set of closely coordinated molecular events concerning a particular protein or a protein complex (see e.g. the Rho GTPase module in Fig. 5a). Furthermore, the regulatory relations (activating/repressing) connecting the modules (Fig. 4a, black and red arrows) are directly derived from the relations encoded in the underlying detailed map (Fig. 2) and the modular view thus represents central decision-making aspects of CCKR signaling. Additional file 3 gives details for each module including the complete list of components and a discussion of their roles in the signaling mechanisms.
Gastrin-regulated cellular processes involve different sets of signaling modules
Depending on cell type and the state of cells, gastrin can induce different cellular outcomes, such as proliferation, migration and apoptosis. In order to investigate what is known about the involvement of each of the signaling pathway modules in these processes, we checked the scientific literature that we have used to build the CCKR map. We chose to focus on gastrin-mediated effects, because molecular mechanism studies of gastrin-regulated cellular responses are more abundantly reported in literature. Specifically, we recorded all instances where experimental evidence was reported indicating that the function of a specific signaling component had an effect on the cellular outcome in response to gastrin treatment. Based on this literature survey, involvement in cellular outcomes was indicated in Fig. 4a. While central modules such as PKC, AKT1, Rho GTPase, MAP3K11, MAPK1/3 and AP1 are reportedly involved in all three cellular outcomes, other signaling mechanisms are more specific in their effect, e.g. the BCL-module signaling is only observed to be involved in apoptosis (Fig. 4a). Below, we provide a brief description of signaling modules and components involved in each of the cellular outcomes to gastrin treatment.
Molecular mechanisms underlying gastrin-mediated proliferation involve regulation of protein synthesis and cell cycle. Protein synthesis is stimulated via the AKT1-module component mTOR triggering p70 S6 kinase [69, 70]. Gastrin-induced transcription of Cyclin D1, a central regulator of cell cycle progression, is mediated by JUN, FOS, CREB1, and TCF7L2 [71–75], which are components of the modules: AP1, ATF2 and β-catenin. The modular representation (Fig. 4a) shows that EGFR-associated signaling enhances gastrin-induced proliferation by feed-forward mechanisms involving SRC module components. The modular map (Fig. 4a) shows that the AKT1 module inhibits RAF1 and consequently MAPK1/3-module pathways. The molecular mechanism for this is AKT1 kinase-mediated phosphorylation of RAF1 . The modular map thus allows us to hypothesize that the activating role of AP1- and ATF2-submodule signaling in proliferation is more likely to proceed via SRC-MAP3K11 pathways, since these are not inhibited by AKT1.
Gastrin promotes cell migration by activating transcription of MMP7 and MMP9 [77, 78] via transcription factors SNAI1, β-catenin and JUN  represented in the β-catenin- (SNAI1, β-catenin)- and AP1-(JUN) modules. Cell adhesion, tightly linked to cell migration, is regulated through components in the FAK1/2 module (FAK1 and FAK2; Paxillin, CAS, and CRK [79, 80]), which is controlled via both PKC and SRC-modules. The module representation (Fig. 4a) indicates that PKC- and SRC-modules exert a positive feedback regulation on the FAK1/2 module.
Apoptosis is inhibited by gastrin via several mechanisms including BCL-mediated repression of pro-apoptotic caspases and AP1-activated expression of Clusterin [64, 65, 81, 82]. The modular representation reveals that these cellular responses are regulated both by PKC independent and PKC dependent mechanisms. This applies to NFκB and its downstream anti-apoptotic BIRC2 and BIRC3 target genes, which can be activated either directly by PKC or independently of PKC through the Rho GTPase module. Likewise, the AKT1-involvement in regulation of the BCL-module can be mediated by PKC-dependent mechanisms or independently of PKC by the CCK2R - Rho GTPase pathway. Activation of AP1, on the other hand, seems to be strictly dependent on PKC that mediates its effect via either RAF1-MAPK1/3 or SRC-MAP3K11 cascades.
We note that the AKT1-module both inhibits the RAF1-MAPK1/3-route to AP1-activation and enhances Rho GTPase activation of the BCL-module. Thus, AKT1 can potentially promote BCL-module apoptosis-regulating mechanisms and at the same time block MAPK1/3-mediated AP1-activation. In the latter configuration the cell would rely on MAP3K11 to bypass the inhibitory effect of AKT1 on AP1-mediated regulation of gene expression.
Extending the CCKR map with PPI data
The comprehensive CCKR signaling map (Fig. 2) has been constructed with a knowledge-driven approach based on molecular reactions and interactions that are well-documented by scientific evidence in the literature. This thereby inevitably leaves significant gaps concerning signaling events and mechanisms that have not yet been studied and evidenced. We have therefore exploited large-scale PPI data to assess how it can complement the CCKR signaling map by a data-driven strategy.
We identified 4119 proteins interacting with at least one CCKR signaling protein (Methods, Additional file 4). Among 199 unique proteins present in our CCKR signal transduction map, 146 were part of this large protein-protein interaction network. We then set out to identify protein interactors that could be of high interest as potential interactors or modulators of the complete CCKR signaling map (Fig. 4a). For this, we analyzed each of the 4119 proteins interacting with the CCKR map proteins using the PathExpand approach . PathExpand checks each interactor to see if it satisfies a number of topological criteria leading to increased compactness of the complete CCKR map. We identified a total of 102 proteins that qualified as PathExpand interactors, 32 of them being components of the original CCKR map. The 70 PathExpand interactors not present in the original CCKR map are listed in Table 3. Interestingly, the set of 70 new CCKR map candidates include 30 proteins that are not known to participate in any signaling pathways listed in the pathway databases KEGG, Reactome, PANTHER and Biocarta (Table 3). A GO term overrepresentation analysis  showed that the set of map candidates is enriched in molecular functions relating to protein kinases, protein phosphatases and GTPase-regulators, indicating that many of them could potentially regulate the CCKR pathway via phosphorylation-dephosphorylation mechanisms and by interfering with small GTPase signaling. Among these are the protein kinases STK4, CSK21, CSK22, ITPKA, FLT4, DDR2, KS6KA4 and MAPK15 and a high number of Dual specificity phosphatases (DUSP1, 2, 4, 5, 7, 9, 22) in addition to phosphatases PHLPP1 and PTPRR.
Inspection of the expanded networks resulting from the combination of PathExpand protein interactors with 7 CCKR modules (Fig. 4b, transcription factor centered modules excluded) shows that 11 of the protein interactors are linked to more than one CCKR module and can contribute both to PKC-independent (e.g. ARHGEF25, ARHGAP31, TBC1D3F) and PKC dependent (e.g. ICMT) signaling routes. In contrast, 21 protein interactors are linked to only one module, suggesting that they may act as preferential regulators of this module. For instance, the Mitogen activated kinase 15 (MAPK15) is a compactness increasing protein only for the SRC module, where it interacts with the two kinases SRC and CSK.
CCKR map and genome-scale PPI generate hypotheses for refinement of Rho GTPase module mechanisms
By taking a closer look at the gastrin regulated Rho GTPase module and its identified PPI extensions, we here discuss putative novel signaling mechanisms involved in gastrin-mediated regulation of the Rho GTPase signaling.
Central to the Rho GTPase module (Fig. 5a) is the activation of small GTPases, RHOA, RAC1, and CDC42 by guanine exchange factor proteins (GEFs), which trigger conversion of the inactive GDP-bound form of small GTPases to the active GTP-bound form. The GEFs in gastrin-mediated Rho GTPase signaling are Leukemia-associated Rho guanine-nucleotide exchange factor (LARG) and the trimeric receptor-associated G-protein complexes Gαq and Gα13. In addition, HRAS may also act to activate RHOA and CDC42. The small GTPases RHOA, RAC1, and CDC42 are documented to activate kinases ROCK1 and PAK1, which are both part of the gastrin-regulated Rho GTPase module, and of kinase AKT1, part of the gastrin-regulated AKT1 signaling module. The inactive form of the small GTPases is restored by GTPase-activating proteins (GAPs) that enhance hydrolyzation of the bound GTP. In gastrin-mediated signaling deactivation of RHOA and RAC1 is effectuated by GAPs, ARHGAP4 and RGS2, respectively.
Eight of the ten PathExpand protein interactors (Fig. 5b, Additional file 5) of Rho GTPase module components interact with one or more of the small GTPases RHOA, RAC1, and CDC42. All of these protein interactors are annotated with GO-terms indicating a known role in regulation of small GTPases, and four of them (IQGAP2, TBC1D3F, ARHGAP31, OPHN1) are known as GTPase activators. We suggest that these eight proteins are strong candidates for potentially novel modulators of gastrin-mediated Rho GTPase module signaling. Three of the proteins (TBC1D3F, ARHGAP25, and ICMT) are not yet present in any of the common pathway databases (see Table 3). The evidence presented here, indicating that the three proteins have tight interaction with multiple proteins involved in gastrin signaling, demonstrates that our approach is well suited to identify testable hypotheses, even for candidates for which knowledge in public databases is still sparse.
The PathExpand interactor BNIP1 (BCL2/adenovirus E1B 19 kDa interacting protein 1) interacts with two downstream BCL module components BCL2L1 and BCL2, both involved in effectuating the gastrin anti-apoptotic signal. Since BNIP1 is known to be involved in SNARE vesicular transport, our findings indicate that it may be an interesting candidate to follow up with investigating a potential role of vesicular transport in modulating BCL-linked signaling mechanisms in the gastrin signaling pathway. The PathExpand interactor programmed cell death 11 (PDCD11) interacting with the two subunits of NFκB transcription factor (NFκB1 and TF65) is known to be involved in maturation of ribosomal RNA but may also play a role in mRNA processing. Our results identify PDCD11 as a potential candidate to follow up in experimental analysis of molecular mechanisms involved in NFKB-mediated gene regulation in cellular responses to gastrin.
In the work presented here, we set out to build a comprehensive and well-annotated molecular interaction map to aid future studies involving gastrin or CCK. A mechanistic understanding of CCK1R and CCK2R signaling networks is essential for experiment design and data interpretation in biological systems involving gastrin- and CCK-regulated processes. Moreover, it may enhance the identification of therapeutic chemicals able to target disease, by using the map as a functional interaction diagram with components that modeling indicates as prime targets for perturbation, as described by Lee et al. , who successful designed a drug-induced rewiring of the ‘state’ of oncogenic signaling networks to maximize the susceptibility to anticancer drugs. The current map depicting molecular signaling mechanisms underlying CCK1R and CCK2R triggered cellular responses may provide a tool to guide further investigations into normo- and pathophysiological processes such as gastrin-regulated stomach mucosa homeostatis and gastrin-linked carcinogenesis [59, 64, 85, 86], and cholecystokinin-induced hypoplasia, cell regeneration and digestive enzyme secretion  as well as to identify potential intervention strategies for associated disease states.
Some pathway proteins and reactions of the CCKR signaling map are also described in the Reactome  and KEGG  databases. However, the CCKR map constitutes a vastly more comprehensive and integrated model, providing detailed signaling reactions linking the receptors CCK1R and CCK2R all the way to regulated genes and cellular responses. For comparison, the recently published gastrin signaling map  comprises approximately half of the molecules and only a quarter of the reactions presented in our CCKR map. Similarly, knowledge currently in Reactome covers less than 5 % of the pathway details presented here. The CCKR map therefore represents a significant increase in curated signaling information. The provided SBML version of the CCKR map can serve as a starting point to generate quantitative mathematical models  for simulation and prediction of cellular outcomes in response to perturbations of the network.
We enhance the applicability of the map for hypothesis generation by two central strategies. First, we provide a computationally modularized version of topologically and functionally connected meta-nodes. This modular view simplifies the complete CCKR map and provides for an improved, higher level comprehension of pathway regulatory aspects concerning cell fate decisions related to proliferation, migration and apoptosis. Secondly, we take advantage of public large scale PPI knowledge to predict new potential regulators of CCKR signaling, including 70 interactors that significantly enhance the compactness of the CCKR map , through tight direct and indirect interactions with map proteins. Although experimental validations are needed to confirm these new CCKR signaling mechanisms, they represent an important source of high quality hypotheses that may be taken as a first step to develop a better comprehension of CCKR pathway functionality.
Compared to other recently published computational approaches for high-throughput hypothesis generation [20, 89, 90], the advantage of our strategy is the combination of i) biological background knowledge encoded in the signaling map, including the modules, manually curated from literature reporting detailed experimental analyses of gastrin- and CCK-signaling, and ii) large-scale PPI information downloaded from available databases of interactions, filtered for binary physical interaction based on selected experimental detection methods.
Our work demonstrates how publicly available tools can aid in the manual curation of an extensive network topology, in order to build a foundation for a systems understanding of CCKR-mediated cellular responses. Furthermore, the integration of a comprehensive model of complex biological networks with genome scale data can provide new knowledge and hypotheses on molecular mechanisms underlying cellular processes. Further development of the resources presented here should be of high interest in translational research aimed at identifying new targets and biomarkers for treatment and diagnostics of gastrin- and/or cholecystokinin-related disease, including cancer.
Rehfeld JF. The New Biology of Gastrointestinal Hormones. Physiol Rev. 1998;78(4):1087–108.
Watson SA, Grabowska AM, El-Zaatari M, Takhar A. Gastrin - active participant or bystander in gastric carcinogenesis? Nat Rev Cancer. 2006;6(12):936–46.
Little TJ, Horowitz M, Feinle-Bisset C. Role of cholecystokinin in appetite control and body weight regulation. Obes Rev. 2005;6(4):297–306.
Saluja AK, Saluja M, Printz H, Zavertnik A, Sengupta A, Steer ML. Experimental pancreatitis is mediated by low-affinity cholecystokinin receptors that inhibit digestive enzyme secretion. Proc Natl Acad Sci U S A. 1989;86:8968–71.
Gukovsky I, Cheng JH, Nam KJ, Lee OT, Lugea A, Fischer L, et al. Phosphatidylinositide 3-kinase gamma regulates key pathologic responses to cholecystokinin in pancreatic acinar cells. Gastroenterology. 2004;126:554–66.
Dabrowski A, Grady T, Logsdon CD, Williams JA. Jun kinases are rapidly activated by cholecystokinin in rat pancreas both in vitro and in vivo. J Biol Chem. 1996;271:5686–90.
Gibbs J, Young RC, Smith GP. Cholecystokinin decreases food intake in rats. J Comp Physiol Psychol. 1973;84:488–95.
Witkamp RF. Current and future drug targets in weight management. Pharm Res. 2011;28:1792–818.
Varga G, Bálint A, Burghardt B, D'Amato M. Involvement of endogenous CCK and CCK1 receptors in colonic motor function. Br J Pharmacol. 2004;141:1275–84.
Dufresne M, Seva C, Fourmy D. Cholecystokinin and gastrin receptors. Physiol Rev. 2006;86(3):805–47.
Cawston EE, Miller LJ. Therapeutic potential for novel drugs targeting the type 1 cholecystokinin receptor. Br J Pharmacol. 2010;159(5):1009–21.
Konturek PC, Konturek SJ, Brzozowski T. Helicobacter pylori infection in gastric cancerogenesis. J Physiol Pharmacol. 2009;60(3):3–21.
Matysiak-Budnik T, Mégraud F. Helicobacter pylori infection and gastric cancer. Eur J Cancer. 2006;42(6):708–16.
Croft D, O'Kelly G, Wu G, Haw R, Gillespie M, Matthews L, et al. Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res. 2011;39(Database issue):D691–7.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(D1):D109–14.
Kaizu K, Ghosh S, Matsuoka Y, Moriya H, Shimizu-Yoshida Y, Kitano H. A comprehensive molecular interaction map of the budding yeast cell cycle. Mol Syst Biol. 2010;6:415.
Gloaguen P, Crépieux P, Heitzler D, Poupon A, Reiter E. Mapping the follicle-stimulating hormone-induced signalling networks. Front Endocrinol. 2011;2:45.
Oda K, Matsuoka Y, Funahashi A, Kitano H. A comprehensive pathway map of epidermal growth factor receptor signaling. Mol Syst Biol. 2005;1:2005.0010.
Calzone L, Gelay A, Zinovyev A, Radvanyi F, Barillot E. A comprehensive modular map of molecular interactions in RB/E2F pathway. Mol Syst Biol. 2008;4:173.
Caron E, Ghosh S, Matsuoka Y, Ashton-Beaucage D, Therrien M, Lemieux S, et al. A comprehensive map of the mTOR signaling network. Mol Syst Biol. 2010;6:453.
Fink MY, Pincas H, Choi SG, Nudelman G, Sealfon SC. Research Resource: Gonadotropin-Releasing Hormone Receptor-Mediated Signaling Network in L beta T2 Cells: A Pathway-Based Web-Accessible Knowledgebase. Mol Endocrinol. 2010;24(9):1863–71.
Mizuno S, Iijima R, Ogishima S, Kikuchi M, Matsuoka Y, Ghosh S, et al. AlzPathway: a comprehensive map of signaling pathways of Alzheimer’s disease. BMC Syst Biol. 2012;6:52.
Patil S, Pincas H, Seto J, Nudelman G, Nudelman I, Sealfon SC. Signaling network of dendritic cells in response to pathogens: a community-input supported knowledgebase. BMC Syst Biol. 2010;4:137.
Raza S, McDerment N, Lacaze PA, Robertson K, Watterson S, Chen Y, et al. Construction of a large scale integrated map of macrophage pathogen recognition and effector systems. BMC Syst Biol. 2010;4:63.
Zhong J, Sharma J, Raju R, Palapetta SM, Prasad TS, Huang TC, et al. TSLP signaling pathway map: a platform for analysis of TSLP-mediated signaling. Database. 2014;2014:bau007.
Subbannayya Y, Anuja K, Advani J, Ojha UK, Nanjappa V, George B, et al. A network map of the gastrin signaling pathway. J Cell Commun Signal. 2014;8(2):165–70.
Zinovyev A, Viara E, Calzone L, Barillot E. BiNoM: a Cytoscape plugin for manipulating and analyzing biological networks. Bioinformatics. 2008;24(6):876–7.
Glaab E, Baudot A, Krasnogor N, Valencia A. Extending pathways and processes using molecular interaction networks to analyse cancer genome data. BMC Bioinformatics. 2010;11:597.
Novere NL, Hucka M, Mi H, Moodie S, Schreiber F, Sorokin A, et al. The systems biology graphical notation. Nat Biotech. 2009;27(8):735–41.
Kitano H, Funahashi A, Matsuoka Y, Oda K. Using process diagrams for the graphical representation of biological networks. Nat Biotechnol. 2005;23(8):961–6.
Novere NL, Finney A, Hucka M, Bhalla US, Campagne F, Collado-Vides J, et al. Minimum information requested in the annotation of biochemical models (MIRIAM). Nat Biotech. 2005;23(12):1509–15.
Frisch M, Klocke B, Haltmeier M, Frech K. LitInspector: literature and signal transduction pathway mining in PubMed abstracts. Nucleic Acids Res. 2009;37 suppl 2:W135–40.
Hoffmann R, Valencia A. A gene network for navigating the literature. Nat Genet. 2004;36(7):664.
Matsuoka Y, Ghosh S, Kikuchi N, Kitano H. Payao: a community platform for SBML pathway model curation. Bioinformatics. 2010;26(10):1381–3.
Mi H, Thomas P. PANTHER pathway: an ontology-based pathway database coupled with data analysis tools. Methods Mol Biol. 2009;563:123–40.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
Assenov Y, Ramirez F, Schelhorn SE, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008;24(2):282–4.
Souiai O, Becker E, Prieto C, Benkahla A, De Las Rivas J, Brun C. Functional integrative levels in the human interactome recapitulate organ organization. PLoS One. 2011;6(7):e22051.
Noble F, Wank SA, Crawley JN, Bradwejn J, Seroogy KB, Hamon M, et al. International Union of Pharmacology. XXI. Structure, Distribution, and Functions of Cholecystokinin Receptors. Pharmacol Rev. 1999;51(4):745–81.
Christophe J. Pancreatic tumoral cell line AR42J: an amphicrine model. Am J Physiol. 1994;266(6):G963–71.
Dabrowski A, VanderKuur JA, CarterSu C, Williams JA. Cholecystokinin stimulates formation of Shc-Grb2 complex in rat pancreatic acinar cells through a protein kinase C-dependent mechanism. J Biol Chem. 1996;271(43):27125–9.
Dabrowski A, Groblewski GE, Schafer C, Guan KL, Williams JA. Cholecystokinin and EGF activate a MAPK cascade by different mechanisms in rat pancreatic acinar cells. Am J Physiol Cell Physiol. 1997;273(5):C1472–9.
Sabbatini ME, Bi Y, Ji B, Ernst SA, Williams JA. CCK activates RhoA and Rac1 differentially through Gα13 and Gαq in mouse pancreatic acini. Am J Physiol Cell Physiol. 2010;298(3):C592–601.
Duan RD, Williams JA. Cholecystokinin rapidly activates mitogen-activated protein kinase in rat pancreatic acini. Am J Physiol Gastrointest Liver Physiol. 1994;267(3):G401–8.
Paillasse MR, de Medina P, Amouroux G, Mhamdi L, Poirot M, Silvente-Poirot S. Signaling through cholesterol esterification: a new pathway for the cholecystokinin 2 receptor involved in cell growth and invasion. J Lipid Res. 2009;50(11):2203–11.
Sancho V, Berna MJ, Thill M, Jensen RT. PKCθ activation in pancreatic acinar cells by gastrointestinal hormones/neurotransmitters and growth factors is needed for stimulation of numerous important cellular signaling cascades. Biochim Biophys Acta. 2011;1813(12):2145–56.
Quattrone A, Dewaele B, Wozniak A, Bauters M, Vanspauwen V, Floris G, et al. Promoting role of cholecystokinin 2 receptor (CCK2R) in gastrointestinal stromal tumours pathogenesis. J Pathol. 2012;228(4):565–74.
Wu V, Yang M, McRoberts JA, Ren J, Seensalu R, Zeng N, et al. First intracellular loop of the human cholecystokinin-A receptor is essential for cyclic AMP signaling in transfected HEK-293 cells. J Biol Chem. 1997;272:9037–42.
Le Page SL, Bi Y, Williams JA. CCK-A receptor activates RhoA through G alpha 12/13 in NIH3T3 cells. Am J Physiol Cell Physiol. 2003;285:C1197–206.
Moustafa A, Sakamoto KQ, Habara Y. A fundamental role for NO-PLC signaling pathway in mediating intracellular Ca2+ oscillation in pancreatic acini. Nitric Oxide Biol Chem. 2011;24:139–50.
Cordelier P, Estève JP, Rivard N, Marletta M, Vaysse N, Susini C, et al. The activation of neuronal NO synthase is mediated by G-protein betagamma subunit and the tyrosine phosphatase SHP-2. FASEB J. 1999;13:2037–50.
Thorn P, Gerasimenko O, Petersen OH. Cyclic ADP-ribose regulation of ryanodine receptors involved in agonist evoked cytosolic Ca2+ oscillations in pancreatic acinar cells. EMBO J. 1994;13:2038–43.
Cosker F, Cheviron N, Yamasaki M, Menteyne A, Lund FE, Moutin M-J, et al. The ecto-enzyme CD38 is a nicotinic acid adenine dinucleotide phosphate (NAADP) synthase that couples receptor activation to Ca2+ mobilization from lysosomes in pancreatic acinar cells. J Biol Chem. 2010;285:38251–9.
Calcraft PJ, Ruas M, Pan Z, Cheng X, Arredouani A, Hao X, et al. NAADP mobilizes calcium from acidic organelles through two-pore channels. Nature. 2009;459:596–600.
Gurda GT, Guo L, Lee SH, Molkentin JD, Williams JA. Cholecystokinin activates pancreatic calcineurin-NFAT signaling in vitro and in vivo. Mol Biol Cell. 2008;19(1):198–206.
Lankisch TO, Nozu F, Owyang C, Tsunoda Y. High-affinity cholecystokinin type A receptor/cytosolic phospholipase A2 pathways mediate Ca2+ oscillations via a positive feedback regulation by calmodulin kinase in pancreatic acini. Eur J Cell Biol. 1999;78(9):632–41.
Yoshida H, Nozu F, Lankisch TO, Mitamura K, Owyang C, Tsunoda Y. A possible role for Ca(2+)/calmodulin-dependent protein kinase IV during pancreatic acinar stimulus-secretion coupling. Biochim Biophys Acta. 2000;1497(1):155–67.
Sinclair NF, Ai W, Raychowdhury R, Bi M, Wang TC, Koh TJ, et al. Gastrin regulates the heparin-binding epidermal-like growth factor promoter via a PKC/EGFR-dependent mechanism. Am J Physiol Gastrointest Liver Physiol. 2004;286(6):G992–9.
Miyazaki Y, Shinomura Y, Tsutsui S, Zushi S, Higashimoto Y, Kanayama S, et al. Gastrin induces heparin-binding epidermal growth factor-like growth factor in rat gastric epithelial cells transfected with gastrin receptor. Gastroenterology. 1999;116(1):78–89.
Sturany S, Van Lint J, Gilchrist A, Vandenheede JR, Adler G, Seufferlein T. Mechanism of Activation of Protein Kinase D2(PKD2) by the CCKB/Gastrin Receptor. J Biol Chem. 2002;277(33):29431–6.
Yassin RR, Little KM. Early signalling mechanism in colonic epithelial cell response to gastrin. Biochem J. 1995;311(Pt 3):945–50.
He H, Shulkes A, Baldwin GS. PAK1 interacts with beta-catenin and is required for the regulation of the beta-catenin signalling pathway by gastrins. Biochim Biophys Acta. 2008;1783(10):1943–54.
Mishra P, Senthivinayagam S, Rana A, Rana B. Glycogen Synthase Kinase-3beta regulates Snail and beta-catenin during gastrin-induced migration of gastric cancer cells. J Mol Signal. 2010;5:9.
He H, Baldwin GS. Rho GTPases and p21-activated kinase in the regulation of proliferation and apoptosis by gastrins. Int J Biochem Cell Biol. 2008;40(10):2018–22.
He H, Yim M, Liu KH, Cody SC, Shulkes A, Baldwin GS. Involvement of G proteins of the Rho family in the regulation of Bcl-2-like protein expression and caspase 3 activation by Gastrins. Cell Signal. 2008;20(1):83–93.
Guo Y-S, Cheng J-Z, Jin G-F, Gutkind JS, Hellmich MR, Townsend CM. Gastrin stimulates cyclooxygenase-2 expression in intestinal epithelial cells through multiple signaling pathways. Evidence for involvement of ERK5 kinase and transactivation of the epidermal growth factor receptor. J Biol Chem. 2002;277:48755–63.
von Blume J, Knippschild U, Dequiedt F, Giamas G, Beck A, Auer A, et al. Phosphorylation at Ser244 by CK1 determines nuclear localization and substrate targeting of PKD2. EMBO J. 2007;26(22):4619–33.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Seva C, Kowalski-Chauvel A, Daulhac L, Barthez C, Vaysse N, Pradayrol L. Wortmannin-Sensitive Activation of p70S6-Kinase and MAP-Kinase by the G Protein-Coupled Receptor, G/CCKB. Biochem Biophys Res Commun. 1997;238(1):202–6.
Kikani CK, Dong LQ, Liu F. “New”-clear functions of PDK1: beyond a master kinase in the cytosol? J Cell Biochem. 2005;96(6):1157–62.
Pradeep A, Sharma C, Sathyanarayana P, Albanese C, Fleming JV, Wang TC, et al. Gastrin-mediated activation of cyclin D1 transcription involves beta-catenin and CREB pathways in gastric cancer cells. Oncogene. 2004;23(20):3689–99.
Steigedal TS, Bruland T, Misund K, Thommesen L, Laegreid A. Inducible cAMP early repressor suppresses gastrin-mediated activation of cyclin D1 and c-fos gene expression. Am J Physiol Gastrointest Liver Physiol. 2007;292(4):G1062–9.
Dehez S, Daulhac L, Kowalski-Chauvel A, Fourmy D, Pradayrol L, Seva C. Gastrin-induced DNA synthesis requires p38-MAPK activation via PKC/Ca2+ and Src-dependent mechanisms. FEBS Lett. 2001;496(1):25–30.
Dehez S, Bierkamp C, Kowalski-Chauvel A, Daulhac L, Escrieut C, Susini C, et al. c-Jun NH2-terminal Kinase Pathway in Growth-promoting Effect of the G Protein-coupled Receptor Cholecystokinin B Receptor: A Protein Kinase C/Src-dependent-Mechanism. Cell Growth Differ. 2002;13(8):375–85.
Stepan VM, Dickinson CJ, del Valle J, Matsushima M, Todisco A. Cell type-specific requirement of the MAPK pathway for the growth factor action of gastrin. Am J Physiol. 1999;276(6 Pt 1):G1363–72.
Zimmermann S, Moelling K. Phosphorylation and Regulation of Raf by Akt (Protein Kinase B). Science. 1999;286(5445):1741–4.
Bierkamp C, Kowalski-Chauvel A, Dehez S, Fourmy D, Pradayrol L, Seva C. Gastrin mediated cholecystokinin-2 receptor activation induces loss of cell adhesion and scattering in epithelial MDCK cells. Oncogene. 2002;21(50):7656–70.
Mishra P, Senthivinayagam S, Rangasamy V, Sondarva G, Rana B. Mixed Lineage Kinase-3/JNK1 Axis Promotes Migration of Human Gastric Cancer Cells following Gastrin Stimulation. Mol Endocrinol. 2010;24(3):598–607.
Yu H-G, Nam J-O, Miller NLG, Tanjoni I, Walsh C, Shi L, et al. p190RhoGEF (Rgnef) Promotes Colon Carcinoma Tumor Progression via Interaction with Focal Adhesion Kinase. Cancer Res. 2011;71(2):360–70.
Ding J, Yu JP, Li D, Yu HG, Luo HS, Wei WZ. Effect of gastrin on invasiveness of human colon cancer cells. Zhonghua zhong liu za zhi. 2005;27(4):213–5.
Fjeldbo CS, Bakke I, Erlandsen SE, Holmseth J, Laegreid A, Sandvik AK, et al. Gastrin upregulates the prosurvival factor secretory clusterin in adenocarcinoma cells and in oxyntic mucosa of hypergastrinemic rats. Am J Physiol Gastrointest Liver Physiol. 2012;302(1):G21–33.
Ramamoorthy S, Stepan V, Todisco A. Intracellular mechanisms mediating the anti-apoptotic action of gastrin. Biochem Biophys Res Commun. 2004;323(1):44–8.
Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in Biological Networks. Bioinformatics. 2005;21(16):3448–9.
Lee MJ, Ye AS, Gardino AK, Heijink AM, Sorger PK, MacBeath G, et al. Sequential application of anticancer drugs enhances cell death by rewiring apoptotic signaling networks. Cell. 2012;149(4):780–94.
Jun Cao J-PY, Chao-Hong L, Lan Z, Hong-Gang Y. Effects of gastrin 17 on β-catenin/Tcf-4 pathway in Colo320WT colon cancer cells. World J Gastroenterol. 2006;12(46):7482–7.
Clarke PA, Dickson JH, Harris JC, Grabowska A, Watson SA. Gastrin enhances the angiogenic potential of endothelial cells via modulation of heparin-binding epidermal-like growth factor. Cancer Res. 2006;66(7):3504–12.
Trulsson LM, Gasslander T, Svanvik J. Cholecystokinin-8-induced hypoplasia of the rat pancreas: influence of nitric oxide on cell proliferation and programmed cell death. Basic Clin Pharmacol Toxicol. 2004;95(4):183–90.
Tiger CF, Krause F, Cedersund G, Palmer R, Klipp E, Hohmann S, et al. A framework for mapping, visualisation and automatic model creation of signal-transduction networks. Mol Syst Biol. 2012;8:578.
Leach SM, Tipney H, Feng W, Baumgartner WA, Kasliwal P, Schuyler RP, et al. Biomedical discovery acceleration, with applications to craniofacial development. PLoS Comput Biol. 2009;5(3):e1000215.
Gitter A, Carmi M, Barkai N, Bar-Joseph Z. Linking the signaling cascades and dynamic regulatory networks controlling stress responses. Genome Res. 2012;23(2):365–76.
We thank Charles E. Chapple for sharing his protein-protein interaction network and Enrico Glaab for “a la carte” PathExpand analyses. This work was supported by The Norwegian Cancer Society, The Liaison Committee between the Central Norway Regional Health Authority (RHA), the Norwegian University of Science and Technology (NTNU) and Sør-Trøndelag University College (HiST).
The authors declare that they have no competing interests.
ST conceived the idea of the CCKR map, participated in map construction, curation and modular decomposition, interpretation of PPI- analyses and manuscript writing. ÅF participated in map construction and curation, and manuscript writing. KC participated in modular decomposition with BiNoM analysis, carried out network topology studies and helped in manuscript writing. AB provided protein-protein interaction network analysis data and manuscript writing. TB participated in pathways curation and manuscript writing. LT raised funding, participated in pathway curation and helped to draft the manuscript, MK raised funding, participated in interpretation of PPI- analyses and helped to draft the manuscript. AL raised funding, participated in pathway curation, interpretation of PPI- and helped to draft the manuscript. All authors read and approved the final manuscript.
The file is the SBML version of the complete CCKR map.
This file is an original Cytoscape session file containing the BiNoM generated modules of the CCKR map.
The file contains a description of individual modules of the comprehensive CCKR map constructed using the BiNoM plugin in Cytoscape.
The file contains a list of CCKR model protein interactors identified from the large scale protein-protein interaction network.
The file contains the PathExpand analysis of BiNoM modules and their protein interactors.
About this article
Cite this article
Tripathi, S., Flobak, Å., Chawla, K. et al. The gastrin and cholecystokinin receptors mediated signaling network: a scaffold for data analysis and new hypotheses on regulatory mechanisms. BMC Syst Biol 9, 40 (2015). https://doi.org/10.1186/s12918-015-0181-z