Abstract: BACKGROUND: Glioblastoma Multiforme, an aggressive primary brain tumor, has a poor prognosis and no effective standard of care treatments. Most patients undergoing radiotherapy, along with Temozolomide chemotherapy, develop resistance to the drug, and recurrence of the tumor is a common issue after the treatment. We propose to model the pathways active in Glioblastoma using Boolean network techniques. The network captures the genetic interactions and possible mutations that are involved in the development of the brain tumor. The model is used to predict the theoretical efficacies of drugs for the treatment of cancer. RESULTS: We use the Boolean network to rank the critical intervention points in the pathway to predict an effective therapeutic strategy for Glioblastoma. Drug repurposing helps to identify non-cancer drugs that could be effective in cancer treatment. We predict the effectiveness of drug combinations of anti-cancer and non-cancer drugs for Glioblastoma. CONCLUSIONS: Given the genetic profile of a GBM tumor, the Boolean model can predict the most effective targets for treatment. We also identified two-drug combinations that could be more effective in killing GBM cells than conventional chemotherapeutic agents. The non-cancer drug Aspirin could potentially increase the cytotoxicity of TMZ in GBM patients.

Abstract: Small molecules play a critical role in modulating biological systems. Knowledge of chemical-protein interactions helps address fundamental and practical questions in biology and medicine. However, with the rapid emergence of newly sequenced genes, the endogenous or surrogate ligands of a vast number of proteins remain unknown. Homology modeling and machine learning are two major methods for assigning new ligands to a protein but mostly fail when sequence homology between an unannotated protein and those with known functions or structures is low. In this study, we develop a new deep learning framework to predict chemical binding to evolutionary divergent unannotated proteins, whose ligand cannot be reliably predicted by existing methods. By incorporating evolutionary information into self-supervised learning of unlabeled protein sequences, we develop a novel method, distilled sequence alignment embedding (DISAE), for the protein sequence representation. DISAE can utilize all protein sequences and their multiple sequence alignment (MSA) to capture functional relationships between proteins without the knowledge of their structure and function. Followed by the DISAE pretraining, we devise a module-based fine-tuning strategy for the supervised learning of chemical-protein interactions. In the benchmark studies, DISAE significantly improves the generalizability of machine learning models and outperforms the state-of-the-art methods by a large margin. Comprehensive ablation studies suggest that the use of MSA, sequence distillation, and triplet pretraining critically contributes to the success of DISAE. The interpretability analysis of DISAE suggests that it learns biologically meaningful information. We further use DISAE to assign ligands to human orphan G-protein coupled receptors (GPCRs) and to cluster the human GPCRome by integrating their phylogenetic and ligand relationships. The promising results of DISAE open an avenue for exploring the chemical landscape of entire sequenced genomes.

Abstract: To reduce time and cost, virtual ligand screening (VLS) often precedes experimental ligand screening in modern drug discovery. Traditionally, high-resolution structure-based docking approaches rely on experimental structures, while ligand-based approaches need known binders to the target protein and only explore their nearby chemical space. In contrast, our structure-based FINDSITE(comb2.0) approach takes advantage of predicted, low-resolution structures and information from ligands that bind distantly related proteins whose binding sites are similar to the target protein. Using a boosted tree regression machine learning framework, we significantly improved FINDSITE(comb2.0) by integrating ligand fragment scores as encoded by molecular fingerprints with the global ligand similarity scores of FINDSITE(comb2.0). The new approach, FRAGSITE, exploits our observation that ligand fragments, e.g., rings, tend to interact with stereochemically conserved protein subpockets that also occur in evolutionarily unrelated proteins. FRAGSITE was benchmarked on the 102 protein DUD-E set, where any template protein whose sequence identify >30% to the target was excluded. Within the top 100 ranked molecules, FRAGSITE improves VLS precision and recall by 14.3 and 18.5%, respectively, relative to FINDSITE(comb2.0). Moreover, the mean top 1% enrichment factor increases from 25.2 to 30.2. On average, both outperform state-of-the-art deep learning-based methods such as AtomNet. On the more challenging unbiased set LIT-PCBA, FRAGSITE also shows better performance than ligand similarity-based and docking approaches such as two-dimensional ECFP4 and Surflex-Dock v.3066. On a subset of 23 targets from DEKOIS 2.0, FRAGSITE shows much better performance than the boosted tree regression-based, vScreenML scoring function. Experimental testing of FRAGSITE's predictions shows that it has more hits and covers a more diverse region of chemical space than FINDSITE(comb2.0). For the two proteins that were experimentally tested, DHFR, a well-studied protein that catalyzes the conversion of dihydrofolate to tetrahydrofolate, and the kinase ACVR1, FRAGSITE identified new small-molecule nanomolar binders. Interestingly, one new binder of DHFR is a kinase inhibitor predicted to bind in a new subpocket. For ACVR1, FRAGSITE identified new molecules that have diverse scaffolds and estimated nanomolar to micromolar affinities. Thus, FRAGSITE shows significant improvement over prior state-of-the-art ligand virtual screening approaches. A web server is freely available for academic users at http:/sites.gatech.edu/cssb/FRAGSITE.

Abstract: There is an urgent need to identify new therapies that prevent SARS-CoV-2 infection and improve the outcome of COVID-19 patients. This pandemic has thus spurred intensive research in most scientific areas and in a short period of time, several vaccines have been developed. But, while the race to find vaccines for COVID-19 has dominated the headlines, other types of therapeutic agents are being developed. In this mini-review, we report several databases and online tools that could assist the discovery of anti-SARS-CoV-2 small chemical compounds and peptides. We then give examples of studies that combined in silico and in vitro screening, either for drug repositioning purposes or to search for novel bioactive compounds. Finally, we question the overall lack of discussion and plan observed in academic research in many countries during this crisis and suggest that there is room for improvement.

Abstract: OBJECTIVE: Drug repurposing, the application of existing therapeutics to new indications, holds promise in achieving rapid clinical effects at a much lower cost than that of de novo drug development. The aim of our study was to perform a more comprehensive drug repurposing prediction of diseases, particularly cancers. METHODS: Here, by targeting 4,096 human diseases, including 384 cancers, we propose a greedy computational model based on a heterogeneous multilayer network for the repurposing of 1,419 existing drugs in DrugBank. We performed additional experimental validation for the dominant repurposed drugs in cancer. RESULTS: The overall performance of the model was well supported by cross-validation and literature mining. Focusing on the top-ranked repurposed drugs in cancers, we verified the anticancer effects of 5 repurposed drugs widely used clinically in drug sensitivity experiments. Because of the distinctive antitumor effects of nifedipine (an antihypertensive agent) and nortriptyline (an antidepressant drug) in prostate cancer, we further explored their underlying mechanisms by using quantitative proteomics. Our analysis revealed that both nifedipine and nortriptyline affected the cancer-related pathways of DNA replication, the cell cycle, and RNA transport. Moreover, in vivo experiments demonstrated that nifedipine and nortriptyline significantly inhibited the growth of prostate tumors in a xenograft model. CONCLUSIONS: Our predicted results, which have been released in a public database named The Predictive Database for Drug Repurposing (PAD), provide an informative resource for discovering and ranking drugs that may potentially be repurposed for cancer treatment and determining new therapeutic effects of existing drugs.

Abstract: The global spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires an urgent need to find effective therapeutics for the treatment of coronavirus disease 2019 (COVID-19). In this study, we developed an integrative drug repositioning framework, which fully takes advantage of machine learning and statistical analysis approaches to systematically integrate and mine large-scale knowledge graph, literature and transcriptome data to discover the potential drug candidates against SARS-CoV-2. Our in silico screening followed by wet-lab validation indicated that a poly-ADP-ribose polymerase 1 (PARP1) inhibitor, CVL218, currently in Phase I clinical trial, may be repurposed to treat COVID-19. Our in vitro assays revealed that CVL218 can exhibit effective inhibitory activity against SARS-CoV-2 replication without obvious cytopathic effect. In addition, we showed that CVL218 can interact with the nucleocapsid (N) protein of SARS-CoV-2 and is able to suppress the LPS-induced production of several inflammatory cytokines that are highly relevant to the prevention of immunopathology induced by SARS-CoV-2 infection.

Abstract: With the advent of deep generative models in computational chemistry, in-silico drug design is undergoing an unprecedented transformation. Although deep learning approaches have shown potential in generating compounds with desired chemical properties, they disregard the cellular environment of target diseases. Bridging systems biology and drug design, we present a reinforcement learning method for de novo molecular design from gene expression profiles. We construct a hybrid Variational Autoencoder that tailors molecules to target-specific transcriptomic profiles, using an anticancer drug sensitivity prediction model (PaccMann) as reward function. Without incorporating information about anticancer drugs, the molecule generation is biased toward compounds with high predicted efficacy against cell lines or cancer types. The generation can be further refined by subsidiary constraints such as toxicity. Our cancer-type-specific candidate drugs are similar to cancer drugs in drug-likeness, synthesizability, and solubility and frequently exhibit the highest structural similarity to compounds with known efficacy against these cancer types.

Abstract: Drug discovery is a cost and time-intensive process that is often assisted by computational methods, such as virtual screening, to speed up and guide the design of new compounds. For many years, machine learning methods have been successfully applied in the context of computer-aided drug discovery. Recently, thanks to the rise of novel technologies as well as the increasing amount of available chemical and bioactivity data, deep learning has gained a tremendous impact in rational active compound discovery. Herein, recent applications and developments of machine learning, with a focus on deep learning, in virtual screening for active compound design are reviewed. This includes introducing different compound and protein encodings, deep learning techniques as well as frequently used bioactivity and benchmark data sets for model training and testing. Finally, the present state-of-the-art, including the current challenges and emerging problems, are examined and discussed.

Abstract: BACKGROUND: Systemic juvenile idiopathic arthritis (sJIA) and adult-onset Still's disease (AOSD) are manifestations of an autoinflammatory disorder with complex pathophysiology and significant morbidity, together also termed Still's disease. The objective of the current study is to set in silico models based on systems biology and investigate the optimal treat-to-target strategy for Still's disease as a proof-of-concept of the modeling approach. METHODS: Molecular characteristics of Still's disease and data on biological inhibitors of interleukin (IL)-1 (anakinra, canakinumab), IL-6 (tocilizumab, sarilumab), and glucocorticoids as well as conventional disease-modifying anti-rheumatic drugs (DMARDs, methotrexate) were used to construct in silico mechanisms of action (MoA) models by means of Therapeutic Performance Mapping System (TPMS) technology. TPMS combines artificial neuronal networks, sampling-based methods, and artificial intelligence. Model outcomes were validated with published expression data from sJIA patients. RESULTS: Biologicals demonstrated more pathophysiology-directed efficiency than non-biological drugs. IL-1 blockade mainly acts on proteins implicated in the innate immune system, while IL-6 signaling blockade has a weaker effect on innate immunity and rather affects adaptive immune mechanisms. The MoA models showed that in the autoinflammatory/systemic phases of Still's disease, in which the innate immunity plays a pivotal role, the IL-1beta-neutralizing antibody canakinumab is more efficient than the IL-6 receptor-inhibiting antibody tocilizumab. MoA models reproduced 67% of the information obtained from expression data. CONCLUSIONS: Systems biology-based modeling supported the preferred use of biologics as an immunomodulatory treatment strategy for Still's disease. Our results reinforce the role for IL-1 blockade on innate immunity regulation, which is critical in systemic autoinflammatory diseases. This further encourages early use on Still's disease IL-1 blockade to prevent the development of disease or drug-related complications. Further analysis at the clinical level will validate the findings and help determining the timeframe of the window of opportunity for canakinumab treatment.

Abstract: Generally, bioequivalence (BE) studies of drug products for pediatric patients are conducted in adults due to ethical reasons. Given the lack of direct BE assessment in pediatric populations, the aim of this work is to develop a database of BE and relative bioavailability (relative BA) studies conducted in pediatric populations and to enable the identification of risk factors associated with certain drug substances or products that may lead to failed BE or different pharmacokinetic (PK) parameters in relative BA studies in pediatrics. A literature search from 1965 to 2020 was conducted in PubMed, Cochrane Library, and Google Scholar to identify BE studies conducted in pediatric populations and relative BA studies conducted in pediatric populations. Overall, 79 studies covering 37 active pharmaceutical ingredients (APIs) were included in the database: 4 bioequivalence studies with data that passed BE evaluations; 2 studies showed bioinequivalence results; 34 relative BA studies showing comparable PK parameters, and 39 relative BA studies showing differences in PK parameters between test and reference products. Based on the above studies, common putative risk factors associated with differences in relative bioavailability (DRBA) in pediatric populations include age-related absorption effects, high inter-individual variability, and poor study design. A database containing 79 clinical studies on BE or relative BA in pediatrics has been developed. Putative risk factors associated with DRBA in pediatric populations are summarized.