Prediction of Drug-Target Interactions with High-Quality Negative Samples and A Network-Based Deep Learning Framework.

Identification of drug-target interactions (DTIs) plays a crucial role in drug discovery. Compared to traditional experimental methods, computer-based methods for predicting DTIs can significantly reduce the time and financial burdens of drug development. In recent years, numerous machine learning-based methods have been proposed for predicting potential DTIs. However, a common limitation among these methods is the absence of high-quality negative samples. Moreover, the effective extraction of multisource information of drugs and proteins for DTI prediction remains a significant challenge. In this paper, we investigated two aspects: the selection of high-quality negative samples and the construction of a high-performance DTI prediction framework. Specifically, we found two types of hidden biases when randomly selecting negative samples from unlabeled drug-protein pairs and proposed a negative sample selection approach based on complex network theory. Furthermore, we proposed a novel DTI prediction method named HNetPa-DTI, which integrates topological information from the drug-protein-disease heterogeneous network and gene ontology (GO) and pathway annotation information of proteins. Specifically, we extracted topological information of the drug-protein-disease heterogeneous network using heterogeneous graph neural networks, and obtained GO and pathway annotation information of proteins from the GO term semantic similarity networks, GO term-protein bipartite networks, and pathway-protein bipartite network using graph neural networks. Experimental results show that HNetPa-DTI outperforms the baseline methods on four types of prediction tasks, demonstrating the superiority of our method. Our code and datasets are available at https://github.com/study-czx/HNetPa-DTI.

GCNA-Cluster: A Gene Co-Expression Network Alignment to Cluster Cancer Patients Algorithm for Identifying Subtypes of Pancreatic Ductal Adenocarcinoma.

Cancer heterogeneity makes it necessary to use different treatment strategies for patients with the same pathological features. Accurate identification of cancer subtypes is a crucial step in this approach. The current studies of pancreatic ductal adenocarcinoma (PDAC) subtypes mainly focus on single genes and ignore the synergistic effects of genes. Here we proposed a network alignment algorithm GCNA-cluster to cluster patients based on gene co-expression networks. We constructed weighted gene co-expression networks for patients and aligned the networks of two patients to estimate the similarity of patients and their cancer subtypes. A scoring function is defined to measure the network alignment result and the score can indicate the similarity between patients. Then, the patients are clustered based on their similarities. We validated the accuracy of the algorithm on the GEO-PDAC dataset with real labels, and the experimental results show that the GCNA-cluster algorithm has better results than classical cancer subtyping algorithms. In addition, the GCNA-cluster algorithm applied to the TCGA-PDAC dataset identified two subtypes based on the Silhouette Coefficient. Biomarkers identified for the PDAC subtypes hint to cell growth, cell cycle or apoptosis as targets for new therapeutic strategies.

Combining non-negative matrix factorization with graph Laplacian regularization for predicting drug-miRNA associations based on multi-source information fusion.

Increasing evidences suggest that miRNAs play a key role in the occurrence and progression of many complex human diseases. Therefore, targeting dysregulated miRNAs with small molecule drugs in the clinical has become a new treatment. Nevertheless, it is high cost and time-consuming for identifying miRNAs-targeted with drugs by biological experiments. Thus, more reliable computational method for identification associations of drugs with miRNAs urgently need to be developed. In this study, we proposed an efficient method, called GNMFDMA, to predict potential associations of drug with miRNA by combining graph Laplacian regularization with non-negative matrix factorization. We first calculated the overall similarity matrices of drugs and miRNAs according to the collected different biological information. Subsequently, the new drug-miRNA association adjacency matrix was reformulated based on the K nearest neighbor profiles so as to put right the false negative associations. Finally, graph Laplacian regularization collaborative non-negative matrix factorization was used to calculate the association scores of drugs with miRNAs. In the cross validation, GNMFDMA obtains AUC of 0.9193, which outperformed the existing methods. In addition, case studies on three common drugs (i.e., 5-Aza-CdR, 5-FU and Gemcitabine), 30, 31 and 34 of the top-50 associations inferred by GNMFDMA were verified. These results reveal that GNMFDMA is a reliable and efficient computational approach for identifying the potential drug-miRNA associations.

A weighted non-negative matrix factorization approach to predict potential associations between drug and disease.

BACKGROUND: Associations of drugs with diseases provide important information for expediting drug development. Due to the number of known drug-disease associations is still insufficient, and considering that inferring associations between them through traditional in vitro experiments is time-consuming and costly. Therefore, more accurate and reliable computational methods urgent need to be developed to predict potential associations of drugs with diseases. METHODS: In this study, we present the model called weighted graph regularized collaborative non-negative matrix factorization for drug-disease association prediction (WNMFDDA). More specifically, we first calculated the drug similarity and disease similarity based on the chemical structures of drugs and medical description information of diseases, respectively. Then, to extend the model to work for new drugs and diseases, weighted [Formula: see text] nearest neighbor was used as a preprocessing step to reconstruct the interaction score profiles of drugs with diseases. Finally, a graph regularized non-negative matrix factorization model was used to identify potential associations between drug and disease. RESULTS: During the cross-validation process, WNMFDDA achieved the AUC values of 0.939 and 0.952 on Fdataset and Cdataset under ten-fold cross validation, respectively, which outperforms other competing prediction methods. Moreover, case studies for several drugs and diseases were carried out to further verify the predictive performance of WNMFDDA. As a result, 13(Doxorubicin), 13(Amiodarone), 12(Obesity) and 12(Asthma) of the top 15 corresponding candidate diseases or drugs were confirmed by existing databases. CONCLUSIONS: The experimental results adequately demonstrated that WNMFDDA is a very effective method for drug-disease association prediction. We believe that WNMFDDA is helpful for relevant biomedical researchers in follow-up studies.

SPCMLMI: A structural perturbation-based matrix completion method to predict lncRNA-miRNA interactions.

Accumulating evidence indicated that the interaction between lncRNA and miRNA is crucial for gene regulation, which can regulate gene transcription, further affecting the occurrence and development of many complex diseases. Accurate identification of interactions between lncRNAs and miRNAs is helpful for the diagnosis and therapeutics of complex diseases. However, the number of known interactions of lncRNA with miRNA is still very limited, and identifying their interactions through biological experiments is time-consuming and expensive. There is an urgent need to develop more accurate and efficient computational methods to infer lncRNA-miRNA interactions. In this work, we developed a matrix completion approach based on structural perturbation to infer lncRNA-miRNA interactions (SPCMLMI). Specifically, we first calculated the similarities of lncRNA and miRNA, including the lncRNA expression profile similarity, miRNA expression profile similarity, lncRNA sequence similarity, and miRNA sequence similarity. Second, a bilayer network was constructed by integrating the known interaction network, lncRNA similarity network, and miRNA similarity network. Finally, a structural perturbation-based matrix completion method was used to predict potential interactions of lncRNA with miRNA. To evaluate the prediction performance of SPCMLMI, five-fold cross validation and a series of comparison experiments were implemented. SPCMLMI achieved AUCs of 0.8984 and 0.9891 on two different datasets, which is superior to other compared methods. Case studies for lncRNA XIST and miRNA hsa-mir-195-5-p further confirmed the effectiveness of our method in inferring lncRNA-miRNA interactions. Furthermore, we found that the structural consistency of the bilayer network was higher than that of other related networks. The results suggest that SPCMLMI can be used as a useful tool to predict interactions between lncRNAs and miRNAs.

Characteristics, influence factors, and health risk assessment of volatile organic compounds through one year of high-resolution measurement at a refinery.

From January 2020 to December 2020, high-resolution data of volatile organic compound (VOC) concentrations were monitored by online instruments at a petroleum refinery. The measurement results showed that the external contaminants, meteorological conditions and photochemical reactions had a great influence on the VOC data measured in the petroleum refineries. Some significant differences were observed in the emission composition of different refineries, while propene (34.2%), propane (10.2%), n-butane (5.6%), i-pentane (5.0%) were the dominant species emitted from the refinery in this study. The correlations between compounds with similar atmospheric lifetimes were strong (R2 > 0.9), which indicated that the diagnostic ratios of these compounds could be used as indicators to identify the refinery emission source. Chronic health effects of non-carcinogenic risk results showed that acrolein had the highest non-carcinogenic risk and other compound-specific health risks may be of less concern in the refining area. Halogenates and aromatics accounted for 97.4% of the total carcinogenic risk values, while 1,2-dibromoethane, chloromethane, benzene, trichloromethane, 1,2-dichloroethane contributed approximately 80% of the total carcinogenic risk assessment values. This research has recorded valuable data about the VOC emission characteristics from the perspective of the high-resolution monitoring of the petroleum refinery. The results of this work will provide a reference to accurately quantify and identify the emission of petroleum refineries and further throw some light on effective VOC abatement strategies.

Condition-Specific Molecular Network Analysis Revealed That Flagellar Proteins Are Involved in Electron Transfer Processes of Shewanella piezotolerans WP3.

Because of the ability to metabolize a large number of electron acceptors such as nitrate, nitrite, fumarate, and metal oxides, Shewanella species have attracted much attention in recent years. Generally, the use of these electron acceptors is mainly achieved through electron transfer proteins and their interactions which will dynamically change across different environmental conditions in cells. Therefore, functional analysis of condition-specific molecular networks can reveal biological information on electron transfer processes. By integrating expression data and molecular networks, we constructed condition-specific molecular networks for Shewanella piezotolerans WP3. We then identified condition-specific key genes and studied their potential functions with an emphasis on their roles in electron transfer processes. Functional module analysis showed that different flagellar assembly modules appeared under these conditions and suggested that flagellar proteins are important for these conditions. We also identified the electron transfer modules underlying these various environmental conditions. The present results could help with screening electron transfer genes and understanding electron transfer processes under various environmental conditions for the Shewanella species.

Genome-scale mutant fitness reveals versatile c-type cytochromes in Shewanella oneidensis MR-1.

Shewanella has been widely investigated for its metabolic versatility and use of a large number of extracellular electron acceptors. Many c-type cytochromes are responsible for this diversity, mainly in condition-specific fashions. By using genome-scale mutant fitness data, we studied which genes (particularly c-type cytochromes) were used to coordinate various electron transfer processes in the present work. First, by integrating fitness profiles with protein-protein interaction (PPI) networks, we showed that the genes with a high total fitness value were generally more important in PPI networks than those with low fitness values. Then, we identified genes that are important across many experiments, and further fitness analysis confirmed five versatile c-type cytochromes: ScyA (SO0264), PetC (SO0610), CcoP (SO2361), CcoO (SO2363) and CytcB (SO4666), which are considered to be crucial in most experimental conditions. Finally, we demonstrated a mediating role in the periplasm for the less-reported CytcB by combining protein structure, subcellular localization and disordered region analysis. Comparative genome analysis further revealed that it is distinctive in Shewanella species. Collectively, these results suggest that periplasmic electron transfer processes are more diverse and flexible than previously reported, giving insight for further experimental studies of Shewanella oneidensis MR-1.

Transcriptional regulatory module analysis reveals that bridge proteins reconcile multiple signals in extracellular electron transfer pathways.

Shewanella oneidensis MR-1 shows remarkable respiratory versatility with a large variety of extracellular electron acceptors (termed extracellular electron transfer, EET). To utilize the various electron acceptors, the bacterium must employ complex regulatory mechanisms to elicit the relevant EET pathways. To investigate the relevant mechanisms, we integrated EET genes and related transcriptional factors (TFs) into transcriptional regulatory modules (TRMs) and showed that many bridge proteins in these modules were signal proteins, which generally contained one or more signal processing domains (eg, GGDEF, EAL, PAS, etc.). Since Shewanella has to respond to diverse environmental conditions despite encoding few EET-relevant TFs, the overabundant signal proteins involved in the TRMs can help decipher the mechanism by which these microbes elicit a wide range of condition-specific responses. By combining proteomic data and protein bioinformatic analysis, we demonstrated that diverse signal proteins reconciled the different EET pathways, and we discussed the functional roles of signal proteins involved in the well-known MtrCAB pathway. Additionally, we showed that the signal proteins SO_2145 and SO_1417 played central roles in triggering EET pathways in anaerobic environments. Taken together, our results suggest that signal proteins have a profound impact on the transcriptional regulation of EET genes and thus have potential applications in microbial fuel cells.

A Comparative Study of Network Motifs in the Integrated Transcriptional Regulation and Protein Interaction Networks of Shewanella.

The Shewanella species shows a remarkable respiratory versatility with a great variety of extracellular electron acceptors (termed Extracellular Electron Transfer, EET). To explore relevant mechanisms from the network motif view, we constructed the integrated networks that combined transcriptional regulation interactions (TRIs) and protein-protein interactions (PPIs) for 13 Shewanella species, identified and compared the network motifs in these integrated networks. We found that the network motifs were evolutionary conserved in these integrated networks. The functional significance of the highly conserved motifs was discussed, especially the important ones that were potentially involved in the Shewanella EET processes. More importantly, we found that: 1) the motif co-regulated PPI took a role in the "standby mode" of protein utilization, which will be helpful for cells to rapidly response to environmental changes; and 2) the type II cofactors, which involved in the motif TRI interacting with a third protein, mainly carried out a signalling role in Shewanella oneidensis MR-1.

Network-Based Methods for Identifying Key Active Proteins in the Extracellular Electron Transfer Process in Shewanella oneidensis MR-1.

Shewanella oneidensis MR-1 can transfer electrons from the intracellular environment to the extracellular space of the cells to reduce the extracellular insoluble electron acceptors (Extracellular Electron Transfer, EET). Benefiting from this EET capability, Shewanella has been widely used in different areas, such as energy production, wastewater treatment, and bioremediation. Genome-wide proteomics data was used to determine the active proteins involved in activating the EET process. We identified 1012 proteins with decreased expression and 811 proteins with increased expression when the EET process changed from inactivation to activation. We then networked these proteins to construct the active protein networks, and identified the top 20 key active proteins by network centralization analysis, including metabolism- and energy-related proteins, signal and transcriptional regulatory proteins, translation-related proteins, and the EET-related proteins. We also constructed the integrated protein interaction and transcriptional regulatory networks for the active proteins, then found three exclusive active network motifs involved in activating the EET process-Bi-feedforward Loop, Regulatory Cascade with a Feedback, and Feedback with a Protein-Protein Interaction (PPI)-and identified the active proteins involved in these motifs. Both enrichment analysis and comparative analysis to the whole-genome data implicated the multiheme c-type cytochromes and multiple signal processing proteins involved in the process. Furthermore, the interactions of these motif-guided active proteins and the involved functional modules were discussed. Collectively, by using network-based methods, this work reported a proteome-wide search for the key active proteins that potentially activate the EET process.

Predicting and Interpreting the Structure of Type IV Pilus of Electricigens by Molecular Dynamics Simulations.

Nanowires that transfer electrons to extracellular acceptors are important in organic matter degradation and nutrient cycling in the environment. Geobacter pili of the group of Type IV pilus are regarded as nanowire-like biological structures. However, determination of the structure of pili remains challenging due to the insolubility of monomers, presence of surface appendages, heterogeneity of the assembly, and low-resolution of electron microscopy techniques. Our previous study provided a method to predict structures for Type IV pili. In this work, we improved on our previous method using molecular dynamics simulations to optimize structures of Neisseria gonorrhoeae (GC), Neisseria meningitidis and Geobacter uraniireducens pilus. Comparison between the predicted structures for GC and Neisseria meningitidis pilus and their native structures revealed that proposed method could predict Type IV pilus successfully. According to the predicted structures, the structural basis for conductivity in G.uraniireducens pili was attributed to the three N-terminal aromatic amino acids. The aromatics were interspersed within the regions of charged amino acids, which may influence the configuration of the aromatic contacts and the rate of electron transfer. These results will supplement experimental research into the mechanism of long-rang electron transport along pili of electricigens.

The accessible chromatin landscape during conversion of human embryonic stem cells to trophoblast by bone morphogenetic protein 4.

Human embryonic stem cells (hESCs) exposed to the growth factor bone morphogenetic protein 4 (BMP4) in the absence of FGF2 have been used as a model to study the development of placental development. However, little is known about the cis-regulatory mechanisms underlying this important process. In this study, we used the public available chromatin accessibility data of hESC H1 cells and BMP4-induced trophoblast (TB) cell lines to identify DNase I hypersensitive sites (DHSs) in the two cell lines, as well as the transcription factor (TF) binding sites within the DHSs. By comparing read profiles in H1 and TB, we identified 17 472 TB-specific DHSs. The TB-specific DHSs are enriched in terms of "blood vessel" and "trophectoderm," consisting of TF motifs family: Leucine Zipper, Helix-Loop-Helix, GATA, and ETS. To validate differential expression of the TFs binding to these motifs, we analyzed public available RNA-seq and microarray data in the same context. Finally, by integrating the protein-protein interaction data, we constructed a TF network for placenta development and identified top 20 key TFs through centrality analysis in the network. Our results indicate BMP4-induced TB system provided an invaluable model for the study of TB development and highlighted novel candidate genes in placenta development in human.

K-shell Analysis Reveals Distinct Functional Parts in an Electron Transfer Network and Its Implications for Extracellular Electron Transfer.

Shewanella oneidensis MR-1 is capable of extracellular electron transfer (EET) and hence has attracted considerable attention. The EET pathways mainly consist of c-type cytochromes, along with some other proteins involved in electron transfer processes. By whole genome study and protein interactions inquisition, we constructed a large-scale electron transfer network containing 2276 interactions among 454 electron transfer related proteins in S. oneidensis MR-1. Using the k-shell decomposition method, we identified and analyzed distinct parts of the electron transfer network. We found that there was a negative correlation between the k s (k-shell values) and the average DR_100 (disordered regions per 100 amino acids) in every shell, which suggested that disordered regions of proteins played an important role during the formation and extension of the electron transfer network. Furthermore, proteins in the top three shells of the network are mainly located in the cytoplasm and inner membrane; these proteins can be responsible for transfer of electrons into the quinone pool in a wide variety of environmental conditions. In most of the other shells, proteins are broadly located throughout the five cellular compartments (cytoplasm, inner membrane, periplasm, outer membrane, and extracellular), which ensures the important EET ability of S. oneidensis MR-1. Specifically, the fourth shell was responsible for EET and the c-type cytochromes in the remaining shells of the electron transfer network were involved in aiding EET. Taken together, these results show that there are distinct functional parts in the electron transfer network of S. oneidensis MR-1, and the EET processes could achieve high efficiency through cooperation through such an electron transfer network.

Identifying the potential extracellular electron transfer pathways from a c-type cytochrome network.

Extracellular electron transfer (EET) is the key feature of some bacteria, such as Geobacter sulfurreducens and Shewanella oneidensis. Via EET processes, these bacteria can grow on electrode surfaces and make current output of microbial fuel cells. c-Type cytochromes can be used as carriers to transfer electrons, which play an important role in EET processes. Typically, from the inner (cytoplasmic) membrane through the periplasm to the outer membrane, they could form EET pathways. Recent studies suggest that a group of c-type cytochromes could form a network which extended the well-known EET pathways. We obtained the protein interaction information for all 41 c-type cytochromes in Shewanella oneidensis MR-1, constructed a large-scale protein interaction network, and studied its structural characteristics and functional significance. Centrality analysis has identified the top 10 key proteins of the network, and 7 of them are associated with electricity production in the bacteria, which suggests that the ability of Shewanella oneidensis MR-1 to produce electricity might be derived from the unique structure of the c-type cytochrome network. By modularity analysis, we obtained 5 modules from the network. The subcellular localization study has shown that the proteins in these modules all have diversiform cellular compartments, which reflects their potential to form EET pathways. In particular, combination of protein subcellular localization and operon analysis, the well-known and new candidate EET pathways are obtained from the Mtr-like module, indicating that potential EET pathways could be obtained from such a c-type cytochrome network.