Prediction of interface residue based on the features of residue interaction network

Fusion Genes Landscape of Lung Cancer Patients From Inner Mongolia, China

Lung cancer is the leading cause of cancer-related deaths globally. Gene fusion, a key driver of tumorigenesis, has led to the identification of numerous driver gene fusions for lung cancer diagnosis and treatment. However, previous studies focused on Western populations, leaving the possibility of unrecognized lung cancer-associated gene fusions specific to Inner Mongolia due to its unique genetic background and dietary habits. To address this, we conducted DNA sequencing analysis on tumor and adjacent nontumor tissues from 1200 individuals with lung cancer in Inner Mongolia. Our analysis established a comprehensive fusion gene landscape specific to lung cancer in Inner Mongolia, shedding light on potential region-specific molecular mechanisms underlying the disease. Compared to Western cohorts, we observed a higher occurrence of ALK and RET fusions in Inner Mongolian patients. Additionally, we discovered eight novel fusion genes in three patients: SLC34A2-EPHB1, CCT6P3-GSTP1, BARHL2-APC, HRAS-MELK, FAM134B-ERBB2, ABCB1-GIPC1, GPR98-ALK, and FAM134B-SALL1. These previously unreported fusion genes suggest potential regional specificity. Furthermore, we characterized the fusion genes' structures based on breakpoints and described their impact on major functional gene domains. Importantly, the identified novel fusion genes exhibited significant clinical and pathological relevance. Notably, patients with SLC34A2-EPHB1, CCT6P3-GSTP1, and BARHL2-APC fusions showed sensitivity to the combination of chemotherapy and immunotherapy. Patients with HRAS-MELK, FAM134B-ERBB2, and ABCB1-GIPC1 fusions showed sensitivity to chemotherapy. In summary, our study provides novel insights into the frequency, distribution, and characteristics of specific fusion genes, offering valuable guidance for the development of effective clinical treatments, particularly in Inner Mongolia.

Prioritization of Trypanosoma brucei editosome protein interactions interfaces at residue resolution through proteome-scale network analysis

BACKGROUND

Trypanosoma brucei is the causative agent for trypanosomiasis in humans and livestock, which presents a growing challenge due to drug resistance. While identifying novel drug targets is vital, the process is delayed due to a lack of functional information on many of the pathogen's proteins. Accordingly, this paper presents a computational framework for prioritizing drug targets within the editosome, a vital molecular machinery responsible for mitochondrial RNA processing in T. brucei. Importantly, this framework may eliminate the need for prior gene or protein characterization, potentially accelerating drug discovery efforts.

RESULTS

By integrating protein-protein interaction (PPI) network analysis, PPI structural modeling, and residue interaction network (RIN) analysis, we quantitatively ranked and identified top hub editosome proteins, their key interaction interfaces, and hotspot residues. Our findings were cross-validated and further prioritized by incorporating them into gene set analysis and differential expression analysis of existing quantitative proteomics data across various life stages of T. brucei. In doing so, we highlighted PPIs such as KREL2-KREPA1, RESC2-RESC1, RESC12A-RESC13, and RESC10-RESC6 as top candidates for further investigation. This includes examining their interfaces and hotspot residues, which could guide drug candidate selection and functional studies.

CONCLUSION

RNA editing offers promise for target-based drug discovery, particularly with proteins and interfaces that play central roles in the pathogen's life cycle. This study introduces an integrative drug target identification workflow combining information from the PPI network, PPI 3D structure, and reside-level information of their interface which can be applicable to diverse pathogens. In the case of T. brucei, via this pipeline, the present study suggested potential drug targets with residue-resolution from RNA editing machinery. However, experimental validation is needed to fully realize its potential in advancing urgently needed antiparasitic drug development.

Discovery of ellagic acid as a competitive inhibitor of Src homology phosphotyrosyl phosphatase 2 (SHP2) for cancer treatment: In vitro and in silico study

Targeting SHP2 has become a potential cancer treatment strategy. In this study, ellagic acid was first reported as a competitive inhibitor of SHP2, with an IC50 value of 0.69 ± 0.07 μM, and its inhibitory potency was 34.86 times higher that of the positive control NSC87877. Ellagic acid also had high inhibitory activity on the SHP2-E76K and SHP2-E76A mutants, with the IC50 values of 1.55 ± 0.17 μM and 0.39 ± 0.05 μM, respectively. Besides, the IC50 values of ellagic acid on homologous proteins SHP1, PTP1B, and TCPTP were 0.93 ± 0.08 μM, 2.04 ± 0.28 μM, and 11.79 ± 0.83 μM, with selectivity of 1.35, 2.96, and 17.09 times, respectively. The CCK8 proliferation experiment exhibited that ellagic acid would inhibit the proliferation of various cancer cells. It was worth noting that the combination of ellagic acid and KRASG12C inhibitor AMG510 would produce a strong synergistic effect in inhibiting NCI-H358 cells. Western blot experiment exhibited that ellagic acid would downregulate the phosphorylation levels of Erk and Akt in NCI-H358 and MDA-MB-468 cells. Molecular docking and molecular dynamics studies revealed the binding information between SHP2 and ellagic acid. In summary, this study provides new ideas for the development of SHP2 inhibitors.

Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding

Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.

Scoring of protein–protein docking models utilizing predicted interface residues

Scoring docking solutions is a difficult task, and many methods have been developed for this purpose. In docking, only a handful of the hundreds of thousands of models generated by docking algorithms are acceptable, causing difficulties when developing scoring functions. Today's best scoring functions can significantly increase the number of top‐ranked models but still fail for most targets. Here, we examine the possibility of utilizing predicted interface residues to score docking models generated during the scan stage of a docking algorithm. Many methods have been developed to infer the regions of a protein surface that interact with another protein, but most have not been benchmarked using docking algorithms. This study systematically tests different interface prediction methods for scoring >300.000 low‐resolution rigid‐body template free docking decoys. Overall we find that contact‐based interface prediction by BIPSPI is the best method to score docking solutions, with >12% of first ranked docking models being acceptable. Additional experiments indicated precision as a high‐importance metric when estimating interface prediction quality, focusing on docking constraints production. Finally, we discussed several limitations for adopting interface predictions as constraints in a docking protocol.

Exploring the cause of the dual allosteric targeted inhibition attaching to allosteric sites enhancing SHP2 inhibition

SHP2 is a protein tyrosine phosphatase (PTP) that can regulate the tyrosine phosphorylation level. Overexpression of SHP2 will promote the development of cancer diseases, so SHP2 has become one of the popular targets for the treatment of cancer. Studies have reported that both SHP099 and SHP844 are inhibitors of SHP2 and bind to different allosteric sites 1 and 2, respectively. Studies have shown that combining SHP099 with SHP844 will enhance pharmacological pathway inhibition in cells. This study uses molecular dynamic simulations to explore the dual allosteric targeted inhibition mechanism. The result shows that the residues THR108-TRP112 (allosteric site 1) move to LEU236-GLN245 (αB-αC link loop in PTP domain) , the residues of GLN79-GLN87 (allosteric site 2) get close to LEU262-GLN269 (αA-αB link loop in PTP domain) and HIS458-ARG465 (P-loop) come near to ARG501-THR507 (Q-loop) in SHP2-SHP099-SHP844 system, which makes the "inactive conformation" more stable and prevents the substrate from entering the catalytic site. Meanwhile, residue GLU110 (allosteric site 1), ARG265 (allosteric site 2), and ARG501 (Q-loop) are speculated to be the key residues that causing the SHP2 protein in auto-inhibition conformation. It is hoped that this study will provide clues for the development of the dual allosteric targeted inhibition of SHP2.

Centrality Measures in Residue Interaction Networks to Highlight Amino Acids in Protein–Protein Binding

Residue interaction networks (RINs) describe a protein structure as a network of interacting residues. Central nodes in these networks, identified by centrality analyses, highlight those residues that play a role in the structure and function of the protein. However, little is known about the capability of such analyses to identify residues involved in the formation of macromolecular complexes. Here, we performed six different centrality measures on the RINs generated from the complexes of the SKEMPI 2 database of changes in protein–protein binding upon mutation in order to evaluate the capability of each of these measures to identify major binding residues. The analyses were performed with and without the crystallographic water molecules, in addition to the protein residues. We also investigated the use of a weight factor based on the inter-residue distances to improve the detection of these residues. We show that for the identification of major binding residues, closeness, degree, and PageRank result in good precision, whereas betweenness, eigenvector, and residue centrality analyses give a higher sensitivity. Including water in the analysis improves the sensitivity of all measures without losing precision. Applying weights only slightly raises the sensitivity of eigenvector centrality analysis. We finally show that a combination of multiple centrality analyses is the optimal approach to identify residues that play a role in protein–protein interaction.

Role of protein-protein interactions in allosteric drug design for DNA methyltransferases

DNA methyltransferases (DNMTs) not only play key roles in epigenetic gene regulation, but also serve as emerging targets for several diseases, especially for cancers. Due to the multi-domains of DNMT structures, targeting allosteric sites of protein-protein interactions (PPIs) is becoming an attractive strategy in epigenetic drug discovery. This chapter aims to review the major contemporary approaches utilized for the drug discovery based on PPIs in different dimensions, from the enumeration of allosteric mechanism to the identification of allosteric pockets. These include the construction of protein structure networks (PSNs) based on molecular dynamics (MD) simulations, performing elastic network models (ENMs) and perturbation response scanning (PRS) calculation, the sequence-based conservation and coupling analysis, and the allosteric pockets identification. Furthermore, we complement this methodology by highlighting the role of computational approaches in promising practical applications for the computer-aided drug design, with special focus on two DNMTs, namely, DNMT1 and DNMT3A.

Identification of novel inhibitor of protein tyrosine phosphatases delta: structure-based pharmacophore modeling, virtual screening, flexible docking, molecular dynamics simulation, and post-molecular dynamics analysis

Owing to their unique functions in regulating the synapse activity of protein tyrosine phosphatases delta (PTPδ) that has drawn special attention for developing drugs to autism spectrum disorders (ASDs). In this study, the PTPδ pharmacophore was first established by the structure-based pharmacophore method. Subsequently, 10 compounds contented Lipinski's rule of five was acquired by the virtual screening of the PTPδ pharmacophore against ZINC and PubChem databases. Then, the 10 identified molecules were discovered that had better binding affinity than a known PTPδ inhibitors compound SCHEMBL16375396. Two compounds SCHEMBL16375408 and ZINC19796658 with high binding score, low toxicity were gained. They were observed by docking analysis and molecular dynamics simulations that the novel potential inhibitors not only possessed the same function as SCHEMBL16375396 did in inhibiting PTPδ, but also had more favorable conformation to bind with the catalytic active regions. This study provides a new method for identify PTPδ inhibitor for the treatment of ASDs disease.Communicated by Ramaswamy H. Sarma.

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.

FusionGDB: fusion gene annotation DataBase

Gene fusion is one of the hallmarks of cancer genome via chromosomal rearrangement initiated by DNA double-strand breakage. To date, many fusion genes (FGs) have been established as important biomarkers and therapeutic targets in multiple cancer types. To better understand the function of FGs in cancer types and to promote the discovery of clinically relevant FGs, we built FusionGDB (Fusion Gene annotation DataBase) available at https://ccsm.uth.edu/FusionGDB. We collected 48 117 FGs across pan-cancer from three representative fusion gene resources: the improved database of chimeric transcripts and RNA-seq data (ChiTaRS 3.1), an integrative resource for cancer-associated transcript fusions (TumorFusions), and The Cancer Genome Atlas (TCGA) fusions by Gao et al. For these ∼48K FGs, we performed functional annotations including gene assessment across pan-cancer fusion genes, open reading frame (ORF) assignment, and retention search of 39 protein features based on gene structures of multiple isoforms with different breakpoints. We also provided the fusion transcript and amino acid sequences according to multiple breakpoints and transcript isoforms. Our analyses identified 331, 303 and 667 in-frame FGs with retaining kinase, DNA-binding, and epigenetic factor domains, respectively, as well as 976 FGs lost protein-protein interaction. FusionGDB provides six categories of annotations: FusionGeneSummary, FusionProtFeature, FusionGeneSequence, FusionGenePPI, RelatedDrug and RelatedDisease.