Charting the NF-κB Pathway Interactome Map

Inflammation is part of a complex physiological response to harmful stimuli and pathogenic stress. The five components of the Nuclear Factor κB (NF-κB) family are prominent mediators of inflammation, acting as key transcriptional regulators of hundreds of genes. Several signaling pathways activated by diverse stimuli converge on NF-κB activation, resulting in a regulatory system characterized by high complexity. It is increasingly recognized that the number of components that impinges upon phenotypic outcomes of signal transduction pathways may be higher than those taken into consideration from canonical pathway representations. Scope of the present analysis is to provide a wider, systemic picture of the NF-κB signaling system. Data from different sources such as literature, functional enrichment web resources, protein-protein interaction and pathway databases have been gathered, curated, integrated and analyzed in order to reconstruct a single, comprehensive picture of the proteins that interact with, and participate to the NF-κB activation system. Such a reconstruction shows that the NF-κB interactome is substantially different in quantity and quality of components with respect to canonical representations. The analysis highlights that several neglected but topologically central proteins may play a role in the activation of NF-κB mediated responses. Moreover the interactome structure fits with the characteristics of a bow tie architecture. This interactome is intended as an open network resource available for further development, refinement and analysis.


Introduction
Inflammation, a complex, pervasive process Inflammation is a complex, systemic, multi-scale physiological process necessary to cope with damaging agents and fundamental for survival, involving a variety of cells, organs and organ systems. The complexity of the inflammatory process has escaped reductionist and linear approaches, since it is characterized, among other features, by non-proportional kinetics as well as numerous and nested feedback loops [1]. Apart from its protective, physiological role, inflammation is an important concomitant cause of many major age-associated pathologies, such as cancer, neurodegeneration, Type II Diabetes and metabolic syndromes [2]. To this regard, it has been observed that a systemic, chronic, low-grade inflammatory status characterizes the aging process, and that markers of inflammation increase with age [3][4][5][6][7]. We proposed to indicate this phenomenon as ''inflammaging'' [8][9][10]. In this perspective, we hypothesised that the phenomenon of inflammaging, likely resulting from the chronic exposure to environmental stressors such as chronic viral infections [6,11], can be an important cause of these age-associated diseases, and therefore the genetic control of inflammation appears to be a crucial determinant of longevity, today a pressing topic given the aging of population all over the world.
A pivotal player of the inflammatory response is the NF-kB (nuclear factor kappa-light-chain-enhancer of activated B cells) transcription factor. Able to be activated by a great number of stimuli and to participate in the regulation of hundreds of genes [12,13], NF-kB has a wide spectrum of actions, such as to induce survival and proliferation, and in particular it is considered the master regulator of inflammation [14]. Thus, it is of great interest to deepen the basic knowledge regarding the topology and the dynamics of the signaling pathways and regulatory network underpinning the activation of such a pleiotropic transcription factor. While the literature dedicated to specific aspects of NF-kB functions or regulation is huge, a wider representation of its complex regulating system is still missing.
The NF-kB system NF-kB is a protein complex that both induces and represses gene expression by binding to discrete DNA sequences, known as kB elements, in gene promoters and enhancers [15,16]. In mammalian cells, there are five NF-kB family members, RelA (p65), RelB, c-Rel, p50/p105 (NF-kB1) and p52/p100 (NF-kB2), and different NF-kB complexes are formed from their homo-and hetero-dimers. All proteins of the NF-kB family share a Rel homology domain. Moreover, RelA, RelB, and c-Rel have a transactivation domain in their C-terminus required for gene activation. NF-kB 1 and NF-kB 2 proteins are synthesized as large precursors, p105 and p100, which undergo processing to generate the mature NF-kB subunits, p50 and p52, respectively. The processing of p105 and p100 is mediated by the ubiquitin/ proteasome pathway and involves selective degradation of their Cterminal region containing ankyrin repeats. Whereas the generation of p52 from p100 is a tightly-regulated process, p50 is produced from constitutive processing of p105 [16]. Additionally, the ubiquitin proteasome system controls NF-kB activation also by the degradation of inhibitory proteins, including IkBs members, but there is evidence for different downstream level of NF-kB regulation that employs several non proteolytic mechanisms, including promoter-specific exchange of dimers and modification of the transactivating p65 subunit by post translational modification, such as phosphorylation, acetylation and ubiquitination [17].
In vertebrates, NF-kB is activated by over 150 different stimuli, such as stress, cytokines, free radicals, ultraviolet irradiation, oxidized LDL and bacterial or viral antigens [12]. In turn, there is evidence that active NF-kB participates in the control of transcription of more than 400 genes [13]. These genes include cytokines, chemokines and their modulators, immunoreceptors, proteins involved in antigen presentation, cell adhesion molecules, acute phase proteins, stress response proteins, cell-surface receptors, regulators of apoptosis, growth factors, ligands and their modulators, early response proteins, transcription factors and regulators, and enzymes, controlling several phenomena such as inflammation as well as innate and adaptive immune response [16]. Considering the central role of NF-kB in maintaining cellular homeostasis, it is not surprising that dysregulation of its finely tuned modulation has often been linked to the development of several disorders, most of which with an inflammatory component such as cancer, autoimmune diseases, and chronic inflammatory disorders [18,19].
Twenty-five years of research efforts [20] have not been enough for a sufficient comprehension of NF-kB dynamics [21], and gaining insights into NF-kB and related processes is still a top priority for the understanding of the onset of associated diseases [20].

Moving from NF-kB pathway to NF-kB interactome network
It is increasingly recognized that the number of components that impinge upon phenotypic outcomes of signal transduction pathways may be higher than that taken into consideration from canonical representations [22][23][24]. Indeed, recent screens suggest the participation of hundreds of components, instead of the dozens classically involved in canonical signaling pathway representations [22], proposing a switch from linear or branched signal cascades to networks with complex interdependencies and feedbacks. Moreover, the interpretation of new findings [23] suggests that cell functionality may be based on a single ''mega-network with limited isolation'' [22] between pathways' elements, modules (i.e. a set of proteins capable of independent functionality) and networks. The difficulty in drawing contours and boundaries of signaling pathways becomes more evident when the study of the system is approached with formal engineering failure analysis [25]. Systems biology approaches may provide crucial clues about the architecture and the logic of the regulation of the NF-kB system and related gene expression [20]. Moreover, comprehensive maps of complex signaling pathways are proving to be important tools to facilitate systems-level study [26,27].
Signaling crosstalk and combinatorial control of NF-kB pathway are very complex [28][29][30][31]. In-depth studies often report fragmented knowledge, since they forcedly focus only on a limited part of the puzzle. In this scenario limiting analyses to canonical NF-kB pathway elements may reveal inadequate in unveiling crucial mechanisms or players in the regulation of this system. The aim of this work is to integrate into one single, more comprehensive and less ''pathway-centric'' [22] picture all the proteins that, on the basis of present knowledge, interact directly or indirectly with the fundamental NF-kB-activation proteins to compose the ''NF-kB pathway interactome''. We propose that such approach may give further insights on the contribution of components that would have been neglected in a more canonical perspective. Such a comprehensive approach may provide insights about how specific stimuli trigger particular subsets of NF-kBtarget genes [21], or about the regulatory hierarchy that rules the selection of target gene expression [32]. To do this, we mined from multiple sources and integrated in an unitary view manually collected data from literature, binary protein-protein interaction (PPI) data, protein annotation data, and NF-kB downstream genes data. We compiled three different sets of proteins that show evidence of involvement in the upstream regulation process of NF-kB, and one set of downstream genes (and related proteins) which expression appears to be regulated by NF-kB. Starting from these sets, we reconstructed three NF-kB protein interaction networks (PINs), or interactomes, by using PPI data from multiple databases, and then checked for the existence of physical interactions between these upstream protein sets and the downstream set. The resulting interactomes also underwent functional enrichment and network analysis to highlight their composition and structure.

Materials and Methods
Part of data retrieval, interactome reconstruction and analysis have been performed with the use of the Cytoscape network analysis platform [33][34][35] and following a general workflow also described in [36]. For PPI data and interactome reconstruction, we used the Agile Protein Interaction Database (APID) [37], a comprehensive resource for protein interaction data, automatically accessed by Cytoscape through the dedicated plugin APID2NET [38]. APID integrates in a single web-based tool all known experimentally validated protein-protein interactions from BIND, BioGRID, DIP, HPRD, IntAct and MINT databases. Annotation and pathway data have been collected mainly from the Universal Protein Resource (UniProt) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), accessed through the suites FaTiGO [39] and Babelomics [40,41]. Additional information on disease pathways have been accessed through DisGeNET [42], and the Pharmacogenomics Knowledge Base PharmGKB [43]. All databases and platforms used here are free and open access. The complete protein lists and relative analysis results are provided in tables S1, S2, S3 and dataset S1. Note that we refer to a ''set'' as to a list of proteins, and to ''interactome'' as to a network whose elements are (all or a part of) the proteins of the respective set, enriched with the available protein interaction data. Interactome acronyms are indicated with a capital ''I'' after the acronym of the respective set.
The network analysis have been performed using the Cytoscape platform and the plugins Networkanalyzer [44] and MCODE [45]. Besides classical parameters (node degree and betweenness centrality, among others), we also calculated the network density (the average number of neighbors for each node), and the network centralization. The degree of a node is the number of edges linked to it. In this case, it is the number of physical interactions that a specific protein shows with other proteins. The betweenness centrality of a node attests the grade of control that this node exerts over the interactions of other nodes in the network, and its biological relevance has been demonstrated [46,47]. The average number of neighbours (connectivity) of a node in the network may also be represented by means of a normalized version of this parameter, i.e. the network density. The density is a value between 0 and 1, and shows how densely the network is populated with links. In calculating density, self-loops are ignored (self-loops are evidence based interactions between two copies of the same protein [44]). Network centrality is a parameter which value ranges between 0 and 1. Networks whose topologies resemble a star (i.e. with a relatively well defined central core) have a centrality close to 1, whereas decentralized networks are characterized by having a centrality close to 0. The interactomes are completely available in Cytoscape format (.cys) in the Dataset S1.

Directly Interacting (DI) NF-kB proteins dataset
The directly-NF-kB-interacting protein set, or hereinafter Direct Interactome, (DI) ( fig. 1a), is composed by all the proteins that show experimental evidence of physical interaction with at least one of the five NF-kB members. Given the procedure used (search for directly NF-kB-interacting proteins), in the case of DI, the protein set coincides with the respective interactome. Data retrieval has been performed with APID2NET and Cytoscape. The query returned a total of 377 proteins (including the five NF-kB members). The DI accounts for 4,119 non-directional interactions (including self-interactions).

Uniprot (U) NF-kB annotated proteins dataset
To populate the set U, we searched the UniProt Knowledge Base with the parameters: «annotation:(type: nonpositional ''nf kappa b'') and organism: ''Homo sapiens (Human) [9606]''» and we retrieved a list of 235 proteins that have been manually screened and checked to obtain a final, validated list of 229 proteins with evidence of implications in the NF-kB functioning (see Table S1 for and MCI (c) and the resulting UNION interactome (d) as from Cytoscape layout, available in the Dataset S1 as .cys file. Datasets and interactomes have been reconstructed following procedures described in materials and methods section. Nodes represent proteins and links are evidence-based physical interactions. Node size and color are proportional to betweenness centrality values (red: high, green: low; not comparable among different interactomes). Isolated proteins show no evidence of physical interactions with any other proteins within the same dataset. doi:10.1371/journal.pone.0032678.g001 description). Starting from this set, experimental protein interaction data are retrieved in order to reconstruct the relative interactome. Two hundred and ten proteins (210), out of 229, are present in the APID database, from which interaction data have been downloaded, and constitute the Uniprot interactome (UI, fig. 1b). One hundred and fifty proteins (out of 210) form a main cluster, accounting for 550 interactions, while other 60 are isolated from the main cluster.

Manually Curated (MC) proteins dataset
We selected and then manually screened a collection of 37 top quality, highly cited literature papers regarding NF-kB to identify proteins that take part with different roles and functions to the signalling cascade leading to NF-kB activation [14,17,. The criteria for protein selection and inclusion in the MC set are based on its presence and function described in each paper, either as directly involved in the cascade dynamics (e.g. interacting protein), or collaterally participating with a function identified and described in the articles considered. One hundred and forty one proteins (141) have been identified. Again, PPI data have been added (through the use of the APID database) to build the Manually Curated interactome (MCI), that accounts for 853 nondirectional interactions (including self-interactions) ( fig. 1c). All proteins in the DI, U and MC sets are considered to work upstream NF-kB, i.e. to participate with a variable role to the modulation and to the signalling cascade leading to NF-kBactivation.

Downstream Genes (DG) dataset
Data extracted from a manually curated list of NF-kBdownstream genes (www.nf-kb.org; [13]) and from the Transcriptional Regulatory Elements Database (TRED; [83,84]) constitute a list of 441 genes that result to be up-or down-regulated in response to the activation of at least one of the NF-kB family members (Table S2). Four hundred forty one (441) valid ENSEMBL gene identifiers mapping to 426 related protein unique identifiers (Uniprot IDs) have been obtained using online ID converter tools [85]. The Babelomics suite for functional analysis identified 16 duplicates leading to a valid list of 410 IDs. Three hundred and eighty four (384) out of 426 protein IDs have been found in the APID database, used for DG interactome reconstruction.

Composition and analysis of the protein sets and interactomes
The three sets DI, U and MC (and corresponding interactomes DI, UI and MCI) show substantial differences with one another in their dimensions and structure as well as in the protein composition. We also considered the interactome resulting from the union of DI, UI and MCI (DI < UI < MCI; hereinafter UNION, Figure 1d), which accounts for 622 proteins and 6,115 interactions (see table 1 for the main parameters) in order to have the widest possible picture of the elements participating in the NF-kB system. The intersection of all three sets (the proteins shared by all three sets: DI > UI > MCI) accounts for only 16 proteins and 89 interactions (see fig. 2 and Dataset S1). The intersection should represent the very ''core'' of the system. Indeed it contains the 5 members of the NF-kB family (p65, RelB, c-Rel, NF-kB1, NF-kB2), the 4 members of the inhibitor of NF-kB (IkB) family (IkBa, IkBb, IkBe, BCL-3) and the three core subunits of the IkB kinase (IKK) complex (NEMO, IKKa, IKKb) as reported by Perkins [14]. Other components of this core subset are IKBKE, MAP3K14, TRAF3 and NKRF, which are very well known elements able to activate or inhibit NF-kB [86][87][88][89]. The MCI is the most dense and less populated by self-loops, while the UNION shows high average neighbours value but also many self-loops ( Table 1). The most centralised interactome is the DI, while the UI is the most decentralized (Table 1).
Network analysis showed that, as expected, being these networks centred on NF-kB, the five NF-kB subunits are among the most central elements for betweenness centrality values ( Table 2). This condition is common to all 3 interactomes. The betweenness centrality of a node attests the amount of control exerted by this node over the capability of interaction of other nodes in the network. In particular, four out of five NF-kB components are in the first five positions for betweenness centrality in the UNION set (table 2; complete list available in Dataset S1), while REL is ranked 21 st . Among the highest ranking proteins, TRAF6, at third position, is the element of connection of crucial membrane receptors. It is a member of the TNF receptor associated factor (TRAF) protein family, and mediates signal transduction from members of the TNF receptor superfamily and of the Toll/IL-1 family [90]. The connection between these two receptor families is interesting because the Toll-Like Receptors (TLRs) and TNF receptors (TNFRs) play the role of external sensors for NF-kB pathway activation. TLRs take a crucial part in early host defense against invading pathogens, recognizing many pathogen-associated molecular patterns (PAMPs), including bacterial cell-surface lipopolysaccharides, peptidoglycans and lipoproteins, viral double-stranded RNA and viral single-stranded RNA, bacterial and viral CpG, as recently reviewed [91]; TNFR is the mediator of TNF effects, one of the universal effectors of innate signaling, involved in host defense and inflammation [92]. Table 2 also shows that three broad activators of gene transcription such as Myc, p53 and STAT3 are found as central elements at position 13, 20 and 25, respectively. This reflects the fact that NF-kB is not only important for linking adaptive and innate immunity but it is also at the very centre of fundamental cell defense mechanisms such as stress response, as well as of other basic cellular functions such as survival and proliferation, thus confirming its primordial role as central actor in the protection of a multicellular organism's integrity [21]. Interestingly, the highly conserved protein UBIQ (ubiquitin) ranks at position 9 (Table 2)  that inhibit NF-kB or target it to proteasomal degradation is a further indication of the centrality of this complex, which, due to its pleiotropic effects, must be tightly controlled [93].
We observed a limited degree of overlapping among the sets (Fig. 2): DI and U share only the 9.1% of their elements, DI and MC the 6.4%, U and MC the 13.6%, while only the 2.6% of the elements are shared by all three interactomes (16 out of 622). These numbers partly reflect the difference and dishomogeneity in sources, databases, data types and relative retrieval method that have been used. Nevertheless, we did not expect such remarkable discrepancy in the composition of the sets, discrepancy that might not be totally explained simply taking into account the wide differences in data sources and the integrative approach used, but, on the contrary, should cast some doubts about the adequacy and the completeness of the classical pathway descriptions.
Another example of such inadequacy is the fact that from the analysis of the datasets TAK1 (MAP3K7) did not result to be high ranking, as it would be expected [94].
By using the Babelomics platform, we analyzed all the sets to identify the metabolic or signalling pathways which elements are shared with the given sets. In Table 1 we show the first ten overrepresented pathways in the DI, U, MC and UNION according to this analysis. The complete lists of KEGG significant terms (pvalue,0.05) are available in the Table S3. We noted that 5 KEGG pathways resulted as ranking within the top 10 in DI, U, and MC, and that 4 out of 5 are pathways linked to the induction of innate immunity response and inflammation (Toll-like receptor signaling pathway; RIG-I-like receptor signalling pathway; NODlike receptor signalling pathway; Adipocytokine signalling pathway). This confirms that, even if the sources of these sets are very different, our analysis was able to capture the role of NF-kB as mediator of inflammation. Interestingly, MAPK pathway ranks always high in terms of percentage of proteins present in all three protein sets, indicating that also signal transduction leading to NF-kB activation is a process whose component are captured by our analysis. The UNION set resumes the characteristics of the three and represents, in number of components, 62% of the whole KEGG-annotated proteasome and 41% of the NOD signalling pathway, among others.

Downstream gene and protein set: composition and characteristics
The same KEGG analysis has been done on NF-kB downstream gene set (410 out of 426 identifiers have been recognized by the Babelomics analysis suite). Among the top 10 identified pathways (Table 3), four are involved in inflammation (Toll-like receptor signalling pathway; NOD-like receptor signalling pathway; cytosolic DNA-sensing pathway; cytokine-cytokine receptor interaction), with two of them also shared by the DI, U and MC sets. The analysis has also pointed out that relevant subsets of DG are related to pathologies with well-known inflammatory aetiology, such as bladder and prostate cancers, among others.
We noted that 49 proteins are present in both UNION and DG sets (Table 4), thus meaning that 13% of the identified NF-kBregulated genes express proteins that play a direct role in the UNION interactome ( Figure 3). These data suggest that the implication of such feedback loop in the dynamics of the interactome and in its regulation and activation of NF-kB may be remarkable. The ability of NF-kB to control the transcription of a non negligible part of its activation pathway interactome surely deserves a deeper attention. Since NF-kB is a constitutively expressed protein complex, and since it has a crucial importance in the regulation of fundamental tasks such as the immune responses and inflammation, it must be tightly self-regulated to prevent malfunctions and inappropriate activation. Several feedback loops are already well-known and studied in the literature, while many others have not been taken under consideration in the perspective of the possible implications in the regulation of the NF-kB pathway dynamics. In fact, as already reported in literature and confirmed by our analysis (Table 4), NF-kB complexes directly control the expression of their own subunits with the exception of TF65 [14]. Thus, NF-kB system appears to be very tightly self-regulated.

Discussion
As mentioned in the introduction, an increase of inflammatory markers has been described as a general feature of the ageing process and proposed as a possible cause of many age-associated diseases. Therefore, the system of NF-kB is likely a crucial player in this process. In this perspective, in order to understand how the changes in composition and abundance of the NF-kB interactome can regulate NF-kB activation, and thus (at least in part) the ageing process, the knowledge about NF-kB pathway as well as the proteins that can interact directly or indirectly with it must be broaden. To this purpose, we integrated data from multiple existing sources in order to chart the NF-kB interactome map. Our map may help in improving the understanding of complex signalling networks and highlights the impact of single elements, their feedback and crosstalk regulations on cellular processes [26]. Beside confirming previous studies and insights [22], the integrative effort proposed here in merging existing, freely available data highlights that the number of elements impinging upon NF-kB-activating pathway outcomes is much higher than that usually taken under consideration in canonical pathway representations.
Furthermore, as other pleiotropic factors that receive many inputs and spread out many outputs, NF-kB core elements are at the centre of a network characterised by a bow tie architecture [95], with a fan in (UNION interactome), a knob (the NF-kB family), a fan out (downstream genes modulated by NF-kB) and feedback loops. This is not surprising, since the bow tie might represent an economical and efficient way to convey on the same target different stimuli to which cells and organs are exposed [95,96]. Moreover, bow tie architecture may be an advantageous solution for the evolution of complex systems, because it would minimize the number of central molecular mediators, ultimately reducing energy expenditure for both the integration of stimuli and the determination of the outcome [97,98]. Bow ties are considered to be fundamental in explaining the coexistence of robustness and evolvability in complex systems due to their ability to facilitate control [95,[99][100][101][102]. Another striking feature of bow ties' functioning is their inner resilience characterised by the presence of regulatory feedback loop, i.e. some products of the fan out are also part of the fan in. In this case the NF-kB system fits perfectly with this architecture, as we checked for proteins present in both the UNION interactome and DG set, thus establishing a feedback loop of the type: interaction with NF-kB pathway R transcription factor activation R transcriptional regulation R interaction with NF-kB pathway (Table 4). Concerning such feedbacks, we report three cases that were correctly captured by our analysis. The Tumor suppressor protein p53 -that appears among the ''feedback'' proteins-can be transcriptionally induced by NF-kB [103][104][105], among others. On its side, p53 can interfere with NF-kB activity [106,107] thus creating a sort of negative feedback loop between the two factors. The two factors not only interact with each other, but also cooperate or compete (depending upon circumstances) to activate downstream genes [108,109]. The  analysis also identified a complex cross-regulatory loop involving NF-kB and PTEN, a ''feedback'' protein, that may serve to balance their mutually antagonistic functions. Indeed, PTEN is down-regulated by NF-kB and is able to negatively regulate NF-kB, establishing a complex regulation which determines cell survival or apoptosis [110,111]. The ''feedback'' ubiquitin-editing protein TNAP3 (TNFAIP3) has been described as a key player in the termination of NF-kB signaling and in controlling NF-kBdependent inflammation [112]. In most cell types basal TNAP3 expression is very low but its transcription is rapidly induced upon NF-kB activation. Once expressed, TNAP3 functions as a negative feedback regulator of NF-kB activation [112]. Based upon these evidences, we propose that a comprehensive and systemic consideration of ''feedback'' genes may complement important findings in autoregulatory feedback loops [113] leading to a deeper understanding of the complex NF-kB regulation and dependent inflammation. In this perspective, the observation of a number of feedback loops deserves a deeper systemic and experimental investigation, since these interactions might be potential critical points for the NF-kB system regulation neglected so far. In this direction, results from this interactome reconstruction can be used as preliminary screening to identify putative key regulators of NF-kB system.
In summary, substantial divergences in the composition of the DI, U and MC sets open questions about the adequacy and comprehensiveness of classical pathway descriptions and representations, and suggest the participation of a number of proteins one order of magnitude higher than that classically taken into consideration, thus transforming the concept of pathway from an isolated entity into an open, unbound one. The map as it is introduced here provides a tool to explore the complexity of the NF-kB system and to make useful qualitative predictions regarding key regulators and mediators of its activating signalling cascade, but it is only a first step in this direction. The next steps will be not only the obvious and continuous updates of existing records and newly produced data, but also the charting of the interactome of the different NF-kB components, and the analysis of interacting proteins taking into account their sub-cellular compartmentalization, which is currently one of most promising perspective for the elucidation of the mechanisms related to ageing [114]. To this regards, it is also plausible that NF-kB interactome can be differently charted considering as a variable the five NF-kB members, as well as the sub-cellular localisation of its components (e.g. nucleus, organelles, membranes, cytoplasm).
A further important step will consist in charting cell typespecific, and time-dependent interactomes, able to show maps of various cellular functions, given that NF-kB has not only a major role in inflammatory status but it has also a crucial role in neuronal cells, for instance in the Schwann cell myelination, at least in mouse models [115]. Dynamic analysis of interactomes is the next fascinating aim, crucial for the better understanding of proteinprotein interactions when data from time-resolved proteomics will be completely available, together with an adequate mathematical and algorithmic contribution [116,117]. This approach could be extremely relevant to definitively disentangle when, where and how protein-protein interactions give their contribution to the different functional levels of cell systems (different types of cells, tissues, different responses to various stimuli etc). To this regard, network inference based on the huge amount of data gathered in the Omics era might represent a successful approach [118].
As a whole, the results of the integrative approach shown and reported here are meant to be used as a starting point to identify new roles of NF-kB in physiology and diseases. The integrated analysis we performed identified elements and pathways that are not immediately linked to NF-kB itself, or not taken into consideration in respect to its regulation and dynamics, and may thus suggest new directions for further studies and analyses.

Supporting Information
Table S1 Uniprot annotations retrieved and used to populate protein dataset U. (XLS) Table S2 gene identifiers extracted from a manually curated list of NF-kB-downstream genes [13] and from the Transcriptional Regulatory Elements Database (TRED; [83,84]), and related protein unique identifiers (Uniprot IDs) obtained using online ID converter tools [85] used to compile DG dataset. (XLS) Dataset S1 complete Cytoscape network session with full datasets and interactomes. (RAR) Figure 3. The UNION interactome is composed by 622 proteins, including the five NF-kB subunits. NF-kB is able to regulate the expression of 426 proteins (DG set, see main text). A subset (384 proteins, present in the APID database, out of 426) has been checked for PPIs (see Dataset S1). Forty nine proteins are shared by both the UNION and the DG sets, establishing a feedback loop: ''interaction with NF-kB pathway R transcription factor activation R transcriptional regulation R interaction with NF-kB pathway'', meaning that 13% of the identified NF-kB-regulated genes express proteins that play a direct role in the UNION interactome. doi:10.1371/journal.pone.0032678.g003