Drug2ways: Reasoning over causal paths in biological networks for drug discovery

Elucidating the causal mechanisms responsible for disease can reveal potential therapeutic targets for pharmacological intervention and, accordingly, guide drug repositioning and discovery. In essence, the topology of a network can reveal the impact a drug candidate may have on a given biological state, leading the way for enhanced disease characterization and the design of advanced therapies. Network-based approaches, in particular, are highly suited for these purposes as they hold the capacity to identify the molecular mechanisms underlying disease. Here, we present drug2ways, a novel methodology that leverages multimodal causal networks for predicting drug candidates. Drug2ways implements an efficient algorithm which reasons over causal paths in large-scale biological networks to propose drug candidates for a given disease. We validate our approach using clinical trial information and demonstrate how drug2ways can be used for multiple applications to identify: i) single-target drug candidates, ii) candidates with polypharmacological properties that can optimize multiple targets, and iii) candidates for combination therapy. Finally, we make drug2ways available to the scientific community as a Python package that enables conducting these applications on multiple standard network formats.     

Challenges and solutions in proteomics

The accelerated growth of proteomics data presents both opportunities and challenges. Large-scale proteomic profiling of biological samples such as cells, organelles or biological fluids has led to discovery of numerous key and novel proteins involved in many biological/disease processes including cancers, as well as to the identification of novel disease biomarkers and potential therapeutic targets. While proteomic data analysis has been greatly assisted by the many bioinformatics tools developed in recent years, a careful analysis of the major steps and flow of data in a typical highthroughput analysis reveals a few gaps that still need to be filled to fully realize the value of the data. To facilitate functional and pathway discovery for large-scale proteomic data, we have developed an integrated proteomic expression analysis system, iProXpress, which facilitates protein identification using a comprehensive sequence library and functional interpretation using integrated data. With its modular design, iProXpress complements and can be integrated with other software in a proteomic data analysis pipeline. This novel approach to complex biological questions involves the interrogation of multiple data sources, thereby facilitating hypothesis generation and knowledge discovery from the genomic-scale studies and fostering disease diagnosis and drug development.     

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

PREDICT: a method for inferring novel drug indications with application to personalized medicine

Inferring potential drug indications, for either novel or approved drugs, is a key step in drug development. Previous computational methods in this domain have focused on either drug repositioning or matching drug and disease gene expression profiles. Here, we present a novel method for the large‐scale prediction of drug indications (PREDICT) that can handle both approved drugs and novel molecules. Our method is based on the observation that similar drugs are indicated for similar diseases, and utilizes multiple drug–drug and disease–disease similarity measures for the prediction task. On cross‐validation, it obtains high specificity and sensitivity (AUC=0.9) in predicting drug indications, surpassing existing methods. We validate our predictions by their overlap with drug indications that are currently under clinical trials, and by their agreement with tissue‐specific expression information on the drug targets. We further show that disease‐specific genetic signatures can be used to accurately predict drug indications for new diseases (AUC=0.92). This lays the computational foundation for future personalized drug treatments, where gene expression signatures from individual patients would replace the disease‐specific signatures.     

Systematic analysis to identify novel disease indications and plausible potential chemical leads of glutamate ionotropic receptor NMDA type subunit 1, GRIN1

Schizophrenia is a mental illness affecting the normal lifestyle of adults and early adolescents incurring major symptoms as jumbled speech, involvement in everyday activities eventually got reduced, patients always struggle with attention and memory, reason being both the genetic and environmental factors responsible for altered brain chemistry and structure, resulting in schizophrenia and associated orphan diseases. The network biology describes the interactions among genes/proteins encoding molecular mechanisms of biological processes, development, and diseases. Besides, all the molecular networks, protein-protein Interaction Networks have been significant in distinguishing the pathogenesis of diseases and thereby drug discovery. The present meta-analysis prioritizes novel disease indications viz. rare and orphan diseases associated with target Glutamate Ionotropic Receptor NMDA Type Subunit 1, GRIN1 using text mining knowledge-based tools. Furthermore, ZINC database was virtually screened, and binding conformation of selected compounds was performed and resulted in the identification of Narciclasine (ZINC04097652) and Alvespimycin (ZINC73138787) as potential inhibitors. Furthermore, docked complexes were subjected to MD simulation studies which suggests that the identified leads could be a better potential drug to recuperate schizophrenia.     

Network based consensus gene signatures for biomarker discovery in breast cancer

Diagnostic and prognostic biomarkers for cancer based on gene expression profiles are viewed as a major step towards a better personalized medicine. Many studies using various computational approaches have been published in this direction during the last decade. However, when comparing different gene signatures for related clinical questions often only a small overlap is observed. This can have various reasons, such as technical differences of platforms, differences in biological samples or their treatment in lab, or statistical reasons because of the high dimensionality of the data combined with small sample size, leading to unstable selection of genes. In conclusion retrieved gene signatures are often hard to interpret from a biological point of view. We here demonstrate that it is possible to construct a consensus signature from a set of seemingly different gene signatures by mapping them on a protein interaction network. Common upstream proteins of close gene products, which we identified via our developed algorithm, show a very clear and significant functional interpretation in terms of overrepresented KEGG pathways, disease associated genes and known drug targets. Moreover, we show that such a consensus signature can serve as prior knowledge for predictive biomarker discovery in breast cancer. Evaluation on different datasets shows that signatures derived from the consensus signature reveal a much higher stability than signatures learned from all probesets on a microarray, while at the same time being at least as predictive. Furthermore, they are clearly interpretable in terms of enriched pathways, disease associated genes and known drug targets. In summary we thus believe that network based consensus signatures are not only a way to relate seemingly different gene signatures to each other in a functional manner, but also to establish prior knowledge for highly stable and interpretable predictive biomarkers.     

Proteomics in the post-genome age.

The genome sequencing effort has helped spawn the burgeoning field of proteomics. This review article examines state-of-the-art proteomics methods that are helping change the discovery paradigm in a variety of biological disciplines and, in particular, protein biochemistry. The review discusses both classical and novel methods to perform high-throughput qualitative and quantitative "global" as well as targeted proteome analysis of complex biological systems. From a drug discovery standpoint, the synergy between genomics and proteomics will help elucidate disease mechanisms, identify novel drug targets, and identify surrogate biomarkers that could be used to conduct clinical trials.     

PhysioSpace: Relating Gene Expression Experiments from Heterogeneous Sources Using Shared Physiological Processes

Relating expression signatures from different sources such as cell lines, in vitro cultures from primary cells and biopsy material is an important task in drug development and translational medicine as well as for tracking of cell fate and disease progression. Especially the comparison of large scale gene expression changes to tissue or cell type specific signatures is of high interest for the tracking of cell fate in (trans-) differentiation experiments and for cancer research, which increasingly focuses on shared processes and the involvement of the microenvironment. These signature relation approaches require robust statistical methods to account for the high biological heterogeneity in clinical data and must cope with small sample sizes in lab experiments and common patterns of co-expression in ubiquitous cellular processes. We describe a novel method, called PhysioSpace, to position dynamics of time series data derived from cellular differentiation and disease progression in a genome-wide expression space. The PhysioSpace is defined by a compendium of publicly available gene expression signatures representing a large set of biological phenotypes. The mapping of gene expression changes onto the PhysioSpace leads to a robust ranking of physiologically relevant signatures, as rigorously evaluated via sample-label permutations. A spherical transformation of the data improves the performance, leading to stable results even in case of small sample sizes. Using PhysioSpace with clinical cancer datasets reveals that such data exhibits large heterogeneity in the number of significant signature associations. This behavior was closely associated with the classification endpoint and cancer type under consideration, indicating shared biological functionalities in disease associated processes. Even though the time series data of cell line differentiation exhibited responses in larger clusters covering several biologically related patterns, top scoring patterns were highly consistent with a priory known biological information and separated from the rest of response patterns.     

Network controllability-based algorithm to target personalized driver genes for discovering combinatorial drugs of individual patients

Multiple driver genes in individual patient samples may cause resistance to individual drugs in precision medicine. However, current computational methods have not studied how to fill the gap between personalized driver gene identification and combinatorial drug discovery for individual patients. Here, we developed a novel structural network controllability-based personalized driver genes and combinatorial drug identification algorithm (CPGD), aiming to identify combinatorial drugs for an individual patient by targeting personalized driver genes from network controllability perspective. On two benchmark disease datasets (i.e. breast cancer and lung cancer datasets), performance of CPGD is superior to that of other state-of-the-art driver gene-focus methods in terms of discovery rate among prior-known clinical efficacious combinatorial drugs. Especially on breast cancer dataset, CPGD evaluated synergistic effect of pairwise drug combinations by measuring synergistic effect of their corresponding personalized driver gene modules, which are affected by a given targeting personalized driver gene set of drugs. The results showed that CPGD performs better than existing synergistic combinatorial strategies in identifying clinical efficacious paired combinatorial drugs. Furthermore, CPGD enhanced cancer subtyping by computationally providing personalized side effect signatures for individual patients. In addition, CPGD identified 90 drug combinations candidates from SARS-COV2 dataset as potential drug repurposing candidates for recently spreading COVID-19.     

生物信息技术加速开发旧药新用途

传统的技术路线研发新药,不仅周期很长而且耗资巨大,开发已获批准药物新的治疗用途,又称为药物重定位,比传统的新药研发具有明显的优势．基于芯片的基因表达谱分析,已常规地广泛用于各种人类疾病的临床研究,提供了在全基因组水平描述疾病状态的特征信号．同时,基因芯片也广泛地用于对比药物处理前后细胞基因表达模式的变化,这也提供了反映药物效应的高质量信号．最近出版的Science Translational Medicine杂志同时发表了一个研究组的两篇论文,为我们展示了如何利用生物信息学手段重新解析和比较全基因组基因表达谱数据,以高效地预测药物的新用途．这两篇论文使用了公共数据库中的100种疾病基因表达谱数据,以及164种药物处理前后细胞基因表达谱数据,通过比较和配对疾病与药物基因表达谱,得到了一些可以逆转疾病异常表达基因的药物,其中证实了一些已知的药物-疾病组合,也预测了一些新的药物-疾病组合．最后通过实验验证了抗溃疡药可用于治疗肺癌,而抗癫痫药可治疗炎症性肠道疾病,进一步证实了他们所采用研究策略的正确性．于是,肺癌和炎性肠道疾病这两种临床上难治的疾病有了新的候选治疗药物,我们也有了一种挖掘已有数据快速发现药物新用途的思路和方法．     

On reliable discovery of molecular signatures

Background  

Molecular signatures are sets of genes, proteins, genetic variants or other variables that can be used as markers for a particular phenotype. Reliable signature discovery methods could yield valuable insight into cell biology and mechanisms of human disease. However, it is currently not clear how to control error rates such as the false discovery rate (FDR) in signature discovery. Moreover, signatures for cancer gene expression have been shown to be unstable, that is, difficult to replicate in independent studies, casting doubts on their reliability.     

Bayes optimal informer sets for early-stage drug discovery

An important experimental design problem in early-stage drug discovery is how to prioritize available compounds for testing when very little is known about the target protein. Informer-based ranking (IBR) methods address the prioritization problem when the compounds have provided bioactivity data on other potentially relevant targets. An IBR method selects an informer set of compounds, and then prioritizes the remaining compounds on the basis of new bioactivity experiments performed with the informer set on the target. We formalize the problem as a two-stage decision problem and introduce the Bayes Optimal Informer SEt (BOISE) method for its solution. BOISE leverages a flexible model of the initial bioactivity data, a relevant loss function, and effective computational schemes to resolve the two-step design problem. We evaluate BOISE and compare it to other IBR strategies in two retrospective studies, one on protein-kinase inhibition and the other on anticancer drug sensitivity. In both empirical settings BOISE exhibits better predictive performance than available methods. It also behaves well with missing data, where methods that use matrix completion show worse predictive performance.     

Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers

Common inflammatome gene signatures as well as disease‐specific signatures were identified by analyzing 12 expression profiling data sets derived from 9 different tissues isolated from 11 rodent inflammatory disease models. The inflammatome signature significantly overlaps with known drug targets and co‐expressed gene modules linked to metabolic disorders and cancer. A large proportion of genes in this signature are tightly connected in tissue‐specific Bayesian networks (BNs) built from multiple independent mouse and human cohorts. Both the inflammatome signature and the corresponding consensus BNs are highly enriched for immune response‐related genes supported as causal for adiposity, adipokine, diabetes, aortic lesion, bone, muscle, and cholesterol traits, suggesting the causal nature of the inflammatome for a variety of diseases. Integration of this inflammatome signature with the BNs uncovered 151 key drivers that appeared to be more biologically important than the non‐drivers in terms of their impact on disease phenotypes. The identification of this inflammatome signature, its network architecture, and key drivers not only highlights the shared etiology but also pinpoints potential targets for intervention of various common diseases.     

Prediction of Drug-Target Interactions for Drug Repositioning Only Based on Genomic Expression Similarity

Small drug molecules usually bind to multiple protein targets or even unintended off-targets. Such drug promiscuity has often led to unwanted or unexplained drug reactions, resulting in side effects or drug repositioning opportunities. So it is always an important issue in pharmacology to identify potential drug-target interactions (DTI). However, DTI discovery by experiment remains a challenging task, due to high expense of time and resources. Many computational methods are therefore developed to predict DTI with high throughput biological and clinical data. Here, we initiatively demonstrate that the on-target and off-target effects could be characterized by drug-induced in vitro genomic expression changes, e.g. the data in Connectivity Map (CMap). Thus, unknown ligands of a certain target can be found from the compounds showing high gene-expression similarity to the known ligands. Then to clarify the detailed practice of CMap based DTI prediction, we objectively evaluate how well each target is characterized by CMap. The results suggest that (1) some targets are better characterized than others, so the prediction models specific to these well characterized targets would be more accurate and reliable; (2) in some cases, a family of ligands for the same target tend to interact with common off-targets, which may help increase the efficiency of DTI discovery and explain the mechanisms of complicated drug actions. In the present study, CMap expression similarity is proposed as a novel indicator of drug-target interactions. The detailed strategies of improving data quality by decreasing the batch effect and building prediction models are also effectively established. We believe the success in CMap can be further translated into other public and commercial data of genomic expression, thus increasing research productivity towards valid drug repositioning and minimal side effects.     

Chemogenomics approaches to novel target discovery

Chemogenomics involves the combination of a compound’s effect on biological targets together with modern genomics technologies. The merger of these two methodologies is creating a new way to screen for compound–target interactions, as well as map chemical and biological space in a parallel fashion. The challenge associated with mining complex databases has initiated the development of many novel in silico tools to profile and analyze data in a systematic way. The ability to analyze the combinatorial effects of chemical libraries on biological systems will aid the discovery of new therapeutic entities. Chemogenomics provides a tool for the rapid validation of novel targeted therapeutics, where a specific molecular target is modulated by a small molecule. Along with targeted therapies comes the ability to discovery pathway nodes where a single molecular target might be an essential component of more than one disease. Several disease areas will benefit directly from the chemogenomics approach, the most advanced being cancer. A genetic loss-of-function screen can be modulated in the presence of a compound to search for genes or pathways involved in the compound’s activity. Several recent papers highlight how chemogenomics is changing with RNA interference-based screening and shaping the discovery of new targeted therapies. Together, chemical and RNA interference-based screens open the door for a new way to discovery disease-associated genes and novel targeted therapies.     

Drug discovery for parasitic diseases: powered by technology,enabled by pharmacology,informed by clinical science

While prevention is a bedrock of public health, innovative therapeutics are needed to complement the armamentarium of interventions required to achieve disease control and elimination targets for neglected diseases. Extraordinary advances in drug discovery technologies have occurred over the past decades, along with accumulation of scientific knowledge and experience in pharmacological and clinical sciences that are transforming many aspects of drug R&D across disciplines. We reflect on how these advances have propelled drug discovery for parasitic infections, focusing on malaria, kinetoplastid diseases, and cryptosporidiosis. We also discuss challenges and research priorities to accelerate discovery and development of urgently needed novel antiparasitic drugs.     

Molecular signatures from omics data: From chaos to consensus

In the past 15 years, new "omics" technologies have made it possible to obtain high-resolution molecular snapshots of organisms, tissues, and even individual cells at various disease states and experimental conditions. It is hoped that these developments will usher in a new era of personalized medicine in which an individual's molecular measurements are used to diagnose disease, guide therapy, and perform other tasks more accurately and effectively than is possible using standard approaches. There now exists a vast literature of reported "molecular signatures". However, despite some notable exceptions, many of these signatures have suffered from limited reproducibility in independent datasets, insufficient sensitivity or specificity to meet clinical needs, or other challenges. In this paper, we discuss the process of molecular signature discovery on the basis of omics data. In particular, we highlight potential pitfalls in the discovery process, as well as strategies that can be used to increase the odds of successful discovery. Despite the difficulties that have plagued the field of molecular signature discovery, we remain optimistic about the potential to harness the vast amounts of available omics data in order to substantially impact clinical practice.