Identifying unexpected therapeutic targets via chemical-protein interactome

Drug-drug interaction prediction: databases, web servers and computational models

In clinical treatment, two or more drugs (i.e. drug combination) are simultaneously or successively used for therapy with the purpose of primarily enhancing the therapeutic efficacy or reducing drug side effects. However, inappropriate drug combination may not only fail to improve efficacy, but even lead to adverse reactions. Therefore, according to the basic principle of improving the efficacy and/or reducing adverse reactions, we should study drug-drug interactions (DDIs) comprehensively and thoroughly so as to reasonably use drug combination. In this review, we first introduced the basic conception and classification of DDIs. Further, some important publicly available databases and web servers about experimentally verified or predicted DDIs were briefly described. As an effective auxiliary tool, computational models for predicting DDIs can not only save the cost of biological experiments, but also provide relevant guidance for combination therapy to some extent. Therefore, we summarized three types of prediction models (including traditional machine learning-based models, deep learning-based models and score function-based models) proposed during recent years and discussed the advantages as well as limitations of them. Besides, we pointed out the problems that need to be solved in the future research of DDIs prediction and provided corresponding suggestions.

Uncovering New Drug Properties in Target-Based Drug-Drug Similarity Networks

Despite recent advances in bioinformatics, systems biology, and machine learning, the accurate prediction of drug properties remains an open problem. Indeed, because the biological environment is a complex system, the traditional approach—based on knowledge about the chemical structures—can not fully explain the nature of interactions between drugs and biological targets. Consequently, in this paper, we propose an unsupervised machine learning approach that uses the information we know about drug–target interactions to infer drug properties. To this end, we define drug similarity based on drug–target interactions and build a weighted Drug–Drug Similarity Network according to the drug–drug similarity relationships. Using an energy-model network layout, we generate drug communities associated with specific, dominant drug properties. DrugBank confirms the properties of 59.52% of the drugs in these communities, and 26.98% are existing drug repositioning hints we reconstruct with our DDSN approach. The remaining 13.49% of the drugs seem not to match the dominant pharmacologic property; thus, we consider them potential drug repurposing hints. The resources required to test all these repurposing hints are considerable. Therefore we introduce a mechanism of prioritization based on the betweenness/degree node centrality. Using betweenness/degree as an indicator of drug repurposing potential, we select Azelaic acid and Meprobamate as a possible antineoplastic and antifungal, respectively. Finally, we use a test procedure based on molecular docking to analyze Azelaic acid and Meprobamate’s repurposing.

DrPOCS: Drug Repositioning Based on Projection Onto Convex Sets

Drug repositioning, i.e., identifying new indications for known drugs, has attracted a lot of attentions recently and is becoming an effective strategy in drug development. In literature, several computational approaches have been proposed to identify potential indications of old drugs based on various types of data sources. In this paper, by formulating the drug-disease associations as a low-rank matrix, we propose a novel method, namely DrPOCS, to identify candidate indications of old drugs based on projection onto convex sets (POCS). With the integration of drug structure and disease phenotype information, DrPOCS predicts potential associations between drugs and diseases with matrix completion. Benchmarking results demonstrate that our proposed approach outperforms popular existing approaches with high accuracy. In addition, a number of novel predicted indications are validated with various types of evidences, indicating the predictive power of our proposed approach.

A benchmark driven guide to binding site comparison: An exhaustive evaluation using tailor-made data sets (ProSPECCTs)

The automated comparison of protein-ligand binding sites provides useful insights into yet unexplored site similarities. Various stages of computational and chemical biology research can benefit from this knowledge. The search for putative off-targets and the establishment of polypharmacological effects by comparing binding sites led to promising results for numerous projects. Although many cavity comparison methods are available, a comprehensive analysis to guide the choice of a tool for a specific application is wanting. Moreover, the broad variety of binding site modeling approaches, comparison algorithms, and scoring metrics impedes this choice. Herein, we aim to elucidate strengths and weaknesses of binding site comparison methodologies. A detailed benchmark study is the only possibility to rationalize the selection of appropriate tools for different scenarios. Specific evaluation data sets were developed to shed light on multiple aspects of binding site comparison. An assembly of all applied benchmark sets (ProSPECCTs–Protein Site Pairs for the Evaluation of Cavity Comparison Tools) is made available for the evaluation and optimization of further and still emerging methods. The results indicate the importance of such analyses to facilitate the choice of a methodology that complies with the requirements of a specific scientific challenge.

Binding site similarities are useful in the context of promiscuity prediction, drug repurposing, the analysis of protein-ligand and protein-protein complexes, function prediction, and further fields of general interest in chemical biology and biochemistry. Many years of research have led to the development of a multitude of methods for binding site analysis and comparison. On the one hand, their availability supports research. On the other hand, the huge number of methods hampers the efficient selection of a specific tool. Our research is dedicated to the analysis of different cavity comparison tools. We use several binding site data sets to establish guidelines which can be applied to ensure a successful application of comparison methods by circumventing potential pitfalls.

Computational Multitarget Drug Design

Designing drugs that can simultaneously interact with multiple targets is a promising approach for treating complicated diseases. Compared to using combinations of single target drugs, multitarget drugs have advantages of higher efficacy, improved safety profile, and simpler administration. Many in silico methods have been developed to approach different aspects of this polypharmacology-guided drug design, particularly for drug repurposing and multitarget drug design. In this review, we summarize recent progress in computational multitarget drug design and discuss their advantages and limitations. Perspectives for future drug development will also be discussed.

DPDR-CPI, a server that predicts Drug Positioning and Drug Repositioning via Chemical-Protein Interactome

The cost of developing a new drug has increased sharply over the past years. To ensure a reasonable return-on-investment, it is useful for drug discovery researchers in both industry and academia to identify all the possible indications for early pipeline molecules. For the first time, we propose the term computational “drug candidate positioning” or “drug positioning”, to describe the above process. It is distinct from drug repositioning, which identifies new uses for existing drugs and maximizes their value. Since many therapeutic effects are mediated by unexpected drug-protein interactions, it is reasonable to analyze the chemical-protein interactome (CPI) profiles to predict indications. Here we introduce the server DPDR-CPI, which can make real-time predictions based only on the structure of the small molecule. When a user submits a molecule, the server will dock it across 611 human proteins, generating a CPI profile of features that can be used for predictions. It can suggest the likelihood of relevance of the input molecule towards ~1,000 human diseases with top predictions listed. DPDR-CPI achieved an overall AUROC of 0.78 during 10-fold cross-validations and AUROC of 0.76 for the independent validation. The server is freely accessible via http://cpi.bio-x.cn/dpdr/.

Toxicological Differences Between NMDA Receptor Antagonists and Cholinesterase Inhibitors

Cholinesterase inhibitors (ChEIs), represented by donepezil, rivastigmine, and galantamine, used to be the only approved class of drugs for the treatment of Alzheimer's disease. After the approval of memantine by the Food and Drug Administration (FDA), N-methyl-d-aspartic acid (NMDA) receptor antagonists have been recognized by authorities and broadly used in the treatment of Alzheimer's disease. Along with complementary mechanisms of action, NMDA antagonists and ChEIs differ not only in therapeutic effects but also in adverse reactions, which is an important consideration in clinical drug use. And the number of patients using NMDA antagonists and ChEIs concomitantly has increased, making the matter more complicated. Here we used the FDA Adverse Event Reporting System for statistical analysis , in order to compare the adverse events of memantine and ChEIs. In general, the clinical evidence confirmed the safety advantages of memantine over ChEIs, reiterating the precautions of clinical drug use and the future direction of antidementia drug development.

Impact of Binding Site Comparisons on Medicinal Chemistry and Rational Molecular Design

Modern rational drug design not only deals with the search for ligands binding to interesting and promising validated targets but also aims to identify the function and ligands of yet uncharacterized proteins having impact on different diseases. Additionally, it contributes to the design of inhibitors with distinct selectivity patterns and the prediction of possible off-target effects. The identification of similarities between binding sites of various proteins is a useful approach to cope with those challenges. The main scope of this perspective is to describe applications of different protein binding site comparison approaches to outline their applicability and impact on molecular design. The article deals with various substantial application domains and provides some outstanding examples to show how various binding site comparison methods can be applied to promote in silico drug design workflows. In addition, we will also briefly introduce the fundamental principles of different protein binding site comparison methods.

Opportunities for Web-based Drug Repositioning: Searching for Potential Antihypertensive Agents with Hypotension Adverse Events

Background

Drug repositioning refers to the process of developing new indications for existing drugs. As a phenotypic indicator of drug response in humans, clinical side effects may provide straightforward signals and unique opportunities for drug repositioning.

Objective

We aimed to identify drugs frequently associated with hypotension adverse reactions (ie, the opposite condition of hypertension), which could be potential candidates as antihypertensive agents.

Methods

We systematically searched the electronic records of the US Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS) through the openFDA platform to assess the association between hypotension incidence and antihypertensive therapeutic effect regarding a list of 683 drugs.

Results

Statistical analysis of FAERS data demonstrated that those drugs frequently co-occurring with hypotension events were more likely to have antihypertensive activity. Ranked by the statistical significance of frequent hypotension reporting, the well-known antihypertensive drugs were effectively distinguished from others (with an area under the receiver operating characteristic curve > 0.80 and a normalized discounted cumulative gain of 0.77). In addition, we found a series of antihypertensive agents (particularly drugs originally developed for treating nervous system diseases) among the drugs with top significant reporting, suggesting the good potential of Web-based and data-driven drug repositioning.

Conclusions

We found several candidate agents among the hypotension-related drugs on our list that may be redirected for lowering blood pressure. More important, we showed that a pharmacovigilance system could alternatively be used to identify antihypertensive agents and sustainably create opportunities for drug repositioning.

Hypothyroidism Side Effect in Patients Treated with Sunitinib or Sorafenib: Clinical and Structural Analyses

Tyrosine kinase inhibitors (TKIs) provide more effective targeted treatments for cancer, but are subject to a variety of adverse effects, such as hypothyroidism. TKI-induced hypothyroidism is a highly complicated issue, because of not only the unrealized toxicological mechanisms, but also different incidences of individual TKI drugs. While sunitinib is suspected for causing thyroid dysfunction more often than other TKIs, sorafenib is believed to be less risky. Here we integrated clinical data and in silico drug-protein interactions to examine the pharmacological distinction between sunitinib and sorafenib. Statistical analysis on the FDA Adverse Event Reporting System (FAERS) confirmed that sunitinib is more concurrent with hypothyroidism than sorafenib, which was observed in both female and male patients. Then, we used docking method and identified 3 proteins specifically binding to sunitinib but not sorafenib, i.e., retinoid X receptor alpha, retinoic acid receptors beta and gamma. As potential off-targets of sunitinib, these proteins are well known to assemble with thyroid hormone receptors, which can explain the profound impact of sunitinib on thyroid function. Taken together, we established a strategy of integrated analysis on clinical records and drug off-targets, which can be applied to explore the molecular basis of various adverse drug reactions.

Target-Based Drug Repositioning Using Large-Scale Chemical-Protein Interactome Data

Drug repositioning, or the identification of new indications for known drugs, is a useful strategy for drug discovery. In this study, we developed novel computational methods to predict potential drug targets and new drug indications for systematic drug repositioning using large-scale chemical-protein interactome data. We explored the target space of drugs (including primary targets and off-targets) based on chemical structure similarity and phenotypic effect similarity by making optimal use of millions of compound-protein interactions. On the basis of the target profiles of drugs, we constructed statistical models to predict new drug indications for a wide range of diseases with various molecular features. The proposed method outperformed previous methods in terms of interpretability, applicability, and accuracy. Finally, we conducted a comprehensive prediction of the drug-target-disease association network for 8270 drugs and 1401 diseases and showed biologically meaningful examples of newly predicted drug targets and drug indications. The predictive model is useful to understand the mechanisms of the predicted drug indications.

A bioinformatic and mechanistic study elicits the antifibrotic effect of ursolic acid through the attenuation of oxidative stress with the involvement of ERK, PI3K/Akt, and p38 MAPK signaling pathways in human hepatic stellate cells and rat liver

NADPH oxidases (NOXs) are a predominant mediator of redox homeostasis in hepatic stellate cells (HSCs), and oxidative stress plays an important role in the pathogenesis of liver fibrosis. Ursolic acid (UA) is a pentacyclic triterpenoid with various pharmacological activities, but the molecular targets and underlying mechanisms for its antifibrotic effect in the liver remain elusive. This study aimed to computationally predict the molecular interactome and mechanistically investigate the antifibrotic effect of UA on oxidative stress, with a focus on NOX4 activity and cross-linked signaling pathways in human HSCs and rat liver. Drug-drug interaction via chemical-protein interactome tool, a server that can predict drug-drug interaction via chemical-protein interactome, was used to predict the molecular targets of UA, and Database for Annotation, Visualization, and Integrated Discovery was employed to analyze the signaling pathways of the predicted targets of UA. The bioinformatic data showed that there were 611 molecular proteins possibly interacting with UA and that there were over 49 functional clusters responding to UA. The subsequential benchmarking data showed that UA significantly reduced the accumulation of type I collagen in HSCs in rat liver, increased the expression level of MMP-1, but decreased the expression level of TIMP-1 in HSC-T6 cells. UA also remarkably reduced the gene expression level of type I collagen in HSC-T6 cells. Furthermore, UA remarkably attenuated oxidative stress via negative regulation of NOX4 activity and expression in HSC-T6 cells. The employment of specific chemical inhibitors, SB203580, LY294002, PD98059, and AG490, demonstrated the involvement of ERK, PI3K/Akt, and p38 MAPK signaling pathways in the regulatory effect of UA on NOX4 activity and expression. Collectively, the antifibrotic effect of UA is partially due to the oxidative stress attenuating effect through manipulating NOX4 activity and expression. The results suggest that UA may act as a promising antifibrotic agent. More studies are warranted to evaluate the safety and efficacy of UA in the treatment of liver fibrosis.

A potential target of Tanshinone IIA for acute promyelocytic leukemia revealed by inverse docking and drug repurposing

Tanshinone IIA is a pharmacologically active ingredient extracted from Danshen, a Chinese traditional medicine. Its molecular mechanisms are still unclear. The present study utilized computational approaches to uncover the potential targets of this compound. In this research, PharmMapper server was used as the inverse docking tool and the results were verified by Autodock vina in PyRx 0.8, and by DRAR-CPI, a server for drug repositioning via the chemical-protein interactome. Results showed that the retinoic acid receptor alpha (RARα), a target protein in acute promyelocytic leukemia (APL), was in the top rank, with a pharmacophore model matching well the molecular features of Tanshinone IIA. Moreover, molecular docking and drug repurposing results showed that the complex was also matched in terms of structure and chemical-protein interactions. These results indicated that RARα may be a potential target of Tanshinone IIA for APL. The study can provide useful information for further biological and biochemical research on natural compounds.

Characterizing and optimizing human anticancer drug targets based on topological properties in the context of biological pathways

One of the challenging problems in drug discovery is to identify the novel targets for drugs. Most of the traditional methods for drug targets optimization focused on identifying the particular families of "druggable targets", but ignored their topological properties based on the biological pathways. In this study, we characterized the topological properties of human anticancer drug targets (ADTs) in the context of biological pathways. We found that the ADTs tended to present the following seven topological properties: influence the number of the pathways related to cancer, be localized at the start or end of the pathways, interact with cancer related genes, exhibit higher connectivity, vulnerability, betweenness, and closeness than other genes. We first ranked ADTs based on their topological property values respectively, then fused them into one global-rank using the joint cumulative distribution of an N-dimensional order statistic to optimize human ADTs. We applied the optimization method to 13 anticancer drugs, respectively. Results demonstrated that over 70% of known ADTs were ranked in the top 20%. Furthermore, the performance for mercaptopurine was significant: 6 known targets (ADSL, GMPR2, GMPR, HPRT1, AMPD3, AMPD2) were ranked in the top 15 and other four out of the top 15 (MAT2A, CDKN1A, AREG, JUN) have the potentialities to become new targets for cancer therapy.

Plumbagin elicits differential proteomic responses mainly involving cell cycle, apoptosis, autophagy, and epithelial-to-mesenchymal transition pathways in human prostate cancer PC-3 and DU145 cells

Plumbagin (PLB) has exhibited a potent anticancer effect in preclinical studies, but the molecular interactome remains elusive. This study aimed to compare the quantitative proteomic responses to PLB treatment in human prostate cancer PC-3 and DU145 cells using the approach of stable-isotope labeling by amino acids in cell culture (SILAC). The data were finally validated using Western blot assay. First, the bioinformatic analysis predicted that PLB could interact with 78 proteins that were involved in cell proliferation and apoptosis, immunity, and signal transduction. Our quantitative proteomic study using SILAC revealed that there were at least 1,225 and 267 proteins interacting with PLB and there were 341 and 107 signaling pathways and cellular functions potentially regulated by PLB in PC-3 and DU145 cells, respectively. These proteins and pathways played a critical role in the regulation of cell cycle, apoptosis, autophagy, epithelial to mesenchymal transition (EMT), and reactive oxygen species generation. The proteomic study showed substantial differences in response to PLB treatment between PC-3 and DU145 cells. PLB treatment significantly modulated the expression of critical proteins that regulate cell cycle, apoptosis, and EMT signaling pathways in PC-3 cells but not in DU145 cells. Consistently, our Western blotting analysis validated the bioinformatic and proteomic data and confirmed the modulating effects of PLB on important proteins that regulated cell cycle, apoptosis, autophagy, and EMT in PC-3 and DU145 cells. The data from the Western blot assay could not display significant differences between PC-3 and DU145 cells. These findings indicate that PLB elicits different proteomic responses in PC-3 and DU145 cells involving proteins and pathways that regulate cell cycle, apoptosis, autophagy, reactive oxygen species production, and antioxidation/oxidation homeostasis. This is the first systematic study with integrated computational, proteomic, and functional analyses revealing the networks of signaling pathways and differential proteomic responses to PLB treatment in prostate cancer cells. Quantitative proteomic analysis using SILAC represents an efficient and highly sensitive approach to identify the target networks of anticancer drugs like PLB, and the data may be used to discriminate the molecular and clinical subtypes, and to identify new therapeutic targets and biomarkers, for prostate cancer. Further studies are warranted to explore the potential of quantitative proteomic analysis in the identification of new targets and biomarkers for prostate cancer.

Xyloketal B, a marine compound, acts on a network of molecular proteins and regulates the activity and expression of rat cytochrome P450 3a: a bioinformatic and animal study

Natural compounds are becoming popular for the treatment of illnesses and health promotion, but the mechanisms of action and safety profiles are often unknown. Xyloketal B (XKB) is a novel marine compound isolated from the mangrove fungus Xylaria sp., with potent antioxidative, neuroprotective, and cardioprotective effects. However, its molecular targets and effects on drug-metabolizing enzymes are unknown. This study aimed to investigate the potential molecular targets of XKB using bioinformatic approaches and to examine the effect of XKB on the expression and activity of rat cytochrome P450 3a (Cyp3a) subfamily members using midazolam as a model probe. DDI-CPI, a server that can predict drug–drug interactions via the chemical–protein interactome, was employed to predict the targets of XKB, and the Database for Annotation, Visualization and Integrated Discovery (DAVID) was used to analyze the pathways of the predicted targets of XKB. Homology modeling was performed using the Discovery Studio program 3.1. The activity and expression of rat hepatic Cyp3a were examined after the rats were treated with XKB at 7 and 14 mg/kg for 8 consecutive days. Rat plasma concentrations of midazolam and its metabolite 1′-OH-midazolam were determined using a validated high-performance liquid chromatographic method. Bioinformatic analysis showed that there were over 324 functional proteins and 61 related signaling pathways that were potentially regulated by XKB. A molecular docking study showed that XKB bound to the active site of human cytochrome P450 3A4 and rat Cyp3a2 homology model via the formation of hydrogen bonds. The in vivo study showed that oral administration of XKB at 14 mg/kg to rats for 8 days significantly increased the area under the plasma concentration-time curve (AUC) of midazolam, with a concomitant decrease in the plasma clearance and AUC ratio of 1′-OH-midazolam over midazolam. Further, oral administration of 14 mg/kg XKB for 8 days markedly reduced the activity and expression of hepatic Cyp3a in rats. Taken together, the results show that XKB could regulate networks of molecular proteins and related signaling pathways and that XKB downregulated hepatic Cyp3a in rats. XKB might cause drug interactions through modulation of the activity and expression of Cyp3a members. More studies are warranted to confirm the mechanisms of action of XKB and to investigate the underlying mechanism for the regulating effect of XKB on Cyp3a subfamily members.

Systematic prediction of drug combinations based on clinical side-effects

Drug co-prescription (or drug combination) is a therapeutic strategy widely used as it may improve efficacy and reduce side-effect (SE). Since it is impractical to screen all possible drug combinations for every indication, computational methods have been developed to predict new combinations. In this study, we describe a novel approach that utilizes clinical SEs from post-marketing surveillance and the drug label to predict 1,508 novel drug-drug combinations. It outperforms other prediction methods, achieving an AUC of 0.92 compared to an AUC of 0.69 in a previous method, on a much larger drug combination set (245 drug combinations in our dataset compared to 75 in previous work.). We further found from the feature selection that three FDA black-box warned serious SEs, namely pneumonia, haemorrhage rectum, and retinal bleeding, contributed mostly to the predictions and a model only using these three SEs can achieve an average area under curve (AUC) at 0.80 and accuracy at 0.91, potentially with its simplicity being recognized as a practical rule-of-three in drug co-prescription or making fixed-dose drug combination. We also demonstrate this performance is less likely to be influenced by confounding factors such as biased disease indications or chemical structures.

A network approach to clinical intervention in neurodegenerative diseases

Network biology has become a powerful tool to dissect the molecular mechanisms triggering neurodegeneration. Recent developments in network biology have led to the discovery of disease-causing genes, diagnostic biomarkers, and therapeutic targets for several neurodegenerative diseases including Alzheimer's, Parkinson's, and Huntington's diseases. Network-based approaches have provided the molecular rationale for the relationship among cancer, diabetes, and neurodegenerative diseases, and have uncovered unexpected links between apparently unrelated diseases. Here, we summarize the recent advances in network biology to untangle the molecular underpinnings giving rise to the most prevalent neurodegenerative diseases. We propose that network analysis provides a feasible and practical tool for identifying biologically meaningful biomarkers and potential therapeutic targets for clinical intervention in neurodegenerative diseases.

DDI-CPI, a server that predicts drug-drug interactions through implementing the chemical-protein interactome

Drug–drug interactions (DDIs) may cause serious side-effects that draw great attention from both academia and industry. Since some DDIs are mediated by unexpected drug–human protein interactions, it is reasonable to analyze the chemical–protein interactome (CPI) profiles of the drugs to predict their DDIs. Here we introduce the DDI-CPI server, which can make real-time DDI predictions based only on molecular structure. When the user submits a molecule, the server will dock user's molecule across 611 human proteins, generating a CPI profile that can be used as a feature vector for the pre-constructed prediction model. It can suggest potential DDIs between the user's molecule and our library of 2515 drug molecules. In cross-validation and independent validation, the server achieved an AUC greater than 0.85. Additionally, by investigating the CPI profiles of predicted DDI, users can explore the PK/PD proteins that might be involved in a particular DDI. A 3D visualization of the drug-protein interaction will be provided as well. The DDI-CPI is freely accessible at http://cpi.bio-x.cn/ddi/.

Distinct DNA-based epigenetic switches trigger transcriptional activation of silent genes in human dermal fibroblasts

The influential role of the epigenome in orchestrating genome-wide transcriptional activation instigates the demand for the artificial genetic switches with distinct DNA sequence recognition. Recently, we developed a novel class of epigenetically active small molecules called SAHA-PIPs by conjugating selective DNA binding pyrrole-imidazole polyamides (PIPs) with the histone deacetylase inhibitor SAHA. Screening studies revealed that certain SAHA-PIPs trigger targeted transcriptional activation of pluripotency and germ cell genes in mouse and human fibroblasts, respectively. Through microarray studies and functional analysis, here we demonstrate for the first time the remarkable ability of thirty-two different SAHA-PIPs to trigger the transcriptional activation of exclusive clusters of genes and noncoding RNAs. QRT-PCR validated the microarray data, and some SAHA-PIPs activated therapeutically significant genes like KSR2. Based on the aforementioned results, we propose the potential use of SAHA-PIPs as reagents capable of targeted transcriptional activation.

Prediction of effective drug combinations by chemical interaction, protein interaction and target enrichment of KEGG pathways

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.

Plumbagin induces the apoptosis of human tongue carcinoma cells through the mitochondria-mediated pathway

BACKGROUND

Plumbagin, a quinonoid constituent isolated from the root of Plumbago zeylanica L., has been proven to possess anti-tumor activity both in vitro and in vivo. However, its anti-tumor properties for human tongue carcinoma have not been reported. This study aimed to investigate the inhibitory effect and the underlying mechanism of plumbagin on the growth of human tongue carcinoma cells.

MATERIAL AND METHODS

Cell proliferation ability was detected by EdU incorporation assay and colony formation assay. Cell-cycle distribution was determined by flow cytometric analysis using propidium iodide (PI) staining. Cellular apoptosis was then evaluated by flow cytometry and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay. Western blotting was applied to assay the expression of Bax and Bcl-2.

RESULTS

Plumbagin inhibited the growth and proliferation of Tca8113 cells in vitro in a concentration- and time-dependent manner. The cell cycles of plumbagin-treated Tca8113 cells were arrested at the G2/M phase. Cells treated with plumbagin presented the characteristic morphological changes of apoptosis. The ratio of Bax/Bcl-2 was raised by plumbagin in a concentration-dependent manner.

CONCLUSIONS

These results indicate that plumbagin induces the apoptosis of Tca8113 cells through mitochondria-mediated pathway.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

Life beyond kinases: structure-based discovery of sorafenib as nanomolar antagonist of 5-HT receptors

Of great interest in recent years has been computationally predicting the novel polypharmacology of drug molecules. Here, we applied an "induced-fit" protocol to improve the homology models of 5-HT(2A) receptor, and we assessed the quality of these models in retrospective virtual screening. Subsequently, we computationally screened the FDA approved drug molecules against the best induced-fit 5-HT(2A) models and chose six top scoring hits for experimental assays. Surprisingly, one well-known kinase inhibitor, sorafenib, has shown unexpected promiscuous 5-HTRs binding affinities, K(i) = 1959, 56, and 417 nM against 5-HT(2A), 5-HT(2B), and 5-HT(2C), respectively. Our preliminary SAR exploration supports the predicted binding mode and further suggests sorafenib to be a novel lead compound for 5HTR ligand discovery. Although it has been well-known that sorafenib produces anticancer effects through targeting multiple kinases, carefully designed experimental studies are desirable to fully understand whether its "off-target" 5-HTR binding activities contribute to its therapeutic efficacy or otherwise undesirable side effects.

Examining the interactome of huperzine A by magnetic biopanning

Huperzine A is a bioactive compound derived from traditional Chinese medicine plant Qian Ceng Ta (Huperzia serrata), and was found to have multiple neuroprotective effects. In addition to being a potent acetylcholinesterase inhibitor, it was thought to act through other mechanisms such as antioxidation, antiapoptosis, etc. However, the molecular targets involved with these mechanisms were not identified. In this study, we attempted to exam the interactome of Huperzine A using a cDNA phage display library and also mammalian brain tissue extracts. The drugs were chemically linked on the surface of magnetic particles and the interactive phages or proteins were collected and analyzed. Among the various cDNA expressing phages selected, one was identified to encode the mitochondria NADH dehydrogenase subunit 1. Specific bindings between the drug and the target phages and target proteins were confirmed. Another enriched phage clone was identified as mitochondria ATP synthase, which was also panned out from the proteome of mouse brain tissue lysate. These data indicated the possible involvement of mitochondrial respiratory chain matrix enzymes in Huperzine A's pharmacological effects. Such involvement had been suggested by previous studies based on enzyme activity changes. Our data supported the new mechanism. Overall we demonstrated the feasibility of using magnetic biopanning as a simple and viable method for investigating the complex molecular mechanisms of bioactive molecules.

Drug repositioning for personalized medicine

Human diseases can be caused by complex mechanisms involving aberrations in numerous proteins and pathways. With recent advances in genomics, elucidating the molecular basis of disease on a personalized level has become an attainable goal. In many cases, relevant molecular targets will be identified for which approved drugs already exist, and the potential repositioning of these drugs to a new indication can be investigated. Repositioning is an accelerated route for drug discovery because existing drugs have established clinical and pharmacokinetic data. Personalized medicine and repositioning both aim to improve the productivity of current drug discovery pipelines, which expend enormous time and cost to develop new drugs, only to have them fail in clinical trials because of lack of efficacy or toxicity. Here, we discuss the current state of research in these two fields, focusing on recent large-scale efforts to systematically find repositioning candidates and elucidate individual disease mechanisms in cancer. We also discuss scenarios in which personalized drug repositioning could be particularly rewarding, such as for diseases that are rare or have specific mutations, as well as current challenges in this field. With an increasing number of drugs being approved for rare cancer subtypes, personalized medicine and repositioning approaches are poised to significantly alter the way we diagnose diseases, infer treatments and develop new drugs.

Prediction of drug combinations by integrating molecular and pharmacological data

Combinatorial therapy is a promising strategy for combating complex disorders due to improved efficacy and reduced side effects. However, screening new drug combinations exhaustively is impractical considering all possible combinations between drugs. Here, we present a novel computational approach to predict drug combinations by integrating molecular and pharmacological data. Specifically, drugs are represented by a set of their properties, such as their targets or indications. By integrating several of these features, we show that feature patterns enriched in approved drug combinations are not only predictive for new drug combinations but also provide insights into mechanisms underlying combinatorial therapy. Further analysis confirmed that among our top ranked predictions of effective combinations, 69% are supported by literature, while the others represent novel potential drug combinations. We believe that our proposed approach can help to limit the search space of drug combinations and provide a new way to effectively utilize existing drugs for new purposes.

Drug repositioning for orphan diseases

The need and opportunity to discover therapeutics for rare or orphan diseases are enormous. Due to limited prevalence and/or commercial potential, of the approximately 6000 orphan diseases (defined by the FDA Orphan Drug Act as <200 000 US prevalence), only a small fraction (5%) is of interest to the biopharmaceutical industry. The fact that drug development is complicated, time-consuming and expensive with extremely low success rates only adds to the low rate of therapeutics available for orphan diseases. An alternative and efficient strategy to boost the discovery of orphan disease therapeutics is to find connections between an existing drug product and orphan disease. Drug Repositioning or Drug Repurposing--finding a new indication for a drug--is one way to maximize the potential of a drug. The advantages of this approach are manifold, but rational drug repositioning for orphan diseases is not trivial and poses several formidable challenges--pharmacologically and computationally. Most of the repositioned drugs currently in the market are the result of serendipity. One reason the connection between drug candidates and their potential new applications are not identified in an earlier or more systematic fashion is that the underlying mechanism 'connecting' them is either very intricate and unknown or indirect or dispersed and buried in an ever-increasing sea of information, much of which is emerging only recently and therefore is not well organized. In this study, we will review some of these issues and the current methodologies adopted or proposed to overcome them and translate chemical and biological discoveries into safe and effective orphan disease therapeutics.

Exploring off-targets and off-systems for adverse drug reactions via chemical-protein interactome--clozapine-induced agranulocytosis as a case study

In the era of personalized medical practice, understanding the genetic basis of patient-specific adverse drug reaction (ADR) is a major challenge. Clozapine provides effective treatments for schizophrenia but its usage is limited because of life-threatening agranulocytosis. A recent high impact study showed the necessity of moving clozapine to a first line drug, thus identifying the biomarkers for drug-induced agranulocytosis has become important. Here we report a methodology termed as antithesis chemical-protein interactome (CPI), which utilizes the docking method to mimic the differences in the drug-protein interactions across a panel of human proteins. Using this method, we identified HSPA1A, a known susceptibility gene for CIA, to be the off-target of clozapine. Furthermore, the mRNA expression of HSPA1A-related genes (off-target associated systems) was also found to be differentially expressed in clozapine treated leukemia cell line. Apart from identifying the CIA causal genes we identified several novel candidate genes which could be responsible for agranulocytosis. Proteins related to reactive oxygen clearance system, such as oxidoreductases and glutathione metabolite enzymes, were significantly enriched in the antithesis CPI. This methodology conducted a multi-dimensional analysis of drugs' perturbation to the biological system, investigating both the off-targets and the associated off-systems to explore the molecular basis of an adverse event or the new uses for old drugs.

Idiosyncratic drug reactions (IDR) generally cannot be identified until after a drug is taken by a large population, but usually result in restricted use or withdrawal. Clozapine provides the most effective treatment for schizophrenia but its use is limited because of a life-threatening IDR, i.e., the agranulocytosis. A high impact clinical study demonstrated the necessity of moving clozapine from 3^rd line to 1^st line drug; therefore, intensive research has aimed at identifying genes responsible for clozapine-induced agranulocytosis (CIA). Olanzapine, an analog of clozapine, has much lower incidence of agranulocytosis. Based on this phenomenon, we proposed an in silico methodology termed as antithesis chemical-protein interactome (CPI), which mimics the differences in the drug-protein interactions of the two drugs across a panel of human proteins. e.g., HSPA1A was identified to be targeted by clozapine not olanzapine. Furthermore, the gene expression of the HSPA1A-related gene system was also found up-regulated after clozapine treatment. This approach can examine the system's perturbation in terms of both the off-target and the off-system's interaction with the drug, providing theoretical basis for decoding the adverse drug reactions or the new uses for old drugs.

Chemical-protein interactome and its application in off-target identification

Drugs exert their therapeutic and adverse effects by interacting with molecular targets. Although designed to interact with specific targets in a desirable manner, drug molecules often bind to unexpected proteins (off-targets). By activating or inhibiting off-targets and the associated biological processes and pathways, the resulting chemical-protein interactions can influence drug reaction directly or indirectly. Exploring the relationship between drug and off-targets and the downstream drug reaction can help understand the polypharmacology of the drug, hence significantly advance the drug repositioning pipeline and the application of personalized medicine in understanding and preventing adverse drug reaction. This review summarizes works on predicting off-targets via chemical-protein interactome (CPI), an interaction strength matrix of drugs across multiple human proteins aiming at exploring the unexpected drug-protein interactions, with a variety of computational strategies, including docking, chemical structure comparison and text-mining etc. Effective recall on previous knowledge, de novo prediction and subsequent experimental validation conferred us strong confidence in these methods. Such studies present prospect of large scale in silico methodologies for off-target discovery with low cost and high efficiency.