Abstract: Herein, the DrugBank database which contains 10,036 approved and investigational drugs was explored deeply for potential drugs that target SARS-CoV-2 main protease (M(pro)). Filtration process of the database was conducted using three levels of accuracy for molecular docking calculations. The top 35 drugs with docking scores > -11.0 kcal/mol were then subjected to 10 ns molecular dynamics (MD) simulations followed by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. The results showed that DB02388 and Cobicistat (DB09065) exhibited potential binding affinities towards M(pro) over 100 ns MD simulations, with binding energy values of -49.67 and -46.60 kcal/mol, respectively. Binding energy and structural analyses demonstrated the higher stability of DB02388 over Cobicistat. The potency of DB02388 and Cobicistat is attributed to their abilities to form several hydrogen bonds with the essential amino acids inside the active site of M(pro). Compared to DB02388 and Cobicistat, Darunavir showed a much lower binding affinity of -34.83 kcal/mol. The present study highlights the potentiality of DB02388 and Cobicistat as anti-COVID-19 drugs for clinical trials. Communicated by Ramaswamy H. Sarma.

Abstract: BACKGROUND: Thyroid-associated ophthalmopathy (TAO) is an autoimmune inflammatory disorder, which lacks effective treatment currently. Spica Prunellae (SP) is popularly used for its anti-inflammatory and immune-regulating properties, indicating SP may have potential therapeutic value in TAO. Therefore, the purpose of this study is to identify the efficiency and potential mechanism of SP in treating TAO. METHODS: A network pharmacology integrated molecular docking strategy was used to predict the underlying molecular mechanism of treating TAO. Firstly, the active compounds of SP were obtained from TCMSP database and literature research. Then we collected the putative targets of SP and TAO based on multi-sources databases to generate networks. Network topology analysis, GO and KEGG pathway enrichment analysis were performed to screen the key targets and mechanism. Furthermore, molecular docking simulation provided an assessment tool for verifying drug and target binding. RESULTS: Our results showed that 8 targets (PTGS2, MAPK3, AKT1, TNF, MAPK1, CASP3, IL6, MMP9) were recognized as key therapeutic targets with excellent binding affinity after network analysis and molecular docking-based virtual screening. The results of enrichment analysis suggested that the underlying mechanism was mainly focused on the biological processes and pathways associated with immune inflammation, proliferation, and apoptosis. Notably, the key pathway was considered as the PI3K-AKT signaling pathway. CONCLUSION: In summary, the present study elucidates that SP may suppress inflammation and proliferation and promote apoptosis through the PI3K-AKT pathway, which makes SP a potential treatment against TAO. And this study offers new reference points for future experimental research and provides a scientific basis for more widespread clinical application.

Abstract: Target-based high-throughput compound screening dominates conventional one-drug-one-gene drug discovery process. However, the readout from the chemical modulation of a single protein is poorly correlated with phenotypic response of organism, leading to high failure rate in drug development. Chemical-induced gene expression profile provides an attractive solution to phenotype-based screening. However, the use of such data is currently limited by their sparseness, unreliability, and relatively low throughput. Several methods have been proposed to impute missing values for gene expression datasets. However, few existing methods can perform de novo chemical compound screening. In this study, we propose a mechanism-driven neural network-based method named DeepCE (Deep Chemical Expression) which utilizes graph convolutional neural network to learn chemical representation and multi-head attention mechanism to model chemical substructure-gene and gene-gene feature associations. In addition, we propose a novel data augmentation method which extracts useful information from unreliable experiments in L1000 dataset. The experimental results show that DeepCE achieves the superior performances not only in de novo chemical setting but also in traditional imputation setting compared to state-of-the-art baselines for the prediction of chemical-induced gene expression. We further verify the effectiveness of gene expression profiles generated from DeepCE by comparing them with gene expression profiles in L1000 dataset for downstream classification tasks including drug-target and disease predictions. To demonstrate the value of DeepCE, we apply it to patient-specific drug repurposing of COVID-19 for the first time, and generate novel lead compounds consistent with clinical evidences. Thus, DeepCE provides a potentially powerful framework for robust predictive modeling by utilizing noisy omics data as well as screening novel chemicals for the modulation of systemic response to disease.

Abstract: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spillover infection in December 2019 has caused an unprecedented pandemic. SARS-CoV-2, as other coronaviruses, binds its target cells through the angiotensin-converting enzyme 2 (ACE2) receptor. Accordingly, this makes ACE2 research essential for understanding the zoonotic nature of coronaviruses and identifying novel drugs. Here we present a systematic analysis of the ACE2 conservation and co-evolution protein network across 1,671 eukaryotes, revealing an unexpected conservation pattern in specific metazoans, plants, fungi, and protists. We identified the co-evolved protein network and pinpointed a list of drugs that target this network by using data integration from different sources. Our computational analysis found widely used drugs such as nonsteroidal anti-inflammatory drugs and vasodilators. These drugs are expected to perturb the ACE2 network affecting infectivity as well as the pathophysiology of the disease.

Abstract: Some drugs can be used to treat multiple diseases, suggesting potential patterns in drug treatment. Determination of drug treatment patterns can improve our understanding of the mechanisms of drug action, enabling drug repurposing. A drug can be associated with a multilayer tissue-specific protein-protein interaction (TSPPI) network for the diseases it is used to treat. Proteins usually interact with other proteins to achieve functions that cause diseases. Hence, studying drug treatment patterns is similar to studying common module structures in multilayer TSPPI networks. Therefore, we propose a network-based model to study the treatment patterns of drugs. The method was designated SDTP (studying drug treatment pattern) and was based on drug effects and a multilayer network model. To demonstrate the application of the SDTP method, we focused on analysis of trichostatin A (TSA) in leukemia, breast cancer, and prostate cancer. We constructed a TSPPI multilayer network and obtained candidate drug-target modules from the network. Gene ontology analysis provided insights into the significance of the drug-target modules and co-expression networks. Finally, two modules were obtained as potential treatment patterns for TSA. Through analysis of the significance, composition, and functions of the selected drug-target modules, we validated the feasibility and rationality of our proposed SDTP method for identifying drug treatment patterns. In summary, our novel approach used a multilayer network model to overcome the shortcomings of single-layer networks and combined the network with information on drug activity. Based on the discovered drug treatment patterns, we can predict the potential diseases that the drug can treat. That is, if a disease-related protein module has a similar structure, then the drug is likely to be a potential drug for the treatment of the disease.

Abstract: BACKGROUND: Drug repurposing aims to detect the new therapeutic benefits of the existing drugs and reduce the spent time and cost of the drug development projects. The synthetic repurposing of drugs may prove to be more useful than the single repurposing in terms of reducing toxicity and enhancing efficacy. However, the researchers have not given it serious consideration. To address the issue, a novel datamining method is introduced and applied to repositioning of drugs for hypertension (HT) which is a serious medical condition and needs some improved treatment plans to help treat it. RESULTS: A novel two-step data mining method, which is based on the If-Then association rules as well as a novel discrete optimization algorithm, was introduced and applied to the synthetic repurposing of drugs for HT. The required data were also extracted from DrugBank, KEGG, and DrugR+ databases. The findings indicated that based on the different statistical criteria, the proposed method outperformed the other state-of-the-art approaches. In contrast to the previously proposed methods which had failed to discover a list on some datasets, our method could find a combination list for all of them. CONCLUSION: Since the proposed synthetic method uses medications in small dosages, it might revive some failed drug development projects and put forward a suitable plan for treating different diseases such as COVID-19 and HT. It is also worth noting that applying efficient computational methods helps to produce better results.

Abstract: Propolis is a natural resinous material produced by bees and has been used in folk medicines since ancient times. Due to it possessing a broad spectrum of biological activities, it has gained significant scientific and commercial interest over the last two decades. As a result of searching 122 publications reported up to the end of 2019, we assembled a unique compound database consisting of 578 components isolated from both honey bee propolis and stingless bee propolis, and analyzed the chemical space and chemical diversity of these compounds. The results demonstrated that both honey bee propolis and stingless bee propolis are valuable sources for pharmaceutical and nutraceutical development.

Abstract: Essential genes are those whose loss of function compromises organism viability or results in profound loss of fitness. Recent gene-editing technologies have provided new opportunities to characterize essential genes. Here, we present an integrated analysis that comprehensively and systematically elucidates the genetic and regulatory characteristics of human essential genes. First, we found that essential genes act as 'hubs' in protein-protein interaction networks, chromatin structure and epigenetic modification. Second, essential genes represent conserved biological processes across species, although gene essentiality changes differently among species. Third, essential genes are important for cell development due to their discriminate transcription activity in embryo development and oncogenesis. In addition, we developed an interactive web server, the Human Essential Genes Interactive Analysis Platform (http://sysomics.com/HEGIAP/), which integrates abundant analytical tools to enable global, multidimensional interpretation of gene essentiality. Our study provides new insights that improve the understanding of human essential genes.

Abstract: OBJECTIVE: Coronavirus disease 2019 (COVID-19) is a fatal and fast-spreading viral infection. To date, the number of COVID-19 patients worldwide has crossed over six million with over three hundred and seventy thousand deaths (according to the data from World Health Organization; updated on 2 June 2020). Although COVID-19 can be rapidly diagnosed, efficient clinical treatment of COVID-19 remains unavailable, resulting in high fatality. Some clinical trials have identified vitamin C (VC) as a potent compound pneumonia management. In addition, glycyrrhizic acid (GA) is clinically as an anti-inflammatory medicine against pneumonia-induced inflammatory stress. We hypothesized that the combination of VC and GA is a potential option for treating COVID-19. METHODS: The aim of this study was to determine pharmacological targets and molecular mechanisms of VC + GA treatment for COVID-19, using bioinformational network pharmacology. RESULTS: We uncovered optimal targets, biological processes and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of VC + GA against COVID-19. Our findings suggested that combinatorial VC and GA treatment for COVID-19 was associated with elevation of immunity and suppression of inflammatory stress, including activation of the T cell receptor signaling pathway, regulation of Fc gamma R-mediated phagocytosis, ErbB signaling pathway and vascular endothelial growth factor signaling pathway. We also identified 17 core targets of VC + GA, which suggest as antimicrobial function. CONCLUSIONS: For the first time, our study uncovered the pharmacological mechanism underlying combined VC and GA treatment for COVID-19. These results should benefit efforts to address the most pressing problem currently facing the world.

Abstract: BACKGROUND: Despite continued efforts using chemical similarity methods in virtual screening, currently developed approaches suffer from time-consuming multistep procedures and low success rates. We recently developed a machine learning-based chemical binding similarity model considering common structural features from molecules binding to the same, or evolutionarily related targets. The chemical binding similarity measures the resemblance of chemical compounds in terms of binding site similarity to better describe functional similarities that arise from target binding. In this study, we have shown how the chemical binding similarity could be used in virtual screening together with the conventional structure-based methods. RESULTS: The chemical binding similarity, receptor-based pharmacophore, chemical structure similarity, and molecular docking methods were evaluated to identify an effective virtual screening procedure for desired target proteins. When we tested the chemical binding similarity method with test sets of 51 kinases, it outperformed the traditional structural similarity-based methods as well as structure-based methods, such as molecular docking and receptor-based pharmacophore modeling, in terms of finding active compounds. We further validated the results by performing virtual screening (using the chemical binding similarity and receptor-based pharmacophore methods) against a completely blind dataset for mitogen-activated protein kinase kinase 1 (MEK1), ephrin type-B receptor 4 (EPHB4) and wee1-like protein kinase (WEE1). The in vitro kinase binding assay confirmed that 6 out of 13 (46.2%) for MEK1 and 2 out of 12 (16.7%) for EPHB4 were newly identified only by the chemical binding similarity model. CONCLUSIONS: We report that the virtual screening results could further be improved by combining the chemical binding similarity model with 3D-QSAR pharmacophore and molecular docking models. Not only the new inhibitors are identified in this study, but also many of the identified molecules have low structural similarity scores against already reported inhibitors and that show the revelation of novel scaffolds.