Virtual identification of essential proteins within the protein interaction network of yeast

Topological analysis of large scale protein-protein interaction networks (PINs) is important for understanding the organizational and functional principles of individual proteins. The number of interactions that a protein has in a PIN has been observed to be correlated with its indispensability. Essential proteins generally have more interactions than the nonessential ones. We show here that the lethality associated with removal of a protein from the yeast proteome correlates with different centrality measures of the nodes in the PIN, such as the closeness of a protein to many other proteins, or the number of pairs of proteins which need a specific protein as an intermediary in their communications, or the participation of a protein in different protein clusters in the PIN. These measures are significantly better than random selection in identifying essential proteins in a PIN. Centrality measures based on graph spectral properties of the network, in particular the subgraph centrality, show the best performance in identifying essential proteins in the yeast PIN. Subgraph centrality gives important structural information about the role of individual proteins, and permits the selection of possible targets for rational drug discovery through the identification of essential proteins in the PIN.     

Mapping the protein interaction network in methicillin-resistant Staphylococcus aureus

Mortality attributable to infection with methicillin-resistant Staphylococcus aureus (MRSA) has now overtaken the death rate for AIDS in the United States, and advances in research are urgently needed to address this challenge. We report the results of the systematic identification of protein-protein interactions for the hospital-acquired strain MRSA-252. Using a high-throughput pull-down strategy combined with quantitative proteomics to distinguish specific from nonspecific interactors, we identified 13,219 interactions involving 608 MRSA proteins. Consecutive analyses revealed that this protein interaction network (PIN) exhibits scale-free organization with the characteristic presence of highly connected hub proteins. When clinical and experimental antimicrobial targets were queried in the network, they were generally found to occupy peripheral positions in the PIN with relatively few interacting partners. In contrast, the hub proteins identified in this MRSA PIN that are essential for network integrity and stability have largely been overlooked as drug targets. Thus, this empirical MRSA-252 PIN provides a rich source for identifying critical proteins essential for network stability, many of which can be considered as prospective antimicrobial drug targets.     

Systematic Differences in Signal Emitting and Receiving Revealed by PageRank Analysis of a Human Protein Interactome

Most protein PageRank studies do not use signal flow direction information in protein interactions because this information was not readily available in large protein databases until recently. Therefore, four questions have yet to be answered: A) What is the general difference between signal emitting and receiving in a protein interactome? B) Which proteins are among the top ranked in directional ranking? C) Are high ranked proteins more evolutionarily conserved than low ranked ones? D) Do proteins with similar ranking tend to have similar subcellular locations? In this study, we address these questions using the forward, reverse, and non-directional PageRank approaches to rank an information-directional network of human proteins and study their evolutionary conservation. The forward ranking gives credit to information receivers, reverse ranking to information emitters, and non-directional ranking mainly to the number of interactions. The protein lists generated by the forward and non-directional rankings are highly correlated, but those by the reverse and non-directional rankings are not. The results suggest that the signal emitting/receiving system is characterized by key-emittings and relatively even receivings in the human protein interactome. Signaling pathway proteins are frequent in top ranked ones. Eight proteins are both informational top emitters and top receivers. Top ranked proteins, except a few species-related novel-function ones, are evolutionarily well conserved. Protein-subunit ranking position reflects subunit function. These results demonstrate the usefulness of different PageRank approaches in characterizing protein networks and provide insights to protein interaction in the cell.     

Combining functional and topological properties to identify core modules in protein interaction networks

Advances in large-scale technologies in proteomics, such as yeast two-hybrid screening and mass spectrometry, have made it possible to generate large Protein Interaction Networks (PINs). Recent methods for identifying dense sub-graphs in such networks have been based solely on graph theoretic properties. Therefore, there is a need for an approach that will allow us to combine domain-specific knowledge with topological properties to generate functionally relevant sub-graphs from large networks. This article describes two alternative network measures for analysis of PINs, which combine functional information with topological properties of the networks. These measures, called weighted clustering coefficient and weighted average nearest-neighbors degree, use weights representing the strengths of interactions between the proteins, calculated according to their semantic similarity, which is based on the Gene Ontology terms of the proteins. We perform a global analysis of the yeast PIN by systematically comparing the weighted measures with their topological counterparts. To show the usefulness of the weighted measures, we develop an algorithm for identification of functional modules, called SWEMODE (Semantic WEights for MODule Elucidation), that identifies dense sub-graphs containing functionally similar proteins. The proposed method is based on the ranking of nodes, i.e., proteins, according to their weighted neighborhood cohesiveness. The highest ranked nodes are considered as seeds for candidate modules. The algorithm then iterates through the neighborhood of each seed protein, to identify densely connected proteins with high functional similarity, according to the chosen parameters. Using a yeast two-hybrid data set of experimentally determined protein-protein interactions, we demonstrate that SWEMODE is able to identify dense clusters containing proteins that are functionally similar. Many of the identified modules correspond to known complexes or subunits of these complexes.     

Characterization of protein hubs by inferring interacting motifs from protein interactions

The characterization of protein interactions is essential for understanding biological systems. While genome-scale methods are available for identifying interacting proteins, they do not pinpoint the interacting motifs (e.g., a domain, sequence segments, a binding site, or a set of residues). Here, we develop and apply a method for delineating the interacting motifs of hub proteins (i.e., highly connected proteins). The method relies on the observation that proteins with common interaction partners tend to interact with these partners through a common interacting motif. The sole input for the method are binary protein interactions; neither sequence nor structure information is needed. The approach is evaluated by comparing the inferred interacting motifs with domain families defined for 368 proteins in the Structural Classification of Proteins (SCOP). The positive predictive value of the method for detecting proteins with common SCOP families is 75% at sensitivity of 10%. Most of the inferred interacting motifs were significantly associated with sequence patterns, which could be responsible for the common interactions. We find that yeast hubs with multiple interacting motifs are more likely to be essential than hubs with one or two interacting motifs, thus rationalizing the previously observed correlation between essentiality and the number of interacting partners of a protein. We also find that yeast hubs with multiple interacting motifs evolve slower than the average protein, contrary to the hubs with one or two interacting motifs. The proposed method will help us discover unknown interacting motifs and provide biological insights about protein hubs and their roles in interaction networks.     

A Novel Scoring Approach for Protein Co-Purification Data Reveals High Interaction Specificity

Large-scale protein interaction networks (PINs) have typically been discerned using affinity purification followed by mass spectrometry (AP/MS) and yeast two-hybrid (Y2H) techniques. It is generally recognized that Y2H screens detect direct binary interactions while the AP/MS method captures co-complex associations; however, the latter technique is known to yield prevalent false positives arising from a number of effects, including abundance. We describe a novel approach to compute the propensity for two proteins to co-purify in an AP/MS data set, thereby allowing us to assess the detected level of interaction specificity by analyzing the corresponding distribution of interaction scores. We find that two recent AP/MS data sets of yeast contain enrichments of specific, or high-scoring, associations as compared to commensurate random profiles, and that curated, direct physical interactions in two prominent data bases have consistently high scores. Our scored interaction data sets are generally more comprehensive than those of previous studies when compared against four diverse, high-quality reference sets. Furthermore, we find that our scored data sets are more enriched with curated, direct physical associations than Y2H sets. A high-confidence protein interaction network (PIN) derived from the AP/MS data is revealed to be highly modular, and we show that this topology is not the result of misrepresenting indirect associations as direct interactions. In fact, we propose that the modularity in Y2H data sets may be underrepresented, as they contain indirect associations that are significantly enriched with false negatives. The AP/MS PIN is also found to contain significant assortative mixing; however, in line with a previous study we confirm that Y2H interaction data show weak disassortativeness, thus revealing more clearly the distinctive natures of the interaction detection methods. We expect that our scored yeast data sets are ideal for further biological discovery and that our scoring system will prove useful for other AP/MS data sets.     

Coevolution is a short-distance force at the protein interaction level and correlates with the modular organization of protein networks

We investigated what roles coevolution plays in shaping yeast protein interaction network (PIN). We found that the extent of coevolution between two proteins decreases rapidly as their interacting distance on the PIN increases, suggesting coevolutionary constraint is a short-distance force at the molecular level. We also found that protein-protein interactions (PPIs) with strong coevolution tend to be enriched in interconnected clusters, whereas PPIs with weak coevolution are more frequently present at inter-cluster region. The findings indicate the close relationship between coevolution and modular organization of PINs, and may provide insights into evolution and modularity of cellular networks.     

Evolutionarily Conserved Network Properties of Intrinsically Disordered Proteins

BackgroundIntrinsically disordered proteins (IDPs) lack a stable tertiary structure in isolation. Remarkably, however, a substantial portion of IDPs undergo disorder-to-order transitions upon binding to their cognate partners. Structural flexibility and binding plasticity enable IDPs to interact with a broad range of partners. However, the broader network properties that could provide additional insights into the functional role of IDPs are not known.ResultsHere, we report the first comprehensive survey of network properties of IDP-induced sub-networks in multiple species from yeast to human. Our results show that IDPs exhibit greater-than-expected modularity and are connected to the rest of the protein interaction network (PIN) via proteins that exhibit the highest betweenness centrality and connect to fewer-than-expected IDP communities, suggesting that they form critical communication links from IDP modules to the rest of the PIN. Moreover, we found that IDPs are enriched at the top level of regulatory hierarchy.ConclusionOverall, our analyses reveal coherent and remarkably conserved IDP-centric network properties, namely, modularity in IDP-induced network and a layer of critical nodes connecting IDPs with the rest of the PIN.     

A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory

Background

The study of biological interaction networks is a central theme of systems biology. Here, we investigate the relationships between two distinct types of interaction networks: the metabolic pathway map and the protein-protein interaction network (PIN). It has long been established that successive enzymatic steps are often catalyzed by physically interacting proteins forming permanent or transient multi-enzymes complexes. Inspecting high-throughput PIN data, it was shown recently that, indeed, enzymes involved in successive reactions are generally more likely to interact than other protein pairs. In our study, we expanded this line of research to include comparisons of the underlying respective network topologies as well as to investigate whether the spatial organization of enzyme interactions correlates with metabolic efficiency.

Results

Analyzing yeast data, we detected long-range correlations between shortest paths between proteins in both network types suggesting a mutual correspondence of both network architectures. We discovered that the organizing principles of physical interactions between metabolic enzymes differ from the general PIN of all proteins. While physical interactions between proteins are generally dissortative, enzyme interactions were observed to be assortative. Thus, enzymes frequently interact with other enzymes of similar rather than different degree. Enzymes carrying high flux loads are more likely to physically interact than enzymes with lower metabolic throughput. In particular, enzymes associated with catabolic pathways as well as enzymes involved in the biosynthesis of complex molecules were found to exhibit high degrees of physical clustering. Single proteins were identified that connect major components of the cellular metabolism and may thus be essential for the structural integrity of several biosynthetic systems.

Conclusion

Our results reveal topological equivalences between the protein interaction network and the metabolic pathway network. Evolved protein interactions may contribute significantly towards increasing the efficiency of metabolic processes by permitting higher metabolic fluxes. Thus, our results shed further light on the unifying principles shaping the evolution of both the functional (metabolic) as well as the physical interaction network.     

Interactions among PIN-FORMED and P-glycoprotein auxin transporters in Arabidopsis

Directional transport of the phytohormone auxin is established primarily at the point of cellular efflux and is required for the establishment and maintenance of plant polarity. Studies in whole plants and heterologous systems indicate that PIN-FORMED (PIN) and P-glycoprotein (PGP) transport proteins mediate the cellular efflux of natural and synthetic auxins. However, aromatic anion transport resulting from PGP and PIN expression in nonplant systems was also found to lack the high level of substrate specificity seen in planta. Furthermore, previous reports that PGP19 stabilizes PIN1 on the plasma membrane suggested that PIN-PGP interactions might regulate polar auxin efflux. Here, we show that PGP1 and PGP19 colocalized with PIN1 in the shoot apex in Arabidopsis thaliana and with PIN1 and PIN2 in root tissues. Specific PGP-PIN interactions were seen in yeast two-hybrid and coimmunoprecipitation assays. PIN-PGP interactions appeared to enhance transport activity and, to a greater extent, substrate/inhibitor specificities when coexpressed in heterologous systems. By contrast, no interactions between PGPs and the AUXIN1 influx carrier were observed. Phenotypes of pin and pgp mutants suggest discrete functional roles in auxin transport, but pin pgp mutants exhibited phenotypes that are both additive and synergistic. These results suggest that PINs and PGPs characterize coordinated, independent auxin transport mechanisms but also function interactively in a tissue-specific manner.     

Destabilizing interactions among [PSI(+)] and [PIN(+)] yeast prion variants

The yeast Sup35 and Rnq1 proteins can exist in either the noninfectious soluble forms, [psi-] or [pin-], respectively, or the multiple infectious amyloid-like forms called [PSI+] or [PIN+] prion variants (or prion strains). It was previously shown that [PSI+] and [PIN+] prions enhance one another's de novo appearance. Here we show that specific prion variants of [PSI+] and [PIN+] disrupt each other's stable inheritance. Acquiring [PSI+] often impedes the inheritance of particular [PIN+] variants. Conversely, the presence of some [PIN+] variants impairs the inheritance of weak [PSI+] but not strong [PSI+] variants. These same [PIN+] variants generate a single-dot fluorescence pattern when a fusion of Rnq1 and green fluorescent protein is expressed. Another [PIN+] variant, which forms a distinctly different multiple-dot fluorescence pattern, does not impair [PSI+] inheritance. Thus, destabilization of prions by heterologous prions depends upon the variants involved. These findings may have implications for understanding interactions among other amyloid-forming proteins, including those associated with certain human diseases.     

Mining Temporal Protein Complex Based on the Dynamic PIN Weighted with Connected Affinity and Gene Co-Expression

The identification of temporal protein complexes would make great contribution to our knowledge of the dynamic organization characteristics in protein interaction networks (PINs). Recent studies have focused on integrating gene expression data into static PIN to construct dynamic PIN which reveals the dynamic evolutionary procedure of protein interactions, but they fail in practice for recognizing the active time points of proteins with low or high expression levels. We construct a Time-Evolving PIN (TEPIN) with a novel method called Deviation Degree, which is designed to identify the active time points of proteins based on the deviation degree of their own expression values. Owing to the differences between protein interactions, moreover, we weight TEPIN with connected affinity and gene co-expression to quantify the degree of these interactions. To validate the efficiencies of our methods, ClusterONE, CAMSE and MCL algorithms are applied on the TEPIN, DPIN (a dynamic PIN constructed with state-of-the-art three-sigma method) and SPIN (the original static PIN) to detect temporal protein complexes. Each algorithm on our TEPIN outperforms that on other networks in terms of match degree, sensitivity, specificity, F-measure and function enrichment etc. In conclusion, our Deviation Degree method successfully eliminates the disadvantages which exist in the previous state-of-the-art dynamic PIN construction methods. Moreover, the biological nature of protein interactions can be well described in our weighted network. Weighted TEPIN is a useful approach for detecting temporal protein complexes and revealing the dynamic protein assembly process for cellular organization.     

A lock-and-key model for protein-protein interactions

MOTIVATION: Protein-protein interaction networks are one of the major post-genomic data sources available to molecular biologists. They provide a comprehensive view of the global interaction structure of an organism's proteome, as well as detailed information on specific interactions. Here we suggest a physical model of protein interactions that can be used to extract additional information at an intermediate level: It enables us to identify proteins which share biological interaction motifs, and also to identify potentially missing or spurious interactions. RESULTS: Our new graph model explains observed interactions between proteins by an underlying interaction of complementary binding domains (lock-and-key model). This leads to a novel graph-theoretical algorithm to identify bipartite subgraphs within protein-protein interaction networks where the underlying data are taken from yeast two-hybrid experimental results. By testing on synthetic data, we demonstrate that under certain modelling assumptions, the algorithm will return correct domain information about each protein in the network. Tests on data from various model organisms show that the local and global patterns predicted by the model are indeed found in experimental data. Using functional and protein structure annotations, we show that bipartite subnetworks can be identified that correspond to biologically relevant interaction motifs. Some of these are novel and we discuss an example involving SH3 domains from the Saccharomyces cerevisiae interactome. AVAILABILITY: The algorithm (in Matlab format) is available (see http://www.maths.strath.ac.uk/~aas96106/lock_key.html).     

Delineation of the interaction domains of Agrobacterium tumefaciens VirB7 and VirB9 by use of the yeast two-hybrid assay.

The Agrobacterium tumefaciens VirB proteins are postulated to form a transport pore for the transfer of T-DNA. Formation of the transport pore will involve interactions among the VirB proteins. A powerful genetic method to study protein-protein interaction is the yeast two-hybrid assay. To test whether this method can be used to study interactions among the VirB membrane proteins, we studied the interaction of VirB7 and VirB9 in yeast. We recently demonstrated that VirB7 and VirB9 form a protein complex linked by a disulfide bond between cysteine 24 of VirB7 and cysteine 262 of VirB9 (L. Anderson, A. Hertzel, and A. Das, Proc. Natl. Acad. Sci. USA 93:8889-8894, 1996). We now demonstrate that VirB7 and VirB9 interact in yeast, and this interaction does not require the cysteine residues essential for the disulfide linkage. By using defined segments in fusion constructions, we mapped the VirB7 interaction domain of VirB9 to residues 173 to 275. In tumor formation assays, both virB7C24S and virB9C262S expressed from a multicopy plasmid complemented the respective deletion mutation, indicating that the cysteine residues may not be essential for DNA transfer.     

Peeling the yeast protein network

A set of highly connected proteins (or hubs) plays an important role for the integrity of the protein interaction network of Saccharomyces cerevisae by connecting the network's intrinsic modules. The importance of the hubs' central placement is further confirmed by their propensity to be lethal. However, although highly emphasized, little is known about the topological coherence among the hubs. Applying a core decomposition method which allows us to identify the inherent layer structure of the protein interaction network, we find that the probability of nodes both being essential and evolutionary conserved successively increases toward the innermost cores. While connectivity alone is often not a sufficient criterion to assess a protein's functional, evolutionary and topological relevance, we classify nodes as globally and locally central depending on their appearance in the inner or outer cores. The observation that globally central proteins participate in a substantial number of protein complexes which display an elevated degree of evolutionary conservation allows us to hypothesize that globally central proteins serve as the evolutionary backbone of the proteome. Even though protein interaction data are extensively flawed, we find that our results are very robust against inaccurately determined protein interactions.     

The PIN domain endonuclease Utp24 cleaves pre-ribosomal RNA at two coupled sites in yeast and humans

During ribosomal RNA (rRNA) maturation, cleavages at defined sites separate the mature rRNAs from spacer regions, but the identities of several enzymes required for 18S rRNA release remain unknown. PilT N-terminus (PIN) domain proteins are frequently endonucleases and the PIN domain protein Utp24 is essential for early cleavages at three pre-rRNA sites in yeast (A0, A1 and A2) and humans (A0, 1 and 2a). In yeast, A1 is cleaved prior to A2 and both cleavages require base-pairing by the U3 snoRNA to the central pseudoknot elements of the 18S rRNA. We found that yeast Utp24 UV-crosslinked in vivo to U3 and the pseudoknot, placing Utp24 close to cleavage at site A1. Yeast and human Utp24 proteins exhibited in vitro endonuclease activity on an RNA substrate containing yeast site A2. Moreover, an intact PIN domain in human UTP24 was required for accurate cleavages at sites 1 and 2a in vivo, whereas mutation of another potential site 2a endonuclease, RCL1, did not affect 18S production. We propose that Utp24 cleaves sites A1/1 and A2/2a in yeast and human cells.     

Dynamical analysis of yeast protein interaction network during the sake brewing process

Proteins interact with each other for performing essential functions of an organism. They change partners to get involved in various processes at different times or locations. Studying variations of protein interactions within a specific process would help better understand the dynamic features of the protein interactions and their functions. We studied the protein interaction network of Saccharomyces cerevisiae (yeast) during the brewing of Japanese sake. In this process, yeast cells are exposed to several stresses. Analysis of protein interaction networks of yeast during this process helps to understand how protein interactions of yeast change during the sake brewing process. We used gene expression profiles of yeast cells for this purpose. Results of our experiments revealed some characteristics and behaviors of yeast hubs and non-hubs and their dynamical changes during the brewing process. We found that just a small portion of the proteins (12.8 to 21.6%) is responsible for the functional changes of the proteins in the sake brewing process. The changes in the number of edges and hubs of the yeast protein interaction networks increase in the first stages of the process and it then decreases at the final stages.     

A secreted protein microarray platform for extracellular protein interaction discovery

Characterization of the extracellular protein interactome has lagged far behind that of intracellular proteins, where mass spectrometry and yeast two-hybrid technologies have excelled. Improved methods for identifying receptor-ligand and extracellular matrix protein interactions will greatly accelerate biological discovery in cell signaling and cellular communication. These technologies must be able to identify low-affinity binding events that are often observed between membrane-bound coreceptor molecules during cell-cell or cell-extracellular matrix contact. Here we demonstrate that functional protein microarrays are particularly well-suited for high-throughput screening of extracellular protein interactions. To evaluate the performance of the platform, we screened a set of 89 immunoglobulin (Ig)-type receptors against a highly diverse extracellular protein microarray with 686 genes represented. To enhance detection of low-affinity interactions, we developed a rapid method to assemble bait Fc fusion proteins into multivalent complexes using protein A microbeads. Based on these screens, we developed a statistical methodology for hit calling and identification of nonspecific interactions on protein microarrays. We found that the Ig receptor interactions identified using our methodology are highly specific and display minimal off-target binding, resulting in a 70% true-positive to false-positive hit ratio. We anticipate that these methods will be useful for a wide variety of functional protein microarray users.     

Proteomic parsimony through bipartite graph analysis improves accuracy and transparency

Assembling peptides identified from LC-MS/MS spectra into a list of proteins is a critical step in analyzing shotgun proteomics data. As one peptide sequence can be mapped to multiple proteins in a database, na?ve protein assembly can substantially overstate the number of proteins found in samples. We model the peptide-protein relationships in a bipartite graph and use efficient graph algorithms to identify protein clusters with shared peptides and to derive the minimal list of proteins. We test the effects of this parsimony analysis approach using MS/MS data sets generated from a defined human protein mixture, a yeast whole cell extract, and a human serum proteome after MARS column depletion. The results demonstrate that the bipartite parsimony technique not only simplifies protein lists but also improves the accuracy of protein identification. We use bipartite graphs for the visualization of the protein assembly results to render the parsimony analysis process transparent to users. Our approach also groups functionally related proteins together and improves the comprehensibility of the results. We have implemented the tool in the IDPicker package. The source code and binaries for this protein assembly pipeline are available under Mozilla Public License at the following URL: http://www.mc.vanderbilt.edu/msrc/bioinformatics/.