Computational analysis of human protein interaction networks

Large amounts of human protein interaction data have been produced by experiments and prediction methods. However, the experimental coverage of the human interactome is still low in contrast to predicted data. To gain insight into the value of publicly available human protein network data, we compared predicted datasets, high-throughput results from yeast two-hybrid screens, and literature-curated protein-protein interactions. This evaluation is not only important for further methodological improvements, but also for increasing the confidence in functional hypotheses derived from predictions. Therefore, we assessed the quality and the potential bias of the different datasets using functional similarity based on the Gene Ontology, structural iPfam domain-domain interactions, likelihood ratios, and topological network parameters. This analysis revealed major differences between predicted datasets, but some of them also scored at least as high as the experimental ones regarding multiple quality measures. Therefore, since only small pair wise overlap between most datasets is observed, they may be combined to enlarge the available human interactome data. For this purpose, we additionally studied the influence of protein length on data quality and the number of disease proteins covered by each dataset. We could further demonstrate that protein interactions predicted by more than one method achieve an elevated reliability.     

A kaleidoscopic view of the Arabidopsis core cell cycle interactome

Although protein-protein interaction (PPI) networks have been shown to offer a systems-wide view of cellular processes, only a few plant PPI maps are available. Recently, the core cell cycle of Arabidopsis thaliana has been analyzed by three independent PPI technologies, including yeast two-hybrid systems, bimolecular fluorescence complementation and tandem affinity purification. Here, we merge the three interactomes with literature-curated and computationally predicted interactions, paving the way for a comprehensive picture of the plant core cell cycle machinery. Platform-specific interactions unveil the strengths and weaknesses of each detection method and give insights into the nature of the interactions among cell cycle proteins. Moreover, comparison of the obtained data reveals that a complete interactome can only be obtained when multiple techniques are applied in parallel.     

Charting plant interactomes: possibilities and challenges

Protein-protein interactions are essential for nearly all cellular processes. Therefore, an important goal of post-genomic research for defining gene function and understanding the function of macromolecular complexes involves creating 'interactome' maps from empirical or inferred datasets. Systematic efforts to conduct high-throughput surveys of protein-protein interactions in plants are needed to chart the complex and dynamic interaction networks that occur throughout plant development. However, no single approach can build a complete map of the interactome. Here, we review the utility and potential of various experimental approaches for creating large-scale protein-protein interaction maps in plants. Bioinformatics approaches for curating and assessing the confidence of these datasets through inter-species comparisons will be crucial in achieving a complete understanding of protein interaction networks in plants.     

Global approaches to protein-protein interactions

The availability of complete, annotated genome sequences for a variety of eukaryotic organisms has paved the way for a paradigm shift in biomedical research from the 'one gene-one hypothesis' approach to more global, systematic strategies that analyse genes or proteins on a genome- and proteome-wide scale. One daunting task in the post-genome era is to determine how the complement of expressed cellular proteins - the proteome - is organised into functional, higher-order networks, by mapping all constitutive and dynamic protein-protein interactions. Traditionally, reductionist approaches have typically focused on a few, selected gene products and their interactions in a particular physiological context. In contrast, more holistic strategies aim at understanding complex biological systems, for example global protein-protein interaction networks on a cellular or organismal level. Several large-scale proteomics technologies have been developed to generate comprehensive, cellular protein-protein interaction maps.     

A Comparative Approach to Characterize the Landscape of Host-Pathogen Protein-Protein Interactions

Significant efforts were gathered to generate large-scale comprehensive protein-protein interaction network maps. This is instrumental to understand the pathogen-host relationships and was essentially performed by genetic screenings in yeast two-hybrid systems. The recent improvement of protein-protein interaction detection by a Gaussia luciferase-based fragment complementation assay now offers the opportunity to develop integrative comparative interactomic approaches necessary to rigorously compare interaction profiles of proteins from different pathogen strain variants against a common set of cellular factors.This paper specifically focuses on the utility of combining two orthogonal methods to generate protein-protein interaction datasets: yeast two-hybrid (Y2H) and a new assay, high-throughput Gaussia princeps protein complementation assay (HT-GPCA) performed in mammalian cells.A large-scale identification of cellular partners of a pathogen protein is performed by mating-based yeast two-hybrid screenings of cDNA libraries using multiple pathogen strain variants. A subset of interacting partners selected on a high-confidence statistical scoring is further validated in mammalian cells for pair-wise interactions with the whole set of pathogen variants proteins using HT-GPCA. This combination of two complementary methods improves the robustness of the interaction dataset, and allows the performance of a stringent comparative interaction analysis. Such comparative interactomics constitute a reliable and powerful strategy to decipher any pathogen-host interplays.     

Protein docking by Rotation-Based Uniform Sampling (RotBUS) with fast computing of intermolecular contact distance and residue desolvation

Background  

Protein-protein interactions are fundamental for the majority of cellular processes and their study is of enormous biotechnological and therapeutic interest. In recent years, a variety of computational approaches to the protein-protein docking problem have been reported, with encouraging results. Most of the currently available protein-protein docking algorithms are composed of two clearly defined parts: the sampling of the rotational and translational space of the interacting molecules, and the scoring and clustering of the resulting orientations. Although this kind of strategy has shown some of the most successful results in the CAPRI blind test , more efforts need to be applied. Thus, the sampling protocol should generate a pool of conformations that include a sufficient number of near-native ones, while the scoring function should discriminate between near-native and non-near-native proposed conformations. On the other hand, protocols to efficiently include full flexibility on the protein structures are increasingly needed.     

A Novel Phosphopeptide Microarray Based Interactome Map in Breast Cancer Cells Reveals Phosphoprotein-GRB2 Cell Signaling Networks

The architecture of cellular proteins connected to form signaling pathways in response to internal and external cues is much more complex than a group of simple protein-protein interactions. Post translational modifications on proteins (e.g., phosphorylation of serine, threonine and tyrosine residues on proteins) initiate many downstream signaling events leading to protein-protein interactions and subsequent activation of signaling cascades leading to cell proliferation, cell differentiation and cell death. As evidenced by a rapidly expanding mass spectrometry database demonstrating protein phosphorylation at specific motifs, there is currently a large gap in understanding the functional significance of phosphoproteins with respect to their specific protein connections in the signaling cascades. A comprehensive map that interconnects phospho-motifs in pathways will enable identification of nodal protein interactions that are sensitive signatures indicating a disease phenotype from the physiological hemostasis and provide clues into control of disease. Using a novel phosphopeptide microarray technology, we have mapped endogenous tyrosine-phosphoproteome interaction networks in breast cancer cells mediated by signaling adaptor protein GRB2, which transduces cellular responses downstream of several RTKs through the Ras-ERK signaling cascade. We have identified several previously reported motif specific interactions and novel interactions. The peptide microarray data indicate that various phospho-motifs on a single protein are differentially regulated in various cell types and shows global downregulation of phosphoprotein interactions specifically in cells with metastatic potential. The study has revealed novel phosphoprotein mediated signaling networks, which warrants further detailed analysis of the nodes of protein-protein interaction to uncover their biomarker or therapeutic potential.     

Quantitative assessment of the structural bias in protein-protein interaction assays

With recent publications of several large-scale protein-protein interaction (PPI) studies, the realization of the full yeast interaction network is getting closer. Here, we have analysed several yeast protein interaction datasets to understand their strengths and weaknesses. In particular, we investigate the effect of experimental biases on some of the protein properties suggested to be enriched in highly connected proteins. Finally, we use support vector machines (SVM) to assess the contribution of these properties to protein interactivity. We find that protein abundance is the most important factor for detecting interactions in tandem affinity purifications (TAP), while it is of less importance for Yeast Two Hybrid (Y2H) screens. Consequently, sequence conservation and/or essentiality of hubs may be related to their high abundance. Further, proteins with disordered structure are over-represented in Y2H screens and in one, but not the other, large-scale TAP assay. Hence, disordered regions may be important both in transient interactions and interactions in complexes. Finally, a few domain families seem to be responsible for a large part of all interactions. Most importantly, we show that there are method-specific biases in PPI experiments. Thus, care should be taken before drawing strong conclusions based on a single dataset.     

SAINT-MS1: protein-protein interaction scoring using label-free intensity data in affinity purification-mass spectrometry experiments

We present a statistical method SAINT-MS1 for scoring protein-protein interactions based on the label-free MS1 intensity data from affinity purification-mass spectrometry (AP-MS) experiments. The method is an extension of Significance Analysis of INTeractome (SAINT), a model-based method previously developed for spectral count data. We reformulated the statistical model for log-transformed intensity data, including adequate treatment of missing observations, that is, interactions identified in some but not all replicate purifications. We demonstrate the performance of SAINT-MS1 using two recently published data sets: a small LTQ-Orbitrap data set with three replicate purifications of single human bait protein and control purifications and a larger drosophila data set targeting insulin receptor/target of rapamycin signaling pathway generated using an LTQ-FT instrument. Using the drosophila data set, we also compare and discuss the performance of SAINT analysis based on spectral count and MS1 intensity data in terms of the recovery of orthologous and literature-curated interactions. Given rapid advances in high mass accuracy instrumentation and intensity-based label-free quantification software, we expect that SAINT-MS1 will become a useful tool allowing improved detection of protein interactions in label-free AP-MS data, especially in the low abundance range.     

Larval morphology and host use confirms ecotypic variation in Cactoblastis cactorum (Berg)

Despite their recognized importance in the literature, the contribution of native-range species interactions to invasion success has been inadequately studied. Previous authors have suggested that biases in the sampling of propagules from the native range might influence invasion success, but most contemporary invasion hypotheses focus on the development of novel interactions or a release from native consumers and competitors. When ecotypic variation exists in native host-consumer associations, the specific pattern of sampling across ecotypes could determine invasion success, especially when the genetic diversity among exotic propagules is low. The South American cactus moth, Cactoblastis cactorum (Berg), is an oligophagous consumer whose larvae feed on prickly pear cacti (subfamily Opuntioideae). The moth was collected from a small geographic area along the Argentina-Uruguay border in 1925 and was introduced to multiple continents as a biological control species, which has subsequently invaded North America. Here we show that groups defined by genetic structure in this species’ native range are concordant with distinct patterns of host association and larval morphology. Furthermore, in Florida populations, morphological traits have diverged from those found in the native range, and patterns of host association suggest that strong biases in host preference also occur in invasive populations. The documented history of C. cactorum introductions confirms that multiple attempts were made to export the moth, but that only a single ecotype was exported successfully. Additional work will be necessary to determine whether the observed host biases in North America reflect a rapid adaptation to naïve hosts or a conservation of traits related to specific aspects of the host-consumer association.     

New insights into protein-protein interaction data lead to increased estimates of the <Emphasis Type="Italic">S. cerevisiae</Emphasis> interactome size

Background  

As protein interactions mediate most cellular mechanisms, protein-protein interaction networks are essential in the study of cellular processes. Consequently, several large-scale interactome mapping projects have been undertaken, and protein-protein interactions are being distilled into databases through literature curation; yet protein-protein interaction data are still far from comprehensive, even in the model organism Saccharomyces cerevisiae. Estimating the interactome size is important for evaluating the completeness of current datasets, in order to measure the remaining efforts that are required.     

Three-dimensional reconstruction of protein networks provides insight into human genetic disease

To better understand the molecular mechanisms and genetic basis of human disease, we systematically examine relationships between 3,949 genes, 62,663 mutations and 3,453 associated disorders by generating a three-dimensional, structurally resolved human interactome. This network consists of 4,222 high-quality binary protein-protein interactions with their atomic-resolution interfaces. We find that in-frame mutations (missense point mutations and in-frame insertions and deletions) are enriched on the interaction interfaces of proteins associated with the corresponding disorders, and that the disease specificity for different mutations of the same gene can be explained by their location within an interface. We also predict 292 candidate genes for 694 unknown disease-to-gene associations with proposed molecular mechanism hypotheses. This work indicates that knowledge of how in-frame disease mutations alter specific interactions is critical to understanding pathogenesis. Structurally resolved interaction networks should be valuable tools for interpreting the wealth of data being generated by large-scale structural genomics and disease association studies.     

Discovering disease-genes by topological features in human protein-protein interaction network

MOTIVATION: Mining the hereditary disease-genes from human genome is one of the most important tasks in bioinformatics research. A variety of sequence features and functional similarities between known human hereditary disease-genes and those not known to be involved in disease have been systematically examined and efficient classifiers have been constructed based on the identified common patterns. The availability of human genome-wide protein-protein interactions (PPIs) provides us with new opportunity for discovering hereditary disease-genes by topological features in PPIs network. RESULTS: This analysis reveals that the hereditary disease-genes ascertained from OMIM in the literature-curated (LC) PPIs network are characterized by a larger degree, tendency to interact with other disease-genes, more common neighbors and quick communication to each other whereas those properties could not be detected from the network identified from high-throughput yeast two-hybrid mapping approach (EXP) and predicted interactions (PDT) PPIs network. KNN classifier based on those features was created and on average gained overall prediction accuracy of 0.76 in cross-validation test. Then the classifier was applied to 5262 genes on human genome and predicted 178 novel disease-genes. Some of the predictions have been validated by biological experiments.     

Evolutionary and physiological importance of hub proteins

It has been claimed that proteins with more interaction partners (hubs) are both physiologically more important (i.e., less dispensable) and, owing to an assumed high density of binding sites, slow evolving. Not all analyses, however, support these results, probably because of biased and less-than reliable global protein interaction data. Here we provide the first examination of these issues using a comprehensive literature-curated dataset of well-substantiated protein interactions in Saccharomyces cerevisiae. Whereas use of less reliable yeast two-hybrid data alone can reject the possibility that local connectivity correlates with measures of dispensability, in higher quality datasets a relatively robust correlation is observed. In contrast, local connectivity does not correlate with the rate of protein evolution even in reliable datasets. This perhaps surprising lack of correlation with evolutionary rate appears in part to arise from the fact that hub proteins do not have a higher density of residues associated with binding. However, hub proteins do have at least one other set of unusual features, namely rapid turnover and regulation, as manifest in high mRNA decay rates and a large number of phosphorylation sites. This, we suggest, is an adaptation to minimize unwanted activation of pathways that might be mediated by adventitious binding to hubs, were they to actively persist longer than required at any given time point. We conclude that hub proteins are more important for cellular growth rate and under tight regulation but are not slow evolving.     

Challenges and dreams: physics of weak interactions essential to life

Biological systems display stunning capacities to self-organize. Moreover, their subcellular architectures are dynamic and responsive to changing needs and conditions. Key to these properties are manifold weak “quinary” interactions that have evolved to create specific spatial networks of macromolecules. These specific arrangements of molecules enable signals to be propagated over distances much greater than molecular dimensions, create phase separations that define functional regions in cells, and amplify cellular responses to changes in their environments. A major challenge is to develop biochemical tools and physical models to describe the panoply of weak interactions operating in cells. We also need better approaches to measure the biases in the spatial distributions of cellular macromolecules that result from the integrated action of multiple weak interactions. Partnerships between cell biologists, biochemists, and physicists are required to deploy these methods. Together these approaches will help us realize the dream of understanding the biological “glue” that sustains life at a molecular and cellular level.     

Normalization of a chromosomal contact map

ABSTRACT: BACKGROUND: Chromatin organization has been increasingly studied in relation with its important influence on DNA-related metabolic processes such as replication or regulation of gene expression. Since its original design ten years ago, capture of chromosome conformation (3C) has become an essential tool to investigate the overall conformation of chromosomes. It relies on the capture of long-range trans and cis interactions of chromosomal segments whose relative proportions in the final bank reflect their frequencies of interactions, hence their average spatial proximity. The recent coupling of 3C with deep sequencing approaches now allows the generation of high resolution genome-wide chromosomal contact maps. Different protocols have been used to generate such maps in various organisms. This includes mammals, drosophila and yeast. The massive amount of raw data generated by the genomic 3C has to be carefully processed to alleviate the various biases and byproducts generated by the experiments. Our study aims at proposing a simple normalization procedure to take into account these unwanted but inevitable events. RESULTS: Careful analysis of the raw data generated previously for budding yeast Saccharomyces cerevisiae led to the identification of three main biases affecting the final datasets, including an original bias resulting from the circularization of DNA molecules exhibiting specific lengths in accordance with laws from polymer physics. We then developed a simple normalization procedure to process the data and allow the generation of a normalized, highly contrasted, chromosomal contact map for S. cerevisiae. The same method was then extended to the first human genome contact map. Using the normalized data, we revisited the preferential interactions originally described between subsets of discrete chromosomal features. Notably, the detection of preferential interactions between tRNA in yeast and CTCF, PolII binding sites in human can vary with the normalization procedure used. CONCLUSIONS: We quantitatively reanalyzed the genomic 3C data obtained for S. cerevisiae, identified some of the biases inherent to the technique and proposed a simple normalization procedure to analyze them. Such an approach can be easily generalized for genomic 3C experiments in other organisms. More experiments and analysis will be necessary to reach optimal resolution and accuracies of the maps generated through these approaches. Working with cell population presenting highest levels of homogeneity will prove useful in this regards.