PIER: protein interface recognition for structural proteomics

Protein-protein interaction prediction with deep learning: A comprehensive review

Most proteins perform their biological function by interacting with themselves or other molecules. Thus, one may obtain biological insights into protein functions, disease prevalence, and therapy development by identifying protein-protein interactions (PPI). However, finding the interacting and non-interacting protein pairs through experimental approaches is labour-intensive and time-consuming, owing to the variety of proteins. Hence, protein-protein interaction and protein-ligand binding problems have drawn attention in the fields of bioinformatics and computer-aided drug discovery. Deep learning methods paved the way for scientists to predict the 3-D structure of proteins from genomes, predict the functions and attributes of a protein, and modify and design new proteins to provide desired functions. This review focuses on recent deep learning methods applied to problems including predicting protein functions, protein-protein interaction and their sites, protein-ligand binding, and protein design.

Hierarchical representation for PPI sites prediction

Background

Protein–protein interactions have pivotal roles in life processes, and aberrant interactions are associated with various disorders. Interaction site identification is key for understanding disease mechanisms and design new drugs. Effective and efficient computational methods for the PPI prediction are of great value due to the overall cost of experimental methods. Promising results have been obtained using machine learning methods and deep learning techniques, but their effectiveness depends on protein representation and feature selection.

Results

We define a new abstraction of the protein structure, called hierarchical representations, considering and quantifying spatial and sequential neighboring among amino acids. We also investigate the effect of molecular abstractions using the Graph Convolutional Networks technique to classify amino acids as interface and no-interface ones. Our study takes into account three abstractions, hierarchical representations, contact map, and the residue sequence, and considers the eight functional classes of proteins extracted from the Protein–Protein Docking Benchmark 5.0. The performance of our method, evaluated using standard metrics, is compared to the ones obtained with some state-of-the-art protein interface predictors. The analysis of the performance values shows that our method outperforms the considered competitors when the considered molecules are structurally similar.

Conclusions

The hierarchical representation can capture the structural properties that promote the interactions and can be used to represent proteins with unknown structures by codifying only their sequential neighboring. Analyzing the results, we conclude that classes should be arranged according to their architectures rather than functions.

Predicting protein–membrane interfaces of peripheral membrane proteins using ensemble machine learning

Abnormal protein–membrane attachment is involved in deregulated cellular pathways and in disease. Therefore, the possibility to modulate protein–membrane interactions represents a new promising therapeutic strategy for peripheral membrane proteins that have been considered so far undruggable. A major obstacle in this drug design strategy is that the membrane-binding domains of peripheral membrane proteins are usually unknown. The development of fast and efficient algorithms predicting the protein–membrane interface would shed light into the accessibility of membrane–protein interfaces by drug-like molecules. Herein, we describe an ensemble machine learning methodology and algorithm for predicting membrane-penetrating amino acids. We utilize available experimental data from the literature for training 21 machine learning classifiers and meta-classifiers. Evaluation of the best ensemble classifier model accuracy yields a macro-averaged F1 score = 0.92 and a Matthews correlation coefficient = 0.84 for predicting correctly membrane-penetrating amino acids on unknown proteins of a validation set. The python code for predicting protein–membrane interfaces of peripheral membrane proteins is available at https://github.com/zoecournia/DREAMM.

Mechanisms of allosteric and mixed mode aromatase inhibitors

Aromatase (CYP19) catalyzes the last biosynthetic step of estrogens in mammals and is a primary drug target for hormone-related breast cancer. However, treatment with aromatase inhibitors is often associated with adverse effects and drug resistance. In this study, we used virtual screening targeting a predicted cytochrome P450 reductase binding site on aromatase to discover four novel non-steroidal aromatase inhibitors. The inhibitors have potencies comparable to the noncompetitive tamoxifen metabolite, endoxifen. Our two most potent inhibitors, AR11 and AR13, exhibit both mixed-type and competitive-type inhibition. The cytochrome P450 reductase-CYP19 coupling interface likely acts as a transient binding site. Our modeling shows that our inhibitors bind better at different sites near the catalytic site. Our results predict the location of multiple ligand binding sites on aromatase. The combination of modeling and experimental results supports the important role of the reductase binding interface as a low affinity, promiscuous ligand binding site. Our new inhibitors may be useful as alternative chemical scaffolds that may show different adverse effects profiles than current clinically used aromatase inhibitors.

In silico Approaches for the Design and Optimization of Interfering Peptides Against Protein-Protein Interactions

Large contact surfaces of protein-protein interactions (PPIs) remain to be an ongoing issue in the discovery and design of small molecule modulators. Peptides are intrinsically capable of exploring larger surfaces, stable, and bioavailable, and therefore bear a high therapeutic value in the treatment of various diseases, including cancer, infectious diseases, and neurodegenerative diseases. Given these promising properties, a long way has been covered in the field of targeting PPIs via peptide design strategies. In silico tools have recently become an inevitable approach for the design and optimization of these interfering peptides. Various algorithms have been developed to scrutinize the PPI interfaces. Moreover, different databases and software tools have been created to predict the peptide structures and their interactions with target protein complexes. High-throughput screening of large peptide libraries against PPIs; "hotspot" identification; structure-based and off-structure approaches of peptide design; 3D peptide modeling; peptide optimization strategies like cyclization; and peptide binding energy evaluation are among the capabilities of in silico tools. In the present study, the most recent advances in the field of in silico approaches for the design of interfering peptides against PPIs will be reviewed. The future perspective of the field and its advantages and limitations will also be pinpointed.

Simplified geometric representations of protein structures identify complementary interaction interfaces

Protein-protein interactions are critical to protein function, but three-dimensional (3D) arrangements of interacting proteins have proven hard to predict, even given the identities and 3D structures of the interacting partners. Specifically, identifying the relevant pairwise interaction surfaces remains difficult, often relying on shape complementarity with molecular docking while accounting for molecular motions to optimize rigid 3D translations and rotations. However, such approaches can be computationally expensive, and faster, less accurate approximations may prove useful for large-scale prediction and assembly of 3D structures of multi-protein complexes. We asked if a reduced representation of protein geometry retains enough information about molecular properties to predict pairwise protein interaction interfaces that are tolerant of limited structural rearrangements. Here, we describe a reduced representation of 3D protein accessible surfaces on which molecular properties such as charge, hydrophobicity, and evolutionary rate can be easily mapped, implemented in the MorphProt package. Pairs of surfaces are compared to rapidly assess partner-specific potential surface complementarity. On two available benchmarks of 185 overall known protein complexes, we observe predictions comparable to other structure-based tools at correctly identifying protein interaction surfaces. Furthermore, we examined the effect of molecular motion through normal mode simulation on a benchmark receptor-ligand pair and observed no marked loss of predictive accuracy for distortions of up to 6 Å Cα-RMSD. Thus, a shape reduction of protein surfaces retains considerable information about surface complementarity, offers enhanced speed of comparison relative to more complex geometric representations, and exhibits tolerance to conformational changes.

Classification and prediction of protein-protein interaction interface using machine learning algorithm

Structural insight of the protein-protein interaction (PPI) interface can provide knowledge about the kinetics, thermodynamics and molecular functions of the complex while elucidating its role in diseases and further enabling it as a potential therapeutic target. However, owing to experimental lag in solving protein-protein complex structures, three-dimensional (3D) knowledge of the PPI interfaces can be gained via computational approaches like molecular docking and post-docking analyses. Despite development of numerous docking tools and techniques, success in identification of native like interfaces based on docking score functions is limited. Hence, we employed an in-depth investigation of the structural features of the interface that might successfully delineate native complexes from non-native ones. We identify interface properties, which show statistically significant difference between native and non-native interfaces belonging to homo and hetero, protein-protein complexes. Utilizing these properties, a support vector machine (SVM) based classification scheme has been implemented to differentiate native and non-native like complexes generated using docking decoys. Benchmarking and comparative analyses suggest very good performance of our SVM classifiers. Further, protein interactions, which are proven via experimental findings but not resolved structurally, were subjected to this approach where 3D-models of the complexes were generated and most likely interfaces were predicted. A web server called Protein Complex Prediction by Interface Properties (PCPIP) is developed to predict whether interface of a given protein-protein dimer complex resembles known protein interfaces. The server is freely available at http://www.hpppi.iicb.res.in/pcpip/ .

Decoding Protein-protein Interactions: An Overview

Drug discovery has focused on the paradigm "one drug, one target" for a long time. However, small molecules can act at multiple macromolecular targets, which serves as the basis for drug repurposing. In an effort to expand the target space, and given advances in X-ray crystallography, protein-protein interactions have become an emerging focus area of drug discovery enterprises. Proteins interact with other biomolecules and it is this intricate network of interactions that determines the behavior of the system and its biological processes. In this review, we briefly discuss networks in disease, followed by computational methods for protein-protein complex prediction. Computational methodologies and techniques employed towards objectives such as protein-protein docking, protein-protein interactions, and interface predictions are described extensively. Docking aims at producing a complex between proteins, while interface predictions identify a subset of residues on one protein that could interact with a partner, and protein-protein interaction sites address whether two proteins interact. In addition, approaches to predict hot spots and binding sites are presented along with a representative example of our internal project on the chemokine CXC receptor 3 B-isoform and predictive modeling with IP10 and PF4.

Leveraging a gain-of-function allele of Caenorhabditis elegans paqr-1 to elucidate membrane homeostasis by PAQR proteins

The C. elegans proteins PAQR-2 (a homolog of the human seven-transmembrane domain AdipoR1 and AdipoR2 proteins) and IGLR-2 (a homolog of the mammalian LRIG proteins characterized by a single transmembrane domain and the presence of immunoglobulin domains and leucine-rich repeats in their extracellular portion) form a complex that protects against plasma membrane rigidification by promoting the expression of fatty acid desaturases and the incorporation of polyunsaturated fatty acids into phospholipids, hence increasing membrane fluidity. In the present study, we leveraged a novel gain-of-function allele of PAQR-1, a PAQR-2 paralog, to carry out structure-function studies. We found that the transmembrane domains of PAQR-2 are responsible for its functional requirement for IGLR-2, that PAQR-1 does not require IGLR-2 but acts via the same pathway as PAQR-2, and that the divergent N-terminal cytoplasmic domains of the PAQR-1 and PAQR-2 proteins serve a regulatory function and may regulate access to the catalytic site of these proteins. We also show that overexpression of human AdipoR1 or AdipoR2 alone is sufficient to confer increased palmitic acid resistance in HEK293 cells, and thus act in a manner analogous to the PAQR-1 gain-of-function allele.

CoRNeA: A Pipeline to Decrypt the Inter-Protein Interfaces from Amino Acid Sequence Information

Decrypting the interface residues of the protein complexes provides insight into the functions of the proteins and, hence, the overall cellular machinery. Computational methods have been devised in the past to predict the interface residues using amino acid sequence information, but all these methods have been majorly applied to predict for prokaryotic protein complexes. Since the composition and rate of evolution of the primary sequence is different between prokaryotes and eukaryotes, it is important to develop a method specifically for eukaryotic complexes. Here, we report a new hybrid pipeline for predicting the protein-protein interaction interfaces in a pairwise manner from the amino acid sequence information of the interacting proteins. It is based on the framework of Co-evolution, machine learning (Random Forest), and Network Analysis named CoRNeA trained specifically on eukaryotic protein complexes. We use Co-evolution, physicochemical properties, and contact potential as major group of features to train the Random Forest classifier. We also incorporate the intra-contact information of the individual proteins to eliminate false positives from the predictions keeping in mind that the amino acid sequence of a protein also holds information for its own folding and not only the interface propensities. Our prediction on example datasets shows that CoRNeA not only enhances the prediction of true interface residues but also reduces false positive rates significantly.

Structural Characterization of Carbonic Anhydrase VIII and Effects of Missense Single Nucleotide Variations to Protein Structure and Function

Human carbonic anhydrase 8 (CA-VIII) is an acatalytic isoform of the α -CA family. Though the protein cannot hydrate CO2, CA-VIII is essential for calcium (Ca2+) homeostasis within the body, and achieves this by allosterically inhibiting the binding of inositol 1,4,5-triphosphate (IP3) to the IP3 receptor type 1 (ITPR1) protein. However, the mechanism of interaction of CA-VIII to ITPR1 is not well understood. In addition, functional defects to CA-VIII due to non-synonymous single nucleotide polymorphisms (nsSNVs) result in Ca2+ dysregulation and the development of the phenotypes such as cerebellar ataxia, mental retardation and disequilibrium syndrome 3 (CAMRQ3). The pathogenesis of CAMRQ3 is also not well understood. The structure and function of CA-VIII was characterised, and pathogenesis of CAMRQ3 investigated. Structural and functional characterisation of CA-VIII was conducted through SiteMap and CPORT to identify potential binding site residues. The effects of four pathogenic nsSNVs, S100A, S100P, G162R and R237Q, and two benign S100L and E109D variants on CA-VIII structure and function was then investigated using molecular dynamics (MD) simulations, dynamic cross correlation (DCC) and dynamic residue network (DRN) analysis. SiteMap and CPORT analyses identified 38 unique CA-VIII residues that could potentially bind to ITPR1. MD analysis revealed less conformational sampling within the variant proteins and highlighted potential increases to variant protein rigidity. Dynamic cross correlation (DCC) showed that wild-type (WT) protein residue motion is predominately anti-correlated, with variant proteins showing no correlation to greater residue correlation. DRN revealed variant-associated increases to the accessibility of the N-terminal binding site residues, which could have implications for associations with ITPR1, and further highlighted differences to the mechanism of benign and pathogenic variants. SNV presence is associated with a reduction to the usage of Trp37 in all variants, which has implications for CA-VIII stability. The differences to variant mechanisms can be further investigated to understand pathogenesis of CAMRQ3, enhancing precision medicine-related studies into CA-VIII.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Ambidextrous helical nanotubes from self-assembly of designed helical hairpin motifs

Tandem repeat proteins exhibit native designability and represent potentially useful scaffolds for the construction of synthetic biomimetic assemblies. We have designed 2 synthetic peptides, HEAT_R1 and LRV_M3Δ1, based on the consensus sequences of single repeats of thermophilic HEAT (PBS_HEAT) and Leucine-Rich Variant (LRV) structural motifs, respectively. Self-assembly of the peptides afforded high-aspect ratio helical nanotubes. Cryo-electron microscopy with direct electron detection was employed to analyze the structures of the solvated filaments. The 3D reconstructions from the cryo-EM maps led to atomic models for the HEAT_R1 and LRV_M3Δ1 filaments at resolutions of 6.0 and 4.4 Å, respectively. Surprisingly, despite sequence similarity at the lateral packing interface, HEAT_R1 and LRV_M3Δ1 filaments adopt the opposite helical hand and differ significantly in helical geometry, while retaining a local conformation similar to previously characterized repeat proteins of the same class. The differences in the 2 filaments could be rationalized on the basis of differences in cohesive interactions at the lateral and axial interfaces. These structural data reinforce previous observations regarding the structural plasticity of helical protein assemblies and the need for high-resolution structural analysis. Despite these observations, the native designability of tandem repeat proteins offers the opportunity to engineer novel helical nanotubes. Moreover, the resultant nanotubes have independently addressable and chemically distinguishable interior and exterior surfaces that would facilitate applications in selective recognition, transport, and release.

Advances in In-silico B-cell Epitope Prediction

Identification of B-cell epitopes in target antigens is one of the most crucial steps for epitopebased vaccine development, immunodiagnostic tests, antibody production, and disease diagnosis and therapy. Experimental methods for B-cell epitope mapping are time consuming, costly and labor intensive; in the meantime, various in-silico methods are proposed to predict both linear and conformational B-cell epitopes. The accurate identification of B-cell epitopes presents major challenges for immunoinformaticians. In this paper, we have comprehensively reviewed in-silico methods for B-cell epitope identification. The aim of this review is to stimulate the development of better tools which could improve the identification of B-cell epitopes, and further for the development of therapeutic antibodies and diagnostic tools.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Exploring the potential of 3D Zernike descriptors and SVM for protein-protein interface prediction

Background

The correct determination of protein–protein interaction interfaces is important for understanding disease mechanisms and for rational drug design. To date, several computational methods for the prediction of protein interfaces have been developed, but the interface prediction problem is still not fully understood. Experimental evidence suggests that the location of binding sites is imprinted in the protein structure, but there are major differences among the interfaces of the various protein types: the characterising properties can vary a lot depending on the interaction type and function. The selection of an optimal set of features characterising the protein interface and the development of an effective method to represent and capture the complex protein recognition patterns are of paramount importance for this task.

Results

In this work we investigate the potential of a novel local surface descriptor based on 3D Zernike moments for the interface prediction task. Descriptors invariant to roto-translations are extracted from circular patches of the protein surface enriched with physico-chemical properties from the HQI8 amino acid index set, and are used as samples for a binary classification problem. Support Vector Machines are used as a classifier to distinguish interface local surface patches from non-interface ones. The proposed method was validated on 16 classes of proteins extracted from the Protein–Protein Docking Benchmark 5.0 and compared to other state-of-the-art protein interface predictors (SPPIDER, PrISE and NPS-HomPPI).

Conclusions

The 3D Zernike descriptors are able to capture the similarity among patterns of physico-chemical and biochemical properties mapped on the protein surface arising from the various spatial arrangements of the underlying residues, and their usage can be easily extended to other sets of amino acid properties. The results suggest that the choice of a proper set of features characterising the protein interface is crucial for the interface prediction task, and that optimality strongly depends on the class of proteins whose interface we want to characterise. We postulate that different protein classes should be treated separately and that it is necessary to identify an optimal set of features for each protein class.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2043-3) contains supplementary material, which is available to authorized users.

Interactome INSIDER: a structural interactome browser for genomic studies

We present Interactome INSIDER, a tool to link genomic variant information with structural protein-protein interactomes. Underlying this tool is the application of machine learning to predict protein interaction interfaces for 185,957 protein interactions with previously unresolved interfaces, in human and 7 model organisms, including the entire experimentally determined human binary interactome. Predicted interfaces exhibit similar functional properties as known interfaces, including enrichment for disease mutations and recurrent cancer mutations. Through 2,164 de novo mutagenesis experiments, we show that mutations of predicted and known interface residues disrupt interactions at a similar rate, and much more frequently than mutations outside of predicted interfaces. To spur functional genomic studies, Interactome INSIDER (http://interactomeinsider.yulab.org) enables users to identify whether variants or disease mutations are enriched in known and predicted interaction interfaces at various resolutions. Users may explore known population variants, disease mutations, and somatic cancer mutations, or upload their own set of mutations for this purpose.

Ofatumumab Monoclonal Antibody Affinity Maturation Through in silico Modeling

Background:

Ofatumumab, an anti-CD20 mAb, was approved in 2009 for the treatment of chronic lymphocytic leukemia. This mAb acts through immune-mediated mechanisms, in particular complement-dependent cytotoxicity and antibody-dependent cellular cytotoxicity by natural killer cells as well as antibody-dependent phagocytosis by macrophages. Apoptosis induction is another mechanism of this antibody. Computational docking is the method of predicting the conformation of an antibody-antigen from its separated elements. Validation of the designed antibodies is carried out by docking tools. Increased affinity enhances the biological action of the antibody, which in turn improves the therapeutic effects. Furthermore, the increased antibody affinity can reduce the therapeutic dose of the antibody, resulting in lower toxicity and handling cost.

Methods:

Considering the importance of this issue, using in silico analysis such as docking and molecular dynamics, we aimed to find the important amino acids of the Ofatumumab antibody and then replaced these amino acids with others to improve antibody-binding affinity. Finally, we examined the binding affinity of antibody variants to antigen.

Results:

Our findings showed that variant 3 mutations have improved the characteristics of antibody binding compared to normal Ofatumumab antibodies.

Conclusion:

The designed anti-CD20 antibodies showed potentiality for improved affinity in comparison to commercial Ofatumumab.

IntPred: a structure-based predictor of protein-protein interaction sites

Motivation

Protein–protein interactions are vital for protein function with the average protein having between three and ten interacting partners. Knowledge of precise protein–protein interfaces comes from crystal structures deposited in the Protein Data Bank (PDB), but only 50% of structures in the PDB are complexes. There is therefore a need to predict protein–protein interfaces in silico and various methods for this purpose. Here we explore the use of a predictor based on structural features and which exploits random forest machine learning, comparing its performance with a number of popular established methods.

Results

On an independent test set of obligate and transient complexes, our IntPred predictor performs well (MCC = 0.370, ACC = 0.811, SPEC = 0.916, SENS = 0.411) and compares favourably with other methods. Overall, IntPred ranks second of six methods tested with SPPIDER having slightly better overall performance (MCC = 0.410, ACC = 0.759, SPEC = 0.783, SENS = 0.676), but considerably worse specificity than IntPred. As with SPPIDER, using an independent test set of obligate complexes enhanced performance (MCC = 0.381) while performance is somewhat reduced on a dataset of transient complexes (MCC = 0.303). The trade-off between sensitivity and specificity compared with SPPIDER suggests that the choice of the appropriate tool is application-dependent.

Availability and implementation

IntPred is implemented in Perl and may be downloaded for local use or run via a web server at www.bioinf.org.uk/intpred/.

Supplementary information

Supplementary data are available at Bioinformatics online.

A structural perspective on the interactions of TRAF6 and Basigin during the onset of melanoma: A molecular dynamics simulation study

Metastatic melanoma is the most fatal type of skin cancer. The roles of matrix metalloproteinases (MMPs) have well been established in the onset of melanoma. Basigin (BSG) belongs to the immunoglobulin superfamily and is critical for induction of extracellular MMPs during the onset of various cancers including melanoma. Tumor necrosis factor receptor-associated factor 6 (TRAF6) is an E3-ligase that interacts with BSG and mediates its membrane localization, which leads to MMP expression in melanoma cells. This makes TRAF6 a potential therapeutic target in melanoma. We here conducted protein-protein interaction studies on TRAF6 and BSG to get molecular level insights of the reactions. The structure of human BSG was constructed by protein threading. Molecular-docking method was applied to develop the TRAF6-BSG complex. The refined docked complex was further optimized by molecular dynamics simulations. Results from binding free energy, surface properties, and electrostatic interaction analysis indicate that Lys340 and Glu417 of TRAF6 play as the anchor residues in the protein interaction interface. The current study will be helpful in designing specific modulators of TRAF6 to control melanoma metastasis.

Structural analysis and insight into Zika virus NS5 mediated interferon inhibition

The Zika virus outbreak in 2015-2016 is the largest of its kind for which WHO declared a Public Health Emergency of International Concerns. No FDA approved drug is available for the treatment of the viral infection. The interaction of flavivirus NS5 protein with SIAH2 ubiquitin ligase has been previously known. NS5 of Zika virus has been implicated in the degradation of STAT2 protein, which activates interferon-stimulated antiviral activity. Based on our proposition that NS5 utilizes SIAH2-mediated proteasomal degradation of STAT2, an in-silico study was carried out to characterize the protein-protein interactions between NS5, SIAH2 and STAT2 proteins. The aim of our study was to identify the amino acid residues of NS5 involved in IFN antagonism as well as to find the association between NS5, SIAH2 and STAT2 to predict the interaction pattern of these proteins. Analysis proposed that NS5 recruits SIAH2 for the ubiquitination-dependent degradation of STAT2. NS5 residues involved in interaction with SIAH2 and/or STAT2 were found to be mostly conserved across related flaviviruses. These are novel findings regarding the Zika virus and require confirmation through experimental approaches.

Network analysis and in silico prediction of protein-protein interactions with applications in drug discovery

Protein-protein interactions (PPIs) are vital to maintaining cellular homeostasis. Several PPI dysregulations have been implicated in the etiology of various diseases and hence PPIs have emerged as promising targets for drug discovery. Surface residues and hotspot residues at the interface of PPIs form the core regions, which play a key role in modulating cellular processes such as signal transduction and are used as starting points for drug design. In this review, we briefly discuss how PPI networks (PPINs) inferred from experimentally characterized PPI data have been utilized for knowledge discovery and how in silico approaches to PPI characterization can contribute to PPIN-based biological research. Next, we describe the principles of in silico PPI prediction and survey the existing PPI and PPI site prediction servers that are useful for drug discovery. Finally, we discuss the potential of in silico PPI prediction in drug discovery.

Computational Approaches for Predicting Binding Partners, Interface Residues, and Binding Affinity of Protein-Protein Complexes

Studying protein-protein interactions leads to a better understanding of the underlying principles of several biological pathways. Cost and labor-intensive experimental techniques suggest the need for computational methods to complement them. Several such state-of-the-art methods have been reported for analyzing diverse aspects such as predicting binding partners, interface residues, and binding affinity for protein-protein complexes with reliable performance. However, there are specific drawbacks for different methods that indicate the need for their improvement. This review highlights various available computational algorithms for analyzing diverse aspects of protein-protein interactions and endorses the necessity for developing new robust methods for gaining deep insights about protein-protein interactions.

Elucidating the druggable interface of protein-protein interactions using fragment docking and coevolutionary analysis

Protein-protein interactions play a central role in cellular function. Improving the understanding of complex formation has many practical applications, including the rational design of new therapeutic agents and the mechanisms governing signal transduction networks. The generally large, flat, and relatively featureless binding sites of protein complexes pose many challenges for drug design. Fragment docking and direct coupling analysis are used in an integrated computational method to estimate druggable protein-protein interfaces. (i) This method explores the binding of fragment-sized molecular probes on the protein surface using a molecular docking-based screen. (ii) The energetically favorable binding sites of the probes, called hot spots, are spatially clustered to map out candidate binding sites on the protein surface. (iii) A coevolution-based interface interaction score is used to discriminate between different candidate binding sites, yielding potential interfacial targets for therapeutic drug design. This approach is validated for important, well-studied disease-related proteins with known pharmaceutical targets, and also identifies targets that have yet to be studied. Moreover, therapeutic agents are proposed by chemically connecting the fragments that are strongly bound to the hot spots.

Identification of a Substrate-selective Exosite within the Metalloproteinase Anthrax Lethal Factor

The metalloproteinase anthrax lethal factor (LF) is secreted by Bacillus anthracis to promote disease virulence through disruption of host signaling pathways. LF is a highly specific protease, exclusively cleaving mitogen-activated protein kinase kinases (MKKs) and rodent NLRP1B (NACHT leucine-rich repeat and pyrin domain-containing protein 1B). How LF achieves such restricted substrate specificity is not understood. Previous studies have suggested the existence of an exosite interaction between LF and MKKs that promotes cleavage efficiency and specificity. Through a combination of in silico prediction and site-directed mutagenesis, we have mapped an exosite to a non-catalytic region of LF. Mutations within this site selectively impair proteolysis of full-length MKKs yet have no impact on cleavage of short peptide substrates. Although this region appears important for cleaving all LF protein substrates, we found that mutation of specific residues within the exosite differentially affects MKK and NLRP1B cleavage in vitro and in cultured cells. One residue in particular, Trp-271, is essential for cleavage of MKK3, MKK4, and MKK6 but dispensable for targeting of MEK1, MEK2, and NLRP1B. Analysis of chimeric substrates suggests that this residue interacts with the MKK catalytic domain. We found that LF-W271A blocked ERK phosphorylation and growth in a melanoma cell line, suggesting that it may provide a highly selective inhibitor of MEK1/2 for use as a cancer therapeutic. These findings provide insight into how a bacterial toxin functions to specifically impair host signaling pathways and suggest a general strategy for mapping protease exosite interactions.

Progress and challenges in predicting protein interfaces

The majority of biological processes are mediated via protein-protein interactions. Determination of residues participating in such interactions improves our understanding of molecular mechanisms and facilitates the development of therapeutics. Experimental approaches to identifying interacting residues, such as mutagenesis, are costly and time-consuming and thus, computational methods for this purpose could streamline conventional pipelines. Here we review the field of computational protein interface prediction. We make a distinction between methods which address proteins in general and those targeted at antibodies, owing to the radically different binding mechanism of antibodies. We organize the multitude of currently available methods hierarchically based on required input and prediction principles to provide an overview of the field.

Computational prediction of protein interfaces: A review of data driven methods

Reliably pinpointing which specific amino acid residues form the interface(s) between a protein and its binding partner(s) is critical for understanding the structural and physicochemical determinants of protein recognition and binding affinity, and has wide applications in modeling and validating protein interactions predicted by high-throughput methods, in engineering proteins, and in prioritizing drug targets. Here, we review the basic concepts, principles and recent advances in computational approaches to the analysis and prediction of protein-protein interfaces. We point out caveats for objectively evaluating interface predictors, and discuss various applications of data-driven interface predictors for improving energy model-driven protein-protein docking. Finally, we stress the importance of exploiting binding partner information in reliably predicting interfaces and highlight recent advances in this emerging direction.

Survey of Natural Language Processing Techniques in Bioinformatics

Informatics methods, such as text mining and natural language processing, are always involved in bioinformatics research. In this study, we discuss text mining and natural language processing methods in bioinformatics from two perspectives. First, we aim to search for knowledge on biology, retrieve references using text mining methods, and reconstruct databases. For example, protein-protein interactions and gene-disease relationship can be mined from PubMed. Then, we analyze the applications of text mining and natural language processing techniques in bioinformatics, including predicting protein structure and function, detecting noncoding RNA. Finally, numerous methods and applications, as well as their contributions to bioinformatics, are discussed for future use by text mining and natural language processing researchers.

Evolutionary and structural analyses of heterodimeric proteins composed of subunits with same fold

Heterodimeric proteins with homologous subunits of same fold are involved in various biological processes. The objective of this study is to understand the evolution of structural and functional features of such heterodimers. Using a non-redundant dataset of 70 such heterodimers of known 3D structure and an independent dataset of 173 heterodimers from yeast, we note that the mean sequence identity between interacting homologous subunits is only 23-24% suggesting that, generally, highly diverged paralogues assemble to form such a heterodimer. We also note that the functional roles of interacting subunits/domains are generally quite different. This suggests that, though the interacting subunits/domains are homologous, the high evolutionary divergence characterize their high functional divergence which contributes to a gross function for the heterodimer considered as a whole. The inverse relationship between sequence identity and RMSD of interacting homologues in heterodimers is not followed. We also addressed the question of formation of homodimers of the subunits of heterodimers by generating models of fictitious homodimers on the basis of the 3D structures of the heterodimers. Interaction energies associated with these homodimers suggests that, in overwhelming majority of the cases, such homodimers are unlikely to be stable. Majority of the homologues of heterodimers of known structures form heterodimers (51.8%) and a small proportion (14.6%) form homodimers. Comparison of 3D structures of heterodimers with homologous homodimers suggests that interfacial nature of residues is not well conserved. In over 90% of the cases we note that the interacting subunits of heterodimers are co-localized in the cell.

Fusion of selected regions of mycobacterial antigens for enhancing sensitivity in serodiagnosis of tuberculosis

Serodiagnosis of tuberculosis requires detection of antibodies against multiple antigens of Mycobacterium tuberculosis, because antibody profiles differ among the patients. Using fusion proteins with epitopes from two or more antigens would facilitate in the detection of multiple antibodies. Fusion constructs tn1FbpC1-tnPstS1 and tn2FbpC1-tnPstS1 were produced by linking truncated regions of variable lengths from FbpC1 to the N-terminus of the truncated PstS1. Similarly a truncated fragment of HSP was linked to the N-terminus of a truncated fragment from FbpC1 to produce tnHSP-tn1FbpC1. ELISA analysis of the plasma samples of TB patients against tn2FbpC1-tnPstS1 showed 72.2% sensitivity which is nearly the same as the expected combined value for the two individual antigens. However, the sensitivity of tn1FbpC1-tnPstS1 was lowered to 60%. tnHSP-tn1FbpC1 showed 67.7% sensitivity which is slightly less than the expected combined value for the two individual antigens, but still significantly higher than that of each of the individual antigen. Data for secondary structure analysis by CD spectrometry was in reasonable agreement with the X-ray crystallographic data of the native proteins and the predicted structure of the fusion proteins. Comparative molecular modeling suggests that the epitopes of the constituent proteins are better exposed in tn2FbpC1-tnPstS1 as compared to those in tn1FbpC1-tnPstS1. Therefore, removal of the N-terminal non-epitopic region of FbpC1 from 34-96 amino acids seems to have unmasked at least some of the epitopes, resulting in greater sensitivity. The high level of sensitivity of tn2FbpC1-tnPstS1 and tnHSP-tn1FbpC1, not reported before, shows that these fusion proteins have great potential for use in serodiagnosis of tuberculosis.

Predicting protein interface residues using easily accessible on-line resources

It has been more than a decade since the completion of the Human Genome Project that provided us with a complete list of human proteins. The next obvious task is to figure out how various parts interact with each other. On that account, we review 10 methods for protein interface prediction, which are freely available as web servers. In addition, we comparatively evaluate their performance on a common data set comprising different quality target structures. We find that using experimental structures and high-quality homology models, structure-based methods outperform those using only protein sequences, with global template-based approaches providing the best performance. For moderate-quality models, sequence-based methods often perform better than those structure-based techniques that rely on fine atomic details. We note that post-processing protocols implemented in several methods quantitatively improve the results only for experimental structures, suggesting that these procedures should be tuned up for computer-generated models. Finally, we anticipate that advanced meta-prediction protocols are likely to enhance interface residue prediction. Notwithstanding further improvements, easily accessible web servers already provide the scientific community with convenient resources for the identification of protein-protein interaction sites.

Molecular modeling and computational analyses suggests that the Sinorhizobium meliloti periplasmic regulator protein ExoR adopts a superhelical fold and is controlled by a unique mechanism of proteolysis

The Sinorhizobium meliloti periplasmic ExoR protein and the ExoS/ChvI two-component system form a regulatory mechanism that directly controls the transformation of free-living to host-invading cells. In the absence of crystal structures, understanding the molecular mechanism of interaction between ExoR and the ExoS sensor, which is believed to drive the key regulatory step in the invasion process, remains a major challenge. In this study, we present a theoretical structural model of the active form of ExoR protein, ExoRm , generated using computational methods. Our model suggests that ExoR possesses a super-helical fold comprising 12 α-helices forming six Sel1-like repeats, including two that were unidentified in previous studies. This fold is highly conducive to mediating protein-protein interactions and this is corroborated by the identification of putative protein binding sites on the surface of the ExoRm protein. Our studies reveal two novel insights: (a) an extended conformation of the third Sel1-like repeat that might be important for ExoR regulatory function and (b) a buried proteolytic site that implies a unique proteolytic mechanism. This study provides new and interesting insights into the structure of S. meliloti ExoR, lays the groundwork for elaborating the molecular mechanism of ExoRm cleavage, ExoRm -ExoS interactions, and studies of ExoR homologs in other bacterial host interactions.

Algorithmic approaches to protein-protein interaction site prediction

Interaction sites on protein surfaces mediate virtually all biological activities, and their identification holds promise for disease treatment and drug design. Novel algorithmic approaches for the prediction of these sites have been produced at a rapid rate, and the field has seen significant advancement over the past decade. However, the most current methods have not yet been reviewed in a systematic and comprehensive fashion. Herein, we describe the intricacies of the biological theory, datasets, and features required for modern protein-protein interaction site (PPIS) prediction, and present an integrative analysis of the state-of-the-art algorithms and their performance. First, the major sources of data used by predictors are reviewed, including training sets, evaluation sets, and methods for their procurement. Then, the features employed and their importance in the biological characterization of PPISs are explored. This is followed by a discussion of the methodologies adopted in contemporary prediction programs, as well as their relative performance on the datasets most recently used for evaluation. In addition, the potential utility that PPIS identification holds for rational drug design, hotspot prediction, and computational molecular docking is described. Finally, an analysis of the most promising areas for future development of the field is presented.

Discovery of novel membrane binding structures and functions

The function of a protein is determined by its intrinsic activity in the context of its subcellular distribution. Membranes localize proteins within cellular compartments and govern their specific activities. Discovering such membrane-protein interactions is important for understanding biological mechanisms and could uncover novel sites for therapeutic intervention. We present a method for detecting membrane interactive proteins and their exposed residues that insert into lipid bilayers. Although the development process involved analysis of how C1b, C2, ENTH, FYVE, Gla, pleckstrin homology (PH), and PX domains bind membranes, the resulting membrane optimal docking area (MODA) method yields predictions for a given protein of known three-dimensional structures without referring to canonical membrane-targeting modules. This approach was tested on the Arf1 GTPase, ATF2 acetyltransferase, von Willebrand factor A3 domain, and Neisseria gonorrhoeae MsrB protein and further refined with membrane interactive and non-interactive FAPP1 and PKD1 pleckstrin homology domains, respectively. Furthermore we demonstrate how this tool can be used to discover unprecedented membrane binding functions as illustrated by the Bro1 domain of Alix, which was revealed to recognize lysobisphosphatidic acid (LBPA). Validation of novel membrane-protein interactions relies on other techniques such as nuclear magnetic resonance spectroscopy (NMR), which was used here to map the sites of micelle interaction. Together this indicates that genome-wide identification of known and novel membrane interactive proteins and sites is now feasible and provides a new tool for functional annotation of the proteome.

CRF-based models of protein surfaces improve protein-protein interaction site predictions

Background

The identification of protein-protein interaction sites is a computationally challenging task and important for understanding the biology of protein complexes. There is a rich literature in this field. A broad class of approaches assign to each candidate residue a real-valued score that measures how likely it is that the residue belongs to the interface. The prediction is obtained by thresholding this score.

Some probabilistic models classify the residues on the basis of the posterior probabilities. In this paper, we introduce pairwise conditional random fields (pCRFs) in which edges are not restricted to the backbone as in the case of linear-chain CRFs utilized by Li et al. (2007). In fact, any 3D-neighborhood relation can be modeled. On grounds of a generalized Viterbi inference algorithm and a piecewise training process for pCRFs, we demonstrate how to utilize pCRFs to enhance a given residue-wise score-based protein-protein interface predictor on the surface of the protein under study. The features of the pCRF are solely based on the interface predictions scores of the predictor the performance of which shall be improved.

Results

We performed three sets of experiments with synthetic scores assigned to the surface residues of proteins taken from the data set PlaneDimers compiled by Zellner et al. (2011), from the list published by Keskin et al. (2004) and from the very recent data set due to Cukuroglu et al. (2014). That way we demonstrated that our pCRF-based enhancer is effective given the interface residue score distribution and the non-interface residue score are unimodal.

Moreover, the pCRF-based enhancer is also successfully applicable, if the distributions are only unimodal over a certain sub-domain. The improvement is then restricted to that domain. Thus we were able to improve the prediction of the PresCont server devised by Zellner et al. (2011) on PlaneDimers.

Conclusions

Our results strongly suggest that pCRFs form a methodological framework to improve residue-wise score-based protein-protein interface predictors given the scores are appropriately distributed. A prototypical implementation of our method is accessible at http://ppicrf.informatik.uni-goettingen.de/index.html.

Structural insights into the activation of the RhoA GTPase by the lymphoid blast crisis (Lbc) oncoprotein

The small GTPase RhoA promotes deregulated signaling upon interaction with lymphoid blast crisis (Lbc), the oncogenic form of A-kinase anchoring protein 13 (AKAP13). The onco-Lbc protein is a hyperactive Rho-specific guanine nucleotide exchange factor (GEF), but its structural mechanism has not been reported despite its involvement in cardiac hypertrophy and cancer causation. The pleckstrin homology (PH) domain of Lbc is located at the C-terminal end of the protein and is shown here to specifically recognize activated RhoA rather than lipids. The isolated dbl homology (DH) domain can function as an independent activator with an enhanced activity. However, the DH domain normally does not act as a solitary Lbc interface with RhoA-GDP. Instead it is negatively controlled by the PH domain. In particular, the DH helical bundle is coupled to the structurally dependent PH domain through a helical linker, which reduces its activity. Together the two domains form a rigid scaffold in solution as evidenced by small angle x-ray scattering and ¹H,¹³C,¹⁵N-based NMR spectroscopy. The two domains assume a “chair” shape with its back possessing independent GEF activity and the PH domain providing a broad seat for RhoA-GTP docking rather than membrane recognition. This provides structural and dynamical insights into how DH and PH domains work together in solution to support regulated RhoA activity. Mutational analysis supports the bifunctional PH domain mediation of DH-RhoA interactions and explains why the tandem domain is required for controlled GEF signaling.

Structural evaluation of BTK and PKCδ mediated phosphorylation of MAL at positions Tyr86 and Tyr106

A number of diseases including sepsis, rheumatoid arthritis, diabetes, cardiovascular diseases and hyperinflammatory immune disorders have been associated with Toll like receptor (TLR) 2 and TLR4. Endogenous adaptor protein known as MyD88 adapter-like protein (MAL) bind exclusively to the cytosolic portions of TLR2 and TLR4 to initiate downstream signalling. Brutons tyrosine kinase (BTK) and protein kinase C delta (PKCδ) have been implicated to phosphorylate MAL and activate it to initiate downstream signalling. BTK has been associated with phosphorylation at positions Tyr86 and Tyr106, necessary for the activation of MAL but definite residual target of PKCδ in MAL is still to be explored. To produce a better understanding of the functional domains involved in the formation of MAL-kinase complexes, computer-aided studies were used to characterize the protein-protein interactions (PPIs) of phosphorylated BTK and PKCδ with MAL. Docking and physicochemical studies indicated that BTK was involved in close contact with Tyr86 and Tyr106 of MAL whereas PKCδ may phosphorylate Tyr106 only. Moreover, the electrostatics charge distribution of binding interfaces of BTK and PKCδ were distinct but compatible with respective regions of MAL. Our results implicate that position of Tyr86 is specifically phosphorylated by BTK whereas Tyr106 can be phosphorylated by competitive action of both BTK and PKCδ. Additionally, the residues of MAL which are necessary for interaction with TLR2, TLR4, MyD88 and SOCS-1 also play their roles in maintaining interaction with kinases and can be targeted in future to reduce TLR2 and TLR4 induced pathological responses.

The interaction of CRM1 and the nuclear pore protein Tpr

While much has been devoted to the study of transport mechanisms through the nuclear pore complex (NPC), the specifics of interactions and binding between export transport receptors and the NPC periphery have remained elusive. Recent work has demonstrated a binding interaction between the exportin CRM1 and the unstructured carboxylic tail of Tpr, on the nuclear basket. Strong evidence suggests that this interaction is vital to the functions of CRM1. Using molecular dynamics simulations and a newly refined method for determining binding regions, we have identified nine candidate binding sites on CRM1 for C-Tpr. These include two adjacent to RanGTP--from which one is blocked in the absence of RanGTP--and three next to the binding region of the cargo Snurportin. We report two additional interaction sites between C-Tpr and Snurportin, suggesting a possible role for Tpr import into the nucleus. Using bioinformatics tools we have conducted conservation analysis and functional residue prediction investigations to identify which parts of the obtained binding sites are inherently more important and should be highlighted. Also, a novel measure based on the ratio of available solvent accessible surface (RASAS) is proposed for monitoring the ligand/receptor binding process.

Transient protein-protein interface prediction: datasets, features, algorithms, and the RAD-T predictor

Background

Transient protein-protein interactions (PPIs), which underly most biological processes, are a prime target for therapeutic development. Immense progress has been made towards computational prediction of PPIs using methods such as protein docking and sequence analysis. However, docking generally requires high resolution structures of both of the binding partners and sequence analysis requires that a significant number of recurrent patterns exist for the identification of a potential binding site. Researchers have turned to machine learning to overcome some of the other methods’ restrictions by generalising interface sites with sets of descriptive features. Best practices for dataset generation, features, and learning algorithms have not yet been identified or agreed upon, and an analysis of the overall efficacy of machine learning based PPI predictors is due, in order to highlight potential areas for improvement.

Results

The presence of unknown interaction sites as a result of limited knowledge about protein interactions in the testing set dramatically reduces prediction accuracy. Greater accuracy in labelling the data by enforcing higher interface site rates per domain resulted in an average 44% improvement across multiple machine learning algorithms. A set of 10 biologically unrelated proteins that were consistently predicted on with high accuracy emerged through our analysis. We identify seven features with the most predictive power over multiple datasets and machine learning algorithms. Through our analysis, we created a new predictor, RAD-T, that outperforms existing non-structurally specializing machine learning protein interface predictors, with an average 59% increase in MCC score on a dataset with a high number of interactions.

Conclusion

Current methods of evaluating machine-learning based PPI predictors tend to undervalue their performance, which may be artificially decreased by the presence of un-identified interaction sites. Changes to predictors’ training sets will be integral to the future progress of interface prediction by machine learning methods. We reveal the need for a larger test set of well studied proteins or domain-specific scoring algorithms to compensate for poor interaction site identification on proteins in general.

A conserved motif flags acyl carrier proteins for β-branching in polyketide synthesis

Type I PKSs often utilise programmed β-branching, via enzymes of an “HMG-CoA synthase (HCS) cassette”, to incorporate various side chains at the second carbon from the terminal carboxylic acid of growing polyketide backbones. We identified a strong sequence motif in Acyl Carrier Proteins (ACPs) where β-branching is known. Substituting ACPs confirmed a correlation of ACP type with β-branching specificity. While these ACPs often occur in tandem, NMR analysis of tandem β-branching ACPs indicated no ACP-ACP synergistic effects and revealed that the conserved sequence motif forms an internal core rather than an exposed patch. Modelling and mutagenesis identified ACP Helix III as a probable anchor point of the ACP-HCS complex whose position is determined by the core. Mutating the core affects ACP functionality while ACP-HCS interface substitutions modulate system specificity. Our method for predicting β-carbon branching expands the potential for engineering novel polyketides and lays a basis for determining specificity rules.

Binding interface prediction by combining protein-protein docking results

We developed a method called residue contact frequency (RCF), which uses the complex structures generated by the protein-protein docking algorithm ZDOCK to predict interface residues. Unlike interface prediction algorithms that are based on monomers alone, RCF is binding partner specific. We evaluated the performance of RCF using the area under the precision-recall (PR) curve (AUC) on a large protein docking Benchmark. RCF (AUC = 0.44) performed as well as meta-PPISP (AUC = 0.43), which is one of the best monomer-based interface prediction methods. In addition, we test a support vector machine (SVM) to combine RCF with meta-PPISP and another monomer-based interface prediction algorithm Evolutionary Trace to further improve the performance. We found that the SVM that combined RCF and meta-PPISP achieved the best performance (AUC = 0.47). We used RCF to predict the binding interfaces of proteins that can bind to multiple partners and RCF was able to correctly predict interface residues that are unique for the respective binding partners. Furthermore, we found that residues that contributed greatly to binding affinity (hotspot residues) had significantly higher RCF than other residues.

Bioinformatics resources and tools for conformational B-cell epitope prediction

Identification of epitopes which invoke strong humoral responses is an essential issue in the field of immunology. Localizing epitopes by experimental methods is expensive in terms of time, cost, and effort; therefore, computational methods feature for its low cost and high speed was employed to predict B-cell epitopes. In this paper, we review the recent advance of bioinformatics resources and tools in conformational B-cell epitope prediction, including databases, algorithms, web servers, and their applications in solving problems in related areas. To stimulate the development of better tools, some promising directions are also extensively discussed.