MOTIVATION: Differential kinase interaction networks (DKINs) are networks containing kinase-substrate links that are differentially active between two conditions. Existing methods are either able to predict condition-agnostic kinase-substrate links or condition-specific differential kinase activity, but do not provide differential kinase-substrate links. Moreover, existing methods for predicting kinase-substrate links usually rely on curated biochemical knowledge. Thus, there is a lack of data-driven DKIN inference methods that are also applicable when prior knowledge is scarce. RESULTS: To address this need, we present KINference. KINference combines computation of a baseline KIN representing the space of all possible kinase-substrate links with filters applied to nodes and edges to identify differentially active subnetworks that are relevant in the context of a specific phosphoproteomics dataset. For the node filters, we rely on functional relevance and differential phosphorylation scores; for the edge filters, we make use of prize-collecting Steiner trees and correlations between phosphorylation sites of kinases and their target proteins. Tests on two phosphoproteomics datasets (kinase inhibition in breast cancer cells, SARS-CoV-2 infection in Calu-3 cells) show that the proposed filters produce significant results in terms of overlap with known interactions between kinases and phosphorylation sites. Furthermore, a case study on the SARS-CoV-2 infection data, suggests a potential host pathway linked to virus replication, showcasing the process of hypothesis generation utilizing DKINs computed by KINference. AVAILABILITY AND IMPLEMENTATION: KINference is available as an R package at https://github.com/bionetslab/KINference and https://doi.org/10.5281/zenodo.15411150. Scripts to reproduce the results are available at https://github.com/bionetslab/KINference-Evaluation-Scripts and https://doi.org/10.5281/zenodo.15424599.

Growing evidence supports other regulatory roles for protein ubiquitination in addition to serving as a tag for proteasomal degradation. In contrast to other common post-translational modifications, such as phosphorylation, little is known about how non-degradative ubiquitination modulates protein structure, dynamics, and function. Due to the wealth of knowledge concerning protein kinase structure and regulation, we examined kinase ubiquitination using ubiquitin remnant immunoaffinity enrichment and quantitative mass spectrometry to identify ubiquitinated kinases and the sites of ubiquitination in Jurkat and HEK293 cells. We find that, unlike phosphorylation, ubiquitination most commonly occurs in structured domains, and on the kinase domain, ubiquitination is concentrated in regions known to be important for regulating activity. We hypothesized that ubiquitination, like other post-translational modifications, may alter the conformational equilibrium of the modified protein. We chose one human kinase, ZAP-70, to simulate using molecular dynamics with and without a monoubiquitin modification. In Jurkat cells, ZAP-70 is ubiquitinated at several sites that are not sensitive to proteasome inhibition and thus may have other regulatory roles. Our simulations show that ubiquitination influences the conformational ensemble of ZAP-70 in a site-dependent manner. When monoubiquitinated at K377, near the C-helix, the active conformation of the ZAP-70 C-helix is disrupted. In contrast, when monoubiquitinated at K476, near the kinase hinge region, an active-like ZAP-70 C-helix conformation is stabilized. These results lead to testable hypotheses that ubiquitination directly modulates kinase activity, and that ubiquitination is likely to alter structure, dynamics, and function in other protein classes as well.

The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.

BACKGROUND: The increasing availability of HIV-host interaction datasets, including both physical and genetic interactions, has created a need for software tools to integrate and visualize the data. Because these host-pathogen interactions are extensive and interactions between human proteins are found within many different databases, it is difficult to generate integrated HIV-human interaction networks. RESULTS: We have developed a web-based platform, termed GPS-Prot http://www.gpsprot.org, that allows for facile integration of different HIV interaction data types as well as inclusion of interactions between human proteins derived from publicly-available databases, including MINT, BioGRID and HPRD. The software has the ability to group proteins into functional modules or protein complexes, generating more intuitive network representations and also allows for the uploading of user-generated data. CONCLUSIONS: GPS-Prot is a software tool that allows users to easily create comprehensive and integrated HIV-host networks. A major advantage of this platform compared to other visualization tools is its web-based format, which requires no software installation or data downloads. GPS-Prot allows novice users to quickly generate networks that combine both genetic and protein-protein interactions between HIV and its human host into a single representation. Ultimately, the platform is extendable to other host-pathogen systems.

MOTIVATION: Pairwise experimental perturbation is increasingly used to probe gene and protein function because these studies offer powerful insight into the activity and regulation of biological systems. Symmetric two-dimensional datasets, such as pairwise genetic interactions are amenable to an optimally designed measurement procedure because of the equivalence of cases and conditions where fewer experimental measurements may be required to extract the underlying structure. RESULTS: We show that optimal experimental design can provide improvements in efficiency when collecting data in an iterative manner. We develop a method built on a statistical clustering model for symmetric data and the Fisher information uncertainty estimates, and we also provide simple heuristic approaches that have comparable performance. Using yeast epistatic miniarrays as an example, we show that correct assignment of the major subnetworks could be achieved with <50% of the measurements in the complete dataset. Optimization is likely to become critical as pairwise functional studies extend to more complex mammalian systems where all by all experiments are currently intractable.

MOTIVATION: When working with large-scale protein interaction data, an important analysis task is the assignment of pairs of proteins to groups that correspond to higher order assemblies. Previously a common approach to this problem has been to apply standard hierarchical clustering methods to identify such a groups. Here we propose a new algorithm for aggregating a diverse collection of matrix factorizations to produce a more informative clustering, which takes the form of a 'soft' hierarchy of clusters. RESULTS: We apply the proposed Ensemble non-negative matrix factorization (NMF) algorithm to a high-quality assembly of binary protein interactions derived from two proteome-wide studies in yeast. Our experimental evaluation demonstrates that the algorithm lends itself to discovering small localized structures in this data, which correspond to known functional groupings of complexes. In addition, we show that the algorithm also supports the assignment of putative functions for previously uncharacterized proteins, for instance the protein YNR024W, which may be an uncharacterized component of the exosome.

Recently, a number of advanced screening technologies have allowed for the comprehensive quantification of aggravating and alleviating genetic interactions among gene pairs. In parallel, TAP-MS studies (tandem affinity purification followed by mass spectroscopy) have been successful at identifying physical protein interactions that can indicate proteins participating in the same molecular complex. Here, we propose a method for the joint learning of protein complexes and their functional relationships by integration of quantitative genetic interactions and TAP-MS data. Using 3 independent benchmark datasets, we demonstrate that this method is >50% more accurate at identifying functionally related protein pairs than previous approaches. Application to genes involved in yeast chromosome organization identifies a functional map of 91 multimeric complexes, a number of which are novel or have been substantially expanded by addition of new subunits. Interestingly, we find that complexes that are enriched for aggravating genetic interactions (i.e., synthetic lethality) are more likely to contain essential genes, linking each of these interactions to an underlying mechanism. These results demonstrate the importance of both large-scale genetic and physical interaction data in mapping pathway architecture and function.

BACKGROUND: Identification of protein interaction networks has received considerable attention in the post-genomic era. The currently available biochemical approaches used to detect protein-protein interactions are all time and labour intensive. Consequently there is a growing need for the development of computational tools that are capable of effectively identifying such interactions. RESULTS: Here we explain the development and implementation of a novel Protein-Protein Interaction Prediction Engine termed PIPE. This tool is capable of predicting protein-protein interactions for any target pair of the yeast Saccharomyces cerevisiae proteins from their primary structure and without the need for any additional information or predictions about the proteins. PIPE showed a sensitivity of 61% for detecting any yeast protein interaction with 89% specificity and an overall accuracy of 75%. This rate of success is comparable to those associated with the most commonly used biochemical techniques. Using PIPE, we identified a novel interaction between YGL227W (vid30) and YMR135C (gid8) yeast proteins. This lead us to the identification of a novel yeast complex that here we term vid30 complex (vid30c). The observed interaction was confirmed by tandem affinity purification (TAP tag), verifying the ability of PIPE to predict novel protein-protein interactions. We then used PIPE analysis to investigate the internal architecture of vid30c. It appeared from PIPE analysis that vid30c may consist of a core and a secondary component. Generation of yeast gene deletion strains combined with TAP tagging analysis indicated that the deletion of a member of the core component interfered with the formation of vid30c, however, deletion of a member of the secondary component had little effect (if any) on the formation of vid30c. Also, PIPE can be used to analyse yeast proteins for which TAP tagging fails, thereby allowing us to predict protein interactions that are not included in genome-wide yeast TAP tagging projects. CONCLUSION: PIPE analysis can predict yeast protein-protein interactions. Also, PIPE analysis can be used to study the internal architecture of yeast protein complexes. The data also suggests that a finite set of short polypeptide signals seem to be responsible for the majority of the yeast protein-protein interactions.