Skip to main content

Advertisement

The structural network of Interleukin-10 and its implications in inflammation and cancer

Article metrics

Abstract

Background

Inflammation has significant roles in all phases of tumor development, including initiation, progression and metastasis. Interleukin-10 (IL-10) is a well-known immuno-modulatory cytokine with an anti-inflammatory activity. Lack of IL-10 allows induction of pro-inflammatory cytokines and hinders anti-tumor immunity, thereby favoring tumor growth. The IL-10 network is among the most important paths linking cancer and inflammation. The simple node-and-edge network representation is useful, but limited, hampering the understanding of the mechanistic details of signaling pathways. Structural networks complete the missing parts, and provide details. The IL-10 structural network may shed light on the mechanisms through which disease-related mutations work and the pathogenesis of malignancies.

Results

Using PRISM (a PRotein Interactions by Structural Matching tool), we constructed the structural network of IL-10, which includes its first and second degree protein neighbor interactions. We predicted the structures of complexes involved in these interactions, thereby enriching the available structural data. In order to reveal the significance of the interactions, we exploited mutations identified in cancer patients, mapping them onto key proteins of this network. We analyzed the effect of these mutations on the interactions, and demonstrated a relation between these and inflammation and cancer. Our results suggest that mutations that disrupt the interactions of IL-10 with its receptors (IL-10RA and IL-10RB) and α2-macroglobulin (A2M) may enhance inflammation and modulate anti-tumor immunity. Likewise, mutations that weaken the A2M-APP (amyloid precursor protein) association may increase the proliferative effect of APP through preventing β-amyloid degradation by the A2M receptor, and mutations that abolish the A2M-Kallikrein-13 (KLK13) interaction may lead to cell proliferation and metastasis through the destructive effect of KLK13 on the extracellular matrix.

Conclusions

Prediction of protein-protein interactions through structural matching can enrich the available cellular pathways. In addition, the structural data of protein complexes suggest how oncogenic mutations influence the interactions and explain their potential impact on IL-10 signaling in cancer and inflammation.

Background

Inflammation by innate immunity is the first line of defense against pathogenic infections [1]. It is also involved in all phases of cancer development, including tumor initiation, promotion and metastatic dissemination [24]. By triggering immunosuppressive mechanisms, inflammation creates a tissue microenvironment which permits tumor growth and metastasis [2]. Inflammatory cells provide growth factors that sustain proliferation, and survival factors that allow escape from apoptosis; it also contributes to extracellular matrix (ECM) modifying enzymes, and to pro-angiogenic factors that facilitate angiogenesis, invasion and ultimately metastasis [3].

Several lines of evidence link cancer and inflammation, emphasizing that chronic inflammation contributes to tumor initiation and progression [5, 6]. Chronic inflammatory bowel disease predisposes individuals to colon cancer [6] and individuals with chronic hepatitis are more prone to develop hepatocellular carcinoma [7]. Chronic Helicobacter pylori infection and the resulting inflammation is associated with gastric cancer [8]; chronic bronchitis with lung cancer; and pancreatitis with pancreas cancer [9]. Additionally, long term use of non-steroidal anti-inflammatory drugs (NSAIDs) [10] which inhibit pro-inflammatory cytokines, like TNF-α and IL-1β, decrease cancer incidence [11].

Identified in 1989 [12], IL-10 is an anti-inflammatory cytokine that modulates the immune response: if IL-10 is not present or functional, inflammation becomes possible. It restricts the immune response to pathogens and prevents damage to the host. Secreted by immune cells, IL-10 diversely affects cell types in the immune system. Although it terminates inflammatory responses by suppressing monocyte/macrophage function, it also acts as an immunostimulator to promote Th2 response. IL-10 regulates growth and/or differentiation of B cells, NK cells, cytotoxic and helper T cells, mast cells, granulocytes, dendritic cells, keratinocytes, and endothelial cells. Additionally, it stimulates immunoglobulin secretion, and promotes antibody class switching [13]. Therefore, IL-10 has both immune suppressive (anti-inflammatory) and immune stimulatory roles (B and T-cell development).

IL-10 deficiency increases the production of IL-1 (a pro-inflammatory cytokine) and in the absence of IL-10, IL-1 promotes tumor growth in mice [14]. IL-10 also prevents development of regulatory T cells (Tregs) and Myeloid-derived suppressor cells (MDSCs) [1417]. IL-10 deficiency leads to an increase in the number of Tregs and MDSCs in tumor tissue. Tregs and MDSCs have suppressive roles against tumor-specific immunity that favor tumor growth [4, 14, 15]. Apart from its anti-inflammatory roles, it is associated with activation of anti-tumor immunity [18]. The presence of IL-10 leads to tumor regression and increase in tumor-specific immunogenicity [19]. In contrast, some studies proposed that blockage of IL-10 signaling promote anti-tumor immunity [20]. These controversial findings stem from the pleiotropic effects of IL-10 and different experimental models (human or animal, in vitro or in vivo, solid or hematological tumors, presence of exogenous or endogenous IL-10 or IL-10 inhibitors, etc.) and the varying site of IL-10 production [19].

IL-10 is a dimeric cytokine [21] that signals through a tetrameric transmembrane receptor complex, consisting of two IL-10RA (also known as IL-10R1) and two IL-10RB (also known as IL-10R2) proteins [22, 23]. Both receptors belong to the class II receptor family, and consist of three domains: an intracellular domain, a transmembrane domain, and an extracellular domain [13]. The receptor complex assembles sequentially: IL-10RA, with higher affinity, binds to IL-10 first and then IL-10RB [21]. IL-10 binding to the extracellular domain of IL-10RA leads to phosphorylation of JAK1 (Janus Kinase-1) and TYK2 (Tyrosine Kinase-2), that interacts with IL-10RA and IL-10RB [24], respectively. Specific tyrosine residues on the intracellular domain of the IL-10RA are then phosphorylated by these kinases. The STAT3 (Signal Transducer and Activator of Transcription-3) transcription factor binds to those tyrosine residues and gets phosphorylated. Activated STAT3 then translocates to the nucleus as a homodimer, and activates transcription of anti-apoptotic and cell-cycle-progression genes [25, 26].

IL-10 has many functional partners, one of which is the α2-macroglobulin (A2M). IL-10 forms a stable complex with activated A2M [27]. A2M is a large homotetrameric glycoprotein [28] in the plasma, and in the extracellular space. It is a proteinase inhibitor, and has a peptide stretch, called the 'bait region'. Cleavage of the bait region by a proteinase leads to a conformational change in the protein that causes proteinase trapping, and receptor-mediated endocytosis of the A2M-proteinase complex [29]. A2M is also a cytokine transporter. Many cytokines, including IL-10, and growth factors bind to A2M non-covalently in vivo [30]. When A2M forms complexes with IL-10, TGF-β (anti-inflammatory cytokines) and IFN-γ, it accelerates the appearance of these cytokines in the blood [31, 32]. A2M in its native form increases the half-life of bound cytokines in the plasma by protecting them from proteolysis [30]. Thus, at sites of inflammation, A2M concentration rises as a response to an increase in proteinase level. [27]. A2M also contributes to the anti-inflammatory response of IL-10 by preventing its destruction.

IL-10 and IL-10Rs are able to interact with many partner molecules in the signaling network; however, their detailed protein structural interactions, as well the corresponding mutational mechanisms have not been well illustrated. In this study, we constructed the structural pathway based on protein-protein interactions (PPIs). The commonly used node-and-edge description of pathways, where nodes represent proteins and edges the interactions between them, are useful, but do not provide structural interaction detail [33, 34]. Further, in many cases, such as IL-10 and the receptors in this study, the available structural interaction data of the proteins are incomplete. However, recently developed computational structural approaches, such as PRISM (PRotein Interactions by Structural Matching tool), are capable of predicting PPI and can help filling in the gaps. PRISM [35, 36] is a motif-based protein-protein interaction prediction tool which uses a knowledge-based strategy to construct and analyze structural PPI networks. PRISM is based on the notion that evolution has exploited favorable structural motifs adapting them to different functions, in protein folds and at protein-protein-interfaces [3739], lending robustness to its predictions. PRISM has predicted protein interactions successfully [40, 41] for different pathways, like apoptosis [42], ubiquitination [43], MAPK [41, 44], the Toll-like receptor pathways, [1] and for identifying drug off-targets [45]. The success of PRISM is very close to %100 (87 out of 88 cases) in rigid-body prediction [40]. Recently we have further enhanced it by introducing ensemble docking, by exploiting different conformations, and PRISM could predict two thirds of the 'difficult' cases of a docking benchmark dataset [41]. Here, we applied this enhanced PRISM protocol to construct the structural protein-protein interaction (PPI) network of IL-10 centered signaling. Importantly, the analysis was able to identify mutations falling in the interfaces and to predict their effects on interactions such as of IL-10 with its receptors, IL-10 with A2M, and A2M with APP and KLK13. This allowed us to enrich the structural interaction data of IL-10 with its partners, and to analyze the mechanisms of mutations that lead to inflammation, and cancer through their impact on predicted interactions.

Methods

Reconstruction of the structural PPI network of IL-10 centered signaling

We used the String server [46] for selecting the first and second-degree neighbors of IL-10. Only interactions with experimental evidence and confidence score larger than 0.4 (the default confidence value) were considered. There were 4 first-degree and 45 second-degree neighbor proteins of IL-10 (Additional file 1). Overall, there were 50 proteins comprising the IL-10 centered protein-protein interaction network.

We used PRISM [35, 36] for modeling protein-protein interactions in the IL-10 centered network. PRISM searches for the motifs on the target protein surfaces similar to known interactions considering both geometrical complementarities and evolutionary conservation of hot spots. It treats proteins with at least 15 residues. To model an interface, PRISM requires the 3D structures of the proteins of interest (for further details of the PRISM protocol, see [36]). 39 of these 50 proteins have structural data in the PDB (corresponding to 958 PDB chains) and we could build homology models for the remaining 10. IGHV3-6's sequence information could not be found, so this protein is not included in our analysis (Additional file 1 Table S1). The I-TASSER server [47] was used for homology modeling and the top 5 models generated by the server were included in the predictions.

We reduced the redundancy of similar interface architectures for each protein, using TM (template modeling)-align [48]. PDB structures having TM-scores larger than 0.5 and RMSD under 2.5Å were classified. Then, a representative PDB structure was assigned for each similar structure group and we ended up with 127 representative structures for 39 target proteins. The final IL-10 centered network is composed of 49 proteins and 70 interactions (Additional file 1 and 2).

Mapping oncogenic mutations onto the interfaces of predicted protein-protein complexes and in silico mutagenesis

The reconstructed structural PPI network of IL-10 centered signaling not only reveals many important details about the mechanism of protein-protein interactions but also offers the possibility of observing the effects of oncogenic mutations. In the case studies presented below, mutational data related to the proteins in the network were taken from the COSMIC (Catalogue of Somatic Mutations in Cancer) database [49, 50] and cBioPortal for Cancer Genomics (The Cancer Genome Atlas, TCGA) [51]. To map oncogenic mutations, we identified the interfaces of the modeled protein-protein complexes using the HotPoint web server, which uses conservation, solvent accessibility and pairwise residue potential data to determine computational hot spots [52]. After the mutational mapping, we performed in silico mutagenesis to observe the effects of the mutations on the interactions (Figure 1). We computationally mutated those key residues using the FoldX plugin [53] for the YASARA molecular viewer [54]. We minimized the energies of the proteins both before and after the mutation and then used the mutant structures to re-run PRISM [35, 36] and model the new interaction between the mutant target and its partner (Figure 1).

Figure 1
figure1

Mapping oncogenic mutations onto the interfaces of predicted protein-protein complexes and in silico mutagenesis. Mutational data related to the proteins in the network were taken from the COSMIC database [49, 50] and The Cancer Genome Atlas, TCGA [51]. To map oncogenic mutations, we identified the interfaces of the modeled protein-protein complexes using the HotPoint web server [52]. After the mutational mapping, we performed in silico mutagenesis to observe the effects of mutations - that are on the interface - on the interactions. We computationally mutated those key residues using the FoldX plugin [53] for the YASARA molecular viewer [54]. We minimized the energies of the proteins both before and after the mutation and then used the mutant structures to re-run PRISM [35, 36] and model the new interaction between the mutant target and its partner.

Results and discussion

We constructed the IL-10 centered human structural protein-protein interaction network with first and second-degree protein neighbors. This network is composed of 49 proteins and 70 interactions between them (Figure 2). Among these 70 interactions, the structure of only two (IL-10 - IL-10RA and APOE - LRP1) are deposited in the PDB in a complex form. By using PRISM, we predicted the structures of the PPI interfaces and 40 additional interactions were structurally modeled (Figure 2). As a result, the available structural data increased from 2 to 42.

Figure 2
figure2

The protein-protein interaction network of Interleukin-10. There are 49 proteins and 70 interactions in this network and only 2 of the interactions have structural data in a complex form in the PDB (edges highlighted in yellow). We modeled the interfaces for 40 additional interactions. Thus there are 42 interactions with interface models (edges highlighted in pink). The remaining 28 edges (out of 70) could not be modeled and are shown in cyan.

In mutation databases, like COSMIC [50, 55], there are thousands of experimentally verified cancer mutations but the precise mechanisms of mutations, how they change protein functions and contribute to cancer pathogenesis is unknown [56]. Proteins function through interactions and mutations that disrupt interactions also change protein functions. Some mutations abolish protein interactions, whereas others make the interactions stronger or change protein folding, favoring aggregation of proteins, as in the case of formation of amyloids in Alzheimer's disease [57]. 4% of all mutations in the databases were computationally predicted to be related with protein interactions [56]. Structural knowledge of protein-protein interfaces is important, allowing mapping mutations and SNPs onto the structures of the complexes, identifying those falling on the interfaces, and predicting their effects on protein interactions [33, 5860]. PRISM is a tool that provides the opportunity to observe the effects of these mutations on the interactions. Overall, 879 "missense" and "coding silent" mutations from COSMIC database [49, 50] are mapped onto 29 target proteins (40 target PDB structures) in the IL-10 centered network (Table 1 and Additional file 3). Each of the interactions of these 29 target proteins can be studied for observing the effects of mutations that are on the interfaces of the predicted models. Below, we present four examples of interactions in the reconstructed structural IL-10 centered network and these interactions are related to oncogenesis (Figure 3). In these case studies, we computationally mutated the wild type residues to their corresponding oncogenic mutation variants and analyzed the effects on the original interaction.

Table 1 The distribution of COSMIC "missense" and "coding silent" mutations mapped onto the target structures in IL-10 centered network.
Figure 3
figure3

The case study interactions highlighted on IL-10 centered network. The thick orange lines indicate four case studies explaining the relation between mutations on the interfaces of protein-protein complex models and cancer (IL10-IL10RB, IL10-A2M, A2M-KLK13 and A2M-APP).

Case studies

The relation between the interaction of IL-10 with its receptors and the implications in inflammation & cancer

The structure of the IL-10 and IL-10RA complex is available in the PDB [21], whereas the IL-10 and IL-10RB complex structure has not been solved yet. We were able to model the interaction of IL-10 (PDB code: 2ilkA) with IL-10RB (3lqmA) using PRISM (Figure 4). The template interface used in the prediction of this interaction is 2a6aAB (a glycoprotein endopeptidase homodimer) and the interaction energy score is -15.98 energy units. In order to confirm our model experimentally, we compared the interface of the predicted structure with the critical residues in the binding of IL-10 - IL-10R2 (IL-10RB) [61], which were determined by surface plasmon resonance and cell-based assays. There is a good agreement between the predicted interface residues and experimentally determined critical residues in binding assays. Five out of six residues, Asn21, Arg24, Arg32, His90 and Ser93, which are determined to be critical in binding, are also predicted to be in the interface. In order to execute its function, IL-10 needs to interact with both receptors IL-10RA and IL-10RB at the same time, as the receptor complex is a ternary structure [61]. In the PRISM-predicted model, the receptors can bind to IL-10 without clashing and can therefore form a ternary complex in agreement with functional data (Figure 4). PRISM also predicts the interaction of IL-10 with IL-10RA with -25.54 energy units, utilizing its own complex structure as the template interface (1j7vLR). IL-10RA has higher affinity for IL-10 (-25.54 energy units) compared to the affinity of IL-10RB for IL-10 (-15.98 energy units). This is also in line with the sequential mechanism of the ternary complex formation with IL-10 binding first to IL-10RA, and then IL-10RB is recruited [61].

Figure 4
figure4

The predicted complex structure of IL-10 with its receptors, IL-10RA and IL-10RB. The complex between IL-10 and IL-10RA is available in PDB (PDB ID: 1j7v, L and R chains, respectively) and the interaction of IL-10 (PDB ID: 2ilkA) with IL-10RB (PDB ID: 3lqmA) is modeled by PRISM. Receptors bind to different surfaces of IL-10 that are close to each other.

Structural details in interactions of IL-10 with its receptors may shed light on the mechanisms of mutations. The role of IL-10 in cancer was shown by Tanikawa et al. that IL-10 deficiency causes a rise in the production of IL-1, a pro-inflammatory cytokine, which in turn leads to increased tumor growth in mice [14]. Mutations in IL-10 or its receptors may disrupt their interactions, thereby preventing IL-10 signaling. For instance, E41* nonsense mutation in IL-10RB, abolishes the interaction of IL-10 with IL-10RB, blocking IL-10 signaling. This mutation causes loss of most of the IL-10RB, including a large portion of the interface between IL-10 and IL-10RB (Figure 5, the part drawn in yellow, Table 2). According to the TCGA data, this mutation is observed in lung adenocarcinoma with 2% frequency [51]. Furthermore, the R198W substitution mutation in IL-10RB (from the COSMIC database [49, 50]), also destroys the interaction between IL-10 and IL-10RB (Table 2), and this mutation is observed in endometrioid carcinoma. Importantly, Figure 6 shows that the R198W mutation in IL-10RB falls right next to the interface. Blockage of IL-10 signaling may lead to enhanced inflammation and increased number of Tregs and MDSCs, which inhibit tumor immunity, allowing tumors to grow [15].

Figure 5
figure5

E41* nonsense mutation on IL-10RB protein. IL-10 (2ilkA) is shown in green, blue protein is IL-10RB (3lqmA) and the yellow segment is the deleted part in E41* nonsense mutant of IL-10RB. This mutation causes the loss of a large part of the protein, including most of the interface between IL-10 and IL-10RB, so that it blocks the interaction of IL-10RB with IL-10. If IL-10 cannot interact with its receptors, downstream signaling does not take place and anti-inflammatory outcomes of IL-10 pathway will be lost, allowing inflammation to develop.

Table 2 Comparison of the binding energy scores for wild type and mutant structures
Figure 6
figure6

R198W substitution mutation on IL-10RB protein. IL-10 (2ilkA) is shown in green, blue protein is IL-10RB (3lqmA). R198, highlighted in yellow, is the substituted residue on the mutant IL-10RB. This mutation is very close to the interface between IL-10 and IL-10RB and it abolishes the interaction of IL-10RB with IL-10.

The association of IL-10 with A2M and its implications for inflammation and cancer

Here we concentrate on the interaction of IL-10 with α2-macroglobulin (A2M). A2M mediates the inflammatory response through acting as a cytokine transporter. A2M binding to IL-10 facilitates the recruitment of this cytokine to the site of inflammation and triggers an anti-inflammatory response [31, 32]. A2M also protects IL-10 from proteolysis [30]. Thus, if the interaction between these two is disrupted, inflammation will occur [31, 32], which may favor cancer development. Normally, IL-10 is predicted to bind to A2M with a binding energy score of -39.2 (Figure 7). We mutated IL-10 computationally based on the oncogenic mutation to see whether the interaction is affected. Q56* nonsense mutation in IL-10, that has been seen in lung adenocarcinoma [50, 51], is observed to abolish the association of IL-10 with A2M. As can be seen from Figure 7, this nonsense mutation causes the complete loss of the interface between IL-10 and A2M (the yellow labeled part).

Figure 7
figure7

IL-10 (2ilk, green) is predicted to bind to A2M (1bv8, purple) with a binding energy score of -39.2. The template interface is 1mwqAB, a hypothetical protein from H. influenzae with a putative active-site phosphohistidine. The yellow segment is the deleted part of IL-10 due to Q56* nonsense mutation. The deleted part includes the interface region, so the interaction between IL-10 and A2M is disrupted.

The interaction between A2M with APP and its implications in cancer

Here, we focus on the interaction between A2M, one of the interaction partners of IL-10, and β-amyloid precursor protein (APP). Extracellular cleavage of APP produces C99, a cell membrane-bound fragment, which is further cleaved by γ-secretase, and releases the intracellular domain of APP to produce amyloid-β (Ab) [62]. Ab forms a complex with native A2M and the complex is internalized by the A2M receptor, a low density lipoprotein receptor-related protein (LRP), and degraded [63].. According to the PRISM results, A2M (PDB ID: 4acqC, residues 24-1474) forms a stable complex with APP (PDB ID: 2llmA, residues 686-726) with a binding energy score of -39.12. The template interface used in the prediction is 1fftAC (ubiquinol oxidase from Escherichia coli). APP is an integral type I transmembrane protein with a single transmembrane domain (residues 700-723), a large extracellular ectodomain (residues 18-699), and a short cytoplasmic tail (residues 724-770) [62].

APP has been found to be up-regulated in many cancers, including pancreatic [65], colon [66], and prostate cancer [67] and squamous cell carcinoma [65]. Gain-of-function studies indicated that overexpression of APP causes increased cell proliferation [67, 68]. We performed in silico mutational analysis in order to observe the effect of the mutations on the A2M - APP interaction, and found a possible association between abrogation of this interaction and cancer. Mutations mapping to the interface of the predicted complex were identified, and a missense mutation, R945Q, which is observed in colorectal cancer was selected for further analysis [50, 51]. R945 in human A2M is a hot spot interface residue (Figure 8) and mutating this residue to glutamine resulted in a noticeable change in the binding energy score (Table 2). The score for the mutated structure decreased to -9.12 indicating that substituting arginine with glutamine greatly weakens the interaction. This observation suggests that the weakened A2M - APP interaction could relate to oncogenesis. This idea is supported by recent studies reporting up-regulation of APP in several cancers due to its growth-promoting function [6568]. Mutations that destabilize this interaction may prevent or reduce the degradation of Ab by LRP and cause a rise in APP level [63], which may eventually lead to increased cancer cellular proliferation.

Figure 8
figure8

The predicted complex structure of A2M (4acqC, purple) and APP (2llmA, blue) with an energy score of -39.12. Template interface: 1fftAC (an ubiquinol oxidase). R945, highlighted in yellow, is the mutated residue. The predicted complex structure of mutated A2M (4acqC, purple) and wt APP (2llmA, blue) has an energy score of -9.12.

The structural interaction of A2M with KLK13, and its implication in cancer

Kallikrein-13 (KLK13, hK13) is a member of the kallikrein family and encodes a secreted trypsin-like serine protease, which is regulated by steroid hormones [69, 70]. Experiments show its ability to cleave the major components of the extracellular matrix [71]. It is known that KLK13 can form a complex with the serine protease inhibitors, including A2M [72]. PRISM was able to predict an interaction between KLK13 and A2M with an energy score of -65.22 (Figure 9 and Table 2). The prediction is between the PDB structure 1bv8A (A2M) and a homology model of KLK13 obtained by the I-TASSER server [47]. The template interface used in the prediction is 1g8tAB, a homodimer interface of nuclease.

Figure 9
figure9

The predicted complex structure of mutated A2M (1bv8A, purple) and wt KLK13 (PSA, model, blue) with an energy score of -65.22. Template interface: 1g8tAB, a homodimer interface of nuclease. R263, highlighted in yellow, is the mutated residue.

Although its function is still unclear, KLK13 is used as a new cancer biomarker in various cancers, including prostate [73], breast [74], ovarian [75], salivary gland [71], testicular [76], and non-small cell lung [77] cancers. This protein may be involved in the promotion of cancer cell growth, angiogenesis, metastasis and invasion. The ramification of the R263L mutation on KLK13 is observed in carcinoma from the COSMIC database [49, 50]. We investigated whether this mutation disrupts the interaction of KLK13 - A2M, created a structurally predicted mutant KLK13 protein and minimized its energy. An interaction between the minimized form of wild type KLK13 and A2M is predicted with an energy score of -87.73 (implying that the effect of minimization is insignificant). However, PRISM could not find a favorable interaction between mutant KLK13 and A2M (the best prediction has an energy score of +19.48, Table 2), suggesting disruption of the interaction. Our data is consistent with the prediction from the HotPoint server. Disruption of this interaction can allow KLK13 to react with other proteins, which may lead to cleavage in the major components of the extracellular matrix [71] and help in the promotion of cancer cell growth, metastasis and invasion.

Conclusions

Structural PPI networks indicate not only which proteins interact, but also how they interact and the location of the interaction sites. Computational techniques allow us to predict PPIs, mutate proteins and investigate the effect of those mutations on the PPIs. Here we constructed the structural PPI network to explore mutational and pathogenesis mechanisms in inflammatory diseases and cancer. We focused on the IL-10 centered network, as IL-10 is a well-known cytokine with an anti-inflammatory activity and relation to cancer. Currently available structural data of the IL-10 pathway are incomplete, with only 2 interactions available in the PDB. We utilized homology modeling to obtain the missing protein structures and a motif-based PPI prediction tool to complete the missing network parts. First we modeled the structures of 10 proteins and then provided models for 40 additional interactions. Although PRISM has a high prediction accuracy, its success to predict interactions is dependent on the conformation of the proteins given. We exploited the structures in the PDB. If the PDB does not include a conformation close to the bound form of the protein, PRISM cannot predict the interaction. That is why we missed some interactions on the network. As the PDB gets richer of structures and different conformations, the success of PRISM to predict interactions will increase. However, the structural PPI network was extended form 2 interactions to 42 interactions via predictions. This allowed us to investigate the effect of clinically observed cancer mutations on our IL-10 centered network. Comparing the interaction models of the wild type and mutant proteins, we observed that specific mutations disrupt the interactions, such as between IL-10 and its receptor, IL-10 and A2M, and A2M and its partners, which may disrupt immune regulation in cancer. We discovered that mutations of the residues, which were clinically observed in cancers as hot spots, change the binding energy and abolished or weakened the interactions. Disruption of the interactions of IL-10 with its receptors (IL-10RA and IL-10RB) and α-2-macroglobulin (A2M) may lead to enhanced inflammation, which could promote tumor growth; blockage of the A2M-APP interaction may lead to cancerous cellular proliferation through free APP; blockage of A2M-KLK13 (hK13) interaction can increase free hK13, which can promote cancer cell growth, metastasis and invasion through damage in the extracellular matrix. Collectively, by merging mutational and structural data - available and predicted using our powerful PRISM tool - and combining it with functional data, we are able to reveal the consequences of weakening or abolishing key interactions, and obtain experimentally-testable mechanisms of oncogenic mutations in the IL-10 network.

References

  1. 1.

    Guven Maiorov E, Keskin O, Gursoy A, Nussinov R: The structural network of inflammation and cancer: Merits and challenges. Seminars in cancer biology. 2013, 23 (4): 243-251.

  2. 2.

    Trinchieri G: Cancer and inflammation: an old intuition with rapidly evolving new concepts. Annu Rev Immunol. 2012, 30: 677-706.

  3. 3.

    Trinchieri G: Inflammation in cancer: a therapeutic target?. Oncology (Williston Park). 2011, 25 (5): 418-420.

  4. 4.

    Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell. 2011, 144 (5): 646-674.

  5. 5.

    Vacchelli E, Galluzzi L, Eggermont A, Galon J, Tartour E, Zitvogel L, Kroemer G: Trial Watch: Immunostimulatory cytokines. Oncoimmunology. 2012, 1 (4): 493-506.

  6. 6.

    Balkwill F, Coussens LM: Cancer: an inflammatory link. Nature. 2004, 431 (7007): 405-406.

  7. 7.

    Coussens LM, Werb Z: Inflammation and cancer. Nature. 2002, 420 (6917): 860-867.

  8. 8.

    Bar D, Apte RN, Voronov E, Dinarello CA, Cohen S: A continuous delivery system of IL-1 receptor antagonist reduces angiogenesis and inhibits tumor development. FASEB journal : official publication of the Federation of American Societies for Experimental Biology. 2004, 18 (1): 161-163.

  9. 9.

    Balkwill F, Mantovani A: Inflammation and cancer: back to Virchow?. Lancet. 2001, 357 (9255): 539-545.

  10. 10.

    Robinson SC, Coussens LM: Soluble mediators of inflammation during tumor development. Advances in cancer research. 2005, 93: 159-187.

  11. 11.

    Mantovani A, Allavena P, Sica A, Balkwill F: Cancer-related inflammation. Nature. 2008, 454 (7203): 436-444.

  12. 12.

    Fiorentino DF, Bond MW, Mosmann TR: Two types of mouse T helper cell. IV. Th2 clones secrete a factor that inhibits cytokine production by Th1 clones. The Journal of experimental medicine. 1989, 170 (6): 2081-2095.

  13. 13.

    Moore KW, de Waal Malefyt R, Coffman RL, O'Garra A: Interleukin-10 and the interleukin-10 receptor. Annual review of immunology. 2001, 19: 683-765.

  14. 14.

    Tanikawa T, Wilke CM, Kryczek I, Chen GY, Kao J, Nunez G, Zou W: Interleukin-10 ablation promotes tumor development, growth, and metastasis. Cancer research. 2012, 72 (2): 420-429.

  15. 15.

    Merlo P, Cecconi F: XIAP: inhibitor of two worlds. The EMBO journal. 2013, 32 (16): 2187-2188.

  16. 16.

    Germano G, Allavena P, Mantovani A: Cytokines as a key component of cancer-related inflammation. Cytokine. 2008, 43 (3): 374-379.

  17. 17.

    Allavena P, Germano G, Marchesi F, Mantovani A: Chemokines in cancer related inflammation. Experimental cell research. 2011, 317 (5): 664-673.

  18. 18.

    Wilke CM, Wei S, Wang L, Kryczek I, Kao J, Zou W: Dual biological effects of the cytokines interleukin-10 and interferon-gamma. Cancer immunology, immunotherapy : CII. 2011, 60 (11): 1529-1541.

  19. 19.

    Mocellin S, Marincola FM, Young HA: Interleukin-10 and the immune response against cancer: a counterpoint. Journal of leukocyte biology. 2005, 78 (5): 1043-1051.

  20. 20.

    Kim BG, Joo HG, Chung IS, Chung HY, Woo HJ, Yun YS: Inhibition of interleukin-10 (IL-10) production from MOPC 315 tumor cells by IL-10 antisense oligodeoxynucleotides enhances cell-mediated immune responses. Cancer immunology, immunotherapy : CII. 2000, 49 (8): 433-440.

  21. 21.

    Sun Z, Luo Q, Ye D, Chen W, Chen F: Role of toll-like receptor 4 on the immune escape of human oral squamous cell carcinoma and resistance of cisplatin-induced apoptosis. Molecular cancer. 2012, 11: 33-

  22. 22.

    Sabat R, Grutz G, Warszawska K, Kirsch S, Witte E, Wolk K, Geginat J: Biology of interleukin-10. Cytokine & growth factor reviews. 2010, 21 (5): 331-344.

  23. 23.

    Asadullah K, Sterry W, Volk HD: Interleukin-10 therapy--review of a new approach. Pharmacological reviews. 2003, 55 (2): 241-269.

  24. 24.

    Kotenko SV, Krause CD, Izotova LS, Pollack BP, Wu W, Pestka S: Identification and functional characterization of a second chain of the interleukin-10 receptor complex. The EMBO journal. 1997, 16 (19): 5894-5903.

  25. 25.

    Riley JK, Takeda K, Akira S, Schreiber RD: Interleukin-10 receptor signaling through the JAK-STAT pathway. Requirement for two distinct receptor-derived signals for anti-inflammatory action. The Journal of biological chemistry. 1999, 274 (23): 16513-16521.

  26. 26.

    Donnelly RP, Dickensheets H, Finbloom DS: The interleukin-10 signal transduction pathway and regulation of gene expression in mononuclear phagocytes. Journal of interferon & cytokine research : the official journal of the International Society for Interferon and Cytokine Research. 1999, 19 (6): 563-573.

  27. 27.

    Garber TR, Gonias SL, Webb DJ: Interleukin-4 and IL-10 bind covalently to activated human alpha2-macroglobulin by a mechanism that requires Cys949. Journal of interferon & cytokine research : the official journal of the International Society for Interferon and Cytokine Research. 2000, 20 (2): 125-131.

  28. 28.

    Sheikh AM, Chauhan V, Tsiouris JA, Mehta PD, Burguess K, Fenko MD, Spivack W, Vaughan M, Malik M: Elevated levels of serum alpha(2) macroglobulin in wild black bears during hibernation. Biochimie. 2003, 85 (10): 1027-1032.

  29. 29.

    Rehman AA, Ahsan H, Khan FH: alpha-2-Macroglobulin: a physiological guardian. Journal of cellular physiology. 2013, 228 (8): 1665-1675.

  30. 30.

    Zhevago NA, Samoilova KA: Pro- and anti-inflammatory cytokine content in human peripheral blood after its transcutaneous (in vivo) and direct (in vitro) irradiation with polychromatic visible and infrared light. Photomedicine and laser surgery. 2006, 24 (2): 129-139.

  31. 31.

    Borth W, Scheer B, Urbansky A, Luger TA, Sottrup-Jensen L: Binding of IL-1 beta to alpha-macroglobulins and release by thioredoxin. Journal of immunology. 1990, 145 (11): 3747-3754.

  32. 32.

    Borth W: Alpha 2-macroglobulin. A multifunctional binding and targeting protein with possible roles in immunity and autoimmunity. Annals of the New York Academy of Sciences. 1994, 737: 267-272.

  33. 33.

    Mosca R, Ceol A, Aloy P: Interactome3D: adding structural details to protein networks. Nature methods. 2013, 10 (1): 47-53.

  34. 34.

    Tuncbag N, Gursoy A, Keskin O: Prediction of protein-protein interactions: unifying evolution and structure at protein interfaces. Physical biology. 2011, 8 (3): 035006-

  35. 35.

    Ogmen U, Keskin O, Aytuna AS, Nussinov R, Gursoy A: PRISM: protein interactions by structural matching. Nucleic Acids Research. 2005, 33: W331-W336.

  36. 36.

    Tuncbag N, Gursoy A, Nussinov R, Keskin O: Predicting protein-protein interactions on a proteome scale by matching evolutionary and structural similarities at interfaces using PRISM. Nature Protocols. 2011, 6 (9): 1341-1354.

  37. 37.

    Tsai CJ, Lin SL, Wolfson HJ, Nussinov R: Protein-protein interfaces: Architectures and interactions in protein-protein interfaces and in protein cores. Their similarities and differences. Critical Reviews in Biochemistry and Molecular Biology. 1996, 31 (2): 127-152.

  38. 38.

    Keskin O, Nussinov R, Favorable scaffolds: proteins with different sequence, structure and function may associate in similar ways. Protein Engineering Design & Selection. 2005, 18 (1): 11-24.

  39. 39.

    Keskin O, Gursoy A, Ma B, Nussinov R: Principles of protein-protein interactions: What are the preferred ways for proteins to interact?. Chemical Reviews. 2008, 108 (4): 1225-1244.

  40. 40.

    Tuncbag N, Keskin O, Nussinov R, Gursoy A: Fast and accurate modeling of protein-protein interactions by combining template-interface-based docking with flexible refinement. Proteins-Structure Function and Bioinformatics. 2012, 80 (4): 1239-1249.

  41. 41.

    Kuzu G, Gursoy A, Nussinov R, Keskin O: Exploiting Conformational Ensembles in Modeling Protein-Protein Interactions on the Proteome Scale. Journal of Proteome Research. 2013, 12 (6): 2641-2653.

  42. 42.

    Ozbabacan SEA, Keskin O, Nussinov R, Gursoy A: Enriching the human apoptosis pathway by predicting the structures of protein-protein complexes. Journal of Structural Biology. 2012, 179 (3): 338-346.

  43. 43.

    Kar G, Keskin O, Nussinov R, Gursoy A: Human proteome-scale structural modeling of E2-E3 interactions exploiting interface motifs. J Proteome Res. 2012, 11 (2): 1196-1207.

  44. 44.

    Kuzu G, Keskin O, Gursoy A, Nussinov R: Constructing structural networks of signaling pathways on the proteome scale. Current Opinion in Structural Biology. 2012, 22 (3): 367-377.

  45. 45.

    Engin HB, Keskin O, Nussinov R, Gursoy A: A strategy based on protein-protein interface motifs may help in identifying drug off-targets. J Chem Inf Model. 2012, 52 (8): 2273-2286.

  46. 46.

    Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, et al: STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013, 41 (Database): D808-815.

  47. 47.

    Roy A, Kucukural A, Zhang Y, I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5 (4): 725-738.

  48. 48.

    Zhang Y, Skolnick J, TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005, 33 (7): 2302-2309.

  49. 49.

    DuPage M, Jacks T: Genetically engineered mouse models of cancer reveal new insights about the antitumor immune response. Current opinion in immunology. 2013, 25 (2): 192-199.

  50. 50.

    Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A, et al: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic acids research. 2011, 39 (Database): D945-950.

  51. 51.

    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al: The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer discovery. 2012, 2 (5): 401-404.

  52. 52.

    Tuncbag N, Keskin O, Gursoy A: HotPoint: hot spot prediction server for protein interfaces. Nucleic acids research. 2010, 38 (Web Server): W402-406.

  53. 53.

    Van Durme J, Delgado J, Stricher F, Serrano L, Schymkowitz J, Rousseau F: A graphical interface for the FoldX forcefield. Bioinformatics. 2011, 27 (12): 1711-1712.

  54. 54.

    Krieger E, Koraimann G, Vriend G: Increasing the precision of comparative models with YASARA NOVA--a self-parameterizing force field. Proteins. 2002, 47 (3): 393-402.

  55. 55.

    Forbes SA, Bhamra G, Bamford S, Dawson E, Kok C, Clements J, Menzies A, Teague JW, Futreal PA, Stratton MR: The Catalogue of Somatic Mutations in Cancer (COSMIC). Curr Protoc Hum Genet. 2008, Chapter 10 (Unit 10 11):

  56. 56.

    Schuster-Bockler B, Bateman A: Protein interactions in human genetic diseases. Genome biology. 2008, 9 (1): R9-

  57. 57.

    Chiti F, Dobson CM: Protein misfolding, functional amyloid, and human disease. Annual review of biochemistry. 2006, 75: 333-366.

  58. 58.

    Wang X, Wei X, Thijssen B, Das J, Lipkin SM, Yu H: Three-dimensional reconstruction of protein networks provides insight into human genetic disease. Nature biotechnology. 2012, 30 (2): 159-164.

  59. 59.

    Engin HB, Guney E, Keskin O, Oliva B, Gursoy A: Integrating structure to protein-protein interaction networks that drive metastasis to brain and lung in breast cancer. PloS one. 2013, 8 (11): e81035-

  60. 60.

    Acuner Ozbabacan SE, Gursoy A, Nussinov R, Keskin O: The structural pathway of interleukin 1 (IL-1) initiated signaling reveals mechanisms of oncogenic mutations and SNPs in inflammation and cancer. PLOS Comp Biol. 10 (2): e1003470-

  61. 61.

    Yoon SI, Logsdon NJ, Sheikh F, Donnelly RP, Walter MR: Conformational changes mediate interleukin-10 receptor 2 (IL-10R2) binding to IL-10 and assembly of the signaling complex. J Biol Chem. 2006, 281 (46): 35088-35096.

  62. 62.

    Muller UC, Zheng H: Physiological functions of APP family proteins. Cold Spring Harbor perspectives in medicine. 2012, 2 (2): a006288-

  63. 63.

    Narita M, Holtzman DM, Schwartz AL, Bu G: Alpha2-macroglobulin complexes with and mediates the endocytosis of beta-amyloid peptide via cell surface low-density lipoprotein receptor-related protein. Journal of neurochemistry. 1997, 69 (5): 1904-1911.

  64. 64.

    Yerbury JJ, Kumita JR, Meehan S, Dobson CM, Wilson MR: alpha2-Macroglobulin and haptoglobin suppress amyloid formation by interacting with prefibrillar protein species. The Journal of biological chemistry. 2009, 284 (7): 4246-4254.

  65. 65.

    Hansel DE, Rahman A, Wehner S, Herzog V, Yeo CJ, Maitra A: Increased expression and processing of the Alzheimer amyloid precursor protein in pancreatic cancer may influence cellular proliferation. Cancer research. 2003, 63 (21): 7032-7037.

  66. 66.

    Meng JY, Kataoka H, Itoh H, Koono M: Amyloid beta protein precursor is involved in the growth of human colon carcinoma cell in vitro and in vivo. International journal of cancer Journal international du cancer. 2001, 92 (1): 31-39.

  67. 67.

    Takayama K, Tsutsumi S, Suzuki T, Horie-Inoue K, Ikeda K, Kaneshiro K, Fujimura T, Kumagai J, Urano T, Sakaki Y, et al: Amyloid precursor protein is a primary androgen target gene that promotes prostate cancer growth. Cancer research. 2009, 69 (1): 137-142.

  68. 68.

    Pietrzik CU, Hoffmann J, Stober K, Chen CY, Bauer C, Otero DA, Roch JM, Herzog V: From differentiation to proliferation: the secretory amyloid precursor protein as a local mediator of growth in thyroid epithelial cells. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (4): 1770-1775.

  69. 69.

    Yousef GM, Chang A, Diamandis EP: Identification and characterization of KLK-L4, a new kallikrein-like gene that appears to be down-regulated in breast cancer tissues. Journal of Biological Chemistry. 2000, 275 (16): 11891-11898.

  70. 70.

    Diamandis EP, Yousef GM: Human tissue kallikreins: A family of new cancer biomarkers. Clinical Chemistry. 2002, 48 (8): 1198-1205.

  71. 71.

    Darling MR, Jackson-Boeters L, Daley TD, Diamandis EP: Human kallikrein 13 expression in salivary gland tumors. International Journal of Biological Markers. 2006, 21 (2): 106-110.

  72. 72.

    Kapadia C, Yousef GM, Mellati AA, Magklara A, Wasney GA, Diamandis EP: Complex formation between human kallikrein 13 and serum protease inhibitors. Clinica Chimica Acta. 2004, 339 (1-2): 157-167.

  73. 73.

    Petraki CD, Gregorakis AK, Papanastasiou PA, Karavana VN, Luo LY, Diamandis EP: Immunohistochemical localization of human kallikreins 6, 10 and 13 in benign and malignant prostatic tissues. Prostate Cancer and Prostatic Diseases. 2003, 6 (3): 223-227.

  74. 74.

    Chang A, Yousef GM, Scorilas A, Grass L, Sismondi P, Ponzone R, Diamandis EP: Human kallikrein gene 13 (KLK13) expression by quantitative RT-PCR: an independent indicator of favourable prognosis in breast cancer. British Journal of Cancer. 2002, 86 (9): 1457-1464.

  75. 75.

    Scorilas A, Borgono CA, Harbeck N, Dorn J, Schmalfeldt B, Schmitt M, Diamandis EP: Human kallikrein 13 protein in ovarian cancer cytosols: A new favorable prognostic marker. Journal of Clinical Oncology. 2004, 22 (4): 678-685.

  76. 76.

    Luo LY, Yousef G, Diamandis EP: Human tissue kallikreins and testicular cancer. Apmis. 2003, 111 (1): 225-233.

  77. 77.

    Planque C, Li L, Zheng YG, Soosaipillai A, Reckamp K, Chia D, Diamandis EP, Goodglick L: A multiparametric serum kallikrein panel for diagnosis of non-small cell lung carcinoma. Clinical Cancer Research. 2008, 14 (5): 1355-1362.

Download references

Acknowledgements

Zhong Chen and Carter Van Waes are funded by NIDCD/NIH Intramural projects Z01-DC-00073 and Z01-DC-00074. This project has been funded in whole or in part with Federal funds from the National Cancer Institute, National Institutes of Health, under contract number HHSN261200800001E. This work is partially supported by the Scientific and Technological Research Council of Turkey (TUBITAK) Grant Number 113E164. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. This research was supported (in part) by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.

Declarations

The publication costs for this article were funded by a grant from NIDCD/NIH Intramural projects Z01-DC-00073 and Z01-DC-00074 and this project has been funded in whole or in part with Federal funds from the National Cancer Institute, National Institutes of Health, under contract number HHSN261200800001E.

This article has been published as part of BMC Genomics Volume 15 Supplement 4, 2014: SNP-SIG 2013: Identification and annotation of genetic variants in the context of structure, function, and disease. The full contents of the supplement are available online at http://0-www.biomedcentral.com.brum.beds.ac.uk/bmcgenomics/supplements/15/S4

Supplementary Files

Table S2.The list of interactions in the IL-10 protein-protein interaction network. 2

Author information

Correspondence to Ruth Nussinov.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SEAO, HBE, EGM, GK, SM devised the study, carried out the modeling and analysis and wrote the paper; AB wrote the new more efficient and capable version of PRISM used here; ZC edited the manuscript and provided IL-10 and related references; ZC and CVW provided the expertise in the IL-10 protein interactions, the oncogenic mutations, and the relation to inflammation and cancer; AG, OK and RN conceived and oversaw the project, and helped in the writing. These authors contributed equally and appear in alphabetical order: Saliha Ece Acuner-Ozbabacan, Hatice Billur Engin, Emine Guven-Maiorov, Guray Kuzu, and Serena Muratcioglu.

Electronic supplementary material

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Acuner-Ozbabacan, E.S., Engin, B.H., Guven-Maiorov, E. et al. The structural network of Interleukin-10 and its implications in inflammation and cancer. BMC Genomics 15, S2 (2014) doi:10.1186/1471-2164-15-S4-S2

Download citation

Keywords

  • Oncogenic Mutation
  • Interface Residue
  • Cosmic Database
  • Bait Region
  • Binding Energy Score