Discovering potential indications of novel or approved drugs is a key step in drug development. Previous computational approaches could be categorized into disease-centric and drug-centric based on the starting point of the issues or small-scaled application and large-scale application according to the diversity of the datasets. Here, a classifier has been constructed to predict the indications of a drug based on the assumption that interactive/associated drugs or drugs with similar structures are more likely to target the same diseases using a large drug indication dataset. To examine the classifier, it was conducted on a dataset with 1,573 drugs retrieved from Comprehensive Medicinal Chemistry database for five times, evaluated by 5-fold cross-validation, yielding five 1st order prediction accuracies that were all approximately 51.48%. Meanwhile, the model yielded an accuracy rate of 50.00% for the 1st order prediction by independent test on a dataset with 32 other drugs in which drug repositioning has been confirmed. Interestingly, some clinically repurposed drug indications that were not included in the datasets are successfully identified by our method. These results suggest that our method may become a useful tool to associate novel molecules with new indications or alternative indications with existing drugs.

Myocardial fibrosis (MF) can cause heart remodeling and it is an independent risk factor for malignant arrhythmias, sudden cardiac death, and other malignant cardiovascular events. It is often characterized by myocardial interstitial collagen deposition and hyperproliferation of cardiac fibroblasts (CFs). The transforming growth factor-β1 (TGF-β1) is the most influential profibrogenic factor. Resveratrol (RSV) is an active polyphenol substance that inhibits myocardial fibrosis. The mechanism of RSV-mediated inhibition of the proliferation of CFs at the microRNA level is not fully understood. We used TGF-β1 to induce CFs proliferation to simulate the pathogenesis of myocardial fibrosis. Neonatal rat CFs were treated with TGF-β1 in the presence or absence of resveratrol. Cell proliferation was measured using the CCK-8 and EdU assay. Collagen secretion was measured using hydroxyproline kit. Further, qPCR analysis was performed to determine microRNA levels after TGF-β1 or resveratrol treatment. To identify the target gene for miR-17, miR-17 was overexpressed or silenced, and the mRNA and protein levels of Smad7 were assessed. The effects of miR-17 silencing or Smad7 overexpression on cell proliferation and collagen secretion were also examined. Resveratrol treatment significantly decreased the TGF-β1-induced CF proliferation and collagen secretion. Resveratrol also decreased the levels of miR-17, miR-34a, and miR-181a in TGF-β1-treated CFs. Overexpression of miR-17 decreased the Smad7 mRNA and protein levels while silencing miR-17 increased them. Additionally, silencing miR-17 or overexpressing Smad7 decreased the TGF-β1-induced CFs proliferation and collagen secretion. In conclusion, resveratrol inhibits TGF-β1-induced CFs proliferation and collagen secretion. This inhibitory effect of resveratrol is orchestrated by the downregulation of miR-17 and the regulation of Smad7.

Drug combinatorial therapy could be more effective in treating some complex diseases than single agents due to better efficacy and reduced side effects. Although some drug combinations are being used, their underlying molecular mechanisms are still poorly understood. Therefore, it is of great interest to deduce a novel drug combination by their molecular mechanisms in a robust and rigorous way. This paper attempts to predict effective drug combinations by a combined consideration of: (1) chemical interaction between drugs, (2) protein interactions between drugs' targets, and (3) target enrichment of KEGG pathways. A benchmark dataset was constructed, consisting of 121 confirmed effective combinations and 605 random combinations. Each drug combination was represented by 465 features derived from the aforementioned three properties. Some feature selection techniques, including Minimum Redundancy Maximum Relevance and Incremental Feature Selection, were adopted to extract the key features. Random forest model was built with its performance evaluated by 5-fold cross-validation. As a result, 55 key features providing the best prediction result were selected. These important features may help to gain insights into the mechanisms of drug combinations, and the proposed prediction model could become a useful tool for screening possible drug combinations.

Accumulating studies demonstrate that dihydromyricetin (DMY), a compound extracted from Chinese traditional herb, Ampelopsis grossedentata, attenuates atherosclerotic process by improvement of endothelial dysfunction. However, the underlying mechanism remains poorly understood. Thus, the aim of this study is to investigate the potential mechanism behind the attenuating effects of DMY on tumor necrosis factor alpha- (TNF-α-) induced endothelial dysfunction. In response to TNF-α, microRNA-21 (miR-21) expression was significantly increased in human umbilical vein endothelial cells (HUVECs), in line with impaired endothelial dysfunction as evidenced by decreased tube formation and migration, endothelial nitric oxide synthase (eNOS) (ser1177) phosphorylation, dimethylarginine dimethylaminohydrolases 1 (DDAH1) expression and metabolic activity, and nitric oxide (NO) concentration as well as increased asymmetric dimethylarginine (ADMA) levels. In contrast, DMY or blockade of miR-21 expression ameliorated endothelial dysfunction in HUVECs treated with TNF-α through downregulation of miR-21 expression, whereas these effects were abolished by overexpression of miR-21. In addition, using a nonspecific NOS inhibitor, L-NAME, also abrogated the attenuating effects of DMY on endothelial dysfunction. Taken together, these data demonstrated that miR-21-mediated DDAH1/ADMA/NO signal pathway plays an important role in TNF-α-induced endothelial dysfunction, and DMY attenuated endothelial dysfunction induced by TNF-α in a miR-21-dependent manner.

Coronaviruses are specific crown-shaped viruses that were first identified in the 1960s, and three typical examples of the most recent coronavirus disease outbreaks include severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), and COVID-19. Particularly, COVID-19 is currently causing a worldwide pandemic, threatening the health of human beings globally. The identification of viral pathogenic mechanisms is important for further developing effective drugs and targeted clinical treatment methods. The delayed revelation of viral infectious mechanisms is currently one of the technical obstacles in the prevention and treatment of infectious diseases. In this study, we proposed a random walk model to identify the potential pathological mechanisms of COVID-19 on a virus-human protein interaction network, and we effectively identified a group of proteins that have already been determined to be potentially important for COVID-19 infection and for similar SARS infections, which help further developing drugs and targeted therapeutic methods against COVID-19. Moreover, we constructed a standard computational workflow for predicting the pathological biomarkers and related pharmacological targets of infectious diseases.

A number of studies have discovered various roles of PAK4 in human tumors, including osteosarcoma. However, the exact role of PAK4 in osteosarcoma and its mechanism have yet to be determined. Therefore, this study focused on interrogating the PAK4 effect on the proliferation and migration ability of osteosarcoma and its underlying mechanisms.

Hepatocellular carcinoma (HCC) is the most common type among primary liver cancers (PLC). With its poor prognosis and survival rate, it is necessary for HCC patients to have a long-term follow-up. We believe that there are currently no relevant reports or literature about nomograms for predicting the cancer-specific mortality of HCC patients. Therefore, the primary goal of this study was to develop and evaluate nomograms to predict cancer-specific mortality and overall mortality. Data of 45,158 cases of HCC patients were collected from the Surveillance, Epidemiology, and End Results (SEER) program database between 2004 and 2013, which were then utilized to develop the nomograms. Finally, the performance of the nomograms was evaluated by the concordance index (C-index) and the area under the time-dependent receiver operating characteristic (ROC) curve (td-AUC). The categories selected to develop a nomogram for predicting cancer-specific mortality included marriage, insurance, radiotherapy, surgery, distant metastasis, lymphatic metastasis, tumor size, grade, sex, and the American Joint Committee on Cancer (AJCC) stage; while the marriage, radiotherapy, surgery, AJCC stage, grade, race, sex, and age were selected to develop a nomogram for predicting overall mortality. The C-indices for predicted 1-, 3-, and 5-year cancer-specific mortality were 0.792, 0.776, and 0.774; the AUC values for 1-, 3-, and 5-year cancer-specific mortality were 0.830, 0.830, and 0.830. The C-indices for predicted 1-, 3-, and 5-year overall mortality were 0.770, 0.755, and 0.752; AUC values for predicted 1-, 3-, and 5-year overall mortality were 0.820, 0.820, and 0.830. The results showed that the nomograms possessed good agreement compared with the observed outcomes. It could provide clinicians with a personalized predicted risk of death information to evaluate the potential changes of the disease-specific condition so that clinicians can adjust therapy options when combined with the actual condition of the patient, which is beneficial to patients.

The mechanisms of lung cancer are highly complex. Not only mRNA gene expression but also microRNAs, DNA methylation, and copy number variation (CNV) play roles in tumorigenesis. It is difficult to incorporate so much information into a single model that can comprehensively reflect all these lung cancer mechanisms. In this study, we analyzed the 129 TCGA (The Cancer Genome Atlas) squamous cell lung carcinoma samples with gene expression, microRNA expression, DNA methylation, and CNV data. First, we used variance inflation factor (VIF) regression to build the whole genome integrative network. Then, we isolated the lung cancer subnetwork by identifying the known lung cancer genes and their direct regulators. This subnetwork was refined by the Bayesian method, and the directed regulations among mRNA genes, microRNAs, methylations, and CNVs were obtained. The novel candidate key drivers in this refined subnetwork, such as the methylation of ARHGDIB and HOXD3, microRNA let-7a and miR-31, and the CNV of AGAP2, were identified and analyzed. On three large public available lung cancer datasets, the key drivers ARHGDIB and HOXD3 demonstrated significant associations with the overall survival of lung cancer patients. Our results provide new insights into lung cancer mechanisms.

Nowadays, pollution levels are rapidly increasing all over the world. One of the most important pollutants is PM2.5. It is known that the pollution environment may cause several problems, such as greenhouse effect and acid rain. Among them, the most important problem is that pollutants can induce a number of serious diseases. Some studies have reported that PM2.5 is an important etiologic factor for lung cancer. In this study, we extensively investigate the associations between PM2.5 and 22 disease classes recommended by Goh et al., such as respiratory diseases, cardiovascular diseases, and gastrointestinal diseases. The protein-protein interactions were used to measure the linkage between disease genes and genes that have been reported to be modulated by PM2.5. The results suggest that some diseases, such as diseases related to ear, nose, and throat and gastrointestinal, nutritional, renal, and cardiovascular diseases, are influenced by PM2.5 and some evidences were provided to confirm our results. For example, a total of 18 genes related to cardiovascular diseases are identified to be closely related to PM2.5, and cardiovascular disease relevant gene DSP is significantly related to PM2.5 gene JUP.

Most drugs have beneficial as well as adverse effects and exert their biological functions by adjusting and altering the functions of their target proteins. Thus, knowledge of drugs target proteins is essential for the improvement of therapeutic effects and mitigation of undesirable side effects. In the study, we proposed a novel prediction method based on drug/compound ontology information extracted from ChEBI to identify drugs target groups from which the kind of functions of a drug may be deduced. By collecting data in KEGG, a benchmark dataset consisting of 876 drugs, categorized into four target groups, was constructed. To evaluate the method more thoroughly, the benchmark dataset was divided into a training dataset and an independent test dataset. It is observed by jackknife test that the overall prediction accuracy on the training dataset was 83.12%, while it was 87.50% on the test dataset-the predictor exhibited an excellent generalization. The good performance of the method indicates that the ontology information of the drugs contains rich information about their target groups, and the study may become an inspiration to solve the problems of this sort and bridge the gap between ChEBI ontology and drugs target groups.

The recently emerging Influenza A/H7N9 virus is reported to be able to infect humans and cause mortality. However, viral and host factors associated with the infection are poorly understood. It is suggested by the "guilt by association" rule that interacting proteins share the same or similar functions and hence may be involved in the same pathway. In this study, we developed a computational method to identify Influenza A/H7N9 virus infection-related human genes based on this rule from the shortest paths in a virus-human protein interaction network. Finally, we screened out the most significant 20 human genes, which could be the potential infection related genes, providing guidelines for further experimental validation. Analysis of the 20 genes showed that they were enriched in protein binding, saccharide or polysaccharide metabolism related pathways and oxidative phosphorylation pathways. We also compared the results with those from human rhinovirus (HRV) and respiratory syncytial virus (RSV) by the same method. It was indicated that saccharide or polysaccharide metabolism related pathways might be especially associated with the H7N9 infection. These results could shed some light on the understanding of the virus infection mechanism, providing basis for future experimental biology studies and for the development of effective strategies for H7N9 clinical therapies.

Identifying disease genes is one of the most important topics in biomedicine and may facilitate studies on the mechanisms underlying disease. Age-related macular degeneration (AMD) is a serious eye disease; it typically affects older adults and results in a loss of vision due to retina damage. In this study, we attempt to develop an effective method for distinguishing AMD-related genes. Gene ontology and KEGG enrichment analyses of known AMD-related genes were performed, and a classification system was established. In detail, each gene was encoded into a vector by extracting enrichment scores of the gene set, including it and its direct neighbors in STRING, and gene ontology terms or KEGG pathways. Then certain feature-selection methods, including minimum redundancy maximum relevance and incremental feature selection, were adopted to extract key features for the classification system. As a result, 720 GO terms and 11 KEGG pathways were deemed the most important factors for predicting AMD-related genes.

Pancreatic cancer is a serious disease that results in more than thirty thousand deaths around the world per year. To design effective treatments, many investigators have devoted themselves to the study of biological processes and mechanisms underlying this disease. However, it is far from complete. In this study, we tried to extract important gene ontology (GO) terms and KEGG pathways for pancreatic cancer by adopting some existing computational methods. Genes that have been validated to be related to pancreatic cancer and have not been validated were represented by features derived from GO terms and KEGG pathways using the enrichment theory. A popular feature selection method, minimum redundancy maximum relevance, was employed to analyze these features and extract important GO terms and KEGG pathways. An extensive analysis of the obtained GO terms and KEGG pathways was provided to confirm the correlations between them and pancreatic cancer.

Human Schlafen 5 (SLFN5) is reported to inhibit or promote the proliferation of several specific types of cancer cells by our lab and other researchers. We are curious about its implications in lung adenocarcinoma (LUAC), a malignant tumor with a high incidence rate and high mortality.

The regulation of transcriptome expression level is a complex process involving multiple-level interactions among molecules such as protein coding RNA (mRNA), long noncoding RNA (lncRNA), and microRNA (miRNA), which are essential for the transcriptome stability and maintenance and regulation of body homeostasis. The availability of multilevel expression data enables a comprehensive view of the regulatory network. In this study, we analyzed the coding and noncoding gene expression profiles of 301 patients with uterine corpus endometrial carcinoma (UCEC). A new method was proposed to construct a genome-wide integrative network based on variance inflation factor (VIF) regression method. The cross-regulation relations of mRNA, lncRNA, and miRNA were then selected based on clique-searching algorithm from the network, when any two molecules of the three were shown as interacting according to the integrative network. Such relation, which we call the mRNA-lncRNA-miRNA triplet, demonstrated the complexity in transcriptome regulation process. Finally, six UCEC-related triplets were selected in which the mRNA participates in endometrial carcinoma pathway, such as CDH1 and TP53. The multi-type RNAs are proved to be cross-regulated as to each of the six triplets according to literature. All the triplets demonstrated the association with the initiation and progression of UCEC. Our method provides a comprehensive strategy for the investigation of transcriptome regulation mechanism.

Pancreatic cancer (PC) is a highly malignant tumor derived from pancreas tissue and is one of the leading causes of death from cancer. Its molecular mechanism has been partially revealed by validating its oncogenes and tumor suppressor genes; however, the available data remain insufficient for medical workers to design effective treatments. Large-scale identification of PC-related genes can promote studies on PC. In this study, we propose a computational method for mining new candidate PC-related genes. A large network was constructed using protein-protein interaction information, and a shortest path approach was applied to mine new candidate genes based on validated PC-related genes. In addition, a permutation test was adopted to further select key candidate genes. Finally, for all discovered candidate genes, the likelihood that the genes are novel PC-related genes is discussed based on their currently known functions.

Isolation and identification of diverse porcine reproductive and respiratory syndrome viruses (PRRSVs) play a fundamental role in PRRSV research and disease management. However, PRRSV has a restricted cell tropism for infection. MARC-145 cells are routinely used for North American genotype PRRSV isolation and vaccine production. But MARC-145 cells have some limitations such as low virus yield. CD163 is a cellular receptor that mediates productive infection of PRRSV in various nonpermissive cell lines. In this study, we established a high and stable porcine CD163- (pCD163-) expressing MARC-145 cell line toward increasing its susceptibility to PRRSV infection. Indirect immunofluorescence assay (IFA) and Western blotting assays showed that pCD163 was expressed higher in pCD163-MARC cell line than MARC-145 cells. Furthermore, the ability of pCD163-MARC cell line to propagate PRRSV was significantly increased as compared with MARC-145 cells. Finally, we found that pCD163-MARC cell line had a higher isolation rate of clinical PRRSV samples and propagated live attenuated PRRS vaccine strains more efficiently than MARC-145 cells. This pCD163-MARC cell line will be a valuable tool for propagation and research of PRRSV.

One of the most important and challenging problems in biomedicine is how to predict the cancer related genes. Retinoblastoma (RB) is the most common primary intraocular malignancy usually occurring in childhood. Early detection of RB could reduce the morbidity and promote the probability of disease-free survival. Therefore, it is of great importance to identify RB genes. In this study, we developed a computational method to predict RB related genes based on Dagging, with the maximum relevance minimum redundancy (mRMR) method followed by incremental feature selection (IFS). 119 RB genes were compiled from two previous RB related studies, while 5,500 non-RB genes were randomly selected from Ensemble genes. Ten datasets were constructed based on all these RB and non-RB genes. Each gene was encoded with a 13,126-dimensional vector including 12,887 Gene Ontology enrichment scores and 239 KEGG enrichment scores. Finally, an optimal feature set including 1061 GO terms and 8 KEGG pathways was obtained. Analysis showed that these features were closely related to RB. It is anticipated that the method can be applied to predict the other cancer related genes as well.

Protein-protein interaction (PPI) plays an extremely remarkable role in the growth, reproduction, and metabolism of all lives. A thorough investigation of PPI can uncover the mechanism of how proteins express their functions. In this study, we used gene ontology (GO) terms and biological pathways to study an extended version of PPI (protein-protein functional associations) and subsequently identify some essential GO terms and pathways that can indicate the difference between two proteins with and without functional associations. The protein-protein functional associations validated by experiments were retrieved from STRING, a well-known database on collected associations between proteins from multiple sources, and they were termed as positive samples. The negative samples were constructed by randomly pairing two proteins. Each sample was represented by several features based on GO and KEGG pathway information of two proteins. Then, the mutual information was adopted to evaluate the importance of all features and some important ones could be accessed, from which a number of essential GO terms or KEGG pathways were identified. The final analysis of some important GO terms and one KEGG pathway can partly uncover the difference between proteins with and without functional associations.

Oncogene is a special type of genes, which can promote the tumor initiation. Good study on oncogenes is helpful for understanding the cause of cancers. Experimental techniques in early time are quite popular in detecting oncogenes. However, their defects become more and more evident in recent years, such as high cost and long time. The newly proposed computational methods provide an alternative way to study oncogenes, which can provide useful clues for further investigations on candidate genes. Considering the limitations of some previous computational methods, such as lack of learning procedures and terming genes as individual subjects, a novel computational method was proposed in this study. The method adopted the features derived from multiple protein networks, viewing proteins in a system level. A classic machine learning algorithm, random forest, was applied on these features to capture the essential characteristic of oncogenes, thereby building the prediction model. All genes except validated oncogenes were ranked with a measurement yielded by the prediction model. Top genes were quite different from potential oncogenes discovered by previous methods, and they can be confirmed to become novel oncogenes. It was indicated that the newly identified genes can be essential supplements for previous results.