Prediction of interface residues in protein-protein complexes by a consensus neural network method: test against NMR data

Encoding the space of protein-protein binding interfaces by artificial intelligence

The physical interactions between proteins are largely determined by the structural properties at their binding interfaces. It was found that the binding interfaces in distinctive protein complexes are highly similar. The structural properties underlying different binding interfaces could be further captured by artificial intelligence. In order to test this hypothesis, we broke protein-protein binding interfaces into pairs of interacting fragments. We employed a generative model to encode these interface fragment pairs in a low-dimensional latent space. After training, new conformations of interface fragment pairs were generated. We found that, by only using a small number of interface fragment pairs that were generated by artificial intelligence, we were able to guide the assembly of protein complexes into their native conformations. These results demonstrate that the conformational space of fragment pairs at protein-protein binding interfaces is highly degenerate. Features in this degenerate space can be well characterized by artificial intelligence. In summary, our machine learning method will be potentially useful to search for and predict the conformations of unknown protein-protein interactions.

Protein-protein interaction site prediction by model ensembling with hybrid feature and self-attention

BACKGROUND

Protein-protein interactions (PPIs) are crucial in various biological functions and cellular processes. Thus, many computational approaches have been proposed to predict PPI sites. Although significant progress has been made, these methods still have limitations in encoding the characteristics of each amino acid in sequences. Many feature extraction methods rely on the sliding window technique, which simply merges all the features of residues into a vector. The importance of some key residues may be weakened in the feature vector, leading to poor performance.

RESULTS

We propose a novel sequence-based method for PPI sites prediction. The new network model, PPINet, contains multiple feature processing paths. For a residue, the PPINet extracts the features of the targeted residue and its context separately. These two types of features are processed by two paths in the network and combined to form a protein representation, where the two types of features are of relatively equal importance. The model ensembling technique is applied to make use of more features. The base models are trained with different features and then ensembled via stacking. In addition, a data balancing strategy is presented, by which our model can get significant improvement on highly unbalanced data.

CONCLUSION

The proposed method is evaluated on a fused dataset constructed from Dset186, Dset_72, and PDBset_164, as well as the public Dset_448 dataset. Compared with current state-of-the-art methods, the performance of our method is better than the others. In the most important metrics, such as AUPRC and recall, it surpasses the second-best programmer on the latter dataset by 6.9% and 4.7%, respectively. We also demonstrated that the improvement is essentially due to using the ensemble model, especially, the hybrid feature. We share our code for reproducibility and future research at https://github.com/CandiceCong/StackingPPINet .

In silico evaluation of the inhibitory potential of nucleocapsid inhibitors of SARS-CoV-2: a binding and energetic perspective

The COVID-19 outbreak brought on by the SARS-CoV-2 virus continued to infect a sizable population worldwide. The SARS-CoV-2 nucleocapsid (N) protein is the most conserved RNA-binding structural protein and is a desirable target because of its involvement in viral transcription and replication. Based on this aspect, this study focused to repurpose antiviral compounds approved or in development for treating COVID-19. The inhibitors chosen are either FDA-approved or are currently being studied in clinical trials against COVID-19. Initially, they were designed to target stress granules and other RNA biology. We have utilized structure-based molecular docking and all-atom molecular dynamics (MD) simulation approach to investigate in detail the binding energy and binding modes of the different anti-N inhibitors to N protein. The result showed that five drugs including Silmitasterib, Ninetanidinb, Ternatin, Luteolin, Fedratinib, PJ34, and Zotatafin were found interacting with RNA binding sites as well as to predicted protein interface with higher binding energy. Overall, drug binding increases the stability of the complex with maximum stability found in the order, Silmitasertib > PJ34 > Zotatatafin. In addition, the frustration changes due to drug binding brings a decrease in local frustration and this decrease is mainly observed in α-helix, β3, β5, and β6 strands and are important for drug binding. Our in-silico data suggest that an effective interaction occurs for some of the tested drugs and prompt their further validation to reduce the rapid outspreading of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

MpbPPI: a multi-task pre-training-based equivariant approach for the prediction of the effect of amino acid mutations on protein-protein interactions

The accurate prediction of the effect of amino acid mutations for protein-protein interactions (PPI $\Delta \Delta G$) is a crucial task in protein engineering, as it provides insight into the relevant biological processes underpinning protein binding and provides a basis for further drug discovery. In this study, we propose MpbPPI, a novel multi-task pre-training-based geometric equivariance-preserving framework to predict PPI  $\Delta \Delta G$. Pre-training on a strictly screened pre-training dataset is employed to address the scarcity of protein-protein complex structures annotated with PPI $\Delta \Delta G$ values. MpbPPI employs a multi-task pre-training technique, forcing the framework to learn comprehensive backbone and side chain geometric regulations of protein-protein complexes at different scales. After pre-training, MpbPPI can generate high-quality representations capturing the effective geometric characteristics of labeled protein-protein complexes for downstream $\Delta \Delta G$ predictions. MpbPPI serves as a scalable framework supporting different sources of mutant-type (MT) protein-protein complexes for flexible application. Experimental results on four benchmark datasets demonstrate that MpbPPI is a state-of-the-art framework for PPI $\Delta \Delta G$ predictions. The data and source code are available at https://github.com/arantir123/MpbPPI.

E(3) equivariant graph neural networks for robust and accurate protein-protein interaction site prediction

Artificial intelligence-powered protein structure prediction methods have led to a paradigm-shift in computational structural biology, yet contemporary approaches for predicting the interfacial residues (i.e., sites) of protein-protein interaction (PPI) still rely on experimental structures. Recent studies have demonstrated benefits of employing graph convolution for PPI site prediction, but ignore symmetries naturally occurring in 3-dimensional space and act only on experimental coordinates. Here we present EquiPPIS, an E(3) equivariant graph neural network approach for PPI site prediction. EquiPPIS employs symmetry-aware graph convolutions that transform equivariantly with translation, rotation, and reflection in 3D space, providing richer representations for molecular data compared to invariant convolutions. EquiPPIS substantially outperforms state-of-the-art approaches based on the same experimental input, and exhibits remarkable robustness by attaining better accuracy with predicted structural models from AlphaFold2 than what existing methods can achieve even with experimental structures. Freely available at https://github.com/Bhattacharya-Lab/EquiPPIS, EquiPPIS enables accurate PPI site prediction at scale.

Two point mutations in protocadherin-1 disrupt hantavirus recognition and afford protection against lethal infection

Andes virus (ANDV) and Sin Nombre virus (SNV) are the etiologic agents of severe hantavirus cardiopulmonary syndrome (HCPS) in the Americas for which no FDA-approved countermeasures are available. Protocadherin-1 (PCDH1), a cadherin-superfamily protein recently identified as a critical host factor for ANDV and SNV, represents a new antiviral target; however, its precise role remains to be elucidated. Here, we use computational and experimental approaches to delineate the binding surface of the hantavirus glycoprotein complex on PCDH1's first extracellular cadherin repeat domain. Strikingly, a single amino acid residue in this PCDH1 surface influences the host species-specificity of SNV glycoprotein-PCDH1 interaction and cell entry. Mutation of this and a neighboring residue substantially protects Syrian hamsters from pulmonary disease and death caused by ANDV. We conclude that PCDH1 is a bona fide entry receptor for ANDV and SNV whose direct interaction with hantavirus glycoproteins could be targeted to develop new interventions against HCPS.

Identification of Protein Complexes by Integrating Protein Abundance and Interaction Features Using a Deep Learning Strategy

Many essential cellular functions are carried out by multi-protein complexes that can be characterized by their protein-protein interactions. The interactions between protein subunits are critically dependent on the strengths of their interactions and their cellular abundances, both of which span orders of magnitude. Despite many efforts devoted to the global discovery of protein complexes by integrating large-scale protein abundance and interaction features, there is still room for improvement. Here, we integrated >7000 quantitative proteomic samples with three published affinity purification/co-fractionation mass spectrometry datasets into a deep learning framework to predict protein-protein interactions (PPIs), followed by the identification of protein complexes using a two-stage clustering strategy. Our deep-learning-technique-based classifier significantly outperformed recently published machine learning prediction models and in the process captured 5010 complexes containing over 9000 unique proteins. The vast majority of proteins in our predicted complexes exhibited low or no tissue specificity, which is an indication that the observed complexes tend to be ubiquitously expressed throughout all cell types and tissues. Interestingly, our combined approach increased the model sensitivity for low abundant proteins, which amongst other things allowed us to detect the interaction of MCM10, which connects to the replicative helicase complex via the MCM6 protein. The integration of protein abundances and their interaction features using a deep learning approach provided a comprehensive map of protein-protein interactions and a unique perspective on possible novel protein complexes.

ProtInteract: A deep learning framework for predicting protein-protein interactions

Proteins mainly perform their functions by interacting with other proteins. Protein-protein interactions underpin various biological activities such as metabolic cycles, signal transduction, and immune response. However, due to the sheer number of proteins, experimental methods for finding interacting and non-interacting protein pairs are time-consuming and costly. We therefore developed the ProtInteract framework to predict protein-protein interaction. ProtInteract comprises two components: first, a novel autoencoder architecture that encodes each protein's primary structure to a lower-dimensional vector while preserving its underlying sequence attributes. This leads to faster training of the second network, a deep convolutional neural network (CNN) that receives encoded proteins and predicts their interaction under three different scenarios. In each scenario, the deep CNN predicts the class of a given encoded protein pair. Each class indicates different ranges of confidence scores corresponding to the probability of whether a predicted interaction occurs or not. The proposed framework features significantly low computational complexity and relatively fast response. The contributions of this work are twofold. First, ProtInteract assimilates the protein's primary structure into a pseudo-time series. Therefore, we leverage the nature of the time series of proteins and their physicochemical properties to encode a protein's amino acid sequence into a lower-dimensional vector space. This approach enables extracting highly informative sequence attributes while reducing computational complexity. Second, the ProtInteract framework utilises this information to identify protein interactions with other proteins based on its amino acid configuration. Our results suggest that the proposed framework performs with high accuracy and efficiency in predicting protein-protein interactions.

Interplay between C1-inhibitor and group IIA secreted phospholipase A2 impairs their respective function

High levels of human group IIA secreted phospholipase A₂ (hGIIA) have been associated with various inflammatory disease conditions. We have recently shown that hGIIA activity and concentration are increased in the plasma of patients with hereditary angioedema due to C1-inhibitor deficiency (C1-INH-HAE) and negatively correlate with C1-INH plasma activity. In this study, we analyzed whether the presence of both hGIIA and C1-INH impairs their respective function on immune cells. hGIIA, but not recombinant and plasma-derived C1-INH, stimulates the production of IL-6, CXCL8, and TNF-α from peripheral blood mononuclear cells (PBMCs). PBMC activation mediated by hGIIA is blocked by RO032107A, a specific hGIIA inhibitor. Interestingly, C1-INH inhibits the hGIIA-induced production of IL-6, TNF-α, and CXCL8, while it does not affect hGIIA enzymatic activity. On the other hand, hGIIA reduces the capacity of C1-INH at inhibiting C1-esterase activity. Spectroscopic and molecular docking studies suggest a possible interaction between hGIIA and C1-INH but further experiments are needed to confirm this hypothesis. Together, these results provide evidence for a new interplay between hGIIA and C1-INH, which may be important in the pathophysiology of hereditary angioedema.

Building Protein-Protein Interaction Graph Database Using Neo4j

A cell's various components interact with each other in a coordinated manner to respond to environmental cues and intracellular signals. Compared to the other biological networks, the protein-protein interaction (PPI) is mostly responsible for maintaining signaling pathways. Increasing numbers of experimentally verified and predicted PPIs in plants demand a scalable platform to deal with large and complex datasets. Network/graph data can be organized and analyzed using different tools. This chapter uses Neo4j, a graph database management system, to store and analyze plant PPI networks. To make the graph database and analyze network centrality, we used Arabidopsis interactome-1 main (AI-1MAIN) PPI network.

Integrated structure-based protein interface prediction

Background

Identifying protein interfaces can inform how proteins interact with their binding partners, uncover the regulatory mechanisms that control biological functions and guide the development of novel therapeutic agents. A variety of computational approaches have been developed for predicting a protein’s interfacial residues from its known sequence and structure. Methods using the known three-dimensional structures of proteins can be template-based or template-free. Template-based methods have limited success in predicting interfaces when homologues with known complex structures are not available to use as templates. The prediction performance of template-free methods that only rely only upon proteins’ intrinsic properties is limited by the amount of biologically relevant features that can be included in an interface prediction model.

Results

We describe the development of an integrated method for protein interface prediction (ISPIP) to explore the hypothesis that the efficacy of a computational prediction method of protein binding sites can be enhanced by using a combination of methods that rely on orthogonal structure-based properties of a query protein, combining and balancing both template-free and template-based features. ISPIP is a method that integrates these approaches through simple linear or logistic regression models and more complex decision tree models. On a diverse test set of 156 query proteins, ISPIP outperforms each of its individual classifiers in identifying protein binding interfaces.

Conclusions

The integrated method captures the best performance of individual classifiers and delivers an improved interface prediction. The method is robust and performs well even when one of the individual classifiers performs poorly on a particular query protein. This work demonstrates that integrating orthogonal methods that depend on different structural properties of proteins performs better at interface prediction than any individual classifier alone.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04852-2.

ProB-Site: Protein Binding Site Prediction Using Local Features

Protein–protein interactions (PPIs) are responsible for various essential biological processes. This information can help develop a new drug against diseases. Various experimental methods have been employed for this purpose; however, their application is limited by their cost and time consumption. Alternatively, computational methods are considered viable means to achieve this crucial task. Various techniques have been explored in the literature using the sequential information of amino acids in a protein sequence, including machine learning and deep learning techniques. The current efficiency of interaction-site prediction still has growth potential. Hence, a deep neural network-based model, ProB-site, is proposed. ProB-site utilizes sequential information of a protein to predict its binding sites. The proposed model uses evolutionary information and predicted structural information extracted from sequential information of proteins, generating three unique feature sets for every amino acid in a protein sequence. Then, these feature sets are fed to their respective sub-CNN architecture to acquire complex features. Finally, the acquired features are concatenated and classified using fully connected layers. This methodology performed better than state-of-the-art techniques because of the selection of the best features and contemplation of local information of each amino acid.

Multi-task learning to leverage partially annotated data for PPI interface prediction

Protein protein interactions (PPI) are crucial for protein functioning, nevertheless predicting residues in PPI interfaces from the protein sequence remains a challenging problem. In addition, structure-based functional annotations, such as the PPI interface annotations, are scarce: only for about one-third of all protein structures residue-based PPI interface annotations are available. If we want to use a deep learning strategy, we have to overcome the problem of limited data availability. Here we use a multi-task learning strategy that can handle missing data. We start with the multi-task model architecture, and adapted it to carefully handle missing data in the cost function. As related learning tasks we include prediction of secondary structure, solvent accessibility, and buried residue. Our results show that the multi-task learning strategy significantly outperforms single task approaches. Moreover, only the multi-task strategy is able to effectively learn over a dataset extended with structural feature data, without additional PPI annotations. The multi-task setup becomes even more important, if the fraction of PPI annotations becomes very small: the multi-task learner trained on only one-eighth of the PPI annotations—with data extension—reaches the same performances as the single-task learner on all PPI annotations. Thus, we show that the multi-task learning strategy can be beneficial for a small training dataset where the protein’s functional properties of interest are only partially annotated.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

In Silico Identification of Potential Inhibitors of the SARS-CoV-2 Nucleocapsid Through Molecular Docking-Based Drug Repurposing

SARS-CoV-2 is the virus responsible for the COVID-19 pandemic, and its effects on people worldwide continue to grow. Protein-targeted therapeutics are currently unavailable for this virus. As with other coronaviruses, the nucleocapsid (N) protein is the most conserved RNA-binding structural protein of SARS-CoV-2. The N protein is an appealing target because of its functional role in viral transcription and replication. Therefore, molecular docking method for structure-based drug design was used to investigate the binding energy and binding modes of various anti-N inhibitors in depth. The inhibitors selected were originally developed to target stress granules and other molecules involved in RNA biology, and were either FDA-approved or in the process of clinical trials for COVID-19. We aimed at targeting the N-terminal RNA binding domain (NTD) for molecular docking-based screening, on the basis of the first resolved crystal structure of SARS-CoV-2 N protein (PDB ID: 6M3M) and C-terminal domain (CTD) dimerization of the nucleocapsid phosphoprotein of SARS-COV-2 (PDB ID: 6WJI). Silmitasertib, nintedanib, ternatin, luteolin, and fedratinib were found to interact with RNA binding sites and to form a predicted protein interface with high binding energy. Similarly, silmitasertib, sirolimus-rapamycin, dovitinib, nintedanib, and fedratinib were found to interact with the SARS-CoV-2 N protein at its CTD dimerization sites, according to previous studies. In addition, we investigated an information gap regarding the relationships among the energetic landscape and stability and drug binding of the SARS-CoV-2 N NTD and CTD. Our in silico results clearly indicated that several tested drugs as potent putative inhibitors for COVID-19 therapeutics, thus indicating that they should be further validated as treatments to slow the spread of SARS-CoV-2.

Protein-protein interaction networks as miners of biological discovery

Protein-protein interactions (PPIs) form the basis of a myriad of biological pathways and mechanism, such as the formation of protein-complexes or the components of signaling cascades. Here, we reviewed experimental methods for identifying PPI pairs, including yeast two-hybrid, mass spectrometry, co-localization, and co-immunoprecipitation. Furthermore, a range of computational methods leveraging biochemical properties, evolution history, protein structures and more have enabled identification of additional PPIs. Given the wealth of known PPIs, we reviewed important network methods to construct and analyze networks of PPIs. These methods aid biological discovery through identifying hub genes and dynamic changes in the network, and have been thoroughly applied in various fields of biological research. Lastly, we discussed the challenges and future direction of research utilizing the power of PPI networks. This article is protected by copyright. All rights reserved.

EGRET: edge aggregated graph attention networks and transfer learning improve protein-protein interaction site prediction

MOTIVATION

Protein-protein interactions (PPIs) are central to most biological processes. However, reliable identification of PPI sites using conventional experimental methods is slow and expensive. Therefore, great efforts are being put into computational methods to identify PPI sites.

RESULTS

We present Edge Aggregated GRaph Attention NETwork (EGRET), a highly accurate deep learning-based method for PPI site prediction, where we have used an edge aggregated graph attention network to effectively leverage the structural information. We, for the first time, have used transfer learning in PPI site prediction. Our proposed edge aggregated network, together with transfer learning, has achieved notable improvement over the best alternate methods. Furthermore, we systematically investigated EGRET's network behavior to provide insights about the causes of its decisions.

AVAILABILITY

EGRET is freely available as an open source project at https://github.com/Sazan-Mahbub/EGRET.

CONTACT

shams_bayzid@cse.buet.ac.bd.

Predicting Diagnostic Potential of Cathepsin in Epithelial Ovarian Cancer: A Design Validated by Computational, Biophysical and Electrochemical Data

Background: Epithelial ovarian cancer remains one of the leading variants of gynecological cancer with a high mortality rate. Feasibility and technical competence for screening and detection of epithelial ovarian cancer remain a major obstacle and the development of point of care diagnostics (POCD) may offer a simple solution for monitoring its progression. Cathepsins have been implicated as biomarkers for cancer progression and metastasis; being a protease, it has an inherent tendency to interact with Cystatin C, a cysteine protease inhibitor. This interaction was assessed for designing a POCD module. Methods: A combinatorial approach encompassing computational, biophysical and electron-transfer kinetics has been used to assess this protease-inhibitor interaction. Results: Calculations predicted two cathepsin candidates, Cathepsin K and Cathepsin L based on their binding energies and structural alignment and both predictions were confirmed experimentally. Differential pulse voltammetry was used to verify the potency of Cathepsin K and Cathepsin L interaction with Cystatin C and assess the selectivity and sensitivity of their electrochemical interactions. Electrochemical measurements indicated selectivity for both the ligands, but with increasing concentrations, there was a marked difference in the sensitivity of the detection. Conclusions: This work validated the utility of dry-lab integration in the wet-lab technique to generate leads for the design of electrochemical diagnostics for epithelial ovarian cancer.

A Deep Learning and XGBoost-Based Method for Predicting Protein-Protein Interaction Sites

Knowledge about protein-protein interactions is beneficial in understanding cellular mechanisms. Protein-protein interactions are usually determined according to their protein-protein interaction sites. Due to the limitations of current techniques, it is still a challenging task to detect protein-protein interaction sites. In this article, we presented a method based on deep learning and XGBoost (called DeepPPISP-XGB) for predicting protein-protein interaction sites. The deep learning model served as a feature extractor to remove redundant information from protein sequences. The Extreme Gradient Boosting algorithm was used to construct a classifier for predicting protein-protein interaction sites. The DeepPPISP-XGB achieved the following results: area under the receiver operating characteristic curve of 0.681, a recall of 0.624, and area under the precision-recall curve of 0.339, being competitive with the state-of-the-art methods. We also validated the positive role of global features in predicting protein-protein interaction sites.

CATH functional families predict functional sites in proteins

MOTIVATION

Identification of functional sites in proteins is essential for functional characterization, variant interpretation and drug design. Several methods are available for predicting either a generic functional site, or specific types of functional site. Here, we present FunSite, a machine learning predictor that identifies catalytic, ligand-binding and protein-protein interaction functional sites using features derived from protein sequence and structure, and evolutionary data from CATH functional families (FunFams).

RESULTS

FunSite's prediction performance was rigorously benchmarked using cross-validation and a holdout dataset. FunSite outperformed other publicly available functional site prediction methods. We show that conserved residues in FunFams are enriched in functional sites. We found FunSite's performance depends greatly on the quality of functional site annotations and the information content of FunFams in the training data. Finally, we analyze which structural and evolutionary features are most predictive for functional sites.

AVAILABILITYAND IMPLEMENTATION

https://github.com/UCL/cath-funsite-predictor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Identification and characterization of Poly(ADP-ribose) polymerase-1 interacting proteins during development of Dictyostelium discoideum

Poly (ADP-ribose) polymerase-1 (PARP-1) is a multifunctional protein that is associated with various biological processes like chromatin remodeling, DNA damage, cell death etc. In Dictyostelium discoideum, PARP-1 has also been implicated in cellular differentiation and development. However, its interacting proteins during multicellular development are not yet explored. Hence, the present study aims to identify PARP-1 interacting proteins during multicellular development of D. discoideum. BRCA1 C-terminus (BRCT) domain of PARP-1, which is mainly involved in protein-protein interactions was cloned in pGEX4T1 vector and developmental interactome of PARP-1 were analyzed by affinity purification-mass spectrometry. These interactions were further confirmed by in-silico protein-protein docking analysis, which led to identification of the proteins that show high affinity for BRCT domain. Initially, the protein structures were modeled on SWISS MODEL and PHYRE2 servers, refined by 3Drefine and validated by PROCHECK. Further, interaction sites of BRCT and the conserved regions in all interacting proteins were predicted using cons-PPISP and ConSurf, respectively. Finally, protein-protein docking analysis was done by HADDOCK. Our results identified 19 possible BRCT interacting proteins during D. discoideum development. Furthermore, interacting residues involved in the interactions and functional regions were explored. This is the first report where PARP-1's developmental interactome in D. discoideum is well established. The current findings demonstrate PARP-1's developmental interactome in D. discoideum and provide the groundwork to understand its regulated functions in developmental biology which would undoubtedly extend our perception towards developmental diseases in higher complex organisms and their treatment.

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.

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.

Deep Learning for Protein-Protein Interaction Site Prediction

Protein-protein interactions (PPIs) are central to cellular functions. Experimental methods for predicting PPIs are well developed but are time and resource expensive and suffer from high false-positive error rates at scale. Computational prediction of PPIs is highly desirable for a mechanistic understanding of cellular processes and offers the potential to identify highly selective drug targets. In this chapter, details of developing a deep learning approach to predicting which residues in a protein are involved in forming a PPI-a task known as PPI site prediction-are outlined. The key decisions to be made in defining a supervised machine learning project in this domain are here highlighted. Alternative training regimes for deep learning models to address shortcomings in existing approaches and provide starting points for further research are discussed. This chapter is written to serve as a companion to developing deep learning approaches to protein-protein interaction site prediction, and an introduction to developing geometric deep learning projects operating on protein structure graphs.

Antitoxin autoregulation of M. tuberculosis toxin-antitoxin expression through negative cooperativity arising from multiple inverted repeat sequences

Toxin-antitoxin systems play key roles in bacterial adaptation, including protection from antibiotic assault and infection by bacteriophages. The type IV toxin-antitoxin system AbiE encodes a DUF1814 nucleotidyltransferase-like toxin, and a two-domain antitoxin. In Streptococcus agalactiae, the antitoxin AbiEi negatively autoregulates abiE expression through positively co-operative binding to inverted repeats within the promoter. The human pathogen Mycobacterium tuberculosis encodes four DUF1814 putative toxins, two of which have antitoxins homologous to AbiEi. One such M. tuberculosis antitoxin, named Rv2827c, is required for growth and whilst the structure has previously been solved, the mode of regulation is unknown. To complete the gaps in our understanding, we first solved the structure of S. agalactiae AbiEi to 1.83 Å resolution for comparison with M. tuberculosis Rv2827c. AbiEi contains an N-terminal DNA binding domain and C-terminal antitoxicity domain, with bilateral faces of opposing charge. The overall AbiEi fold is similar to Rv2827c, though smaller, and with a 65° difference in C-terminal domain orientation. We further demonstrate that, like AbiEi, Rv2827c can autoregulate toxin-antitoxin operon expression. In contrast with AbiEi, the Prv2827c promoter contains two sets of inverted repeats, which bind Rv2827c with differing affinities depending on the sequence consensus. Surprisingly, Rv2827c bound with negative co-operativity to the full Prv2827c promoter, demonstrating an unexpectedly complex form of transcriptional regulation.

Targeting SARS-CoV-2 nucleocapsid oligomerization: Insights from molecular docking and molecular dynamics simulations

The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.

Structure of fish Toll-like receptors (TLR) and NOD-like receptors (NLR)

Innate immunity driven by pattern recognition receptor (PRR) protects the host from invading pathogens. Aquatic animals like fish where the adaptive immunity is poorly developed majorly rely on their innate immunity modulated by PRRs like toll-like receptors (TLR) and NOD-like receptors (NLR). However, current development to improve the fish immunity via TLR/NLR signaling is affected by a poor understanding of its mechanistic and structural features. This review discusses the structure of fish TLRs/NLRs and its interaction with pathogen associated molecular patterns (PAMPs) and downstream signaling molecules. Over the past one decade, significant progress has been done in studying the structure of TLRs/NLRs in higher eukaryotes; however, structural studies on fish innate immune receptors are undermined. Several novel TLR genes are identified in fish that are absent in higher eukaryotes, but the function is still poorly understood. Unlike the fundamental progress achieved in developing antagonist/agonist to modulate human innate immunity, analogous studies in fish are nearly lacking due to structural inadequacy. This underlies the importance of exploring the structural and mechanistic details of fish TLRs/NLRs at an atomic and molecular level. This review outlined the mechanistic and structural basis of fish TLR and NLR activation.

Protein-protein interaction site prediction through combining local and global features with deep neural networks

MOTIVATION

Protein-protein interactions (PPIs) play important roles in many biological processes. Conventional biological experiments for identifying PPI sites are costly and time-consuming. Thus, many computational approaches have been proposed to predict PPI sites. Existing computational methods usually use local contextual features to predict PPI sites. Actually, global features of protein sequences are critical for PPI site prediction.

RESULTS

A new end-to-end deep learning framework, named DeepPPISP, through combining local contextual and global sequence features, is proposed for PPI site prediction. For local contextual features, we use a sliding window to capture features of neighbors of a target amino acid as in previous studies. For global sequence features, a text convolutional neural network is applied to extract features from the whole protein sequence. Then the local contextual and global sequence features are combined to predict PPI sites. By integrating local contextual and global sequence features, DeepPPISP achieves the state-of-the-art performance, which is better than the other competing methods. In order to investigate if global sequence features are helpful in our deep learning model, we remove or change some components in DeepPPISP. Detailed analyses show that global sequence features play important roles in DeepPPISP.

AVAILABILITY AND IMPLEMENTATION

The DeepPPISP web server is available at http://bioinformatics.csu.edu.cn/PPISP/. The source code can be obtained from https://github.com/CSUBioGroup/DeepPPISP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein–Protein Interaction Methods and Protein Phase Separation

In the last decade, newly developed experimental methods have made it possible to highlight that macromolecules in the cell milieu physically interact to support physiology. This has shifted the problem of protein–protein interaction from a microscopic, electron-density scale to a mesoscopic one. Further, nowadays there is increasing evidence that proteins in the nucleus and in the cytoplasm can aggregate in membraneless organelles for different physiological reasons. In this scenario, it is urgent to face the problem of biomolecule functional annotation with efficient computational methods, suited to extract knowledge from reliable data and transfer information across different domains of investigation. Here, we revise the present state of the art of our knowledge of protein–protein interaction and the computational methods that differently implement it. Furthermore, we explore experimental and computational features of a set of proteins involved in phase separation.

Genome-wide and structural analyses of pseudokinases encoded in the genome of Arabidopsis thaliana provide functional insights

Protein Kinase-Like Non-Kinases (PKLNKs), commonly known as "pseudokinases", are homologous to eukaryotic Ser/Thr/Tyr protein kinases (PKs) but lack the crucial aspartate residue in the catalytic loop, indispensable for phosphotransferase activity. Therefore, they are predicted to be "catalytically inactive" enzyme homologs. Analysis of protein-kinase like sequences from Arabidopsis thaliana led to the identification of more than 120 pseudokinases lacking catalytic aspartate, majority of which are closely related to the plant-specific receptor-like kinase family. These pseudokinases engage in different biological processes, enabled by their diverse domain architectures and specific subcellular localizations. Structural comparison of pseudokinases with active and inactive conformations of canonical PKs, belonging to both plant and animal origin, revealed unique structural differences. The currently available crystal structures of pseudokinases show that the loop topologically equivalent to activation segment of PKs adopts a distinct-folded conformation, packing against the pseudoenzyme core, in contrast to the extended and inhibitory geometries observed for active and inactive states, respectively, of catalytic PKs. Salt-bridge between ATP-binding Lys and DFG-Asp as well as hydrophobic interactions between the conserved nonpolar residue C-terminal to the equivalent DFG motif and nonpolar residues in C-helix mediate such a conformation in pseudokinases. This results in enhanced solvent accessibility of the pseudocatalytic loop in pseudokinases that can possibly serve as an interacting surface while associating with other proteins. Specifically, our analysis identified several residues that may be involved in pseudokinase regulation and hints at the repurposing of pseudocatalytic residues to achieve mechanistic control over noncatalytic functions of pseudoenzymes.

In-silico analysis to identify the role of MEN1 missense mutations in breast cancer

Background: The multiple endocrine neoplasia type 1 (MEN1) gene located on chromosome 11q13 encodes menin protein. Previously reported mutations were thought to result in loss of function of menin protein and that they are associated with multiple endocrine neoplasia 1 disorder. However, recently menin has also been characterized as an oncosuppressor protein and it was suggested that mutations in it are associated with various other tumors. Studies indicate that the menin protein stimulates the estrogen receptor (ER) that in turn increases the predisposition for inherited breast cancer.

Methods: Here, we used our supervised in-house combinatory in-silico predictor method to investigate the impact of unclassified missense mutations in MEN1 gene found in breast cancer tissue. We also examined the biophysical and biochemical properties to predict the effects of these missense variants on the menin protein stability and interactions. The results are compared with the effects of known pathogenic mutations in menin causing neoplasia.

Results: Our analysis indicates that some of the variants found in breast cancer tissue show similar pattern of destabilizing the menin protein and its interactions as the pathogenic variants associated with neoplasia. Taking together with the results of our in-silico consensus predictor, we classify missense mutations in menin protein found in breast cancer tissue into pathogenic and benign, and thus, suggesting as an indicator for early detection of elevated breast cancer risk.

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.

ConvsPPIS: Identifying Protein-protein Interaction Sites by an Ensemble Convolutional Neural Network with Feature Graph

Background/Objective:

Protein-protein interactions are essentials for most cellular processes and thus, unveiling how proteins interact with is a crucial question that can be better understood by recognizing which residues participate in the interaction. Although many computational approaches have been proposed to predict interface residues, their feature perspective and model learning ability are not enough to achieve ideal results. So, our objective is to improve the predictive performance under considering feature perspective and new learning algorithm.

Method:

In this study, we proposed an ensemble deep convolutional neural network, which explores the context and positional context of consecutive residues within a protein sub-sequence. Specifically, unlike the feature view of previous methods, ConvsPPIS uses evolutionary, physicochemical, and structural protein characteristics to construct their own feature graph respectively. After that, three independent deep convolutional neural networks are trained on each type of feature graph for learning the underlying pattern in sub-sequence. Lastly, we integrated those three deep networks into an ensemble predictor with leveraging complementary information of those features to predict potential interface residues.

Results:

Some comparative experiments have conducted through 10-fold cross-validation. The results indicated that ConvsPPIS achieved superior performance on DBv5-Sel dataset with an accuracy of 88%. Additional experiments on CAPRI-Alone dataset demonstrated ConvsPPIS has also better prediction performance.

Conclusion:

The ConvsPPIS method provided a new perspective to capture protein feature expression for identifying protein-protein interaction sites. The results proved the superiority of this method.

A Novel Index of Contact Frequency from Noise Protein-Protein Interaction Data Help for Accurate Interface Residue Pair Prediction

Protein-protein interactions are important for most biological processes and have been studied for decades. However, the detailed formation mechanism of protein-protein interaction interface is still ambiguous, which makes it difficult to accurately predict the protein-protein interaction interface residue pairs. Here, we extract the interface residue-residue contacts from the decoys in the ZDOCK protein-protein complex decoy set with RMSD mostly larger than 3 Å. To accurately compute the interface residue-residue contacts, we define a new constant called interface residue pairs frequency, which counts the atom contact numbers between two interface residues. We normalize interface residue pairs frequency to pick out the top residue-residue pairs from all the possible pairs preferential to be on correct protein-protein interaction interface. When tested on 37 protein dimers from the decoy set where most decoys are incorrect, our method successfully predicts 30 protein dimers with a success rate of up to 81.1%. Higher accuracy than some other state-of-the-art methods confirmed the performance of our method.

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.

Analysis of protein missense alterations by combining sequence- and structure-based methods

BACKGROUND

Different types of in silico approaches can be used to predict the phenotypic consequence of missense variants. Such algorithms are often categorized as sequence based or structure based, when they necessitate 3D structural information. In addition, many other in silico tools, not dedicated to the analysis of variants, can be used to gain additional insights about the possible mechanisms at play.

METHODS

Here we applied different computational approaches to a set of 20 known missense variants present on different proteins (CYP, complement factor B, antithrombin and blood coagulation factor VIII). The tools that were used include fast computational approaches and web servers such as PolyPhen-2, PopMusic, DUET, MaestroWeb, SAAFEC, Missense3D, VarSite, FlexPred, PredyFlexy, Clustal Omega, meta-PPISP, FTMap, ClusPro, pyDock, PPM, RING, Cytoscape, and ChannelsDB.

RESULTS

We observe some conflicting results among the methods but, most of the time, the combination of several engines helped to clarify the potential impacts of the amino acid substitutions.

CONCLUSION

Combining different computational approaches including some that were not developed to investigate missense variants help to predict the possible impact of the amino acid substitutions. Yet, when the modified residues are involved in a salt-bridge, the tools tend to fail, even when the analysis is performed in 3D. Thus, interactive structural analysis with molecular graphics packages such as Chimera or PyMol or others are still needed to clarify automatic prediction.

Prediction of Protein-Protein Interaction Sites Using Convolutional Neural Network and Improved Data Sets

Protein-protein interaction (PPI) sites play a key role in the formation of protein complexes, which is the basis of a variety of biological processes. Experimental methods to solve PPI sites are expensive and time-consuming, which has led to the development of different kinds of prediction algorithms. We propose a convolutional neural network for PPI site prediction and use residue binding propensity to improve the positive samples. Our method obtains a remarkable result of the area under the curve (AUC) = 0.912 on the improved data set. In addition, it yields much better results on samples with high binding propensity than on randomly selected samples. This suggests that there are considerable false-positive PPI sites in the positive samples defined by the distance between residue atoms.

The Cap-Snatching SFTSV Endonuclease Domain Is an Antiviral Target

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a tick-borne virus with 12%-30% case mortality rates and is related to the Heartland virus (HRTV) identified in the United States. Together, SFTSV and HRTV are emerging segmented, negative-sense RNA viral (sNSV) pathogens with potential global health impact. Here, we characterize the amino-terminal cap-snatching endonuclease domain of SFTSV polymerase (L) and solve a 2.4-Å X-ray crystal structure. While the overall structure is similar to those of other cap-snatching sNSV endonucleases, differences near the C terminus of the SFTSV endonuclease suggest divergence in regulation. Influenza virus endonuclease inhibitors, including the US Food and Drug Administration (FDA) approved Baloxavir (BXA), inhibit the endonuclease activity in in vitro enzymatic assays and in cell-based studies. BXA displays potent activity with a half maximal inhibitory concentration (IC50) of ∼100 nM in enzyme inhibition and an EC50 value of ∼250 nM against SFTSV and HRTV in plaque assays. Together, our data support sNSV endonucleases as an antiviral target.

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.

Identification and visualization of protein binding regions with the ArDock server

ArDock (ardock.ibcp.fr) is a structural bioinformatics web server for the prediction and the visualization of potential interaction regions at protein surfaces. ArDock ranks the surface residues of a protein according to their tendency to form interfaces in a set of predefined docking experiments between the query protein and a set of arbitrary protein probes. The ArDock methodology is derived from large scale cross-docking studies where it was observed that randomly chosen proteins tend to dock in a non-random way at protein surfaces. The method predicts interaction site of the protein, or alternate interfaces in the case of proteins with multiple interaction modes. The server takes a protein structure as input and computes a score for each surface residue. Its output focuses on the interactive visualization of results and on interoperability with other services.

Protein-protein interaction sites prediction by ensemble random forests with synthetic minority oversampling technique

MOTIVATION

The prediction of protein-protein interaction (PPI) sites is a key to mutation design, catalytic reaction and the reconstruction of PPI networks. It is a challenging task considering the significant abundant sequences and the imbalance issue in samples.

RESULTS

A new ensemble learning-based method, Ensemble Learning of synthetic minority oversampling technique (SMOTE) for Unbalancing samples and RF algorithm (EL-SMURF), was proposed for PPI sites prediction in this study. The sequence profile feature and the residue evolution rates were combined for feature extraction of neighboring residues using a sliding window, and the SMOTE was applied to oversample interface residues in the feature space for the imbalance problem. The Multi-dimensional Scaling feature selection method was implemented to reduce feature redundancy and subset selection. Finally, the Random Forest classifiers were applied to build the ensemble learning model, and the optimal feature vectors were inserted into EL-SMURF to predict PPI sites. The performance validation of EL-SMURF on two independent validation datasets showed 77.1% and 77.7% accuracy, which were 6.2-15.7% and 6.1-18.9% higher than the other existing tools, respectively.

AVAILABILITY AND IMPLEMENTATION

The source codes and data used in this study are publicly available at http://github.com/QUST-AIBBDRC/EL-SMURF/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gain of native conformation of Aurora A S155R mutant by small molecules

Aurora A is a mitotic serine/threonine kinase protein that is a proposed target of the first-line anticancer drug design. It has been found to be overexpressed in many human cancer cells, including hematological, breast, and colorectal. Here, we focus on a particular somatic mutant S155R of Aurora kinase A protein, whose activity decreases because of loss of interaction with a TPX2 protein that results in ectopic expression of the Aurora kinase A protein, which contributes chromosome instability, centrosome amplification, and oncogenic transformation. The primary target of this study is to select a drug molecule whose binding results in gaining S155R mutant interaction with TPX2. The computational methodology applied in this study involves mapping of hotspots (for uncompetitive binding), virtual screening, protein-ligand docking, postdocking optimization, and protein-protein docking approach. In this study, we screen and validate ZINC968264, which acts as a potential molecule that can improve the loss of function occurred because of mutation (S155R) in Aurora A. Our approaches pave a suitable path to design a potential drug against physiological condition manifested because of S155R mutant in Aurora A.

Predicting Protein-Protein Interfaces that Bind Intrinsically Disordered Protein Regions

A long-standing goal in biology is the complete annotation of function and structure on all protein-protein interactions, a large fraction of which is mediated by intrinsically disordered protein regions (IDRs). However, knowledge derived from experimental structures of such protein complexes is disproportionately small due, in part, to challenges in studying interactions of IDRs. Here, we introduce IDRBind, a computational method that by combining gradient boosted trees and conditional random field models predicts binding sites of IDRs with performance approaching state-of-the-art globular interface predictions, making it suitable for proteome-wide applications. Although designed and trained with a focus on molecular recognition features, which are long interaction-mediating-elements in IDRs, IDRBind also predicts the binding sites of short peptides more accurately than existing specialized predictors. Consistent with IDRBind's specificity, a comparison of protein interface categories uncovered uniform trends in multiple physicochemical properties, positioning molecular recognition feature interfaces between peptide and globular interfaces.

The Ferredoxin-Like Protein FerR Regulates PrbP Activity in Liberibacter asiaticus

In Liberibacter asiaticus, PrbP is an important transcriptional accessory protein that regulates gene expression through interactions with the RNA polymerase β-subunit and a specific sequence on the promoter region. The constitutive expression of prbP observed upon chemical inactivation of PrbP-DNA interactions in vivo indicated that the expression of prbP was not autoregulated at the level of transcription. This observation suggested that a modulatory mechanism via protein-protein interactions may be involved. In silico genome association analysis identified FerR (CLIBASIA_01505), a putative ferredoxin-like protein, as a PrbP-interacting protein. Using a bacterial two-hybrid system and immunoprecipitation assays, interactions between PrbP and FerR were confirmed. In vitro transcription assays were used to show that FerR can increase the activity of PrbP by 16-fold when present in the PrbP-RNA polymerase reaction mixture. The FerR protein-protein interaction surface was predicted by structural modeling and followed by site-directed mutagenesis. Amino acids V20, V23, and C40 were identified as the most important residues in FerR involved in the modulation of PrbP activity in vitro The regulatory mechanism of FerR abundance was examined at the transcription level. In contrast to prbP of L. asiaticus (prbP _Las), mRNA levels of ferR of L. asiaticus (ferR _Las) are induced by an increase in osmotic pressure. The results of this study revealed that the activity of the transcriptional activator PrbP_Las is modulated via interactions with FerR_Las The induction of ferR _Las expression by osmolarity provides insight into the mechanisms of adjusting gene expression in response to host environmental signals in L. asiaticus IMPORTANCE The rapid spread and aggressive progression of huanglongbing (HLB) in the major citrus-producing areas have raised global recognition of and vigilance to this disease. As a result, the causative agent, Liberibacter asiaticus, has been investigated from various perspectives. However, gene expression regulatory mechanisms that are important for the survival and persistence of this intracellular pathogen remain largely unexplored. PrbP is a transcriptional accessory protein important for L. asiaticus survival in the plant host. In this study, we investigated the interactions between PrbP in L. asiaticus (PrbP_Las) and a ferredoxin-like protein (FerR) in L. asiaticus, FerR_Las We show that the presence of FerR stabilizes and augments the activity of PrbP_Las In addition, we demonstrate that the expression of ferR is induced by increases in osmolarity in Liberibacter crescens Altogether, these results suggest that FerR_Las and PrbP_Las may play important roles in the regulation of gene expression in response to changing environmental signals during L. asiaticus infection in the citrus host.

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.