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

Predicting human and viral protein variants affecting COVID-19 susceptibility and repurposing therapeutics

The COVID-19 disease is an ongoing global health concern. Although vaccination provides some protection, people are still susceptible to re-infection. Ostensibly, certain populations or clinical groups may be more vulnerable. Factors causing these differences are unclear and whilst socioeconomic and cultural differences are likely to be important, human genetic factors could influence susceptibility. Experimental studies indicate SARS-CoV-2 uses innate immune suppression as a strategy to speed-up entry and replication into the host cell. Therefore, it is necessary to understand the impact of variants in immunity-associated human proteins on susceptibility to COVID-19. In this work, we analysed missense coding variants in several SARS-CoV-2 proteins and their human protein interactors that could enhance binding affinity to SARS-CoV-2. We curated a dataset of 19 SARS-CoV-2: human protein 3D-complexes, from the experimentally determined structures in the Protein Data Bank and models built using AlphaFold2-multimer, and analysed the impact of missense variants occurring in the protein-protein interface region. We analysed 468 missense variants from human proteins and 212 variants from SARS-CoV-2 proteins and computationally predicted their impacts on binding affinities for the human viral protein complexes. We predicted a total of 26 affinity-enhancing variants from 13 human proteins implicated in increased binding affinity to SARS-CoV-2. These include key-immunity associated genes (TOMM70, ISG15, IFIH1, IFIT2, RPS3, PALS1, NUP98, AXL, ARF6, TRIMM, TRIM25) as well as important spike receptors (KREMEN1, AXL and ACE2). We report both common (e.g., Y13N in IFIH1) and rare variants in these proteins and discuss their likely structural and functional impact, using information on known and predicted functional sites. Potential mechanisms associated with immune suppression implicated by these variants are discussed. Occurrence of certain predicted affinity-enhancing variants should be monitored as they could lead to increased susceptibility and reduced immune response to SARS-CoV-2 infection in individuals/populations carrying them. Our analyses aid in understanding the potential impact of genetic variation in immunity-associated proteins on COVID-19 susceptibility and help guide drug-repurposing strategies.

Assessing the functional impact of protein binding site definition

Many biomedical applications, such as classification of binding specificities or bioengineering, depend on the accurate definition of protein binding interfaces. Depending on the choice of method used, substantially different sets of residues can be classified as belonging to the interface of a protein. A typical approach used to verify these definitions is to mutate residues and measure the impact of these changes on binding. Besides the lack of exhaustive data, this approach also suffers from the fundamental problem that a mutation introduces an unknown amount of alteration into an interface, which potentially alters the binding characteristics of the interface. In this study we explore the impact of alternative binding site definitions on the ability of a protein to recognize its cognate ligand using a pharmacophore approach, which does not affect the interface. The study also shows that methods for protein binding interface predictions should perform above approximately F-score = 0.7 accuracy level to capture the biological function of a protein.

SurfPro-NN: A 3D point cloud neural network for the scoring of protein-protein docking models based on surfaces features and protein language models

Protein-protein interactions (PPI) play a crucial role in numerous key biological processes, and the structure of protein complexes provides valuable clues for in-depth exploration of molecular-level biological processes. Protein-protein docking technology is widely used to simulate the spatial structure of proteins. However, there are still challenges in selecting candidate decoys that closely resemble the native structure from protein-protein docking simulations. In this study, we introduce a docking evaluation method based on three-dimensional point cloud neural networks named SurfPro-NN, which represents protein structures as point clouds and learns interaction information from protein interfaces by applying a point cloud neural network. With the continuous advancement of deep learning in the field of biology, a series of knowledge-rich pre-trained models have emerged. We incorporate protein surface representation models and language models into our approach, greatly enhancing feature representation capabilities and achieving superior performance in protein docking model scoring tasks. Through comprehensive testing on public datasets, we find that our method outperforms state-of-the-art deep learning approaches in protein-protein docking model scoring. Not only does it significantly improve performance, but it also greatly accelerates training speed. This study demonstrates the potential of our approach in addressing protein interaction assessment problems, providing strong support for future research and applications in the field of biology.

Rosmarinic acid alleviates septic acute respiratory distress syndrome in mice by suppressing the bronchial epithelial RAS-mediated ferroptosis

Activating angiotensin-converting enzyme 2 (ACE2) is an important player in the pathogenesis of septic-related acute respiratory distress syndrome (ARDS). Rosmarinic acid (RA) as a prominent polyphenolic secondary metabolite derived from Rosmarinus officinalis modulates ACE2 in sepsis remains unclear, although its impact on ACE inhibition and septic-associated lung injury has been explored. The study investigated the ACE2 expression in lipopolysaccharide (LPS)-induced lungs in mice and BEAS2B cells. Additionally, molecular docking, protein-protein interaction (PPI) network analysis, and western blotting were employed to predict and evaluate the molecular mechanism of RA on LPS-induced ferroptosis in vivo and in vitro. LPS-induced glutathione peroxidase 4 (GPX4) downregulation, ACE/ACE2 imbalance, and alteration of frequency of breathing (BPM), minute volume (MV), and the expiratory flow at 50% expired volume (EF50) were reversed by captopril pretreatment in vitro and in vivo. RA notably inhibited the infiltration into the lungs of neutrophils and monocytes with increased amounts of GPX4 and ACE2 proteins, lung function improvement, and decreased inflammatory cytokines levels and ER stress in LPS-induced ARDS in mice. Molecular docking showed RA was able to interact with ACE and ACE2. Moreover, combined with different pharmacological inhibitors to block ACE and ferroptosis, RA still significantly inhibited inflammatory cytokines Interleukin-1ß (IL-1ß), tumor necrosis factor-α (TNF-α), and C-X-C motif chemokine 2 (CXCL2) levels, as well as improved lung function, and enhanced GPX4 expression. Particularly, the anti-ferroptosis effect of RA in LPS-induced septic ARDS is RAS-dependent.

The Role of AI in Drug Discovery

The emergence of Artificial Intelligence (AI) in drug discovery marks a pivotal shift in pharmaceutical research, blending sophisticated computational techniques with conventional scientific exploration to break through enduring obstacles. This review paper elucidates the multifaceted applications of AI across various stages of drug development, highlighting significant advancements and methodologies. It delves into AI's instrumental role in drug design, polypharmacology, chemical synthesis, drug repurposing, and the prediction of drug properties such as toxicity, bioactivity, and physicochemical characteristics. Despite AI's promising advancements, the paper also addresses the challenges and limitations encountered in the field, including data quality, generalizability, computational demands, and ethical considerations. By offering a comprehensive overview of AI's role in drug discovery, this paper underscores the technology's potential to significantly enhance drug development, while also acknowledging the hurdles that must be overcome to fully realize its benefits.

Database-Driven Identification of Structurally Similar Protein-Protein Interfaces

Analyzing the similarity of protein interfaces in protein-protein interactions gives new insights into protein function and assists in discovering new drugs. Usually, tools that assess the similarity focus on the interactions between two protein interfaces, while sometimes we only have one predicted interface. Herein, we present PiMine, a database-driven protein interface similarity search. It compares interface residues of one or two interacting chains by calculating and searching tetrahedral geometric patterns of α-carbon atoms and calculating physicochemical and shape-based similarity. On a dedicated, tailor-made dataset, we show that PiMine outperforms commonly used comparison tools in terms of early enrichment when considering interfaces of sequentially and structurally unrelated proteins. In an application example, we demonstrate its usability for protein interaction partner prediction by comparing predicted interfaces to known protein-protein interfaces.

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

Tribuloside: Mechanisms and Efficacy in Treating Acute Lung Injury Revealed by Network Pharmacology and Experimental Validation

Background

Acute lung injury (ALI) is a serious illness that has few treatment options available. Tribuloside, a natural flavonoid extracted from the Tribulus Terrestris plant in China, is potent in addressing many health issues such as headaches, dizziness, itching, and vitiligo.

Objective

This study intends to explore the mechanisms of action of Tribuloside in treating ALI through a combination of network pharmacology and experimental validation.

Methods

We obtained the 2D structure and SMILES number of Tribuloside from the PubChem database. We used the SwissTargetPrediction database to identify pharmacological targets. We found 1215 targets linked to ALI by examining the GeneCards database. We used the String database and Cytoscape software to create the "drug or disease-target" network as well as the protein-protein interactions (PPI). Key targets were identified by evaluating associated biological processes and pathway enrichment. A Venny Diagram showed 49 intersection points between Tribuloside and ALI. Molecular docking with AutoDockTools found that Tribuloside had a high affinity for IL6, BCL2, TNF, STAT3, IL1B, and MAPK3, the top 6 targets in the PPI network by Degree values. To test Tribuloside's therapeutic efficacy in ALI, an acute lung damage model in mice was constructed using lipopolysaccharide. Tribuloside treatment reduced inflammatory cell infiltration, decreased fibrotic area, repaired damaged alveoli, and suppressed inflammatory factors IL-6, TNF-α, and IL-1β in the lungs through many pathways and targets.

Conclusion

This study reveals that Tribuloside has the potential to treat ALI by targeting various pathways and targets, according to network pharmacology predictions and experimental confirmation.

Structurally-informed human interactome reveals proteome-wide perturbations by disease mutations

Human genome sequencing studies have identified numerous loci associated with complex diseases. However, translating human genetic and genomic findings to disease pathobiology and therapeutic discovery remains a major challenge at multiscale interactome network levels. Here, we present a deep-learning-based ensemble framework, termed PIONEER (Protein-protein InteractiOn iNtErfacE pRediction), that accurately predicts protein binding partner-specific interfaces for all known protein interactions in humans and seven other common model organisms, generating comprehensive structurally-informed protein interactomes. We demonstrate that PIONEER outperforms existing state-of-the-art methods. We further systematically validated PIONEER predictions experimentally through generating 2,395 mutations and testing their impact on 6,754 mutation-interaction pairs, confirming the high quality and validity of PIONEER predictions. We show that disease-associated mutations are enriched in PIONEER-predicted protein-protein interfaces after mapping mutations from ~60,000 germline exomes and ~36,000 somatic genomes. We identify 586 significant protein-protein interactions (PPIs) enriched with PIONEER-predicted interface somatic mutations (termed oncoPPIs) from pan-cancer analysis of ~11,000 tumor whole-exomes across 33 cancer types. We show that PIONEER-predicted oncoPPIs are significantly associated with patient survival and drug responses from both cancer cell lines and patient-derived xenograft mouse models. We identify a landscape of PPI-perturbing tumor alleles upon ubiquitination by E3 ligases, and we experimentally validate the tumorigenic KEAP1-NRF2 interface mutation p.Thr80Lys in non-small cell lung cancer. We show that PIONEER-predicted PPI-perturbing alleles alter protein abundance and correlates with drug responses and patient survival in colon and uterine cancers as demonstrated by proteogenomic data from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium. PIONEER, implemented as both a web server platform and a software package, identifies functional consequences of disease-associated alleles and offers a deep learning tool for precision medicine at multiscale interactome network levels.

ARCTIC-3D: automatic retrieval and clustering of interfaces in complexes from 3D structural information

The formation of a stable complex between proteins lies at the core of a wide variety of biological processes and has been the focus of countless experiments. The huge amount of information contained in the protein structural interactome in the Protein Data Bank can now be used to characterise and classify the existing biological interfaces. We here introduce ARCTIC-3D, a fast and user-friendly data mining and clustering software to retrieve data and rationalise the interface information associated with the protein input data. We demonstrate its use by various examples ranging from showing the increased interaction complexity of eukaryotic proteins, 20% of which on average have more than 3 different interfaces compared to only 10% for prokaryotes, to associating different functions to different interfaces. In the context of modelling biomolecular assemblies, we introduce the concept of "recognition entropy", related to the number of possible interfaces of the components of a protein-protein complex, which we demonstrate to correlate with the modelling difficulty in classical docking approaches. The identified interface clusters can also be used to generate various combinations of interface-specific restraints for integrative modelling. The ARCTIC-3D software is freely available at github.com/haddocking/arctic3d and can be accessed as a web-service at wenmr.science.uu.nl/arctic3d.

DeepRank-GNN-esm: a graph neural network for scoring protein-protein models using protein language model

Motivation

Protein-Protein interactions (PPIs) play critical roles in numerous cellular processes. By modelling the 3D structures of the correspond protein complexes valuable insights can be obtained, providing, e.g. starting points for drug and protein design. One challenge in the modelling process is however the identification of near-native models from the large pool of generated models. To this end we have previously developed DeepRank-GNN, a graph neural network that integrates structural and sequence information to enable effective pattern learning at PPI interfaces. Its main features are related to the Position Specific Scoring Matrices (PSSMs), which are computationally expensive to generate, significantly limits the algorithm's usability.

Results

We introduce here DeepRank-GNN-esm that includes as additional features protein language model embeddings from the ESM-2 model. We show that the ESM-2 embeddings can actually replace the PSSM features at no cost in-, or even better performance on two PPI-related tasks: scoring docking poses and detecting crystal artifacts. This new DeepRank version bypasses thus the need of generating PSSM, greatly improving the usability of the software and opening new application opportunities for systems for which PSSM profiles cannot be obtained or are irrelevant (e.g. antibody-antigen complexes).

Availability and implementation

DeepRank-GNN-esm is freely available from https://github.com/DeepRank/DeepRank-GNN-esm.

Exploration of the Potential Mechanism of Yujin Powder Treating Dampness-heat Diarrhea by Integrating UPLC-MS/MS and Network Pharmacology Prediction

BACKGROUND

Yujin powder (YJP) is a classic prescription for treating dampness-heat diarrhea (DHD) in Traditional Chinese Medicine (TCM), but the main functional active ingredients and the exact mechanisms have not been systematically studied.

OBJECTIVES

This study aimed to preliminarily explore the potential mechanisms of YJP for treating DHD by integrating UPLC-MS/MS and network pharmacology methods.

METHODS

Ultra performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) technology was used to determine the ingredients of YJP. And then, the targets of these components were predicted and screened from TCMSP, SwissTargetPrediction databases. The disease targets related to DHD were obtained by using the databases of GeneCards, OMIM, DisGeNET, TTD, and DrugBank. The protein-protein interaction networks (PPI) of YJP-DHD were constructed using the STRING database and Origin 2022 software to identify the cross-targets by screening the core-acting targets and a network diagram by Cytoscape 3.8.2 software was also constructed. Metascape database was used for performing GO and KEGG enrichment anlysis on the core genes. Finally, molecular docking was used to verify the results with AutoDock 4.2.6, AutoDock Tools 1.5.6, PyMOL 2.4.0, and Open Babel 2.3.2 software.

RESULTS

597 components in YJP were detected, and 153 active components were obtained through database screening, among them the key active ingredients include coptisine, berberine, baicalein, etc. There were 362 targets treating DHD, among them the core targets included TNF, IL-6, ALB, etc. The enriched KEGG pathways mainly involve PI3K-Akt, TNF, MAPK, etc. Molecular docking results showed that coptisine, berberine, baicalein, etc., had a strong affinity with TNF, IL-6, and MAPK14. Therefore, TNF, IL-6, MAPK14, ALB, etc., are the key targets of the active ingredients of YJP coptisine, baicalein, and berberine, etc. They have the potential to regulate PI3K-Akt, MAPK, and TNF signalling pathways. The component-target-disease network diagram revealed that YJP treated DHD through the effects of anti-inflammation, anti-diarrhea, immunoregulation, and improving intestinal mucosal injury.

CONCLUSION

It is demonstrated that YJP treats DHD mainly through the main active ingredients coptisine, berberine, baicalein, etc. comprehensively exerting the effects of anti-inflammation, anti-diarrhea, immunoregulation, and improving intestinal mucosal injury, which will provide evidence for further in-depth studying the mechanism of YJP treating DHD.

Exploring the Molecular Mechanism of HongTeng Decoction against Inflammation based on Network Analysis and Experiments Validation

BACKGROUND

HongTeng Decoction (HTD) is a traditional Chinese medicine that is widely used to treat bacterial infections and chronic inflammation. However, its pharmacological mechanism is not clear. Here, network pharmacology and experimental verification were applied to investigate the drug targets and potential mechanisms of HTD in inflammation treatment.

METHODS

The active ingredients of HTD were collected from the multi-source databases and clarified by Q Exactive Orbitrap analysis in the treatment of inflammation. Then, molecular docking technology was used to explore the binding ability of key active ingredients and targets in HTD. In vitro experiments, the inflammatory factors and MAPK signaling pathways are detected to verify the anti-inflammatory effect of HTD on the RAW264.7 cells. Finally, the anti-inflammatory effect of HTD was evaluated in LPS induced mice model.

RESULTS

A total of 236 active compounds and 492 targets of HTD were obtained through database screening, and 954 potential targets of inflammation were identified. Finally, 164 possible targets of HTD acting on inflammation were obtained. The PPI analysis and KEGG enrichment analyses showed that the targets of HTD in inflammation were mostly related to the MAPK signaling pathway, the IL-17 signaling pathway, and the TNF signaling pathway. By integrating the results of the network analysis, the core targets of HTD in inflammation mainly include MAPK3, TNF, MMP9, IL6, EGFR, and NFKBIA. The molecular docking results indicated solid binding activity between MAPK3-naringenin and MAPK3-paeonol. It has been shown that HTD could inhibit the levels of inflammatory factors, IL6 and TNF-α, as well as the splenic index in the LPS-stimulated mice. Moreover, HTD could regulate protein expression levels of p-JNK1/2, and p-ERK1/2, which reflects the inhibitory effect of HTD on the MAPKS signaling pathway.

CONCLUSION

Our study is expected to provide the pharmacological mechanisms by which HTD may be a promising anti-inflammatory drug for future clinical trials.

Challenges of Protein-Protein Docking of the Membrane Proteins

Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.

Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

PICKLUSTER: a protein-interface clustering and analysis plug-in for UCSF ChimeraX

SUMMARY

The identification and characterization of interfaces in protein complexes is crucial for understanding the mechanisms of molecular recognition. These interfaces are also attractive targets for protein inhibition. However, targeting protein interfaces can be challenging for large interfaces that consist of multiple interacting regions. We present PICKLUSTER [Protein Interface C(K)luster]-a program for identifying "sub-interfaces" in protein-protein complexes using distance clustering. The division of the interface into smaller "sub-interfaces" offers a more focused approach for targeting protein-protein interfaces.

AVAILABILITY AND IMPLEMENTATION

PICKLUSTER is implemented as a plug-in for the molecular visualization program UCSF ChimeraX 1.4 and subsequent versions. It is freely available for download in the ChimeraX Toolshed and https://gitlab.com/topf-lab/pickluster.git.

Inferring the Effects of Protein Variants on Protein-Protein Interactions with Interpretable Transformer Representations

Identifying pathogenetic variants and inferring their impact on protein-protein interactions sheds light on their functional consequences on diseases. Limited by the availability of experimental data on the consequences of protein interaction, most existing methods focus on building models to predict changes in protein binding affinity. Here, we introduced MIPPI, an end-to-end, interpretable transformer-based deep learning model that learns features directly from sequences by leveraging the interaction data from IMEx. MIPPI was specifically trained to determine the types of variant impact (increasing, decreasing, disrupting, and no effect) on protein-protein interactions. We demonstrate the accuracy of MIPPI and provide interpretation through the analysis of learned attention weights, which exhibit correlations with the amino acids interacting with the variant. Moreover, we showed the practicality of MIPPI in prioritizing de novo mutations associated with complex neurodevelopmental disorders and the potential to determine the pathogenic and driving mutations. Finally, we experimentally validated the functional impact of several variants identified in patients with such disorders. Overall, MIPPI emerges as a versatile, robust, and interpretable model, capable of effectively predicting mutation impacts on protein-protein interactions and facilitating the discovery of clinically actionable variants.

New insights on the differential interaction of sulfiredoxin with members of the peroxiredoxin family revealed by protein-protein docking and experimental studies

Sulfiredoxin (Srx) is the enzyme that restores the peroxidase activity of peroxiredoxins (Prxs) through catalyzing the reduction of hyperoxidized Prxs back to their active forms. This process involves protein-protein interaction in an enzyme-substrate binding manner. The integrity of the Srx-Prx axis contributes to the pathogenesis of various oxidative stress related human disorders including cancer, inflammation, cardiovascular and neurological diseases. The purpose of this study is to understand the structural and molecular biology of the Srx-Prx interaction, which may be of significance for prediction of target site for the novel drug-discovery. Homology modeling and protein-protein docking approaches were applied to examine the Srx-Prx interaction using online platforms including ITASSER, Phyre2, Swissmodel, AlphaFold, MZDOCK and ZDOCK. By in-silico studies, A 26-amino acid motif at the C-terminus of Prx1 was predicted to cause a steric hindrance for the kinetics of the Srx-Prx1 interaction. These predictions were tested in-vitro using purified recombinant proteins including Srx, full-length Prxs, and C-terminus deleted Prxs. We confirmed that deletion of the C-terminus of Prxs significantly enhanced its rate of association with Srx (i.e. >1000 fold increase in the ka of the Srx-Prx1 interaction) with minimal effect on the rate of dissociation (kd). Differential interaction of Srx with individual members of the Prx family was further examined in cultured cells. Taken together, these data add novel molecular and structural insights critical for the understanding of the biology of the Srx-Prx interaction that may be of value for the development of targeted therapy for human disorders.

xProtCAS: A Toolkit for Extracting Conserved Accessible Surfaces from Protein Structures

The identification of protein surfaces required for interaction with other biomolecules broadens our understanding of protein function, their regulation by post-translational modification, and the deleterious effect of disease mutations. Protein interaction interfaces are often identifiable as patches of conserved residues on a protein's surface. However, finding conserved accessible surfaces on folded regions requires an understanding of the protein structure to discriminate between functional and structural constraints on residue conservation. With the emergence of deep learning methods for protein structure prediction, high-quality structural models are now available for any protein. In this study, we introduce tools to identify conserved surfaces on AlphaFold2 structural models. We define autonomous structural modules from the structural models and convert these modules to a graph encoding residue topology, accessibility, and conservation. Conserved surfaces are then extracted using a novel eigenvector centrality-based approach. We apply the tool to the human proteome identifying hundreds of uncharacterised yet highly conserved surfaces, many of which contain clinically significant mutations. The xProtCAS tool is available as open-source Python software and an interactive web server.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

Exploring the molecular mechanism of Gan Shuang granules for the treatment of non-alcoholic steatohepatitis using network pharmacology, molecular docking, and experimental verification

Background: With the gradual increase in prevalence in recent years, non-alcoholic steatohepatitis (NASH) has become one of the significant health problems that urgently needs to be addressed worldwide. GanShuang Granules (GSG) is derived from the classical Chinese formula Xiaoyao San and mainly used in the clinical treatment of chronic liver diseases. Objective: In this study, we aim to gain a deeper insight into the inhibiting effects of GSG on non-alcoholic fatty liver disease (NAFLD) rats and preliminarily elucidate the underlying intervention mechanisms. Methods: First, High performance liquid chromatography (UHPLC-Q/Orbitrap-MS/MS) was used for the active compounds prediction in GSG. Then the data was mapped to mzCloud database. The targets corresponding to GSG compounds were collected from public databases, along with disease genes for NAFLD. The core targets and molecular mechanisms of GSG for NAFLD treatment were predicted by protein-protein interaction (PPI) network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses. Molecular docking of the core target-component interactions was simulated using AutoDock Vina software. The effect of GSG on NASH rats was evaluated by pathological staining and analysis of various index results. Finally, the candidate targets were further validated by ELISA and western blot (WB) analyses. Results: Combining UHPLC-Q/Orbitrap-MS/MS data analysis and public database data, a total of 346 cross-targets were obtained, corresponding to 81 compounds. The subnetwork with an MCODE score of 53.623 is a potential core target group for this study. GO and KEGG enrichment analyses showed that the targets of GSG in NAFLD were mostly related to oxidative stress, the NF-κB signaling pathway, and the apoptosis signaling pathway. By integrating the results of network pharmacology analysis, the core objectives of this study mainly include AKT1, CASP9, TNF, and CASP8. The core ingredients are related to resveratrol and fisetin. The molecular docking results indicated key binding activity between AKT1-fisetin, AKT1-Resveratrol, and CASP8-fisetin. Moreover, GSG could improve the inflammatory status and restore the abnormal lipid accumulation of NAFLD/NASH liver, and these levels are further verified by pathological staining and detection of related indicators. Mechanistically, GSG could regulate protein expression levels in the liver for P65, p-P65, IKB, p-IKB, IKK, caspase-3, -8, -9, and cytochrome C, etc. It reflects the inhibitory effect of GSG on the NF-κB/IκB signaling pathway. Conclusion: Our results suggested that GSG demonstrated therapeutic effects on NAFLD/NASH rats, and these may be mainly reflected in the inhibitory effects on the NF-κB/IκB signaling pathway and its downstream inflammation and apoptosis signals.

Growing ecosystem of deep learning methods for modeling protein-protein interactions

Numerous cellular functions rely on protein-protein interactions. Efforts to comprehensively characterize them remain challenged however by the diversity of molecular recognition mechanisms employed within the proteome. Deep learning has emerged as a promising approach for tackling this problem by exploiting both experimental data and basic biophysical knowledge about protein interactions. Here, we review the growing ecosystem of deep learning methods for modeling protein interactions, highlighting the diversity of these biophysically informed models and their respective trade-offs. We discuss recent successes in using representation learning to capture complex features pertinent to predicting protein interactions and interaction sites, geometric deep learning to reason over protein structures and predict complex structures, and generative modeling to design de novo protein assemblies. We also outline some of the outstanding challenges and promising new directions. Opportunities abound to discover novel interactions, elucidate their physical mechanisms, and engineer binders to modulate their functions using deep learning and, ultimately, unravel how protein interactions orchestrate complex cellular behaviors.

Artificial Intelligence and Machine Learning Technology Driven Modern Drug Discovery and Development

The discovery and advances of medicines may be considered as the ultimate relevant translational science effort that adds to human invulnerability and happiness. But advancing a fresh medication is a quite convoluted, costly, and protracted operation, normally costing USD ~2.6 billion and consuming a mean time span of 12 years. Methods to cut back expenditure and hasten new drug discovery have prompted an arduous and compelling brainstorming exercise in the pharmaceutical industry. The engagement of Artificial Intelligence (AI), including the deep-learning (DL) component in particular, has been facilitated by the employment of classified big data, in concert with strikingly reinforced computing prowess and cloud storage, across all fields. AI has energized computer-facilitated drug discovery. An unrestricted espousing of machine learning (ML), especially DL, in many scientific specialties, and the technological refinements in computing hardware and software, in concert with various aspects of the problem, sustain this progress. ML algorithms have been extensively engaged for computer-facilitated drug discovery. DL methods, such as artificial neural networks (ANNs) comprising multiple buried processing layers, have of late seen a resurgence due to their capability to power automatic attribute elicitations from the input data, coupled with their ability to obtain nonlinear input-output pertinencies. Such features of DL methods augment classical ML techniques which bank on human-contrived molecular descriptors. A major part of the early reluctance concerning utility of AI in pharmaceutical discovery has begun to melt, thereby advancing medicinal chemistry. AI, along with modern experimental technical knowledge, is anticipated to invigorate the quest for new and improved pharmaceuticals in an expeditious, economical, and increasingly compelling manner. DL-facilitated methods have just initiated kickstarting for some integral issues in drug discovery. Many technological advances, such as "message-passing paradigms", "spatial-symmetry-preserving networks", "hybrid de novo design", and other ingenious ML exemplars, will definitely come to be pervasively widespread and help dissect many of the biggest, and most intriguing inquiries. Open data allocation and model augmentation will exert a decisive hold during the progress of drug discovery employing AI. This review will address the impending utilizations of AI to refine and bolster the drug discovery operation.

MetaScore: A Novel Machine-Learning-Based Approach to Improve Traditional Scoring Functions for Scoring Protein-Protein Docking Conformations

Protein-protein interactions play a ubiquitous role in biological function. Knowledge of the three-dimensional (3D) structures of the complexes they form is essential for understanding the structural basis of those interactions and how they orchestrate key cellular processes. Computational docking has become an indispensable alternative to the expensive and time-consuming experimental approaches for determining the 3D structures of protein complexes. Despite recent progress, identifying near-native models from a large set of conformations sampled by docking-the so-called scoring problem-still has considerable room for improvement. We present MetaScore, a new machine-learning-based approach to improve the scoring of docked conformations. MetaScore utilizes a random forest (RF) classifier trained to distinguish near-native from non-native conformations using their protein-protein interfacial features. The features include physicochemical properties, energy terms, interaction-propensity-based features, geometric properties, interface topology features, evolutionary conservation, and also scores produced by traditional scoring functions (SFs). MetaScore scores docked conformations by simply averaging the score produced by the RF classifier with that produced by any traditional SF. We demonstrate that (i) MetaScore consistently outperforms each of the nine traditional SFs included in this work in terms of success rate and hit rate evaluated over conformations ranked among the top 10; (ii) an ensemble method, MetaScore-Ensemble, that combines 10 variants of MetaScore obtained by combining the RF score with each of the traditional SFs outperforms each of the MetaScore variants. We conclude that the performance of traditional SFs can be improved upon by using machine learning to judiciously leverage protein-protein interfacial features and by using ensemble methods to combine multiple scoring functions.

Pathogen-driven cancers from a structural perspective: Targeting host-pathogen protein-protein interactions

Host-pathogen interactions (HPIs) affect and involve multiple mechanisms in both the pathogen and the host. Pathogen interactions disrupt homeostasis in host cells, with their toxins interfering with host mechanisms, resulting in infections, diseases, and disorders, extending from AIDS and COVID-19, to cancer. Studies of the three-dimensional (3D) structures of host-pathogen complexes aim to understand how pathogens interact with their hosts. They also aim to contribute to the development of rational therapeutics, as well as preventive measures. However, structural studies are fraught with challenges toward these aims. This review describes the state-of-the-art in protein-protein interactions (PPIs) between the host and pathogens from the structural standpoint. It discusses computational aspects of predicting these PPIs, including machine learning (ML) and artificial intelligence (AI)-driven, and overviews available computational methods and their challenges. It concludes with examples of how theoretical computational approaches can result in a therapeutic agent with a potential of being used in the clinics, as well as future directions.

Mechanisms of Xiaochaihu Decoction on Treating Hepatic Fibrosis Explored by Network Pharmacology

Purpose. To explore the material basis and pharmacological mechanism of Xiaochaihu Decoction (XCHD), the classic Traditional Chinese Medicine (TCM) formula in inhibiting hepatic fibrosis (HF). Methods. The main components in XCHD were screened from the TCMSP database, ETCM database, and literature, and their potential targets were detected and predicted using the Swiss Target Prediction platform. The HF-related targets were retrieved and screened through GeneCard database and OMIM database, combined with GEO gene chips. The XCHD targets and HF targets were mapped to search common targets. The protein-protein interaction (PPI) network was acquired via the STRING11.0 database and analyzed visually using Cytoscape 3.8.0 software. The potential mechanisms of the common targets identified through GO and KEGG pathway enrichment analysis were analyzed by using Metascape database. The results were visualized through OmicShare Tools. The “XCHD compound-HF target” network was visually constructed by Cytoscape 3.8.0 software. AutoDockVina1.1.2 and PyMoL software were used to verify the molecular docking of XCHD main active compounds and HF key targets. Results. A total of 164 potential active compounds from XCHD were screened to act on 95 HF-related targets. Bioinformatics analysis revealed that quercetin, β-sitosterol, and kaempferol may be candidate agents, which acted on multiple targets like PTGS2, HSP90AA1, and PTGS1 and regulate multiple key biological pathways like IL-17 signaling pathway, TNF signaling pathway and PI3K-Akt signaling pathway to relieve HF. Moreover, molecular docking suggested that quercetin and PTGS2 could statically bind and interact with each other through amino acid residues val-349, LEU-352, PHE-381, etc. Conclusion. This work provides a systems perspective to study the relationship between Chinese medicines and diseases. The therapeutic efficacy of XCHD on HF was the sum of multitarget and multi-approach effects from the bioactive ingredients. This study could be one of the cornerstones for further research.

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.

A Transendothelial Leukocyte Transmigration Model Based on Computational Fluid Dynamics and BP Neural Network

The mechanism of immune infiltration involving immune cells is closely related to various diseases. A key issue in immune infiltration is the transendothelial transmigration of leukocytes. Previous studies have primarily interpreted the leukocyte infiltration of from biomedical perspective. The physical mechanism of leukocyte infiltration remains to be explored. By integrating the immune cell transmigration computational fluid dynamics (CFD) data, the paper builds a time-dependent leukocyte transmigration prediction model based on the bio-inspired methods, namely back propagation neural networks (BPNN) model. The model can efficiently predict the immune cell transmigration in a special microvascular environment, and obtain good prediction accuracy. The model accurately predicted the cell movement and flow field changes during the transmigration. In the test data set, it has high prediction accuracy for cell deformation, motion velocity and flow lift forces during downstream motion, and maintains a good prediction accuracy for drag force. The two prediction models achieved the prediction of leukocyte transmigration in a specific microvascular environment and maintained a high prediction accuracy, indicating the feasibility and robustness of the BPNN model applied to the prediction of immune cell infiltration. Compared with traditional CFD simulations, BPNN models avoid complex and time-dependent physical modeling and computational processes.

Characterizing and explaining the impact of disease-associated mutations in proteins without known structures or structural homologs

Mutations in human proteins lead to diseases. The structure of these proteins can help understand the mechanism of such diseases and develop therapeutics against them. With improved deep learning techniques, such as RoseTTAFold and AlphaFold, we can predict the structure of proteins even in the absence of structural homologs. We modeled and extracted the domains from 553 disease-associated human proteins without known protein structures or close homologs in the Protein Databank. We noticed that the model quality was higher and the Root mean square deviation (RMSD) lower between AlphaFold and RoseTTAFold models for domains that could be assigned to CATH families as compared to those which could only be assigned to Pfam families of unknown structure or could not be assigned to either. We predicted ligand-binding sites, protein–protein interfaces and conserved residues in these predicted structures. We then explored whether the disease-associated missense mutations were in the proximity of these predicted functional sites, whether they destabilized the protein structure based on ddG calculations or whether they were predicted to be pathogenic. We could explain 80% of these disease-associated mutations based on proximity to functional sites, structural destabilization or pathogenicity. When compared to polymorphisms, a larger percentage of disease-associated missense mutations were buried, closer to predicted functional sites, predicted as destabilizing and pathogenic. Usage of models from the two state-of-the-art techniques provide better confidence in our predictions, and we explain 93 additional mutations based on RoseTTAFold models which could not be explained based solely on AlphaFold models.

Identifying vaccine escape sites via statistical comparisons of short-term molecular dynamics

The identification of viral mutations that confer escape from antibodies is crucial for understanding the interplay between immunity and viral evolution. We describe a molecular dynamics (MD)-based approach that goes beyond contact mapping, scales well to a desktop computer with a modern graphics processor, and enables the user to identify functional protein sites that are prone to vaccine escape in a viral antigen. We first implement our MD pipeline to employ site-wise calculation of Kullback-Leibler divergence in atom fluctuation over replicate sets of short-term MD production runs thus enabling a statistical comparison of the rapid motion of influenza hemagglutinin (HA) in both the presence and absence of three well-known neutralizing antibodies. Using this simple comparative method applied to motions of viral proteins, we successfully identified in silico all previously empirically confirmed sites of escape in influenza HA, predetermined via selection experiments and neutralization assays. Upon the validation of our computational approach, we then surveyed potential hotspot residues in the receptor binding domain of the SARS-CoV-2 virus in the presence of COVOX-222 and S2H97 antibodies. We identified many single sites in the antigen-antibody interface that are similarly prone to potential antibody escape and that match many of the known sites of mutations arising in the SARS-CoV-2 variants of concern. In the Omicron variant, we find only minimal adaptive evolutionary shifts in the functional binding profiles of both antibodies. In summary, we provide an inexpensive and accurate computational method to monitor hotspots of functional evolution in antibody binding footprints.

CSM-Potential: mapping protein interactions and biological ligands in 3D space using geometric deep learning

Recent advances in protein structural modelling have enabled the accurate prediction of the holo 3D structures of almost any protein, however protein function is intrinsically linked to the interactions it makes. While a number of computational approaches have been proposed to explore potential biological interactions, they have been limited to specific interactions, and have not been readily accessible for non-experts or use in bioinformatics pipelines. Here we present CSM-Potential, a geometric deep learning approach to identify regions of a protein surface that are likely to mediate protein-protein and protein-ligand interactions in order to provide a link between 3D structure and biological function. Our method has shown robust performance, outperforming existing methods for both predictive tasks. By assessing the performance of CSM-Potential on independent blind tests, we show that our method was able to achieve ROC AUC values of up to 0.81 for the identification of potential protein-protein binding sites, and up to 0.96 accuracy on biological ligand classification. Our method is freely available as a user-friendly and easy-to-use web server and API at http://biosig.unimelb.edu.au/csm_potential.

Complementarity of the residue-level protein function and structure predictions in human proteins

Sequence-based predictors of the residue-level protein function and structure cover a broad spectrum of characteristics including intrinsic disorder, secondary structure, solvent accessibility and binding to nucleic acids. They were catalogued and evaluated in numerous surveys and assessments. However, methods focusing on a given characteristic are studied separately from predictors of other characteristics, while they are typically used on the same proteins. We fill this void by studying complementarity of a representative collection of methods that target different predictions using a large, taxonomically consistent, and low similarity dataset of human proteins. First, we bridge the gap between the communities that develop structure-trained vs. disorder-trained predictors of binding residues. Motivated by a recent study of the protein-binding residue predictions, we empirically find that combining the structure-trained and disorder-trained predictors of the DNA-binding and RNA-binding residues leads to substantial improvements in predictive quality. Second, we investigate whether diverse predictors generate results that accurately reproduce relations between secondary structure, solvent accessibility, interaction sites, and intrinsic disorder that are present in the experimental data. Our empirical analysis concludes that predictions accurately reflect all combinations of these relations. Altogether, this study provides unique insights that support combining results produced by diverse residue-level predictors of protein function and structure.

Network Pharmacology Analysis of Huangqi Jianzhong Tang Targets in Gastric Cancer

Background: The Chinese medicine, Huangqi Jianzhong Tang (HJT), is widely used to treat gastric cancer (GC). In this study, network pharmacological methods were used to analyze the potential therapeutic targets and pharmacological mechanisms of HJT in GC.

Methods: Bioactive components and targets of HJT and GC-related targets were identified using public databases. The protein-protein interaction network of potential targets of HJT in GC was constructed using the Cytoscape plug-in (v3.8.0), CytoHubba. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed, in addition to molecular docking and animal experiments to verify the results of network pharmacology analysis.

Results: A total of 538 GC-related targets were identified. The bioactive components of HJT were selected for drug-likeness evaluation and binomial statistical model screening, which revealed 63 bioactive components and 72 targets. Based on GO enrichment analysis, all targets in the protein-protein interaction network were mainly involved in the response to oxidative stress and neuronal death. Further, KEGG enrichment analysis suggested that the treatment of GC with HJT mainly involved the Wnt signaling pathway, PI3K-Akt signaling pathway, TGF-β signaling pathway, and MAPK signaling pathway, thereby providing insights into the mechanism of the effects of HJT on GC.

Conclusion: This study revealed the potential bioactive components and molecular mechanisms of HJT, which may be useful for the treatment of GC, and provided insights into the development of new drugs for GC.

A novel computational strategy for defining the minimal protein molecular surface representation

Most proteins perform their biological function by interacting with one or more molecular partners. In this respect, characterizing local features of the molecular surface, that can potentially be involved in the interaction with other molecules, represents a step forward in the investigation of the mechanisms of recognition and binding between molecules. Predictive methods often rely on extensive samplings of molecular patches with the aim to identify hot spots on the surface. In this framework, analysis of large proteins and/or many molecular dynamics frames is often unfeasible due to the high computational cost. Thus, finding optimal ways to reduce the number of points to be sampled maintaining the biological information (including the surface shape) carried by the molecular surface is pivotal. In this perspective, we here present a new theoretical and computational algorithm with the aim of defining a set of molecular surfaces composed of points not uniformly distributed in space, in such a way as to maximize the information of the overall shape of the molecule by minimizing the number of total points. We test our procedure’s ability in recognizing hot-spots by describing the local shape properties of portions of molecular surfaces through a recently developed method based on the formalism of 2D Zernike polynomials. The results of this work show the ability of the proposed algorithm to preserve the key information of the molecular surface using a reduced number of points compared to the complete surface, where all points of the surface are used for the description. In fact, the methodology shows a significant gain of the information stored in the sampling procedure compared to uniform random sampling.

Molecular docking analysis of two bioactive molecules KLUF10 and KLUF13 isolated from the marine bacteria Micrococcus sp. OUS9 with TNF alpha

Tumor necrosis factor-alpha (TNF-α) is known to be linked with tumor. Therefore, it is of interest to document the Molecular docking analysis of two bioactive molecules KLUF10 and KLUF13 isolated from the marine bacteria Micrococcus sp. OUS9 with TNF alpha. We report the molecular interactions of KLUF10 and KLUF13 with TNF alpha.

DeepDISOBind: accurate prediction of RNA-, DNA- and protein-binding intrinsically disordered residues with deep multi-task learning

Proteins with intrinsically disordered regions (IDRs) are common among eukaryotes. Many IDRs interact with nucleic acids and proteins. Annotation of these interactions is supported by computational predictors, but to date, only one tool that predicts interactions with nucleic acids was released, and recent assessments demonstrate that current predictors offer modest levels of accuracy. We have developed DeepDISOBind, an innovative deep multi-task architecture that accurately predicts deoxyribonucleic acid (DNA)-, ribonucleic acid (RNA)- and protein-binding IDRs from protein sequences. DeepDISOBind relies on an information-rich sequence profile that is processed by an innovative multi-task deep neural network, where subsequent layers are gradually specialized to predict interactions with specific partner types. The common input layer links to a layer that differentiates protein- and nucleic acid-binding, which further links to layers that discriminate between DNA and RNA interactions. Empirical tests show that this multi-task design provides statistically significant gains in predictive quality across the three partner types when compared to a single-task design and a representative selection of the existing methods that cover both disorder- and structure-trained tools. Analysis of the predictions on the human proteome reveals that DeepDISOBind predictions can be encoded into protein-level propensities that accurately predict DNA- and RNA-binding proteins and protein hubs. DeepDISOBind is available at https://www.csuligroup.com/DeepDISOBind/.

DeepRank: a deep learning framework for data mining 3D protein-protein interfaces

Three-dimensional (3D) structures of protein complexes provide fundamental information to decipher biological processes at the molecular scale. The vast amount of experimentally and computationally resolved protein-protein interfaces (PPIs) offers the possibility of training deep learning models to aid the predictions of their biological relevance. We present here DeepRank, a general, configurable deep learning framework for data mining PPIs using 3D convolutional neural networks (CNNs). DeepRank maps features of PPIs onto 3D grids and trains a user-specified CNN on these 3D grids. DeepRank allows for efficient training of 3D CNNs with data sets containing millions of PPIs and supports both classification and regression. We demonstrate the performance of DeepRank on two distinct challenges: The classification of biological versus crystallographic PPIs, and the ranking of docking models. For both problems DeepRank is competitive with, or outperforms, state-of-the-art methods, demonstrating the versatility of the framework for research in structural biology.

The role of complementary shape in protein dimerization

Shape guides colloidal nanoparticles to form complex assemblies, but its role in defining interfaces in biomolecular complexes is less clear. In this work, we isolate the role of shape in protein complexes by studying the reversible binding processes of 46 protein dimer pairs, and investigate when entropic effects from shape complementarity alone are sufficient to predict the native protein binding interface. We employ depletants using a generic, implicit depletion model to amplify the magnitude of the entropic forces arising from lock-and-key binding and isolate the effect of shape complementarity in protein dimerization. For 13% of the complexes studied here, protein shape is sufficient to predict native complexes as equilibrium assemblies. We elucidate the results by analyzing the importance of competing binding configurations and how it affects the assembly. A machine learning classifier, with a precision of 89.14% and a recall of 77.11%, is able to identify the cases where shape alone predicts the native protein interface.

Machine Learning in Drug Discovery: A Review

This review provides the feasible literature on drug discovery through ML tools and techniques that are enforced in every phase of drug development to accelerate the research process and deduce the risk and expenditure in clinical trials. Machine learning techniques improve the decision-making in pharmaceutical data across various applications like QSAR analysis, hit discoveries, de novo drug architectures to retrieve accurate outcomes. Target validation, prognostic biomarkers, digital pathology are considered under problem statements in this review. ML challenges must be applicable for the main cause of inadequacy in interpretability outcomes that may restrict the applications in drug discovery. In clinical trials, absolute and methodological data must be generated to tackle many puzzles in validating ML techniques, improving decision-making, promoting awareness in ML approaches, and deducing risk failures in drug discovery.

Network Pharmacology Prediction and Molecular Docking-Based Strategy to Discover the Potential Pharmacological Mechanism of Huai Hua San Against Ulcerative Colitis

Background

Huai Hua San (HHS), a famous Traditional Chinese Medicine (TCM) formula, has been widely applied in treating ulcerative colitis (UC). However, the interaction of bioactives from HHS with the targets involved in UC has not been elucidated yet.

Aim

A network pharmacology-based approach combined with molecular docking and in vitro validation was performed to determine the bioactives, key targets, and potential pharmacological mechanism of HHS against UC.

Materials and Methods

Bioactives and potential targets of HHS, as well as UC-related targets, were retrieved from public databases. Crucial bioactive ingredients, potential targets, and signaling pathways were acquired through bioinformatics analysis, including protein-protein interaction (PPI), as well as the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Subsequently, molecular docking was carried out to predict the combination of active compounds with core targets. Lastly, in vitro experiments were conducted to further verify the findings.

Results

A total of 28 bioactive ingredients of HHS and 421 HHS-UC-related targets were screened. Bioinformatics analysis revealed that quercetin, luteolin, and nobiletin may be potential candidate agents. JUN, TP53, and ESR1 could become potential therapeutic targets. PI3K-AKT signaling pathway might play an important role in HHS against UC. Moreover, molecular docking suggested that quercetin, luteolin, and nobiletin combined well with JUN, TP53, and ESR1, respectively. Cell experiments showed that the most important ingredient of HHS, quercetin, could inhibit the levels of inflammatory factors and phosphorylated c-Jun, as well as PI3K-Akt signaling pathway in LPS-induced RAW264.7 cells, which further confirmed the prediction by network pharmacology strategy and molecular docking.

Conclusion

Our results comprehensively illustrated the bioactives, potential targets, and molecular mechanism of HHS against UC. It also provided a promising strategy to uncover the scientific basis and therapeutic mechanism of TCM formulae in treating diseases.

Identifying hydrophobic protein patches to inform protein interaction interfaces

Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.

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.

Marker residue types at the structural regions of transmembrane alpha-helical and beta-barrel interfaces

Membrane proteins play a variety of biological functions to the survival of organisms and functionalities of these proteins are often due to their homo- or hetero-complexation. Encoded by ~30% of the genome in most organisms, they represent the target of over half of nowadays drugs. Spanning the entirety of the cell membrane, transmembrane proteins are the most common type of membrane proteins and can be classified by secondary structures: alpha-helical and beta-barrel structures. Protein-protein interaction (PPI) have been widely studied for globular proteins and many computational tools are available for predicting PPI sites and construct models of complexes. Here, the structural regions of a non-redundant set of 232 alpha-helical and 37 beta-barrel transmembrane complexes and their interfaces are analyzed. Using the residue composition, frequency and propensity, this study brings the light on the marker residue types located at the structural regions of alpha-helical and beta-barrel transmembrane homomeric protein complexes and of their interfaces. This study also shows the necessity to relate the frequency to the composition into a ratio for immediately figuring out residue types presenting high frequencies at the interface and/or at one of its structural regions despite being a minor contributor compared to other residue types to that location's residue composition.

Computational investigation of FDA approved drugs as selective PARP-1 inhibitors by targeting BRCT domain for cancer therapy

Poly(ADP-ribose) polymerase-1 is a promising target for the treatment of cancer due to its involvement in base excision repair pathways for repairing DNA single-strand breaks. However, available PARP-1 inhibitors target a highly conserved PARPs catalytic domain, which causes toxicity due to the off-target activity. Therefore, the present study was hypothesized to identify selective inhibitors by targeting specific protein-protein interacting (PPI) PARP-1 BRCT domain. Moreover, PPI hotspot residues (Gly399, Lys400, Leu401, Lys441 & Lys442) and a druggable pocket was detected to screen small molecule inhibitors. Hence, two FDA approved drug molecules (levoleucovorin and balsalazide) were recognized to fit in the druggable pocket. Since they are already under investigation for anti-cancer activity, thus could be further explored in PARP-1 sensitive cancer cells to expand their selectivity and develop as effective anti-cancer agents. Besides, the study also provides detailed structural insight of PARP-1 and XRCC1 complex through their BRCT domains.