A Computational Investigation of Small-Molecule Engagement of Hot Spots at Protein-Protein Interaction Interfaces

Identification of potent anti-immunogenic agents through virtual screening, 3D-QSAR studies, and in vitro experiments

A wealth of literature has highlighted the discovery of various immune modulators, frequently used in clinical practice, yet associated with numerous drawbacks. In light of this pharmacological deficiency, medical scientists are motivated to develop new immune modulators with minimized adverse effects yet retaining the improved therapeutic potential. T-cell differentiation and growth are central to human defense and are regulated by interleukin-2 (IL-2), an immune-modulatory cytokine. However, scientific investigation is hindered due to its flat binding site and widespread hotspot residues. In this regard, a prompt and logical investigation guided by integrated computational techniques was undertaken to unravel new and potential leads against IL-2. In particular, the combination of score-based and pharmacophore-based virtual screening approaches were employed, reducing the data from millions of small molecules to a manageable number. Subsequent docking and 3D-QSAR prediction via CoMFA further helped remove false positives from the data. The reliability of the model was assessed via standard metrics, which explain the model's fitness and the robustness of the model in predicting the activity of new compounds. The extensive virtual screening herein led to the identification of a total of 24 leads with potential anti-IL-2 activity. Furthermore, the theoretical findings were corroborated with in vitro testing, further endorsing the anti-inflammatory potential of the identified leads.

Hot-Spot Residue-Based Virtual Screening of Novel Selective Estrogen-Receptor Degraders for Breast Cancer Treatment

The estrogen-receptor alfa (ERα) is considered pivotal for breast cancer treatment. Although selective estrogen-receptor degraders (SERDs) have been developed to induce ERα degradation and antagonism, their agonistic effect on the uterine tissue and poor pharmacokinetic properties limit further application of ERα; thus, discovering novel SERDs is necessary. The ligand preferentially interacts with several key residues of the protein (defined as hot-spot residues). Improving the interaction with hot-spot residues of ERα offers a promising avenue for obtaining novel SERDs. In this study, pharmacophore modeling, molecular mechanics/generalized Born surface area (MM/GBSA), and amino-acid mutation were combined to determine several hot-spot residues. Focusing on the interaction with these hot-spot residues, hit fragments A1-A3 and A9 were virtually screened from two fragment libraries. Finally, these hit fragments were linked to generate compounds B1-B3, and their biological activities were evaluated. Remarkably, compound B1 exhibited potent antitumor activity against MCF-7 cells (IC50 = 4.21 nM), favorable ERα binding affinity (Ki = 14.6 nM), and excellent ERα degradative ability (DC50 = 9.7 nM), which indicated its potential to evolve as a promising SERD for breast cancer treatment.

In Silico molecular docking and dynamic analysis of natural compounds against major non-structural proteins of SARS-COV-2

COVID-19 has infected millions and significantly affected the global economy and healthcare systems. Despite continuous lockdowns, symptomatic management with currently available medications, and numerous vaccination drives, it is still far more difficult to control. Against COVID-19 infection, the pressure to develop vaccines and drugs has led to using some currently available medications like remdesivir, azithromycin, hydroxychloroquine and ritonavir. Understanding the importance and potential of harmless molecules to tackle SARS-COV-2, we designed the present study to identify potential natural phytocompounds. In the present study, we docked natural compounds and standard drugs against SARS-COV-2 proteins: papain-like protease, main protease and helicase. ADME/T and ProTox-II analyses were used to determine the toxicity of phytocompounds and drugs. The docking analysis revealed that podophyllotoxin gave the highest binding affinity scores of -8.1, -7.1 and -7.4 kcal/mol against PLpro, Mpro and helicase, respectively. Among the control drugs, doxycycline hydrochloride showed the highest binding affinity of -10.5, -8.4 and -8.8 kcal/mol against PLpro, Mpro and helicase. The results of this study revealed that podophyllotoxin and doxycycline hydrochloride could be promising inhibitors against SARS-Cov-2. Molecular dynamic simulations were executed for the best docked (PLpro-podophyllotoxin) complex, and the results displayed stable conformation and convergence. Energy plot results predicted a global minima average energy of -95 kcal/mol and indicated podophyllotoxin's role in stabilizing protein and making it compact and complex. FarPPI server used MM/GBSA approach to determine free binding affinity, and helicase-gallic acid complex showed the highest affinity, respectively. Therefore, it can be concluded that there is still a need for in vitro and in vivo studies to support further and validate these findings and validate these findings.Communicated by Ramaswamy H. Sarma.

Analysis of physicochemical properties of protein-protein interaction modulators suggests stronger alignment with the "rule of five"

Despite the important roles played by protein-protein interactions (PPIs) in disease, they have been long considered as 'undruggable'. However, recent advances have suggested that PPIs are druggable but may not follow conventional rules of 'drug ability'. Here we explore which physicochemical parameters are essential for a PPI modulator to be a clinical drug by analysing the physicochemical properties of small-molecule PPI modulators in the market, in clinical trials, and published. Our analysis reveals that those compounds currently on the market have a larger range of values for most of the physicochemical parameters, whereas those in clinical trials fit much more stringently to standard drug-like parameters. This observation was particularly true for molecular weight, clog P and topological polar surface area, where aside from a few outliers, most of the compounds in clinical trials fit within standard drug-like parameters. This implies that the newer PPI modulators are more drug-like than those currently on the market, suggesting that designing new PPI-specific screening libraries should remain within standard drug-like parameters in order to obtain a clinical candidate. Taken together, our analysis has important implications for designing future drug discovery campaigns aimed at targeting PPIs.

The protein-protein interaction network alignment using recurrent neural network

The main challenge of biological network alignment is that the problem of finding the alignments in two graphs is NP-hard. The discovery of protein-protein interaction (PPI) networks is of great importance in bioinformatics due to their utilization in identifying the cellular pathways, finding new medicines, and disease recognition. In this regard, we describe the network alignment method in the form of a classification problem for the very first time and introduce a deep network that finds the alignment of nodes present in the two networks. We call this method RENA, which means Network Alignment using REcurrent neural network. The proposed solution consists of three steps; in the first phase, we obtain the sequence and topological similarities from the networks' structure. For the second phase, the dataset needed for the transformation of the problem into a classification problem is created from obtained features. In the third phase, we predict the nodes' alignment between two networks using deep learning. We used Biogrid dataset for RENA evaluation. The RENA method is compared with three classification approaches of support vector machine, K-nearest neighbors, and linear discriminant analysis. The experimental results demonstrate the efficiency of the RENA method and 100% accuracy in PPI network alignment prediction.

Small-Molecule Inhibition of the uPAR ⋅ uPA Interaction by Conformational Selection

The urokinase receptor (uPAR) is a cell surface receptor that binds to the serine protease urokinase-type plasminogen activator (uPA) with high affinity. This interaction is beneficial for extravascular fibrin clearance, but it has also been associated with a broad range of pathological conditions including cancer, atherosclerosis, and kidney disease. Here, starting with a small molecule that we previously discovered by virtual screening and cheminformatics analysis, we design and synthesize several derivatives that were tested for binding and inhibition of the uPAR ⋅ uPA interaction. To confirm the binding site and establish a binding mode of the compounds, we carried out biophysical studies using uPAR mutants, among them uPAR^H47C-N259C , a mutant previously developed to mimic the structure of uPA-bound uPAR. Remarkably, a substantial increase in potency is observed for inhibition of uPAR^H47C-N259C binding to uPA compared to wild-type uPAR, consistent with our use of the structure of uPAR in its uPA-bound state to design small-molecule uPAR ⋅ uPA antagonists. Combined with the biophysical studies, molecular docking followed by extensive explicit-solvent molecular dynamics simulations and MM-GBSA free energy calculations yielded the most favorable binding pose of the compound. Collectively, these results suggest that potent inhibition of uPAR binding to uPA with small molecules will likely only be achieved by developing small molecules that exhibit high-affinity to solution apo structures of uPAR, rather than uPA-bound structures of the receptor.

HawkDock: a web server to predict and analyze the protein-protein complex based on computational docking and MM/GBSA

Protein–protein interactions (PPIs) play an important role in the different functions of cells, but accurate prediction of the three-dimensional structures for PPIs is still a notoriously difficult task. In this study, HawkDock, a free and open accessed web server, was developed to predict and analyze the structures of PPIs. In the HawkDock server, the ATTRACT docking algorithm, the HawkRank scoring function developed in our group and the MM/GBSA free energy decomposition analysis were seamlessly integrated into a multi-functional platform. The structures of PPIs were predicted by combining the ATTRACT docking and the HawkRank re-scoring, and the key residues for PPIs were highlighted by the MM/GBSA free energy decomposition. The molecular visualization was supported by 3Dmol.js. For the structural modeling of PPIs, HawkDock could achieve a better performance than ZDOCK 3.0.2 in the benchmark testing. For the prediction of key residues, the important residues that play an essential role in PPIs could be identified in the top 10 residues for ∼81.4% predicted models and ∼95.4% crystal structures in the benchmark dataset. To sum up, the HawkDock server is a powerful tool to predict the binding structures and identify the key residues of PPIs. The HawkDock server is accessible free of charge at http://cadd.zju.edu.cn/hawkdock/.

Identifying protein-protein interface via a novel multi-scale local sequence and structural representation

Background

Protein-protein interaction plays a key role in a multitude of biological processes, such as signal transduction, de novo drug design, immune responses, and enzymatic activities. Gaining insights of various binding abilities can deepen our understanding of the interaction. It is of great interest to understand how proteins in a complex interact with each other. Many efficient methods have been developed for identifying protein-protein interface.

Results

In this paper, we obtain the local information on protein-protein interface, through multi-scale local average block and hexagon structure construction. Given a pair of proteins, we use a trained support vector regression (SVR) model to select best configurations. On Benchmark v4.0, our method achieves average I_rmsd value of 3.28Å and overall F_nat value of 63%, which improves upon I_rmsd of 3.89Å and F_nat of 49% for ZRANK, and I_rmsd of 3.99Å and F_nat of 46% for ClusPro. On CAPRI targets, our method achieves average I_rmsd value of 3.45Å and overall F_nat value of 46%, which improves upon I_rmsd of 4.18Å and F_nat of 40% for ZRANK, and I_rmsd of 5.12Å and F_nat of 32% for ClusPro. The success rates by our method, FRODOCK 2.0, InterEvDock and SnapDock on Benchmark v4.0 are 41.5%, 29.0%, 29.4% and 37.0%, respectively.

Conclusion

Experiments show that our method performs better than some state-of-the-art methods, based on the prediction quality improved in terms of CAPRI evaluation criteria. All these results demonstrate that our method is a valuable technological tool for identifying protein-protein interface.

Chemical Space Overlap with Critical Protein-Protein Interface Residues in Commercial and Specialized Small-Molecule Libraries

There is growing interest in the use of structure-based virtual screening to identify small molecules that inhibit challenging protein-protein interactions (PPIs). In this study, we investigated how effectively chemical library members docked at the PPI interface mimic the position of critical side-chain residues known as "hot spots". Three compound collections were considered, a commercially available screening collection (ChemDiv), a collection of diversity-oriented synthesis (DOS) compounds that contains natural-product-like small molecules, and a library constructed using established reactions (the "screenable chemical universe based on intuitive data organization", SCUBIDOO). Three different tight PPIs for which hot-spot residues have been identified were selected for analysis: uPAR⋅uPA, TEAD4⋅Yap1, and CaV α⋅CaV β. Analysis of library physicochemical properties was followed by docking to the PPI receptors. A pharmacophore method was used to measure overlap between small-molecule substituents and hot-spot side chains. Fragment-like conformationally restricted small molecules showed better hot-spot overlap for interfaces with well-defined pockets such as uPAR⋅uPA, whereas better overlap was observed for more complex DOS compounds in interfaces lacking a well-defined binding site such as TEAD4⋅Yap1. Virtual screening of conformationally restricted compounds targeting uPAR⋅uPA and TEAD4⋅Yap1 followed by experimental validation reinforce these findings, as the best hits were fragment-like and had few rotatable bonds for the former, while no hits were identified for the latter. Overall, such studies provide a framework for understanding PPIs in the context of additional chemical matter and new PPI definitions.

Hot spot prediction in protein-protein interactions by an ensemble system

BACKGROUND

Hot spot residues are functional sites in protein interaction interfaces. The identification of hot spot residues is time-consuming and laborious using experimental methods. In order to address the issue, many computational methods have been developed to predict hot spot residues. Moreover, most prediction methods are based on structural features, sequence characteristics, and/or other protein features.

RESULTS

This paper proposed an ensemble learning method to predict hot spot residues that only uses sequence features and the relative accessible surface area of amino acid sequences. In this work, a novel feature selection technique was developed, an auto-correlation function combined with a sliding window technique was applied to obtain the characteristics of amino acid residues in protein sequence, and an ensemble classifier with SVM and KNN base classifiers was built to achieve the best classification performance.

CONCLUSION

The experimental results showed that our model yields the highest F1 score of 0.92 and an MCC value of 0.87 on ASEdb dataset. Compared with other machine learning methods, our model achieves a big improvement in hot spot prediction.

AVAILABILITY

http://deeplearner.ahu.edu.cn/web/HotspotEL.htm .

Singular value decomposition for the correlation of atomic fluctuations with arbitrary angle

Many proteins exhibit a critical property called allostery, which enables intra-molecular transmission of information between distal sites. Microscopically, allosteric response is closely related to correlated atomic fluctuations. Conventional correlation analysis correlates the atomic fluctuations at two sites by taking the dot product (DP) between the fluctuations, which accounts only for the parallel and antiparallel components. Here, we present a singular value decomposition (SVD) method that analyzes the correlation coefficient of fluctuation dynamics with an arbitrary angle between the correlated directions. In a model allosteric system, the second PDZ domain (PDZ2) in the human PTP1E protein, approximately one third of the strong correlations have near-perpendicular directions, which are underestimated in the conventional method. The discrimination becomes more prominent for residue pairs with larger separation. The results of the proposed SVD method are more consistent with the experimentally determined PDZ2 dynamics than those of conventional method. In addition, the SVD method improved the prediction accuracy of the allosteric sites in a dataset of 23 known allosteric monomer proteins. The proposed method may inspire extended investigation not only into allostery, but also into protein dynamics and drug design.