Alignment of virus-host protein-protein interaction networks by integer linear programming: SARS-CoV-2

Deciphering the similarities and disparities of molecular mechanisms behind respiratory epithelium response to HCoV-229E and SARS-CoV-2 and drug repurposing, a systems biology approach

PURPOSE

Identifying the molecular mechanisms behind SARS-CoV-2 disparities and similarities will help find new treatments. The present study determines networks' shared and non-shared (specific) crucial elements in response to HCoV-229E and SARS-CoV-2 viruses to recommend candidate medications.

METHODS

We retrieved the omics data on respiratory cells infected with HCoV-229E and SARS-CoV-2, constructed PPIN and GRN, and detected clusters and motifs. Using a drug-gene interaction network, we determined the similarities and disparities of mechanisms behind their host response and drug-repurposed.

RESULTS

CXCL1, KLHL21, SMAD3, HIF1A, and STAT1 were the shared DEGs between both viruses' protein-protein interaction network (PPIN) and gene regulatory network (GRN). The NPM1 was a specific critical node for HCoV-229E and was a Hub-Bottleneck shared between PPI and GRN in HCoV-229E. The HLA-F, ADCY5, TRIM14, RPF1, and FGA were the seed proteins in subnetworks of the SARS-CoV-2 PPI network, and HSPA1A and RPL26 proteins were the seed in subnetworks of the PPI network of HCOV-229E. TRIM14, STAT2, and HLA-F played the same role for SARS-CoV-2. Top enriched KEGG pathways included cell cycle and proteasome in HCoV-229E and RIG-I-like receptor, Chemokine, Cytokine-cytokine, NOD-like receptor, and TNF signaling pathways in SARS-CoV-2. We suggest some candidate medications for COVID-19 patient lungs, including Noscapine, Isoetharine mesylate, Cycloserine, Ethamsylate, Cetylpyridinium, Tretinoin, Ixazomib, Vorinostat, Venetoclax, Vorinostat, Ixazomib, Venetoclax, and epoetin alfa for further in-vitro and in-vivo investigations.

CONCLUSION

We suggested CXCL1, KLHL21, SMAD3, HIF1A, and STAT1, ADCY5, TRIM14, RPF1, and FGA, STAT2, and HLA-F as critical genes and Cetylpyridinium, Cycloserine, Noscapine, Ethamsylate, Epoetin alfa, Isoetharine mesylate, Ribavirin, and Tretinoin drugs to study further their importance in treating COVID-19 lung complications.

Prediction of viral-host interactions of COVID-19 by computational methods

Experimental approaches are currently used to determine viral-host interactions, but these approaches are both time-consuming and costly. For these reasons, computational-based approaches are recommended. In this study, using computational-based approaches, viral-host interactions of SARS-CoV-2 virus and human proteins were predicted. The study consists of four different stages; in the first stage viral and host protein sequences were obtained. In the second stage, protein sequences were converted into numerical expressions by various protein mapping methods. These methods are entropy-based, AVL-tree, FIBHASH, binary encoding, CPNR, PAM250, BLOSUM62, Atchley factors, Meiler parameters, EIIP, AESNN1, Miyazawa energies, Micheletti potentials, Z-scale, and hydrophobicity. In the third stage, a deep learning model was designed and BiLSTM was used for this. In the last stage, the protein sequences were classified, and the viral-host interactions were predicted. The performances of protein mapping methods were determined by accuracy, F1-score, specificity, sensitivity, and AUC scores. According to the classification results, the best classification process was obtained by the entropy-based method. With this method, 94.74% accuracy, and 0.95 AUC score were calculated. Then, the most successful classification process was performed with the Z-scale and 91.23% accuracy, and 0.96 AUC score were obtained. Although other protein mapping methods are not as efficient as Z-scale and entropy-based methods, they have achieved successful classification. AVL-tree, FIBHASH, binary encoding, CPNR, PAM250, BLOSUM62, Atchley factors, Meiler parameters and AESNN1 methods showed over 80% accuracy, F1-score, and AUC score. Accuracy scores of EIIP, Miyazawa energies, Micheletti potentials and hydrophobicity methods remained below 80%. When the results were examined in general, it was observed that the computational approaches were successful in predicting viral-host interactions between SARS-CoV-2 virus and human proteins.