Novel peptide-mediated interactions derived from high-resolution 3-dimensional structures

A dual-population multi-objective evolutionary algorithm driven by generative adversarial networks for benchmarking and protein-peptide docking

Multi-objective optimization problems (MOPs) are characterized as optimization problems in which multiple conflicting objective functions are optimized simultaneously. To solve MOPs, some algorithms used machine learning models to drive the evolutionary algorithms, leading to the design of a variety of model-based evolutionary algorithms. However, model collapse occurs during the generation of candidate solutions, which results in local optima and poor diversity in model-based evolutionary algorithms. To address this problem, we propose a dual-population multi-objective evolutionary algorithm driven by Wasserstein generative adversarial network with gradient penalty (DGMOEA), where the dual-populations coordinate and cooperate to generate high-quality solutions, thus improving the performance of the evolutionary algorithm. We compare the proposed algorithm with the 7 state-of-the-art algorithms on 20 multi-objective benchmark functions. Experimental results indicate that DGMOEA achieves significant results in solving MOPs, where the metrics IGD and HV outperform the other comparative algorithms on 15 and 18 out of 20 benchmarks, respectively. Our algorithm is evaluated on the LEADS-PEP dataset containing 53 protein-peptide complexes, and the experimental results on solving the protein-peptide docking problem indicated that DGMOEA can effectively reduce the RMSD between the generated and the original peptide's 3D poses and achieve more competitive results.

Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-Liquid/Liquid Phase Separation, Oligomerization, GO Analysis: A Toolkit for Interaction Data-Based Domain Analysis

Although a large number of databases are available for regulatory elements, a bottleneck has been created by the lack of bioinformatics tools to predict the interaction modes of regulatory elements. To reduce this gap, we developed the Arabidopsis Transcription Regulatory Factor Domain/Domain Interaction Analysis Tool-liquid/liquid phase separation (LLPS), oligomerization, GO analysis (ART FOUNDATION-LOG), a useful toolkit for protein-nucleic acid interaction (PNI) and protein-protein interaction (PPI) analysis based on domain-domain interactions (DDIs). LLPS, protein oligomerization, the structural properties of protein domains, and protein modifications are major components in the orchestration of the spatiotemporal dynamics of PPIs and PNIs. Our goal is to integrate PPI/PNI information into the development of a prediction model for identifying important genetic variants in peaches. Our program unified interdatabase relational keys based on protein domains to facilitate inference from the model species. A key advantage of this program lies in the integrated information of related features, such as protein oligomerization, LOG analysis, structural characterizations of domains (e.g., domain linkers, intrinsically disordered regions, DDIs, domain-motif (peptide) interactions, beta sheets, and transmembrane helices), and post-translational modification. We provided simple tests to demonstrate how to use this program, which can be applied to other eukaryotic organisms.

Targeting Protein-Protein Interfaces with Peptides: The Contribution of Chemical Combinatorial Peptide Library Approaches

Protein-protein interfaces play fundamental roles in the molecular mechanisms underlying pathophysiological pathways and are important targets for the design of compounds of therapeutic interest. However, the identification of binding sites on protein surfaces and the development of modulators of protein-protein interactions still represent a major challenge due to their highly dynamic and extensive interfacial areas. Over the years, multiple strategies including structural, computational, and combinatorial approaches have been developed to characterize PPI and to date, several successful examples of small molecules, antibodies, peptides, and aptamers able to modulate these interfaces have been determined. Notably, peptides are a particularly useful tool for inhibiting PPIs due to their exquisite potency, specificity, and selectivity. Here, after an overview of PPIs and of the commonly used approaches to identify and characterize them, we describe and evaluate the impact of chemical peptide libraries in medicinal chemistry with a special focus on the results achieved through recent applications of this methodology. Finally, we also discuss the role that this methodology can have in the framework of the opportunities, and challenges that the application of new predictive approaches based on artificial intelligence is generating in structural biology.

Targeting peptide-mediated interactions in omics

Peptide-mediated interactions (PMIs) play a crucial role in cell signaling network, which are responsible for about half of cellular protein-protein associations in the human interactome and have recently been recognized as a new kind of promising druggable target for drug development and disease therapy. In this article, we give a systematic review regarding the proteome-wide discovery of PMIs and targeting druggable PMIs (dPMIs) with chemical drugs, self-inhibitory peptides (SIPs) and protein agents, particularly focusing on their implications and applications for therapeutic purpose in omics. We also introduce computational peptidology strategies used to model, analyze, and design PMI-targeted molecular entities and further extend the concepts of protein context, direct/indirect readout, and enthalpy/entropy effect involved in PMIs. Current issues and future perspective on this topic are discussed. There is still a long way to go before establishment of efficient therapeutic strategies to target PMIs on the omics scale.

Protein Receptors Evolved from Homologous Cohesion Modules That Self-Associated and Are Encoded by Interactive Networked Genes

Previously, it was proposed that protein receptors evolved from self-binding peptides that were encoded by self-interacting gene segments (inverted repeats) widely dispersed in the genome. In addition, self-association of the peptides was thought to be mediated by regions of amino acid sequence similarity. To extend these ideas, special features of receptors have been explored, such as their degree of homology to other proteins, and the arrangement of their genes for clues about their evolutionary origins and dynamics in the genome. As predicted, BLASTP searches for homologous proteins detected a greater number of unique hits for queries with receptor sequences than for sequences of randomly-selected, non-receptor proteins. This suggested that the building blocks (cohesion modules) for receptors were duplicated, dispersed, and maintained in the genome, due to structure/function relationships discussed here. Furthermore, the genes coding for a representative panel of receptors participated in a larger number of gene-gene interactions than for randomly-selected genes. This could conceivably reflect a greater evolutionary conservation of the receptor genes, with their more extensive integration into networks along with inherent properties of the genes themselves. In support of the latter possibility, some receptor genes were located in active areas of adaptive gene relocation/amalgamation to form functional blocks of related genes. It is suggested that adaptive relocation might allow for their joint regulation by common promoters and enhancers, and affect local chromatin structural domains to facilitate or repress gene expression. Speculation is included about the nature of the coordinated communication between receptors and the genes that encode them.

Large-scale survey and database of high affinity ligands for peptide recognition modules

Many proteins involved in signal transduction contain peptide recognition modules (PRMs) that recognize short linear motifs (SLiMs) within their interaction partners. Here, we used large-scale peptide-phage display methods to derive optimal ligands for 163 unique PRMs representing 79 distinct structural families. We combined the new data with previous data that we collected for the large SH3, PDZ, and WW domain families to assemble a database containing 7,984 unique peptide ligands for 500 PRMs representing 82 structural families. For 74 PRMs, we acquired enough new data to map the specificity profiles in detail and derived position weight matrices and binding specificity logos based on multiple peptide ligands. These analyses showed that optimal peptide ligands resembled peptides observed in existing structures of PRM-ligand complexes, indicating that a large majority of the phage-derived peptides are likely to target natural peptide-binding sites and could thus act as inhibitors of natural protein-protein interactions. The complete dataset has been assembled in an online database (http://www.prm-db.org) that will enable many structural, functional, and biological studies of PRMs and SLiMs.

Exploring Proteomic Drug Targets, Therapeutic Strategies and Protein - Protein Interactions in Cancer: Mechanistic View

Protein-Protein Interactions (PPIs) drive major signalling cascades and play critical role in cell proliferation, apoptosis, angiogenesis and trafficking. Deregulated PPIs are implicated in multiple malignancies and represent the critical targets for treating cancer. Herein, we discuss the key protein-protein interacting domains implicated in cancer notably PDZ, SH2, SH3, LIM, PTB, SAM and PH. These domains are present in numerous enzymes/kinases, growth factors, transcription factors, adaptor proteins, receptors and scaffolding proteins and thus represent essential sites for targeting cancer. This review explores the candidature of various proteins involved in cellular trafficking (small GTPases, molecular motors, matrix-degrading enzymes, integrin), transcription (p53, cMyc), signalling (membrane receptor proteins), angiogenesis (VEGFs) and apoptosis (BCL-2family), which could possibly serve as targets for developing effective anti-cancer regimen. Interactions between Ras/Raf; X-linked inhibitor of apoptosis protein (XIAP)/second mitochondria-derived activator of caspases (Smac/DIABLO); Frizzled (FRZ)/Dishevelled (DVL) protein; beta-catenin/T Cell Factor (TCF) have also been studied as prospective anticancer targets. Efficacy of diverse molecules/ drugs targeting such PPIs although evaluated in various animal models/cell lines, there is an essential need for human-based clinical trials. Therapeutic strategies like the use of biologicals, high throughput screening (HTS) and fragment-based technology could play an imperative role in designing cancer therapeutics. Moreover, bioinformatic/computational strategies based on genome sequence, protein sequence/structure and domain data could serve as competent tools for predicting PPIs. Exploring hot spots in proteomic networks represents another approach for developing targetspecific therapeutics. Overall, this review lays emphasis on a productive amalgamation of proteomics, genomics, biochemistry, and molecular dynamics for successful treatment of cancer.

Is protein context responsible for peptide-mediated interactions?

Many cell signaling pathways are orchestrated by the weak, transient, and reversible peptide-mediated interactions (PMIs). Here, the role of protein context in contributing to the stability and specificity of PMIs is investigated systematically.

Design, cyclization, and optimization of MMP13-TIMP1 interaction-derived self-inhibitory peptides against chondrocyte senescence in osteoarthritis

The matrix metallopeptidase 13 (MMP13) is a central regulator of chondrocyte senescence that contributes to the development and progression of osteoarthritis (OA). In the present study, the native inhibitory structure of MMP13 in complex with its natural cognate inhibitor, the tissue inhibitor of metalloproteinases 1 (TIMP1), was modeled at atomic level using a grafting-based structural bioinformatics method with existing crystal structures. The modeled complex structure was then examined in detail, from which a TIMP1 inhibitory site that directly inserts into the active site of MMP13 enzyme was identified. The inhibitory site contains a coiled inhibitory loop (ILP) and a stretched N-terminal tail (NTT); they are highly structured in the intact MMP13-TIMP1 complex interface, but exhibit a large flexibility and intrinsic disorder when split from the interface context. In vitro binding assays demonstrated that the isolated ILP and NTT peptides cannot effectively rebind at the MMP13 active site (K_d > ~100 μM or = n.d.), although they have all key interacting residues in the enzyme inhibition. In silico simulations revealed that splitting of the peptide segments from TIMP1 inhibitory site does not influence the direct intermolecular interaction between MMP13 and the peptides substantially; instead, the large conformational flexibility of these isolated peptides in absence of interface context is primarily responsible for the affinity impairment, which would incur a considerable entropy penalty upon the peptide binding to MMP13. An extended version of ILP peptide, namely eILP (⁶³TPAMESVCGY⁷²), was redesigned with a rational strategy to derive a number of its cyclized counterparts by introducing a disulfide bridge across the peptide two-termini; the redesign reduces the peptide flexibility in free state and constrains the peptide pre-folding to a native-like conformation, which would help the peptide binding with minimized entropy penalty. Binding assays substantiated that the affinity K_d values of four designed cyclic peptides (, , and ) were improved to 23, 67, 42 and 18 μM, respectively, from the 96 μM of linear eILP peptide.

Exhaustive search of linear information encoding protein-peptide recognition

High-throughput in vitro methods have been extensively applied to identify linear information that encodes peptide recognition. However, these methods are limited in number of peptides, sequence variation, and length of peptides that can be explored, and often produce solutions that are not found in the cell. Despite the large number of methods developed to attempt addressing these issues, the exhaustive search of linear information encoding protein-peptide recognition has been so far physically unfeasible. Here, we describe a strategy, called DALEL, for the exhaustive search of linear sequence information encoded in proteins that bind to a common partner. We applied DALEL to explore binding specificity of SH3 domains in the budding yeast Saccharomyces cerevisiae. Using only the polypeptide sequences of SH3 domain binding proteins, we succeeded in identifying the majority of known SH3 binding sites previously discovered either in vitro or in vivo. Moreover, we discovered a number of sites with both non-canonical sequences and distinct properties that may serve ancillary roles in peptide recognition. We compared DALEL to a variety of state-of-the-art algorithms in the blind identification of known binding sites of the human Grb2 SH3 domain. We also benchmarked DALEL on curated biological motifs derived from the ELM database to evaluate the effect of increasing/decreasing the enrichment of the motifs. Our strategy can be applied in conjunction with experimental data of proteins interacting with a common partner to identify binding sites among them. Yet, our strategy can also be applied to any group of proteins of interest to identify enriched linear motifs or to exhaustively explore the space of linear information encoded in a polypeptide sequence. Finally, we have developed a webserver located at http://michnick.bcm.umontreal.ca/dalel, offering user-friendly interface and providing different scenarios utilizing DALEL.

Here we describe the first strategy for the exhaustive search of the linear information encoding protein-peptide recognition; an approach that has previously been physically unfeasible because the combinatorial space of polypeptide sequences is too vast. The search covers the entire space of sequences with no restriction on motif length or composition, and includes all possible combinations of amino acids at distinct positions of each sequence, as well as positions with correlated preferences for amino acids.

Prediction of GCRV virus-host protein interactome based on structural motif-domain interactions

BACKGROUND

Grass carp hemorrhagic disease, caused by grass carp reovirus (GCRV), is the most fatal causative agent in grass carp aquaculture. Protein-protein interactions between virus and host are one avenue through which GCRV can trigger infection and induce disease. Experimental approaches for the detection of host-virus interactome have many inherent limitations, and studies on protein-protein interactions between GCRV and its host remain rare.

RESULTS

In this study, based on known motif-domain interaction information, we systematically predicted the GCRV virus-host protein interactome by using motif-domain interaction pair searching strategy. These proteins derived from different domain families and were predicted to interact with different motif patterns in GCRV. JAM-A protein was successfully predicted to interact with motifs of GCRV Sigma1-like protein, and shared the similar binding mode compared with orthoreovirus. Differentially expressed genes during GCRV infection process were extracted and mapped to our predicted interactome, the overlapped genes displayed different tissue expression distributions on the whole, the overall expression level in intestinal is higher than that of other three tissues, which may suggest that the functions of these genes are more active in intestinal. Function annotation and pathway enrichment analysis revealed that the host targets were largely involved in signaling pathway and immune pathway, such as interferon-gamma signaling pathway, VEGF signaling pathway, EGF receptor signaling pathway, B cell activation, and T cell activation.

CONCLUSIONS

Although the predicted PPIs may contain some false positives due to limited data resource and poor research background in non-model species, the computational method still provide reasonable amount of interactions, which can be further validated by high throughput experiments. The findings of this work will contribute to the development of system biology for GCRV infectious diseases, and help guide the identification of novel receptors of GCRV in its host.

Evolution of domain-peptide interactions to coadapt specificity and affinity to functional diversity

Evolution of complexity in eukaryotic proteomes has arisen, in part, through emergence of modular independently folded domains mediating protein interactions via binding to short linear peptides in proteins. Over 30 years, structural properties and sequence preferences of these peptides have been extensively characterized. Less successful, however, were efforts to establish relationships between physicochemical properties and functions of domain-peptide interactions. To our knowledge, we have devised the first strategy to exhaustively explore the binding specificity of protein domain-peptide interactions. We applied the strategy to SH3 domains to determine the properties of their binding peptides starting from various experimental data. The strategy identified the majority (∼70%) of experimentally determined SH3 binding sites. We discovered mutual relationships among binding specificity, binding affinity, and structural properties and evolution of linear peptides. Remarkably, we found that these properties are also related to functional diversity, defined by depth of proteins within hierarchies of gene ontologies. Our results revealed that linear peptides evolved to coadapt specificity and affinity to functional diversity of domain-peptide interactions. Thus, domain-peptide interactions follow human-constructed gene ontologies, which suggest that our understanding of biological process hierarchies reflect the way chemical and thermodynamic properties of linear peptides and their interaction networks, in general, have evolved.

Motif mediated protein-protein interactions as drug targets

Protein-protein interactions (PPI) are involved in virtually every cellular process and thus represent an attractive target for therapeutic interventions. A significant number of protein interactions are frequently formed between globular domains and short linear peptide motifs (DMI). Targeting these DMIs has proven challenging and classical approaches to inhibiting such interactions with small molecules have had limited success. However, recent new approaches have led to the discovery of potent inhibitors, some of them, such as Obatoclax, ABT-199, AEG-40826 and SAH-p53-8 are likely to become approved drugs. These novel inhibitors belong to a wide range of different molecule classes, ranging from small molecules to peptidomimetics and biologicals. This article reviews the main reasons for limited success in targeting PPIs, discusses how successful approaches overcome these obstacles to discovery promising inhibitors for human protein double minute 2 (HDM2), B-cell lymphoma 2 (Bcl-2), X-linked inhibitor of apoptosis protein (XIAP), and provides a summary of the promising approaches currently in development that indicate the future potential of PPI inhibitors in drug discovery.

High-throughput methods for identification of protein-protein interactions involving short linear motifs

Interactions between modular domains and short linear motifs (3–10 amino acids peptide stretches) are crucial for cell signaling. The motifs typically reside in the disordered regions of the proteome and the interactions are often transient, allowing for rapid changes in response to changing stimuli. The properties that make domain-motif interactions suitable for cell signaling also make them difficult to capture experimentally and they are therefore largely underrepresented in the known protein-protein interaction networks. Most of the knowledge on domain-motif interactions is derived from low-throughput studies, although there exist dedicated high-throughput methods for the identification of domain-motif interactions. The methods include arrays of peptides or proteins, display of peptides on phage or yeast, and yeast-two-hybrid experiments. We here provide a survey of scalable methods for domain-motif interaction profiling. These methods have frequently been applied to a limited number of ubiquitous domain families. It is now time to apply them to a broader set of peptide binding proteins, to provide a comprehensive picture of the linear motifs in the human proteome and to link them to their potential binding partners. Despite the plethora of methods, it is still a challenge for most approaches to identify interactions that rely on post-translational modification or context dependent or conditional interactions, suggesting directions for further method development.

Fully Blind Peptide-Protein Docking with pepATTRACT

Peptide-protein interactions are ubiquitous in the cell and form an important part of the interactome. Computational docking methods can complement experimental characterization of these complexes, but current protocols are not applicable on the proteome scale. Here, we present a new fully blind flexible peptide-protein docking protocol, pepATTRACT, which combines a rapid coarse-grained global peptide docking search of the entire protein surface with a two-stage atomistic flexible refinement. Global unbound-unbound docking yielded near-native models for 70% of the docking cases when testing the protocol on the largest benchmark of peptide-protein complexes available to date. This performance is similar to that of state-of-the-art local docking protocols that rely on information about the binding site. Upon restricting the search to the peptide binding region, the resulting pepATTRACT-local approach outperformed existing methods. Docking scripts for pepATTRACT and pepATTRACT-local can be generated via a web interface at www.attract.ph.tum.de/peptide.html.

Predicting peptide-mediated interactions on a genome-wide scale

We describe a method to predict protein-protein interactions (PPIs) formed between structured domains and short peptide motifs. We take an integrative approach based on consensus patterns of known motifs in databases, structures of domain-motif complexes from the PDB and various sources of non-structural evidence. We combine this set of clues using a Bayesian classifier that reports the likelihood of an interaction and obtain significantly improved prediction performance when compared to individual sources of evidence and to previously reported algorithms. Our Bayesian approach was integrated into PrePPI, a structure-based PPI prediction method that, so far, has been limited to interactions formed between two structured domains. Around 80,000 new domain-motif mediated interactions were predicted, thus enhancing PrePPI’s coverage of the human protein interactome.

Complexes formed between a structured domain on one protein and an unstructured peptide on another are ubiquitous. However, they are often quite difficult to detect experimentally. The development of computational approaches to predict domain-motif interactions is therefore an important goal. We report a method to predict domain-motif interactions using a Bayesian approach to integrate evidence from a variety of sources, including three-dimensional structural and non-structural information. The method was applied to the entire human proteome and showed significant improvement over existing methods. The method was incorporated into PrePPI, a computational pipeline for the prediction of protein-protein interactions that relies heavily on structural information. Approximately 80,000 new interactions were detected. The new PrePPI database provides easy access to about 400,000 human protein-protein interactions and should thus constitute a valuable resource in a variety of biological applications including the characterization of molecular interaction networks and, more generally, in the study of interactions mediated by proteins in families that may not be extensively studied experimentally.

QSLiMFinder: improved short linear motif prediction using specific query protein data

Motivation: The sensitivity of de novo short linear motif (SLiM) prediction is limited by the number of patterns (the motif space) being assessed for enrichment. QSLiMFinder uses specific query protein information to restrict the motif space and thereby increase the sensitivity and specificity of predictions.

Results: QSLiMFinder was extensively benchmarked using known SLiM-containing proteins and simulated protein interaction datasets of real human proteins. Exploiting prior knowledge of a query protein likely to be involved in a SLiM-mediated interaction increased the proportion of true positives correctly returned and reduced the proportion of datasets returning a false positive prediction. The biggest improvement was seen if a short region of the query protein flanking the interaction site was known.

Availability and implementation: All the tools and data used in this study, including QSLiMFinder and the SLiMBench benchmarking software, are freely available under a GNU license as part of SLiMSuite, at: http://bioware.soton.ac.uk.

Contact:richard.edwards@unsw.edu.au

Supplementary information:Supplementary data are available at Bioinformatics online.

Structural basis of protein-protein interactions

Regulated interactions between proteins govern signaling pathways within and between cells. Structural studies on protein complexes formed reversibly and/or transiently illustrate the remarkable diversity of interactions, both in terms of interfacial size and nature. In recent years, "domain-peptide" interactions have gained much greater recognition and may be viewed as both pre-translational and posttranslational-dependent functional switches. Our understanding of the multistep regulation of auto-inhibited multidomain proteins has also grown. Their activity may be understood as the "combinatorial" output of multiple input signals, including phosphorylation, location, and mechanical force. The prospects for bridging the gap between the new "systems biology" data and the traditional "reductionist" data are also discussed.

Efficient prediction of human protein-protein interactions at a global scale

BACKGROUND

Our knowledge of global protein-protein interaction (PPI) networks in complex organisms such as humans is hindered by technical limitations of current methods.

RESULTS

On the basis of short co-occurring polypeptide regions, we developed a tool called MP-PIPE capable of predicting a global human PPI network within 3 months. With a recall of 23% at a precision of 82.1%, we predicted 172,132 putative PPIs. We demonstrate the usefulness of these predictions through a range of experiments.

CONCLUSIONS

The speed and accuracy associated with MP-PIPE can make this a potential tool to study individual human PPI networks (from genomic sequences alone) for personalized medicine.

Non-redundant unique interface structures as templates for modeling protein interactions

Improvements in experimental techniques increasingly provide structural data relating to protein-protein interactions. Classification of structural details of protein-protein interactions can provide valuable insights for modeling and abstracting design principles. Here, we aim to cluster protein-protein interactions by their interface structures, and to exploit these clusters to obtain and study shared and distinct protein binding sites. We find that there are 22604 unique interface structures in the PDB. These unique interfaces, which provide a rich resource of structural data of protein-protein interactions, can be used for template-based docking. We test the specificity of these non-redundant unique interface structures by finding protein pairs which have multiple binding sites. We suggest that residues with more than 40% relative accessible surface area should be considered as surface residues in template-based docking studies. This comprehensive study of protein interface structures can serve as a resource for the community. The dataset can be accessed at http://prism.ccbb.ku.edu.tr/piface.

Peptide docking and structure-based characterization of peptide binding: from knowledge to know-how

Peptide-mediated interactions are gaining increased attention due to their predominant roles in the many regulatory processes that involve dynamic interactions between proteins. The structures of such interactions provide an excellent starting point for their characterization and manipulation, and can provide leads for targeted inhibitor design. The relatively few experimentally determined structures of peptide-protein complexes can be complemented with an outburst of modeling approaches that have been introduced in recent years, with increasing accuracy and applicability to ever more systems. We review different methods to address the considerable challenges in modeling the binding of a short yet highly flexible peptide to its partner. These methods apply an array of sampling strategies and draw from a recent amassing of knowledge about the biophysical nature of peptide-protein interactions. We elaborate on applications of these structure-based approaches and in particular on the characterization of peptide binding specificity to different peptide-binding domains and enzymes. Such applications can identify new biological targets and thus complement our current view of protein-protein interactions in living organisms. Accurate peptide-protein docking is of particular importance in the light of increased appreciation of the crucial functional roles of disordered regions and the many linear binding motifs embedded within.

3did: a catalog of domain-based interactions of known three-dimensional structure

The database of 3D interacting domains (3did, available online for browsing and bulk download at http://3did.irbbarcelona.org) is a catalog of protein–protein interactions for which a high-resolution 3D structure is known. 3did collects and classifies all structural templates of domain–domain interactions in the Protein Data Bank, providing molecular details for such interactions. The current version also includes a pipeline for the discovery and annotation of novel domain–motif interactions. For every interaction, 3did identifies and groups different binding modes by clustering similar interfaces into ‘interaction topologies’. By maintaining a constantly updated collection of domain-based structural interaction templates, 3did is a reference source of information for the structural characterization of protein interaction networks. 3did is updated every 6 months.

Functional characterization of human genes from exon expression and RNA interference results

Complex biological systems comprise a large number of interacting molecules. The identification and detailed characterization of the functions of the involved genes and proteins are crucial for modeling and understanding such systems. To interrogate the various cellular processes, high-throughput techniques such as the Affymetrix Exon Array or RNA interference (RNAi) screens are powerful experimental approaches for functional genomics. However, they typically yield long gene lists that require computational methods to further analyze and functionally annotate the experimental results and to gain more insight into important molecular interactions. Here, we focus on bioinformatics software tools for the functional interpretation of exon expression data to discover alternative splicing events and their impact on gene and protein architecture, molecular networks, and pathways. We additionally demonstrate how to explore large lists of candidate genes as they also result from RNAi screens. In particular, our exemplary application studies show how to analyze the function of human genes that play a major role in human stem cells or viral infections.

Predicting binding within disordered protein regions to structurally characterised peptide-binding domains

Disordered regions of proteins often bind to structured domains, mediating interactions within and between proteins. However, it is difficult to identify a priori the short disordered regions involved in binding. We set out to determine if docking such peptide regions to peptide binding domains would assist in these predictions.We assembled a redundancy reduced dataset of SLiM (Short Linear Motif) containing proteins from the ELM database. We selected 84 sequences which had an associated PDB structures showing the SLiM bound to a protein receptor, where the SLiM was found within a 50 residue region of the protein sequence which was predicted to be disordered. First, we investigated the Vina docking scores of overlapping tripeptides from the 50 residue SLiM containing disordered regions of the protein sequence to the corresponding PDB domain. We found only weak discrimination of docking scores between peptides involved in binding and adjacent non-binding peptides in this context (AUC 0.58).Next, we trained a bidirectional recurrent neural network (BRNN) using as input the protein sequence, predicted secondary structure, Vina docking score and predicted disorder score. The results were very promising (AUC 0.72) showing that multiple sources of information can be combined to produce results which are clearly superior to any single source.We conclude that the Vina docking score alone has only modest power to define the location of a peptide within a larger protein region known to contain it. However, combining this information with other knowledge (using machine learning methods) clearly improves the identification of peptide binding regions within a protein sequence. This approach combining docking with machine learning is primarily a predictor of binding to peptide-binding sites, and is not intended as a predictor of specificity of binding to particular receptors.

A viral-human interactome based on structural motif-domain interactions captures the human infectome

Protein interactions between a pathogen and its host are fundamental in the establishment of the pathogen and underline the infection mechanism. In the present work, we developed a single predictive model for building a host-viral interactome based on the identification of structural descriptors from motif-domain interactions of protein complexes deposited in the Protein Data Bank (PDB). The structural descriptors were used for searching, in a database of protein sequences of human and five clinically important viruses; therefore, viral and human proteins sharing a descriptor were predicted as interacting proteins. The analysis of the host-viral interactome allowed to identify a set of new interactions that further explain molecular mechanism associated with viral infections and showed that it was able to capture human proteins already associated to viral infections (human infectome) and non-infectious diseases (human diseasome). The analysis of human proteins targeted by viral proteins in the context of a human interactome showed that their neighbors are enriched in proteins reported with differential expression under infection and disease conditions. It is expected that the findings of this work will contribute to the development of systems biology for infectious diseases, and help guide the rational identification and prioritization of novel drug targets.

Discovering interacting domains and motifs in protein-protein interactions

Many important biological processes, such as the signaling pathways, require protein-protein interactions (PPIs) that are designed for fast response to stimuli. These interactions are usually transient, easily formed, and disrupted, yet specific. Many of these transient interactions involve the binding of a protein domain to a short stretch (3-10) of amino acid residues, which can be characterized by a sequence pattern, i.e., a short linear motif (SLiM). We call these interacting domains and motifs domain-SLiM interactions. Existing methods have focused on discovering SLiMs in the interacting proteins' sequence data. With the recent increase in protein structures, we have a new opportunity to detect SLiMs directly from the proteins' 3D structures instead of their linear sequences. In this chapter, we describe a computational method called SLiMDIet to directly detect SLiMs on domain interfaces extracted from 3D structures of PPIs. SLiMDIet comprises two steps: (1) interaction interfaces belonging to the same domain are extracted and grouped together using structural clustering and (2) the extracted interaction interfaces in each cluster are structurally aligned to extract the corresponding SLiM. Using SLiMDIet, de novo SLiMs interacting with protein domains can be computationally detected from structurally clustered domain-SLiM interactions for PFAM domains which have available 3D structures in the PDB database.

The proteomics and interactomics of human erythrocytes

In this minireview, we focus on advances in our knowledge of the human erythrocyte proteome and interactome that have occurred since our seminal review on the topic published in 2007. As will be explained, the number of unique proteins has grown from 751 in 2007 to 2289 as of today. We describe how proteomics and interactomics tools have been used to probe critical protein changes in disorders impacting the blood. The primary example used is the work done on sickle cell disease where biomarkers of severity have been identified, protein changes in the erythrocyte membranes identified, pharmacoproteomic impact of hydroxyurea studied and interactomics used to identify erythrocyte protein changes that are predicted to have the greatest impact on protein interaction networks.

Protein-peptide complex prediction through fragment interaction patterns

The number of protein-peptide interactions in a cell is so large that experimental determination of all these complex structures would be a daunting task. Although homology modeling and refinement protocols have vastly improved the number and quality of predicted structural models, ab initio methods are still challenged by both the large number of possible docking sites and the conformational space accessible to flexible peptides. We present a method that addresses these challenges by sampling the entire accessible surface of a protein with a reduced conformational space of interacting backbone fragment pairs from unrelated structures. We demonstrate its potential by predicting ab initio the bound structure for a variety of protein-peptide complexes. In addition, we show the potential of our method for the discovery of domain interaction sites and domain-domain docking.

Disordered binding regions and linear motifs--bridging the gap between two models of molecular recognition

Intrinsically disordered proteins (IDPs) exist without the presence of a stable tertiary structure in isolation. These proteins are often involved in molecular recognition processes via their disordered binding regions that can recognize partner molecules by undergoing a coupled folding and binding process. The specific properties of disordered binding regions give way to specific, yet transient interactions that enable IDPs to play central roles in signaling pathways and act as hubs of protein interaction networks. An alternative model of protein-protein interactions with largely overlapping functional properties is offered by the concept of linear interaction motifs. This approach focuses on distilling a short consensus sequence pattern from proteins with a common interaction partner. These motifs often reside in disordered regions and are considered to mediate the interaction roughly independent from the rest of the protein. Although a connection between linear motifs and disordered binding regions has been established through common examples, the complementary nature of the two concepts has yet to be fully explored. In many cases the sequence based definition of linear motifs and the structural context based definition of disordered binding regions describe two aspects of the same phenomenon. To gain insight into the connection between the two models, prediction methods were utilized. We combined the regular expression based prediction of linear motifs with the disordered binding region prediction method ANCHOR, each specialized for either model to get the best of both worlds. The thorough analysis of the overlap of the two methods offers a bioinformatics tool for more efficient binding site prediction that can serve a wide range of practical implications. At the same time it can also shed light on the theoretical connection between the two co-existing interaction models.

Constructing structural networks of signaling pathways on the proteome scale

Proteins function through their interactions, and the availability of protein interaction networks could help in understanding cellular processes. However, the known structural data are limited and the classical network node-and-edge representation, where proteins are nodes and interactions are edges, shows only which proteins interact; not how they interact. Structural networks provide this information. Protein-protein interface structures can also indicate which binding partners can interact simultaneously and which are competitive, and can help forecasting potentially harmful drug side effects. Here, we use a powerful protein-protein interactions prediction tool which is able to carry out accurate predictions on the proteome scale to construct the structural network of the extracellular signal-regulated kinases (ERK) in the mitogen-activated protein kinase (MAPK) signaling pathway. This knowledge-based method, PRISM, is motif-based, and is combined with flexible refinement and energy scoring. PRISM predicts protein interactions based on structural and evolutionary similarity to known protein interfaces.

The identification of short linear motif-mediated interfaces within the human interactome

MOTIVATION

Eukaryotic proteins are highly modular, containing multiple interaction interfaces that mediate binding to a network of regulators and effectors. Recent advances in high-throughput proteomics have rapidly expanded the number of known protein-protein interactions (PPIs); however, the molecular basis for the majority of these interactions remains to be elucidated. There has been a growing appreciation of the importance of a subset of these PPIs, namely those mediated by short linear motifs (SLiMs), particularly the canonical and ubiquitous SH2, SH3 and PDZ domain-binding motifs. However, these motif classes represent only a small fraction of known SLiMs and outside these examples little effort has been made, either bioinformatically or experimentally, to discover the full complement of motif instances.

RESULTS

In this article, interaction data are analysed to identify and characterize an important subset of PPIs, those involving SLiMs binding to globular domains. To do this, we introduce iELM, a method to identify interactions mediated by SLiMs and add molecular details of the interaction interfaces to both interacting proteins. The method identifies SLiM-mediated interfaces from PPI data by searching for known SLiM-domain pairs. This approach was applied to the human interactome to identify a set of high-confidence putative SLiM-mediated PPIs.

AVAILABILITY

iELM is freely available at http://elmint.embl.de

CONTACT

toby.gibson@embl.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Prediction of short linear protein binding regions

Short linear motifs in proteins (typically 3-12 residues in length) play key roles in protein-protein interactions by frequently binding specifically to peptide binding domains within interacting proteins. Their tendency to be found in disordered segments of proteins has meant that they have often been overlooked. Here we present SLiMPred (short linear motif predictor), the first general de novo method designed to computationally predict such regions in protein primary sequences independent of experimentally defined homologs and interactors. The method applies machine learning techniques to predict new motifs based on annotated instances from the Eukaryotic Linear Motif database, as well as structural, biophysical, and biochemical features derived from the protein primary sequence. We have integrated these data sources and benchmarked the predictive accuracy of the method, and found that it performs equivalently to a predictor of protein binding regions in disordered regions, in addition to having predictive power for other classes of motif sites such as polyproline II helix motifs and short linear motifs lying in ordered regions. It will be useful in predicting peptides involved in potential protein associations and will aid in the functional characterization of proteins, especially of proteins lacking experimental information on structures and interactions. We conclude that, despite the diversity of motif sequences and structures, SLiMPred is a valuable tool for prioritizing potential interaction motifs in proteins.

Sequence- and interactome-based prediction of viral protein hotspots targeting host proteins: a case study for HIV Nef

Virus proteins alter protein pathways of the host toward the synthesis of viral particles by breaking and making edges via binding to host proteins. In this study, we developed a computational approach to predict viral sequence hotspots for binding to host proteins based on sequences of viral and host proteins and literature-curated virus-host protein interactome data. We use a motif discovery algorithm repeatedly on collections of sequences of viral proteins and immediate binding partners of their host targets and choose only those motifs that are conserved on viral sequences and highly statistically enriched among binding partners of virus protein targeted host proteins. Our results match experimental data on binding sites of Nef to host proteins such as MAPK1, VAV1, LCK, HCK, HLA-A, CD4, FYN, and GNB2L1 with high statistical significance but is a poor predictor of Nef binding sites on highly flexible, hoop-like regions. Predicted hotspots recapture CD8 cell epitopes of HIV Nef highlighting their importance in modulating virus-host interactions. Host proteins potentially targeted or outcompeted by Nef appear crowding the T cell receptor, natural killer cell mediated cytotoxicity, and neurotrophin signaling pathways. Scanning of HIV Nef motifs on multiple alignments of hepatitis C protein NS5A produces results consistent with literature, indicating the potential value of the hotspot discovery in advancing our understanding of virus-host crosstalk.

Protein design with fragment databases

Structure-based computational methods are popular tools for designing proteins and interactions between proteins because they provide the necessary insight and details required for rational engineering. Here, we first argue that large-scale databases of fragments contain a discrete but complete set of building blocks that can be used to design structures. We show that these structural alphabets can be saturated to provide conformational ensembles that sample the native structure space around energetic minima. Second, we show that catalogs of interaction patterns hold the key to overcome the lack of scaffolds when computationally designing protein interactions. Finally, we illustrate the power of database-driven computational protein design methods by recent successful applications and discuss what challenges remain to push this field forward.

Computational design of peptide ligands

Peptides possess several attractive features when compared to small molecule and protein therapeutics, such as high structural compatibility with target proteins, the ability to disrupt protein-protein interfaces, and small size. Efficient design of high-affinity peptide ligands via rational methods has been a major obstacle to the development of this potential drug class. However, structural insights into the architecture of protein-peptide interfaces have recently culminated in several computational approaches for the rational design of peptides that target proteins. These methods provide a valuable alternative to experimental high-resolution structures of target protein-peptide complexes, bringing closer the dream of in silico designed peptides for therapeutic applications.

Three-dimensional modeling of protein interactions and complexes is going 'omics

High-throughput interaction discovery initiatives have revealed the existence of hundreds of multiprotein complexes whose functions are regulated through thousands of protein-protein interactions (PPIs). However, the structural details of these interactions, often necessary to understand their function, are only available for a tiny fraction, and the experimental difficulties surrounding complex structure determination make computational modeling techniques paramount. In this manuscript, we critically review some of the most recent developments in the field of structural bioinformatics applied to the modeling of protein interactions and complexes, from large macromolecular machines to domain-domain and peptide-mediated interactions. In particular, we place a special emphasis on those methods that can be applied in a proteome-wide manner, and discuss how they will help in the ultimate objective of building 3D interactome networks.

A regression framework incorporating quantitative and negative interaction data improves quantitative prediction of PDZ domain-peptide interaction from primary sequence

Motivation: Predicting protein interactions involving peptide recognition domains is essential for understanding the many important biological processes they mediate. It is important to consider the binding strength of these interactions to help us construct more biologically relevant protein interaction networks that consider cellular context and competition between potential binders.

Results: We developed a novel regression framework that considers both positive (quantitative) and negative (qualitative) interaction data available for mouse PDZ domains to quantitatively predict interactions between PDZ domains, a large peptide recognition domain family, and their peptide ligands using primary sequence information. First, we show that it is possible to learn from existing quantitative and negative interaction data to infer the relative binding strength of interactions involving previously unseen PDZ domains and/or peptides given their primary sequence. Performance was measured using cross-validated hold out testing and testing with previously unseen PDZ domain–peptide interactions. Second, we find that incorporating negative data improves quantitative interaction prediction. Third, we show that sequence similarity is an important prediction performance determinant, which suggests that experimentally collecting additional quantitative interaction data for underrepresented PDZ domain subfamilies will improve prediction.

Availability and Implementation: The Matlab code for our SemiSVR predictor and all data used here are available at http://baderlab.org/Data/PDZAffinity.

Contact:gary.bader@utoronto.ca; dengnaiyang@cau.edu.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

3did: identification and classification of domain-based interactions of known three-dimensional structure

The database of three-dimensional interacting domains (3did) is a collection of protein interactions for which high-resolution three-dimensional structures are known. 3did exploits the availability of structural data to provide molecular details on interactions between two globular domains as well as novel domain–peptide interactions, derived using a recently published method from our lab. The interface residues are presented for each interaction type individually, plus global domain interfaces at which one or more partners (domains or peptides) bind. The 3did web server at http://3did.irbbarcelona.org visualizes these interfaces along with atomic details of individual interactions using Jmol. The complete contents are also available for download.