Alterations in membrane proteins (MPs) and their regulated pathways have been established as cancer hallmarks and extensively targeted in clinical applications. However, the analysis of MP-interacting proteins and downstream pathways across human malignancies remains challenging. Here, we present a systematically integrated method to generate a resource of cancer membrane protein-regulated networks (CaMPNets), containing 63,746 high-confidence protein–protein interactions (PPIs) for 1962 MPs, using expression profiles from 5922 tumors with overall survival outcomes across 15 human cancers. Comprehensive analysis of CaMPNets links MP partner communities and regulated pathways to provide MP-based gene sets for identifying prognostic biomarkers and druggable targets. For example, we identify CHRNA9 with 12 PPIs (e.g., ERBB2) can be a therapeutic target and find its anti-metastasis agent, bupropion, for treatment in nicotine-induced breast cancer. This resource is a study to systematically integrate MP interactions, genomics, and clinical outcomes for helping illuminate cancer-wide atlas and prognostic landscapes in tumor homo/heterogeneity.

The primary type of liver cancer, hepatocellular carcinoma (HCC), has been associated with nonalcoholic steatohepatitis, diabetes, and obesity. Previous studies have identified some genetic risk factors, such as hepatitis B virus X antigens, overexpression of SRC oncogene, and mutation of the p53 tumor suppressor gene; however, the synergism between diet and genetic risk factors is still unclear. To investigate the synergism between diet and genetic risk factors in hepatocarcinogenesis, we used zebrafish with four genetic backgrounds and overfeeding or high-fat-diet-induced obesity with an omics-based expression of genes and histopathological changes. The results show that overfeeding and high-fat diet can induce obesity and nonalcoholic steatohepatitis in wild-type fish. In HBx, Src (p53-) triple transgenic zebrafish, diet-induced obesity accelerated HCC formation at five months of age and increased the cancer incidence threefold. We developed a global omics data analysis method to investigate genes, pathways, and biological systems based on microarray and next-generation sequencing (NGS, RNA-seq) omics data of zebrafish with four diet and genetic risk factors. The results show that two Kyoto Encyclopedia of Genes and Genomes (KEGG) systems, metabolism and genetic information processing, as well as the pathways of fatty acid metabolism, steroid biosynthesis, and ribosome biogenesis, are activated during hepatocarcinogenesis. This study provides a systematic view of the synergism between genetic and diet factors in the dynamic liver cancer formation process, and indicate that overfeeding or a high-fat diet and the risk genes have a synergistic effect in causing liver cancer by affecting fatty acid metabolism and ribosome biogenesis.

One of the crucial steps toward understanding the biological functions of a cellular system is to investigate protein–protein interaction (PPI) networks. As an increasing number of reliable PPIs become available, there is a growing need for discovering PPIs to reconstruct PPI networks of interesting organisms. Some interolog-based methods and homologous PPI families have been proposed for predicting PPIs from the known PPIs of source organisms. Here, we propose a multiple-strategy scoring method to identify reliable PPIs for reconstructing the mouse PPI network from two well-known organisms: human and fly. We firstly identified the PPI candidates of target organisms based on homologous PPIs, sharing significant sequence similarities (joint E-value ≤ 1 × 10⁻⁴⁰), from source organisms using generalized interolog mapping. These PPI candidates were evaluated by our multiple-strategy scoring method, combining sequence similarities, normalized ranks, and conservation scores across multiple organisms. According to 106,825 PPI candidates in yeast derived from human and fly, our scoring method can achieve high prediction accuracy and outperform generalized interolog mapping. Experiment results show that our multiple-strategy score can avoid the influence of the protein family size and length to significantly improve PPI prediction accuracy and reflect the biological functions. In addition, the top-ranked and conserved PPIs are often orthologous/essential interactions and share the functional similarity. Based on these reliable predicted PPIs, we reconstructed a comprehensive mouse PPI network, which is a scale-free network and can reflect the biological functions and high connectivity of 292 KEGG modules, including 216 pathways and 76 structural complexes. Experimental results show that our scoring method can improve the predicting accuracy based on the normalized rank and evolutionary conservation from multiple organisms. Our predicted PPIs share similar biological processes and cellular components, and the reconstructed genome-wide PPI network can reflect network topology and modularity. We believe that our method is useful for inferring reliable PPIs and reconstructing a comprehensive PPI network of an interesting organism.

Inhibiting a specific target in cancer cells and reducing unwanted side effects has become a promising strategy in pancreatic cancer treatment. MAP4K4 is associated with pancreatic cancer development and correlates with poor clinical outcomes. By phosphorylating MKK4, proteins associated with cell apoptosis and survival are translated. Therefore, inhibiting MAP4K4 activity in pancreatic tumours is a new therapeutic strategy. Herein, we performed a structure-based virtual screening to identify MAP4K4 inhibitors and discovered the compound F389-0746 with a potent inhibition (IC₅₀ 120.7 nM). The results of kinase profiling revealed that F389-0746 was highly selective to MAP4K4 and less likely to cause side effects. Results of in vitro experiments showed that F389-0746 significantly suppressed cancer cell growth and viability. Results of in vivo experiments showed that F389-0746 displayed comparable tumour growth inhibition with the group treated with gemcitabine. These findings suggest that F389-0746 has promising potential to be further developed as a novel pancreatic cancer treatment. Graphical abstract

Although many studies have shown the association between smoking and the increased incidence and adverse prognosis of head and neck squamous cell carcinoma (HNSCC), the mechanisms and pharmaceutical targets involved remain unclear. Here, we integrated gene expression signatures, genetic alterations, and survival analyses to identify prognostic indicators and therapeutic targets for smoking HNSCC patients, and we discovered that the FDA-approved drug varenicline inhibits the target for cancer cell migration/invasion. We first identified 18 smoking-related and prognostic genes for HNSCC by using RNA-Seq and clinical follow-up data. One of these genes, CHRNB4 (neuronal acetylcholine receptor subunit beta-4), increased the risk of death by approximately threefold in CHRNB4-high expression smokers compared to CHRNB4-low expression smokers (log rank, p = 0.00042; hazard ratio, 2.82; 95% CI, 1.55–5.14), former smokers, and non-smokers. Furthermore, we examined the functional enrichment of co-regulated genes of CHRNB4 and its 246 frequently occurring copy number alterations (CNAs). We found that these genes were involved in promoting angiogenesis, resisting cell death, and sustaining proliferation, and contributed to much worse outcomes for CHRNB4-high patients. Finally, we performed CHRNB4 gene editing and drug inhibition assays, and the results validate these observations. In summary, our study suggests that CHRNB4 is a prognostic indicator for smoking HNSCC patients and provides a potential new therapeutic drug to prevent recurrence or distant metastasis.

The COVID-19 pandemic has had a widespread impact on a global scale, and the evolution of considerable dominants has already taken place. Some variants contained certain key mutations located on the receptor binding domain (RBD) of spike protein, such as E484K and N501Y. It is increasingly worrying that these variants could impair the efficacy of current vaccines or therapies. Therefore, analyzing and predicting the high-risk mutations of SARS-CoV-2 spike glycoprotein is crucial to design future vaccines against the different variants. In this work, we proposed an in silico approach, immune-escaping score (IES), to predict high-risk immune-escaping hot spots on the receptor-binding domain (RBD), implemented through integrated delta binding free energy measured by computational mutagenesis of spike-antibody complexes and mutation frequency calculated from viral genome sequencing data. We identified 23 potentially immune-escaping mutations on the RBD by using IES, nine of which occurred in omicron variants (R346K, K417N, N440K, L452Q, L452R, S477N, T478K, F490S, and N501Y), despite our dataset being curated before the omicron first appeared. The highest immune-escaping score (IES = 1) was found for E484K, which agrees with recent studies stating that the mutation significantly reduced the efficacy of neutralization antibodies. Furthermore, our predicted delta binding free energy and IES show a high correlation with high-throughput deep mutational scanning data (Pearson’s r = 0.70) and experimentally measured neutralization titers data (mean Pearson’s r = −0.80). In summary, our work presents a new method to identify the potentially immune-escaping mutations on the RBD and provides valuable insights into future COVID-19 vaccine design.