Protein interaction mapping in C. elegans using proteins involved in vulval development

Integrative structure modeling of macromolecular assemblies from proteomics data

Proteomics techniques have been used to generate comprehensive lists of protein interactions in a number of species. However, relatively little is known about how these interactions result in functional multiprotein complexes. This gap can be bridged by combining data from proteomics experiments with data from established structure determination techniques. Correspondingly, integrative computational methods are being developed to provide descriptions of protein complexes at varying levels of accuracy and resolution, ranging from complex compositions to detailed atomic structures.

From Experimental Approaches to Computational Techniques: A Review on the Prediction of Protein-Protein Interactions

A crucial step towards understanding the properties of cellular systems in organisms is to map their network of protein-protein interactions (PPIs) on a proteomic-wide scale completely and as accurately as possible. Uncovering the diverse function of proteins and their interactions within the cell may improve our understanding of disease and provide a basis for the development of novel therapeutic approaches. The development of large-scale high-throughput experiments has resulted in the production of a large volume of data which has aided in the uncovering of PPIs. However, these data are often erroneous and limited in interactome coverage. Therefore, additional experimental and computational methods are required to accelerate the discovery of PPIs. This paper provides a review on the prediction of PPIs addressing key prediction principles and highlighting the common experimental and computational techniques currently employed to infer PPI networks along with relevant studies in the area.

Functional modules in the Arabidopsis core cell cycle binary protein-protein interaction network

As in other eukaryotes, cell division in plants is highly conserved and regulated by cyclin-dependent kinases (CDKs) that are themselves predominantly regulated at the posttranscriptional level by their association with proteins such as cyclins. Although over the last years the knowledge of the plant cell cycle has considerably increased, little is known on the assembly and regulation of the different CDK complexes. To map protein-protein interactions between core cell cycle proteins of Arabidopsis thaliana, a binary protein-protein interactome network was generated using two complementary high-throughput interaction assays, yeast two-hybrid and bimolecular fluorescence complementation. Pairwise interactions among 58 core cell cycle proteins were tested, resulting in 357 interactions, of which 293 have not been reported before. Integration of the binary interaction results with cell cycle phase-dependent expression information and localization data allowed the construction of a dynamic interaction network. The obtained interaction map constitutes a framework for further in-depth analysis of the cell cycle machinery.

Building and analyzing protein interactome networks by cross-species comparisons

Background

A genomic catalogue of protein-protein interactions is a rich source of information, particularly for exploring the relationships between proteins. Numerous systems-wide and small-scale experiments have been conducted to identify interactions; however, our knowledge of all interactions for any one species is incomplete, and alternative means to expand these network maps is needed. We therefore took a comparative biology approach to predict protein-protein interactions across five species (human, mouse, fly, worm, and yeast) and developed InterologFinder for research biologists to easily navigate this data. We also developed a confidence score for interactions based on available experimental evidence and conservation across species.

Results

The connectivity of the resultant networks was determined to have scale-free distribution, small-world properties, and increased local modularity, indicating that the added interactions do not disrupt our current understanding of protein network structures. We show examples of how these improved interactomes can be used to analyze a genome-scale dataset (RNAi screen) and to assign new function to proteins. Predicted interactions within this dataset were tested by co-immunoprecipitation, resulting in a high rate of validation, suggesting the high quality of networks produced.

Conclusions

Protein-protein interactions were predicted in five species, based on orthology. An InteroScore, a score accounting for homology, number of orthologues with evidence of interactions, and number of unique observations of interactions, is given to each known and predicted interaction. Our website http://www.interologfinder.org provides research biologists intuitive access to this data.

A unifying view of 21st century systems biology

The idea that multi-scale dynamic complex systems formed by interacting macromolecules and metabolites, cells, organs and organisms underlie some of the most fundamental aspects of life was proposed by a few visionaries half a century ago. We are witnessing a powerful resurgence of this idea made possible by the availability of nearly complete genome sequences, ever improving gene annotations and interactome network maps, the development of sophisticated informatic and imaging tools, and importantly, the use of engineering and physics concepts such as control and graph theory. Alongside four other fundamental "great ideas" as suggested by Sir Paul Nurse, namely, the gene, the cell, the role of chemistry in biological processes, and evolution by natural selection, systems-level understanding of "What is Life" may materialize as one of the major ideas of biology.

Gene-centered regulatory networks

Differential gene expression plays a critical role in the development and physiology of multicellular organisms. At a 'systems level' (e.g. at the level of a tissue, organ or whole organism), this process can be studied using gene regulatory network (GRN) models that capture physical and regulatory interactions between genes and their regulators. In the past years, significant progress has been made toward the mapping of GRNs using a variety of experimental and computational approaches. Here, we will discuss gene-centered approaches that we employed to characterize GRNs and describe insights that we have obtained into the global design principles of gene regulation in complex metazoan systems.

PutidaNET: interactome database service and network analysis of Pseudomonas putida KT2440

Background

Pseudomonas putida KT2440 (P. putida KT2440) is a highly versatile saprophytic soil bacterium. It is a certified bio-safety host for transferring foreign genes. Therefore, the bacterium is used as a model organism for genetic and physiological studies and for the development of biotechnological applications. In order to provide a more systematic application of the organism, we have constructed a protein-protein interaction (PPI) network analysis system of P. putida KT2440.

Results

PutidaNET is a comprehensive interaction database and server of P. putida KT2440 which is generated from three protein-protein interaction (PPI) methods. We used PSIMAP (Protein Structural Interactome MAP), PEIMAP (Protein Experimental Interactome MAP), and Domain-domain interactions using iPfam. PutidaNET contains 3,254 proteins, and 82,019 possible interactions consisting of 61,011 (PSIMAP), 4,293 (PEIMAP), and 30,043 (iPfam) interaction pairs except for self interaction. Also, we performed a case study by integrating a protein interaction network and experimental 1-DE/MS-MS analysis data P. putida. We found that 1) major functional modules are involved in various metabolic pathways and ribosomes, and 2) existing PPI sub-networks that are specific to succinate or benzoate metabolism are not in the center as predicted.

Conclusion

We introduce the PutidaNET which provides predicted interaction partners and functional analyses such as physicochemical properties, KEGG pathway assignment, and Gene Ontology mapping of P. putida KT2440 PutidaNET is freely available at http://sequenceome.kobic.kr/PutidaNET.

'Edgetic' perturbation of a C. elegans BCL2 ortholog

Genes and gene products do not function in isolation but within highly interconnected “interactome” networks, modeled as graphs of nodes and edges representing macromolecules and interactions between them, respectively. We propose to investigate genotype-phenotype associations by methodical use of alleles that lack single interactions, while retaining all others, in contrast to genetic approaches designed to eliminate gene products completely. We describe an integrated strategy based on the reverse yeast two-hybrid system to isolate and characterize such edge-specific, or “edgetic” alleles. We establish a proof-of-concept with CED-9, a C. elegans BCL2 ortholog involved in apoptosis. Using ced-9 edgetic alleles, we uncover a new potential functional link between apoptosis and a centrosomal protein, demonstrating both the interest and efficiency of our strategy. This approach is amenable to higher throughput and is particularly applicable to interactome network analysis in organisms for which transgenesis is straightforward.

Accounting for redundancy when integrating gene interaction databases

During the last years gene interaction networks are increasingly being used for the assessment and interpretation of biological measurements. Knowledge of the interaction partners of an unknown protein allows scientists to understand the complex relationships between genetic products, helps to reveal unknown biological functions and pathways, and get a more detailed picture of an organism's complexity. Being able to measure all protein interactions under all relevant conditions is virtually impossible. Hence, computational methods integrating different datasets for predicting gene interactions are needed. However, when integrating different sources one has to account for the fact that some parts of the information may be redundant, which may lead to an overestimation of the true likelihood of an interaction. Our method integrates information derived from three different databases (Bioverse, HiMAP and STRING) for predicting human gene interactions. A Bayesian approach was implemented in order to integrate the different data sources on a common quantitative scale. An important assumption of the Bayesian integration is independence of the input data (features). Our study shows that the conditional dependency cannot be ignored when combining gene interaction databases that rely on partially overlapping input data. In addition, we show how the correlation structure between the databases can be detected and we propose a linear model to correct for this bias. Benchmarking the results against two independent reference data sets shows that the integrated model outperforms the individual datasets. Our method provides an intuitive strategy for weighting the different features while accounting for their conditional dependencies.

Large-scale RACE approach for proactive experimental definition of C. elegans ORFeome

Although a highly accurate sequence of the Caenorhabditis elegans genome has been available for 10 years, the exact transcript structures of many of its protein-coding genes remain unsettled. Approximately two-thirds of the ORFeome has been verified reactively by amplifying and cloning computationally predicted transcript models; still a full third of the ORFeome remains experimentally unverified. To fully identify the protein-coding potential of the worm genome including transcripts that may not satisfy existing heuristics for gene prediction, we developed a computational and experimental platform adapting rapid amplification of cDNA ends (RACE) for large-scale structural transcript annotation. We interrogated 2000 unverified protein-coding genes using this platform. We obtained RACE data for approximately two-thirds of the examined transcripts and reconstructed ORF and transcript models for close to 1000 of these. We defined untranslated regions, identified new exons, and redefined previously annotated exons. Our results show that as much as 20% of the C. elegans genome may be incorrectly annotated. Many annotation errors could be corrected proactively with our large-scale RACE platform.

Predicting protein-protein interactions in the context of protein evolution

Here we review the methods for the prediction of protein interactions and the ideas in protein evolution that relate to them. The evolutionary assumptions implicit in many of the protein interaction prediction methods are elucidated. We draw attention to the caution needed in deploying certain evolutionary assumptions, in particular cross-organism transfer of interactions by sequence homology, and discuss the known issues in deriving interaction predictions from evidence of co-evolution. We also conject that there is evolutionary knowledge yet to be exploited in the prediction of interactions, in particular the heterogeneity of interactions, the increasing availability of interaction data from multiple species, and the models of protein interaction network growth.

Regulation of endosomal clathrin and retromer-mediated endosome to Golgi retrograde transport by the J-domain protein RME-8

After endocytosis, most cargo enters the pleiomorphic early endosomes in which sorting occurs. As endosomes mature, transmembrane cargo can be sequestered into inwardly budding vesicles for degradation, or can exit the endosome in membrane tubules for recycling to the plasma membrane, the recycling endosome, or the Golgi apparatus. Endosome to Golgi transport requires the retromer complex. Without retromer, recycling cargo such as the MIG-14/Wntless protein aberrantly enters the degradative pathway and is depleted from the Golgi. Endosome-associated clathrin also affects the recycling of retrograde cargo and has been shown to function in the formation of endosomal subdomains. Here, we find that the Caemorhabditis elegans endosomal J-domain protein RME-8 associates with the retromer component SNX-1. Loss of SNX-1, RME-8, or the clathrin chaperone Hsc70/HSP-1 leads to over-accumulation of endosomal clathrin, reduced clathrin dynamics, and missorting of MIG-14 to the lysosome. Our results indicate a mechanism, whereby retromer can regulate endosomal clathrin dynamics through RME-8 and Hsc70, promoting the sorting of recycling cargo into the retrograde pathway.

The human protein coevolution network

Coevolution maintains interactions between phenotypic traits through the process of reciprocal natural selection. Detecting molecular coevolution can expose functional interactions between molecules in the cell, generating insights into biological processes, pathways, and the networks of interactions important for cellular function. Prediction of interaction partners from different protein families exploits the property that interacting proteins can follow similar patterns and relative rates of evolution. Current methods for detecting coevolution based on the similarity of phylogenetic trees or evolutionary distance matrices have, however, been limited by requiring coevolution over the entire evolutionary history considered and are inaccurate in the presence of paralogous copies. We present a novel method for determining coevolving protein partners by finding the largest common submatrix in a given pair of distance matrices, with the size of the largest common submatrix measuring the strength of coevolution. This approach permits us to consider matrices of different size and scale, to find lineage-specific coevolution, and to predict multiple interaction partners. We used MatrixMatchMaker to predict protein-protein interactions in the human genome. We show that proteins that are known to interact physically are more strongly coevolving than proteins that simply belong to the same biochemical pathway. The human coevolution network is highly connected, suggesting many more protein-protein interactions than are currently known from high-throughput and other experimental evidence. These most strongly coevolving proteins suggest interactions that have been maintained over long periods of evolutionary time, and that are thus likely to be of fundamental importance to cellular function.

Simultaneous inference of biological networks of multiple species from genome-wide data and evolutionary information: a semi-supervised approach

MOTIVATION

The existing supervised methods for biological network inference work on each of the networks individually based only on intra-species information such as gene expression data. We believe that it will be more effective to use genomic data and cross-species evolutionary information from different species simultaneously, rather than to use the genomic data alone.

RESULTS

We created a new semi-supervised learning method called Link Propagation for inferring biological networks of multiple species based on genome-wide data and evolutionary information. The new method was applied to simultaneous reconstruction of three metabolic networks of Caenorhabditis elegans, Helicobacter pylori and Saccharomyces cerevisiae, based on gene expression similarities and amino acid sequence similarities. The experimental results proved that the new simultaneous network inference method consistently improves the predictive performance over the individual network inferences, and it also outperforms in accuracy and speed other established methods such as the pairwise support vector machine.

AVAILABILITY

The software and data are available at http://cbio.ensmp.fr/~yyamanishi/LinkPropagation/.

Polo-like kinase 1 (PLK1) regulates interferon (IFN) induction by MAVS

The mitochondria-bound adapter MAVS participates in IFN induction by recruitment of downstream partners such as members of the TRAF family, leading to activation of NF-kappaB, and the IRF3 pathways. A yeast two-hybrid search for MAVS-interacting proteins yielded the Polo-box domain (PBD) of the mitotic Polo-like kinase PLK1. We showed that PBD associates with two different domains of MAVS in both dependent and independent phosphorylation events. The phosphodependent association requires the phosphopeptide binding ability of PBD. It takes place downstream of the proline-rich domain of MAVS, within an STP motif, characteristic of the binding of PLK1 to its targets, where the central Thr234 residue is phosphorylated. Its phosphoindependent association takes place at the C terminus of MAVS. PLK1 strongly inhibits the ability of MAVS to activate the IRF3 and NF-kappaB pathways and to induce IFN. Reciprocally, depletion of PLK1 can increase IFN induction in response to RIG-I/SeV or RIG-I/poly(I)-poly(C) treatments. This inhibition is dependent on the phosphoindependent association of PBD at the C terminus of MAVS where it disrupts the association of MAVS with its downstream partner TRAF3. IFN induction was strongly inhibited in cells arrested in G2/M by nocodazole, which provokes increased expression of endogenous PLK1. Interestingly, depletion of PLK1 from these nocodazole-treated cells could restore, at least partially, IFN induction. Altogether, these data demonstrate a new function for PLK1 as a regulator of IFN induction and provide the basis for the development of inhibitors preventing the PLK1/MAVS association to sustain innate immunity.

Local network topology in human protein interaction data predicts functional association

The use of high-throughput techniques to generate large volumes of protein-protein interaction (PPI) data has increased the need for methods that systematically and automatically suggest functional relationships among proteins. In a yeast PPI network, previous work has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional association. In this study we improved the prediction scheme by developing a new algorithm and applied it on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting function-associated protein pairs. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as benchmarks to compare and evaluate the function relevance. The application of our algorithms to human PPI data yielded 4,233 significant functional associations among 1,754 proteins. Further functional comparisons between them allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made functional inferences from detailed analysis on one subcluster highly enriched in the TGF-β signaling pathway (P<10⁻⁵⁰). Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotation in this post-genomic era.

A multiparameter network reveals extensive divergence between C. elegans bHLH transcription factors

Differences in expression, protein interactions, and DNA binding of paralogous transcription factors ("TF parameters") are thought to be important determinants of regulatory and biological specificity. However, both the extent of TF divergence and the relative contribution of individual TF parameters remain undetermined. We comprehensively identify dimerization partners, spatiotemporal expression patterns, and DNA-binding specificities for the C. elegans bHLH family of TFs, and model these data into an integrated network. This network displays both specificity and promiscuity, as some bHLH proteins, DNA sequences, and tissues are highly connected, whereas others are not. By comparing all bHLH TFs, we find extensive divergence and that all three parameters contribute equally to bHLH divergence. Our approach provides a framework for examining divergence for other protein families in C. elegans and in other complex multicellular organisms, including humans. Cross-species comparisons of integrated networks may provide further insights into molecular features underlying protein family evolution. For a video summary of this article, see the PaperFlick file available with the online Supplemental Data.

Marian Walhout: transcriptional mapmaker. Interviewed by Ben Short

Walhout uses the genome as a base camp for exploring transcriptional regulation.

Evolutionary flexibility of protein complexes

Background

Proteins play a key role in cellular life. They do not act alone but are organised in complexes. Throughout the life of a cell, complexes are dynamic in their composition due to attachments and shared components. Experimental and computational evidence indicate that consecutive addition and secondary losses of components played a major role in the evolution of some complexes, mostly without affecting the core function. Here, we analysed in a large scale approach whether this flexibility in evolution is only limited to a distinct number of complexes or represents a more general trend.

Results

Focussing on human protein complexes, we based our analysis on a manually curated dataset from HPRD. In total, 1,060 complexes with 6,136 proteins from 2,187 unique genes were considered. We computed interologs in 25 different species and predicted the composition of complexes. Over the analysed species, the composition of most complexes was highly flexible and only 25% of all genes were never lost. Even if one component was lost at a particular point in time, the fraction of observed second, independent losses of additional components was high (75% of all complexes affected). Still, loss of whole complexes happened rarely. This biological signal deviated significantly from random models. We exemplified this trend on the anaphase promoting complex (APC) where a core is highly conserved throughout all metazoans, but flexibility in certain components is observable.

Conclusion

Consecutive additions and losses of distinct units is a fundamental process in the evolution of protein complexes. These evolutionary events affecting genes coding for units in human protein complexes showed a significantly different phylogenetic pattern compared to randomly selected genes. Determination of taxon specific attachments or losses might be linked to specific cellular or morphological features. Thus, protein complexes contain not only structural and functional, but also evolutionary cores.

Predicting protein-protein interactions in Arabidopsis thaliana through integration of orthology, gene ontology and co-expression

Background

Large-scale identification of the interrelationships between different components of the cell, such as the interactions between proteins, has recently gained great interest. However, unraveling large-scale protein-protein interaction maps is laborious and expensive. Moreover, assessing the reliability of the interactions can be cumbersome.

Results

In this study, we have developed a computational method that exploits the existing knowledge on protein-protein interactions in diverse species through orthologous relations on the one hand, and functional association data on the other hand to predict and filter protein-protein interactions in Arabidopsis thaliana. A highly reliable set of protein-protein interactions is predicted through this integrative approach making use of existing protein-protein interaction data from yeast, human, C. elegans and D. melanogaster. Localization, biological process, and co-expression data are used as powerful indicators for protein-protein interactions. The functional repertoire of the identified interactome reveals interactions between proteins functioning in well-conserved as well as plant-specific biological processes. We observe that although common mechanisms (e.g. actin polymerization) and components (e.g. ARPs, actin-related proteins) exist between different lineages, they are active in specific processes such as growth, cancer metastasis and trichome development in yeast, human and Arabidopsis, respectively.

Conclusion

We conclude that the integration of orthology with functional association data is adequate to predict protein-protein interactions. Through this approach, a high number of novel protein-protein interactions with diverse biological roles is discovered. Overall, we have predicted a reliable set of protein-protein interactions suitable for further computational as well as experimental analyses.

Detecting hierarchical modularity in biological networks

Spatially or chemically isolated modules that carry out discrete functions are considered fundamental building blocks of cellular organization. However, detecting them in highly integrated biological networks requires a thorough understanding of the organization of these networks. In this chapter I argue that many biological networks are organized into many small, highly connected topologic modules that combine in a hierarchical manner into larger, less cohesive units. On top of a scale-free degree distribution, these networks show a power law scaling of the clustering coefficient with the node degree, a property that can be used as a signature of hierarchical organization. As a case study, I identify the hierarchical modules within the Escherichia coli metabolic network, and show that the uncovered hierarchical modularity closely overlaps with known metabolic functions.

Tissue specificity and the human protein interaction network

A protein interaction network describes a set of physical associations that can occur between proteins. However, within any particular cell or tissue only a subset of proteins is expressed and so only a subset of interactions can occur. Integrating interaction and expression data, we analyze here this interplay between protein expression and physical interactions in humans. Proteins only expressed in restricted cell types, like recently evolved proteins, make few physical interactions. Most tissue-specific proteins do, however, bind to universally expressed proteins, and so can function by recruiting or modifying core cellular processes. Conversely, most ‘housekeeping' proteins that are expressed in all cells also make highly tissue-specific protein interactions. These results suggest a model for the evolution of tissue-specific biology, and show that most, and possibly all, ‘housekeeping' proteins actually have important tissue-specific molecular interactions.

Uncovering new signaling proteins and potential drug targets through the interactome analysis of Mycobacterium tuberculosis

BACKGROUND

Analysis of the pathogen interactome is a powerful approach for dissecting potential signal transduction and virulence pathways. It also offers opportunities for exploring new drug targets.

RESULTS

In this study, a protein-protein interaction (PPI) network of Mycobacterium tuberculosis H37Rv was constructed using a homogenous protein mapping method, which has shown molecular chaperones, ribosomal proteins and ABC transporters to be highly interconnected proteins. A further analysis of this network unraveled the function of hypothetical proteins as well as a potential signaling pathway. A hypothetical protein, Rv2752c, which was linked to a metal cation-transporting ATPase, was characterized as a metal-beta-lactamase, through domain analysis in combination with an in vitro activity experiment. A second hypothetical protein, Rv1354c, and an unknown protein kinase, PknK, interacted with a similar group of inner membrane-associated ABC transporters in the PPI network. The interactions of Rv1354 with these proteins were also confirmed by a further bacterial two-hybrid analysis. According to protein domain structures, the unique M. tuberculosis Rv1354c gene was proposed, for the first time, to be responsible for the turnover of cyclic-di-GMP, a second messenger molecule in this bacterium. A further structure-based inhibitors screening for Rv1354c was also performed in silicon.

CONCLUSION

We constructed a comprehensive protein-protein interaction network for M. tuberculosis consisting of 738 proteins and 5639 interaction pairs. Our analysis unraveled the function of hypothetical proteins as well as a potential signaling pathway. The group of ABC transporters, PknK, and Rv1354c were proposed to constitute a potential membrane-associated signaling pathway that cooperatively responds to environmental stresses in M. tuberculosis. The study therefore provides valuable clues in exploring new signaling proteins, virulence pathways, and drug targets.

A human protein interaction network shows conservation of aging processes between human and invertebrate species

We have mapped a protein interaction network of human homologs of proteins that modify longevity in invertebrate species. This network is derived from a proteome-scale human protein interaction Core Network generated through unbiased high-throughput yeast two-hybrid searches. The longevity network is composed of 175 human homologs of proteins known to confer increased longevity through loss of function in yeast, nematode, or fly, and 2,163 additional human proteins that interact with these homologs. Overall, the network consists of 3,271 binary interactions among 2,338 unique proteins. A comparison of the average node degree of the human longevity homologs with random sets of proteins in the Core Network indicates that human homologs of longevity proteins are highly connected hubs with a mean node degree of 18.8 partners. Shortest path length analysis shows that proteins in this network are significantly more connected than would be expected by chance. To examine the relationship of this network to human aging phenotypes, we compared the genes encoding longevity network proteins to genes known to be changed transcriptionally during aging in human muscle. In the case of both the longevity protein homologs and their interactors, we observed enrichments for differentially expressed genes in the network. To determine whether homologs of human longevity interacting proteins can modulate life span in invertebrates, homologs of 18 human FRAP1 interacting proteins showing significant changes in human aging muscle were tested for effects on nematode life span using RNAi. Of 18 genes tested, 33% extended life span when knocked-down in Caenorhabditis elegans. These observations indicate that a broad class of longevity genes identified in invertebrate models of aging have relevance to human aging. They also indicate that the longevity protein interaction network presented here is enriched for novel conserved longevity proteins.

Studies of longevity in model organisms such as baker's yeast, roundworm, and fruit fly have clearly demonstrated that a diverse array of genetic mutations can result in increased life span. In fact, large-scale genetic screens have identified hundreds of genes that when mutated, knocked down, or deleted will significantly enhance longevity in these organisms. Despite great progress in understanding genetic and genomic determinants of life span in model organisms, the general relevance of invertebrate longevity genes to human aging and longevity has yet to be fully established. In this study, we show that human homologs of invertebrate longevity genes change in their expression levels during aging in human tissue. We also show that human genes encoding proteins that interact with human longevity homolog proteins are also changed in expression during human aging. These observations taken together indicate that the broad patterns underlying genetic control of life span in invertebrates is highly relevant to human aging and longevity. We also present a collection of novel candidate genes and proteins that may influence human life span.

Protein arrays as tools for serum autoantibody marker discovery in cancer

Protein array technology has begun to play a significant role in the study of protein-protein interactions and in the identification of antigenic targets of serum autoantibodies in a variety of autoimmune disorders. More recently, this technology has been applied to the identification of autoantibody signatures in cancer. The identification of tumour-associated antigens (TAAs) recognised by the patient's immune response represents an exciting approach to identify novel diagnostic cancer biomarkers and may contribute towards a better understanding of the molecular mechanisms involved. Circulating autoantibodies have not only been used to identify TAAs as diagnostic/prognostic markers and potential therapeutic targets, they also represent excellent biomarkers for the early detection of tumours and potential markers for monitoring the efficacy of treatment. Protein array technology offers the ability to screen the humoral immune response in cancer against thousands of proteins in a high throughput technique, thus readily identifying new panels of TAAs. Such an approach should not only aid in improved diagnostics, but has already contributed to the identification of complex autoantibody signatures that may represent disease subgroups, early diagnostics and facilitated the analysis of vaccine trials.

Protein-protein interaction as a predictor of subcellular location

Background

Many biological processes are mediated by dynamic interactions between and among proteins. In order to interact, two proteins must co-occur spatially and temporally. As protein-protein interactions (PPIs) and subcellular location (SCL) are discovered via separate empirical approaches, PPI and SCL annotations are independent and might complement each other in helping us to understand the role of individual proteins in cellular networks. We expect reliable PPI annotations to show that proteins interacting in vivo are co-located in the same cellular compartment. Our goal here is to evaluate the potential of using PPI annotation in determining SCL of proteins in human, mouse, fly and yeast, and to identify and quantify the factors that contribute to this complementarity.

Results

Using publicly available data, we evaluate the hypothesis that interacting proteins must be co-located within the same subcellular compartment. Based on a large, manually curated PPI dataset, we demonstrate that a substantial proportion of interacting proteins are in fact co-located. We develop an approach to predict the SCL of a protein based on the SCL of its interaction partners, given sufficient confidence in the interaction itself. The frequency of false positive PPIs can be reduced by use of six lines of supporting evidence, three based on type of recorded evidence (empirical approach, multiplicity of databases, and multiplicity of literature citations) and three based on type of biological evidence (inferred biological process, domain-domain interactions, and orthology relationships), with biological evidence more-effective than recorded evidence. Our approach performs better than four existing prediction methods in identifying the SCL of membrane proteins, and as well as or better for soluble proteins.

Conclusion

Understanding cellular systems requires knowledge of the SCL of interacting proteins. We show how PPI data can be used more effectively to yield reliable SCL predictions for both soluble and membrane proteins. Scope exists for further improvement in our understanding of cellular function through consideration of the biological context of molecular interactions.

Literature-curated protein interaction datasets

High quality datasets are needed to understand how global and local properties of protein-protein interaction, or “interactome”, networks relate to biological mechanisms, and to guide research on individual proteins. Evaluations of existing curation of protein interaction experiments reported in the literature find that curation can be error prone and possibly of lower quality than commonly assumed.

Empirically controlled mapping of the Caenorhabditis elegans protein-protein interactome network

To provide accurate biological hypotheses and elucidate global properties of cellular networks, systematic identification of protein-protein interactions must meet high quality standards.We present an expanded C. elegans protein-protein interaction network, or 'interactome' map, derived from testing a matrix of approximately 10,000 x approximately 10,000 proteins using a highly specific, high-throughput yeast two-hybrid system. Through a new empirical quality control framework, we show that the resulting data set (Worm Interactome 2007, or WI-2007) was similar in quality to low-throughput data curated from the literature. We filtered previous interaction data sets and integrated them with WI-2007 to generate a high-confidence consolidated map (Worm Interactome version 8, or WI8). This work allowed us to estimate the size of the worm interactome at approximately 116,000 interactions. Comparison with other types of functional genomic data shows the complementarity of distinct experimental approaches in predicting different functional relationships between genes or proteins

An integrated workflow for charting the human interaction proteome: insights into the PP2A system

Protein complexes represent major functional units for the execution of biological processes. Systematic affinity purification coupled with mass spectrometry (AP-MS) yielded a wealth of information on the compendium of protein complexes expressed in Saccharomyces cerevisiae. However, global AP-MS analysis of human protein complexes is hampered by the low throughput, sensitivity and data robustness of existing procedures, which limit its application for systems biology research. Here, we address these limitations by a novel integrated method, which we applied and benchmarked for the human protein phosphatase 2A system. We identified a total of 197 protein interactions with high reproducibility, showing the coexistence of distinct classes of phosphatase complexes that are linked to proteins implicated in mitosis, cell signalling, DNA damage control and more. These results show that the presented analytical process will substantially advance throughput and reproducibility in future systematic AP-MS studies on human protein complexes.

Improved homology-driven computational validation of protein-protein interactions motivated by the evolutionary gene duplication and divergence hypothesis

Background

Protein-protein interaction (PPI) data sets generated by high-throughput experiments are contaminated by large numbers of erroneous PPIs. Therefore, computational methods for PPI validation are necessary to improve the quality of such data sets. Against the background of the theory that most extant PPIs arose as a consequence of gene duplication, the sensitive search for homologous PPIs, i.e. for PPIs descending from a common ancestral PPI, should be a successful strategy for PPI validation.

Results

To validate an experimentally observed PPI, we combine FASTA and PSI-BLAST to perform a sensitive sequence-based search for pairs of interacting homologous proteins within a large, integrated PPI database. A novel scoring scheme that incorporates both quality and quantity of all observed matches allows us (1) to consider also tentative paralogs and orthologs in this analysis and (2) to combine search results from more than one homology detection method. ROC curves illustrate the high efficacy of this approach and its improvement over other homology-based validation methods.

Conclusion

New PPIs are primarily derived from preexisting PPIs and not invented de novo. Thus, the hallmark of true PPIs is the existence of homologous PPIs. The sensitive search for homologous PPIs within a large body of known PPIs is an efficient strategy to separate biologically relevant PPIs from the many spurious PPIs reported by high-throughput experiments.

Prediction and integration of regulatory and protein-protein interactions

Knowledge of transcriptional regulatory interactions (TRIs) is essential for exploring functional genomics and systems biology in any organism. While several results from genome-wide analysis of transcriptional regulatory networks are available, they are limited to model organisms such as yeast ( 1 ) and worm ( 2 ). Beyond these networks, experiments on TRIs study only individual genes and proteins of specific interest. In this chapter, we present a method for the integration of various data sets to predict TRIs for 54 organisms in the Bioverse ( 3 ). We describe how to compile and handle various formats and identifiers of data sets from different sources and how to predict TRIs using a homology-based approach, utilizing the compiled data sets. Integrated data sets include experimentally verified TRIs, binding sites of transcription factors, promoter sequences, protein subcellular localization, and protein families. Predicted TRIs expand the networks of gene regulation for a large number of organisms. The integration of experimentally verified and predicted TRIs with other known protein-protein interactions (PPIs) gives insight into specific pathways, network motifs, and the topological dynamics of an integrated network with gene expression under different conditions, essential for exploring functional genomics and systems biology.

High throughput interaction data reveals degree conservation of hub proteins

Research in model organisms relies on unspoken assumptions about the conservation of protein-protein interactions across species, yet several analyses suggest such conservation is limited. Fortunately, for many purposes the crucial issue is not global conservation of interactions, but preferential conservation of functionally important ones. An observed bias towards essentiality in highly-connected proteins implies the functional importance of such "hubs". We therefore define the notion of degree-conservation and demonstrate that hubs are preferentially degree-conserved. We show that a protein is more likely to be a hub if it has a high-degree ortholog, and that once a protein becomes a hub, it tends to remain so. We also identify a positive correlation between the average degree of a protein and the conservation of its interaction partners, and we find that the conservation of individual hub interactions is surprisingly high. Our work has important implications for prediction of protein function, computational inference of PPIs, and interpretation of data from model organisms.

Inferring protein-protein interactions from multiple protein domain combinations

The ever accumulating wealth of knowledge about protein interactions and the domain architecture of involved proteins in different organisms offers ways to understand the intricate interplay between interactome and proteome. Ultimately, the combination of these sources of information will allow the prediction of interactions among proteins where only domain composition is known. Based on the currently available protein-protein interaction and domain data of Saccharomyces cerevisiae and Drosophila melanogaster we introduce a novel method, Maximum Specificity Set Cover (MSSC), to predict potential protein-protein interactions. Utilizing interactions and domain architectures of domains as training sets, this algorithm employs a set cover approach to partition domain pairs, which allows the explanation of the underlying protein interaction to the largest degree of specificity. While MSSC in its basic version only considers domain pairs as the driving force between interactions, we also modified the algorithm to account for combinations of more than two domains that govern a protein-protein interaction. This approach allows us to predict the previously unknown protein-protein interactions in S. cerevisiae and D. melanogaster, with a degree of sensitivity and specificity that clearly outscores other approaches. As a proof of concept we also observe high levels of co-expression and decreasing GO distances between interacting proteins. Although our results are very encouraging, we observe that the quality of predictions significantly depends on the quality of interactions, which were utilized as the training set of the algorithm. The algorithm is part of a Web portal available at http://ppi.cse.nd.edu .

Interaction networks: from protein functions to drug discovery. A review

Most genes, proteins and other components carry out their functions within a complex network of interactions and a single molecule can affect a wide range of other cell components. A global, integrative, approach has been developed for several years, including protein-protein interaction networks (interactomes). In this review, we describe the high-throughput methods used to identify new interactions and to build large interaction datasets. The minimum information required for reporting a molecular interaction experiment (MIMIx) has been defined as a standard for storing data in publicly available interaction databases. Several examples of interaction networks from molecular machines (proteasome) or organelles (phagosome, mitochondrion) to whole organisms (viruses, bacteria, yeast, fly, and worm) are given and attempts to cover the entire human interaction network are discussed. The methods used to perform the topological analysis of interaction networks and to extract biological information from them are presented. These investigations have provided clues on protein functions, signalling and metabolic pathways, and physiological processes, unraveled the molecular basis of some diseases (cancer, infectious diseases), and will be very useful to identify new therapeutic targets and for drug discovery. A major challenge is now to integrate data from different sources (interactome, transcriptome, phenome, localization) to switch from static to dynamic interaction networks. The merging of a viral interactome and the human interactome has been used to simulate viral infection, paving the way for future studies aiming at providing molecular basis of human diseases.

Ortholog-based protein-protein interaction prediction and its application to inter-species interactions

BACKGROUND

The rapid growth of protein-protein interaction (PPI) data has led to the emergence of PPI network analysis. Despite advances in high-throughput techniques, the interactomes of several model organisms are still far from complete. Therefore, it is desirable to expand these interactomes with ortholog-based and other methods.

RESULTS

Orthologous pairs of 18 eukaryotic species were expanded and merged with experimental PPI datasets. The contributions of interologs from each species were evaluated. The expanded orthologous pairs enable the inference of interologs for various species. For example, more than 32,000 human interactions can be predicted. The same dataset has also been applied to the prediction of host-pathogen interactions. PPIs between P. falciparum calmodulin and several H. sapiens proteins are predicted, and these interactions may contribute to the maintenance of host cell Ca2+ concentration. Using comparisons with Bayesian and structure-based approaches, interactions between putative HSP40 homologs of P. falciparum and the H. sapiens TNF receptor associated factor family are revealed, suggesting a role for these interactions in the interference of the human immune response to P. falciparum.

CONCLUSION

The PPI datasets are available from POINT http://point.bioinformatics.tw/ and POINeT http://poinet.bioinformatics.tw/. Further development of methods to predict host-pathogen interactions should incorporate multiple approaches in order to improve sensitivity, and should facilitate the identification of targets for drug discovery and design.

Cell biology and evolution: molecular modules link it all?

Classical studies comparing developing embryos have suggested the importance of modified cell biological processes in the evolution of new phenotypes. Here, I revisit this connection focusing on embryonic development, in particular nematode embryogenesis. I compare phenotypic differences in nematode embryogenesis in two basic cell biological processes, the cell cycle and the localization of the first division axis. The analysis of these and other processes shows that, at the cell biological level, exhaustive variation is found that does not necessarily translate into morphological differences. Modern molecular analyses have led to a view in which molecular complexes, made up of groups of proteins, or modules, that are working together, are responsible for the proper execution of cell biological programs. I discuss how this modular architecture could facilitate the phenotypic changes observed in cell biological processes. Ultimately, understanding the connection between cellular behavior and phenotypic outcome will further elucidate the mechanisms responsible for phenotypic evolution.

Hybrids of the bHLH and bZIP protein motifs display different DNA-binding activities in vivo vs. in vitro

Minimalist hybrids comprising the DNA-binding domain of bHLH/PAS (basic-helix-loop-helix/Per-Arnt-Sim) protein Arnt fused to the leucine zipper (LZ) dimerization domain from bZIP (basic region-leucine zipper) protein C/EBP were designed to bind the E-box DNA site, CACGTG, targeted by bHLHZ (basic-helix-loop-helix-zipper) proteins Myc and Max, as well as the Arnt homodimer. The bHLHZ-like structure of ArntbHLH-C/EBP comprises the Arnt bHLH domain fused to the C/EBP LZ: i.e. swap of the 330 aa PAS domain for the 29 aa LZ. In the yeast one-hybrid assay (Y1H), transcriptional activation from the E-box was strong by ArntbHLH-C/EBP, and undetectable for the truncated ArntbHLH (PAS removed), as detected via readout from the HIS3 and lacZ reporters. In contrast, fluorescence anisotropy titrations showed affinities for the E-box with ArntbHLH-C/EBP and ArntbHLH comparable to other transcription factors (K(d) 148.9 nM and 40.2 nM, respectively), but only under select conditions that maintained folded protein. Although in vivo yeast results and in vitro spectroscopic studies for ArntbHLH-C/EBP targeting the E-box correlate well, the same does not hold for ArntbHLH. As circular dichroism confirms that ArntbHLH-C/EBP is a much more strongly alpha-helical structure than ArntbHLH, we conclude that the nonfunctional ArntbHLH in the Y1H must be due to misfolding, leading to the false negative that this protein is incapable of targeting the E-box. Many experiments, including protein design and selections from large libraries, depend on protein domains remaining well-behaved in the nonnative experimental environment, especially small motifs like the bHLH (60-70 aa). Interestingly, a short helical LZ can serve as a folding- and/or solubility-enhancing tag, an important device given the focus of current research on exploration of vast networks of biomolecular interactions.

A C. elegans genome-scale microRNA network contains composite feedback motifs with high flux capacity

MicroRNAs (miRNAs) and transcription factors (TFs) are primary metazoan gene regulators. Whereas much attention has focused on finding the targets of both miRNAs and TFs, the transcriptional networks that regulate miRNA expression remain largely unexplored. Here, we present the first genome-scale Caenorhabditis elegans miRNA regulatory network that contains experimentally mapped transcriptional TF --> miRNA interactions, as well as computationally predicted post-transcriptional miRNA --> TF interactions. We find that this integrated miRNA network contains 23 miRNA <--> TF composite feedback loops in which a TF that controls a miRNA is itself regulated by that same miRNA. By rigorous network randomizations, we show that such loops occur more frequently than expected by chance and, hence, constitute a genuine network motif. Interestingly, miRNAs and TFs in such loops are heavily regulated and regulate many targets. This "high flux capacity" suggests that loops provide a mechanism of high information flow for the coordinate and adaptable control of miRNA and TF target regulons.

A protein domain-based interactome network for C. elegans early embryogenesis

Many protein-protein interactions are mediated through independently folding modular domains. Proteome-wide efforts to model protein-protein interaction or "interactome" networks have largely ignored this modular organization of proteins. We developed an experimental strategy to efficiently identify interaction domains and generated a domain-based interactome network for proteins involved in C. elegans early-embryonic cell divisions. Minimal interacting regions were identified for over 200 proteins, providing important information on their domain organization. Furthermore, our approach increased the sensitivity of the two-hybrid system, resulting in a more complete interactome network. This interactome modeling strategy revealed insights into C. elegans centrosome function and is applicable to other biological processes in this and other organisms.

High-quality binary protein interaction map of the yeast interactome network

Current yeast interactome network maps contain several hundred molecular complexes with limited and somewhat controversial representation of direct binary interactions. We carried out a comparative quality assessment of current yeast interactome data sets, demonstrating that high-throughput yeast two-hybrid (Y2H) screening provides high-quality binary interaction information. Because a large fraction of the yeast binary interactome remains to be mapped, we developed an empirically controlled mapping framework to produce a "second-generation" high-quality, high-throughput Y2H data set covering approximately 20% of all yeast binary interactions. Both Y2H and affinity purification followed by mass spectrometry (AP/MS) data are of equally high quality but of a fundamentally different and complementary nature, resulting in networks with different topological and biological properties. Compared to co-complex interactome models, this binary map is enriched for transient signaling interactions and intercomplex connections with a highly significant clustering between essential proteins. Rather than correlating with essentiality, protein connectivity correlates with genetic pleiotropy.