Network analyses in systems pharmacology

Systems-Pharmacology Dissection of Traditional Chinese Medicine Compound Saffron Formula Reveals Multi-scale Treatment Strategy for Cardiovascular Diseases

Cardiovascular diseases (CVDs) have been regarding as “the world’s first killer” of human beings in recent years owing to the striking morbidity and mortality, the involved molecular mechanisms are extremely complex and remain unclear. Traditional Chinese medicine (TCM) adheres to the aim of combating complex diseases from an integrative and holistic point of view, which has shown effectiveness in CVDs therapy. However, system-level understanding of such a mechanism of multi-scale treatment strategy for CVDs is still difficult. Here, we developed a system pharmacology approach with the purpose of revealing the underlying molecular mechanisms exemplified by a famous compound saffron formula (CSF) in treating CVDs. First, by systems ADME analysis combined with drug targeting process, 103 potential active components and their corresponding 219 direct targets were retrieved and some key interactions were further experimentally validated. Based on this, the network relationships among active components, targets and diseases were further built to uncover the pharmacological actions of the drug. Finally, a “CVDs pathway” consisted of several regulatory modules was incorporated to dissect the therapeutic effects of CSF in different pathological features-relevant biological processes. All this demonstrates CSF has multi-scale curative activity in regulating CVD-related biological processes, which provides a new potential way for modern medicine in the treatment of complex diseases.

BindingDB in 2015: A public database for medicinal chemistry, computational chemistry and systems pharmacology

BindingDB, www.bindingdb.org, is a publicly accessible database of experimental protein-small molecule interaction data. Its collection of over a million data entries derives primarily from scientific articles and, increasingly, US patents. BindingDB provides many ways to browse and search for data of interest, including an advanced search tool, which can cross searches of multiple query types, including text, chemical structure, protein sequence and numerical affinities. The PDB and PubMed provide links to data in BindingDB, and vice versa; and BindingDB provides links to pathway information, the ZINC catalog of available compounds, and other resources. The BindingDB website offers specialized tools that take advantage of its large data collection, including ones to generate hypotheses for the protein targets bound by a bioactive compound, and for the compounds bound by a new protein of known sequence; and virtual compound screening by maximal chemical similarity, binary kernel discrimination, and support vector machine methods. Specialized data sets are also available, such as binding data for hundreds of congeneric series of ligands, drawn from BindingDB and organized for use in validating drug design methods. BindingDB offers several forms of programmatic access, and comes with extensive background material and documentation. Here, we provide the first update of BindingDB since 2007, focusing on new and unique features and highlighting directions of importance to the field as a whole.

A Dormant Microbial Component in the Development of Preeclampsia

Preeclampsia (PE) is a complex, multisystem disorder that remains a leading cause of morbidity and mortality in pregnancy. Four main classes of dysregulation accompany PE and are widely considered to contribute to its severity. These are abnormal trophoblast invasion of the placenta, anti-angiogenic responses, oxidative stress, and inflammation. What is lacking, however, is an explanation of how these themselves are caused. We here develop the unifying idea, and the considerable evidence for it, that the originating cause of PE (and of the four classes of dysregulation) is, in fact, microbial infection, that most such microbes are dormant and hence resist detection by conventional (replication-dependent) microbiology, and that by occasional resuscitation and growth it is they that are responsible for all the observable sequelae, including the continuing, chronic inflammation. In particular, bacterial products such as lipopolysaccharide (LPS), also known as endotoxin, are well known as highly inflammagenic and stimulate an innate (and possibly trained) immune response that exacerbates the inflammation further. The known need of microbes for free iron can explain the iron dysregulation that accompanies PE. We describe the main routes of infection (gut, oral, and urinary tract infection) and the regularly observed presence of microbes in placental and other tissues in PE. Every known proteomic biomarker of "preeclampsia" that we assessed has, in fact, also been shown to be raised in response to infection. An infectious component to PE fulfills the Bradford Hill criteria for ascribing a disease to an environmental cause and suggests a number of treatments, some of which have, in fact, been shown to be successful. PE was classically referred to as endotoxemia or toxemia of pregnancy, and it is ironic that it seems that LPS and other microbial endotoxins really are involved. Overall, the recognition of an infectious component in the etiology of PE mirrors that for ulcers and other diseases that were previously considered to lack one.

Neuropharmacology beyond reductionism - A likely prospect

Neuropharmacology had several major past successes, but the last few decades did not witness any leap forward in the drug treatment of brain disorders. Moreover, current drugs used in neurology and psychiatry alleviate the symptoms, while hardly curing any cause of disease, basically because the etiology of most neuro-psychic syndromes is but poorly known. This review argues that this largely derives from the unbalanced prevalence in neuroscience of the analytic reductionist approach, focused on the cellular and molecular level, while the understanding of integrated brain activities remains flimsier. The decline of drug discovery output in the last decades, quite obvious in neuropharmacology, coincided with the advent of the single target-focused search of potent ligands selective for a well-defined protein, deemed critical in a given pathology. However, all the widespread neuro-psychic troubles are multi-mechanistic and polygenic, their complex etiology making unsuited the single-target drug discovery. An evolving approach, based on systems biology considers that a disease expresses a disturbance of the network of interactions underlying organismic functions, rather than alteration of single molecular components. Accordingly, systems pharmacology seeks to restore a disturbed network via multi-targeted drugs. This review notices that neuropharmacology in fact relies on drugs which are multi-target, this feature having occurred just because those drugs were selected by phenotypic screening in vivo, or emerged from serendipitous clinical observations. The novel systems pharmacology aims, however, to devise ab initio multi-target drugs that will appropriately act on multiple molecular entities. Though this is a task much more complex than the single-target strategy, major informatics resources and computational tools for the systemic approach of drug discovery are already set forth and their rapid progress forecasts promising outcomes for neuropharmacology.

A novel algorithm for analyzing drug-drug interactions from MEDLINE literature

Drug–drug interaction (DDI) is becoming a serious clinical safety issue as the use of multiple medications becomes more common. Searching the MEDLINE database for journal articles related to DDI produces over 330,000 results. It is impossible to read and summarize these references manually. As the volume of biomedical reference in the MEDLINE database continues to expand at a rapid pace, automatic identification of DDIs from literature is becoming increasingly important. In this article, we present a random-sampling-based statistical algorithm to identify possible DDIs and the underlying mechanism from the substances field of MEDLINE records. The substances terms are essentially carriers of compound (including protein) information in a MEDLINE record. Four case studies on warfarin, ibuprofen, furosemide and sertraline implied that our method was able to rank possible DDIs with high accuracy (90.0% for warfarin, 83.3% for ibuprofen, 70.0% for furosemide and 100% for sertraline in the top 10% of a list of compounds ranked by p-value). A social network analysis of substance terms was also performed to construct networks between proteins and drug pairs to elucidate how the two drugs could interact.

Applicability of gene expression and systems biology to develop pharmacogenetic predictors; antipsychotic-induced extrapyramidal symptoms as an example

Pharmacogenetics has been driven by a candidate gene approach. The disadvantage of this approach is that is limited by our current understanding of the mechanisms by which drugs act. Gene expression could help to elucidate the molecular signatures of antipsychotic treatments searching for dysregulated molecular pathways and the relationships between gene products, especially protein-protein interactions. To embrace the complexity of drug response, machine learning methods could help to identify gene-gene interactions and develop pharmacogenetic predictors of drug response. The present review summarizes the applicability of the topics presented here (gene expression, network analysis and gene-gene interactions) in pharmacogenetics. In order to achieve this, we present an example of identifying genetic predictors of extrapyramidal symptoms induced by antipsychotic.

Systems pharmacology to decipher the combinational anti-migraine effects of Tianshu formula

Migraine is the most common neurovascular disorder that imparts a considerable burden to health care system around the world. However, currently there are still no effective and widely applicable pharmacotherapies for migraine patients. Herbal formulae, characterized as multiple herbs, constituents and targets, have been acknowledged with clinical effects in treating migraine, which attract more and more researchers' attention although their exact molecular mechanisms are still unclear. In this work, a novel systems pharmacology-based method which integrates pharmacokinetic filtering, target fishing and network analysis was developed and exemplified by a probe, i.e. Tianshu formula, a widely clinically used anti-migraine herbal formula in China which comprises of Rhizoma chuanxiong and Gastrodia elata. The results exhibit that 20 active ingredients of Tianshu formula possess favorable pharmacokinetic profiles, which have interactions with 48 migraine-related targets to provide potential synergistic therapeutic effects. Additionally, from systematic analysis, we speculate that R. chuanxiong as the monarch herb mediates the major targets like PTGS2, ESR1, NOS2, HTR1B and NOS3 to regulate the vascular and nervous systems, as well as the inflammation and pain-related pathways to benefit migraine patients. Meanwhile, as an adjuvant herb, G. elata may not only assist the monarch herb to improve the outcome of migraine patients, but also regulate multiple targets like ABAT, HTR1D, ALOX15 and KCND3 to modify migraine accompanying symptoms like vomiting, vertigo and gastrointestinal disorders.

Updates on drug-target network; facilitating polypharmacology and data integration by growth of DrugBank database

Network pharmacology elucidates the relationship between drugs and targets. As the identified targets for each drug increases, the corresponding drug-target network (DTN) evolves from solely reflection of the pharmaceutical industry trend to a portrait of polypharmacology. The aim of this study was to evaluate the potentials of DrugBank database in advancing systems pharmacology. We constructed and analyzed DTN from drugs and targets associations in the DrugBank 4.0 database. Our results showed that in bipartite DTN, increased ratio of identified targets for drugs augmented density and connectivity of drugs and targets and decreased modular structure. To clear up the details in the network structure, the DTNs were projected into two networks namely, drug similarity network (DSN) and target similarity network (TSN). In DSN, various classes of Food and Drug Administration-approved drugs with distinct therapeutic categories were linked together based on shared targets. Projected TSN also showed complexity because of promiscuity of the drugs. By including investigational drugs that are currently being tested in clinical trials, the networks manifested more connectivity and pictured the upcoming pharmacological space in the future years. Diverse biological processes and protein-protein interactions were manipulated by new drugs, which can extend possible target combinations. We conclude that network-based organization of DrugBank 4.0 data not only reveals the potential for repurposing of existing drugs, also allows generating novel predictions about drugs off-targets, drug-drug interactions and their side effects. Our results also encourage further effort for high-throughput identification of targets to build networks that can be integrated into disease networks.

An efficient and unbiased method for sensitivity analysis of stochastic reaction networks

We consider the problem of estimating parameter sensitivity for Markovian models of reaction networks. Sensitivity values measure the responsiveness of an output with respect to the model parameters. They help in analysing the network, understanding its robustness properties and identifying the important reactions for a specific output. Sensitivity values are commonly estimated using methods that perform finite-difference computations along with Monte Carlo simulations of the reaction dynamics. These methods are computationally efficient and easy to implement, but they produce a biased estimate which can be unreliable for certain applications. Moreover, the size of the bias is generally unknown and hence the accuracy of these methods cannot be easily determined. There also exist unbiased schemes for sensitivity estimation but these schemes can be computationally infeasible, even for very simple networks. Our goal in this paper is to present a new method for sensitivity estimation, which combines the computational efficiency of finite-difference methods with the accuracy of unbiased schemes. Our method is easy to implement and it relies on an exact representation of parameter sensitivity that we recently proved elsewhere. Through examples, we demonstrate that the proposed method can outperform the existing methods, both biased and unbiased, in many situations.

Evaluating the antidiabetic effects of Chinese herbal medicine: Xiao-Ke-An in 3T3-L1 cells and KKAy mice using both conventional and holistic omics approaches

BACKGROUND

Xiao-Ke-An (XKA) is a Chinese medicine widely used for treating type 2 diabetes mellitus (T2D). It is composed of eight herbal medicines traditionally used for T2D, including Rehmannia glutinosa Libosch, Anemarrhena asphodeloides Bunge, Coptis chinensis Franch, etc. The aim of the present study was to investigate the antidiabetic effects of XKA with both conventional and holistic omics approaches.

METHODS

The antidiabetic effect of XKA was first investigated in 3T3-L1 cells to study the effect of XKA on adipogenesis in vitro. Oil Red O staining was performed to determine the lipid accumulation. The intracellular total cholesterol (TC) and triglyceride (TG) contents in XKA treated 3T3-L1 cells were also evaluated. The therapeutic effects of XKA was further evaluated in KKAy mice with both conventional and holistic omics approaches. Body weight, fasting and non-fasting blood glucose, and oral glucose tolerance were measured during the experiment. At the time of sacrifice, serum was collected for the measurement of TG, TC, high-density lipoprotein cholesterol (HDL-c) and low-density lipoprotein cholesterol (LDL-c). The liver, kidney, spleen, pancreas, heart and adipose tissues were harvested and weighted. The liver was used for further microarray experiment. Omics approaches were adopted to evaluate the holistic rebalancing effect of XKA at molecular network level.

RESULTS

XKA significantly inhibited adipogenic differentiation, lowered the intracellular TC and TG contents in 3T3-L1 cells. XKA improved the glucose homeostasis and lipid metabolism, ameliorated insulin resistance in KKAy mice. Furthermore, XKA also exhibited effective therapeutic effects by reversing the molecular T2D disease network from an unbalanced state.

CONCLUSIONS

This study investigated the antidiabetic effects of XKA with both conventional and holistic omics approaches, providing both phenotypic evidence and underlying action mechanisms for the clinical use of XKA treating T2D.

CancerHSP: anticancer herbs database of systems pharmacology

The numerous natural products and their bioactivity potentially afford an extraordinary resource for new drug discovery and have been employed in cancer treatment. However, the underlying pharmacological mechanisms of most natural anticancer compounds remain elusive, which has become one of the major obstacles in developing novel effective anticancer agents. Here, to address these unmet needs, we developed an anticancer herbs database of systems pharmacology (CancerHSP), which records anticancer herbs related information through manual curation. Currently, CancerHSP contains 2439 anticancer herbal medicines with 3575 anticancer ingredients. For each ingredient, the molecular structure and nine key ADME parameters are provided. Moreover, we also provide the anticancer activities of these compounds based on 492 different cancer cell lines. Further, the protein targets of the compounds are predicted by state-of-art methods or collected from literatures. CancerHSP will help reveal the molecular mechanisms of natural anticancer products and accelerate anticancer drug development, especially facilitate future investigations on drug repositioning and drug discovery. CancerHSP is freely available on the web at http://lsp.nwsuaf.edu.cn/CancerHSP.php.

A weighted and integrated drug-target interactome: drug repurposing for schizophrenia as a use case

BACKGROUND

Computational pharmacology can uniquely address some issues in the process of drug development by providing a macroscopic view and a deeper understanding of drug action. Specifically, network-assisted approach is promising for the inference of drug repurposing. However, the drug-target associations coming from different sources and various assays have much noise, leading to an inflation of the inference errors. To reduce the inference errors, it is necessary and critical to create a comprehensive and weighted data set of drug-target associations.

RESULTS

In this study, we created a weighted and integrated drug-target interactome (WinDTome) to provide a comprehensive resource of drug-target associations for computational pharmacology. We first collected drug-target interactions from six commonly used drug-target centered data sources including DrugBank, KEGG, TTD, MATADOR, PDSP K(i) Database, and BindingDB. Then, we employed the record linkage method to normalize drugs and targets to the unique identifiers by utilizing the public data sources including PubChem, Entrez Gene, and UniProt. To assess the reliability of the drug-target associations, we assigned two scores (Score_S and Score_R) to each drug-target association based on their data sources and publication references. Consequently, the WinDTome contains 546,196 drug-target associations among 303,018 compounds and 4,113 genes. To assess the application of the WinDTome, we designed a network-based approach for drug repurposing using mental disorder schizophrenia (SCZ) as a case. Starting from 41 known SCZ drugs and their targets, we inferred a total of 264 potential SCZ drugs through the associations of drug-target with Score_S higher than two in WinDTome and human protein-protein interactions. Among the 264 SCZ-related drugs, 39 drugs have been investigated in clinical trials for SCZ treatment and 74 drugs for the treatment of other mental disorders, respectively. Compared with the results using other Score_S cutoff values, single data source, or the data from STITCH, the inference of 264 SCZ-related drugs had the highest performance.

CONCLUSIONS

The WinDTome generated in this study contains comprehensive drug-target associations with confidence scores. Its application to the SCZ drug repurposing demonstrated that the WinDTome is promising to serve as a useful resource for drug repurposing.

Network pharmacology study on the mechanism of traditional Chinese medicine for upper respiratory tract infection

Traditional Chinese medicine (TCM) is a multi-component and multi-target agent and could treat complex diseases in a holistic way, especially infection diseases. However, the underlying pharmacology remains unclear. Fortunately, network pharmacology by integrating system biology and polypharmacology provides a strategy to address this issue. In this work, Reduning Injection (RDN), a well-used TCM treatment in the clinic for upper respiratory tract infections (URTIs), was investigated to interpret the molecular mechanism and predict new clinical directions by integrating molecular docking, network analysis and cell-based assays. 32 active ingredients and 38 potential targets were identified. In vitro experiments confirmed the bioactivities of the compounds against lipopolysaccharide (LPS)-stimulated PGE2 and NO production in RAW264.7 cells. Moreover, network analysis showed that RDN could not only inhibit viral replication but also alleviate the sickness symptoms of URTIs through directly targeting the key proteins in the respiratory viral life cycle and indirectly regulating host immune systems. In addition, other clinical applications of RDN such as neoplasms, cardiovascular diseases and immune system diseases were predicted on the basis of the relationships between targets and diseases.

Detection and characterization of nonspecific, sparsely populated binding modes in the early stages of complexation

A method is proposed to study protein-ligand binding in a system governed by specific and nonspecific interactions. Strong associations lead to narrow distributions in the proteins configuration space; weak and ultraweak associations lead instead to broader distributions, a manifestation of nonspecific, sparsely populated binding modes with multiple interfaces. The method is based on the notion that a discrete set of preferential first-encounter modes are metastable states from which stable (prerelaxation) complexes at equilibrium evolve. The method can be used to explore alternative pathways of complexation with statistical significance and can be integrated into a general algorithm to study protein interaction networks. The method is applied to a peptide-protein complex. The peptide adopts several low-population conformers and binds in a variety of modes with a broad range of affinities. The system is thus well suited to analyze general features of binding, including conformational selection, multiplicity of binding modes, and nonspecific interactions, and to illustrate how the method can be applied to study these problems systematically. The equilibrium distributions can be used to generate biasing functions for simulations of multiprotein systems from which bulk thermodynamic quantities can be calculated.

Systems pharmacology-based drug discovery for marine resources: an example using sea cucumber (Holothurians)

ETHNOPHARMACOLOGICAL RELEVANCE

Sea cucumber, a kind of marine animal, have long been utilized as tonic and traditional remedies in the Middle East and Asia because of its effectiveness against hypertension, asthma, rheumatism, cuts and burns, impotence, and constipation. In this study, an overall study performed on sea cucumber was used as an example to show drug discovery from marine resource by using systems pharmacology model. The value of marine natural resources has been extensively considered because these resources can be potentially used to treat and prevent human diseases. However, the discovery of drugs from oceans is difficult, because of complex environments in terms of composition and active mechanisms. Thus, a comprehensive systems approach which could discover active constituents and their targets from marine resource, understand the biological basis for their pharmacological properties is necessary.

MATERIALS AND METHODS

In this study, a feasible pharmacological model based on systems pharmacology was established to investigate marine medicine by incorporating active compound screening, target identification, and network and pathway analysis.

RESULTS

As a result, 106 candidate components of sea cucumber and 26 potential targets were identified. Furthermore, the functions of sea cucumber in health improvement and disease treatment were elucidated in a holistic way based on the established compound-target and target-disease networks, and incorporated pathways.

CONCLUSIONS

This study established a novel strategy that could be used to explore specific active mechanisms and discover new drugs from marine sources.

Cytoprotection of human endothelial cells from oxidant stress with CDDO derivatives: network analysis of genes responsible for cytoprotection

AIM

To identify drugs that may reduce the impact of oxidant stress on cell viability.

METHODS

Human umbilical vein endothelial cells were treated with 200 nmol/l CDDO-Im (imidazole) and CDDO-Me (methyl) after exposure to menadione and compared to vehicle-treated cells. Cell viability and cytotoxicity were assessed, and gene expression profiling was performed.

RESULTS

CDDO-Im was significantly more cytoprotective and less cytotoxic (p < 0.001) than CDDO-Me. Although both provided cytoprotection by induction of gene transcription, CDDO-Im induced more genes. In addition to a higher induction of the key cytoprotective gene heme oxygenase-1 (38.9-fold increase for CDDO-Im and 26.5-fold increase for CDDO-Me), CDDO-Im also induced greater expression of heat shock proteins (5.5-fold increase compared to 2.8-fold for CDDO-Me).

CONCLUSIONS

Both compounds showed good induction of heme oxygenase, which largely accounted for their cytoprotective effect. Differences were detected in cytotoxicity at higher doses, indicating that CDDO-Me was more cytotoxic than CDDO-Im. Significant differences were detected in the ability of CDDO-Im and CDDO-Me to affect differential gene transcription. CDDO-Im induced more genes than did CDDO-Me. The source of the differences in gene expression patterns between CDDO-Im and CDDO-Me was not determined but may be important in long-term use of this class of drugs.

Integrative analysis of kinase networks in TRAIL-induced apoptosis provides a source of potential targets for combination therapy

Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) is an endogenous secreted peptide and, in preclinical studies, preferentially induces apoptosis in tumor cells rather than in normal cells. The acquisition of resistance in cells exposed to TRAIL or its mimics limits their clinical efficacy. Because kinases are intimately involved in the regulation of apoptosis, we systematically characterized kinases involved in TRAIL signaling. Using RNA interference (RNAi) loss-of-function and cDNA overexpression screens, we identified 169 protein kinases that influenced the dynamics of TRAIL-induced apoptosis in the colon adenocarcinoma cell line DLD-1. We classified the kinases as sensitizers or resistors or modulators, depending on the effect that knockdown and overexpression had on TRAIL-induced apoptosis. Two of these kinases that were classified as resistors were PX domain-containing serine/threonine kinase (PXK) and AP2-associated kinase 1 (AAK1), which promote receptor endocytosis and may enable cells to resist TRAIL-induced apoptosis by enhancing endocytosis of the TRAIL receptors. We assembled protein interaction maps using mass spectrometry-based protein interaction analysis and quantitative phosphoproteomics. With these protein interaction maps, we modeled information flow through the networks and identified apoptosis-modifying kinases that are highly connected to regulated substrates downstream of TRAIL. The results of this analysis provide a resource of potential targets for the development of TRAIL combination therapies to selectively kill cancer cells.

Colorectal cancer drug target prediction using ontology-based inference and network analysis

Identification of novel drug targets is a critical step in drug development. Many recent studies have produced multiple types of data, which provides an opportunity to mine the relationships among them to predict drug targets. In this study, we present a novel integrative approach that combines ontology reasoning with network-assisted gene ranking to predict new drug targets. We utilized colorectal cancer (CRC) as a proof-of-concept use case to illustrate the approach. Starting from FDA-approved CRC drugs and the relationships among disease, drug, gene, pathway, and SNP in an ontology representing PharmGKB data, we inferred 113 potential CRC drug targets. We further prioritized these genes based on their relationships with CRC disease genes in the context of human protein-protein interaction networks. Thus, among the 113 potential drug targets, 15 were selected as the promising drug targets, including some genes that are supported by previous studies. Among them, EGFR, TOP1 and VEGFA are known targets of FDA-approved drugs. Additionally, CCND1 (cyclin D1), and PTGS2 (prostaglandin-endoperoxide synthase 2) have reported to be relevant to CRC or as potential drug targets based on the literature search. These results indicate that our approach is promising for drug target prediction for CRC treatment, which might be useful for other cancer therapeutics.

Systems Pharmacology Links GPCRs with Retinal Degenerative Disorders

In most biological systems, second messengers and their key regulatory and effector proteins form links between multiple cellular signaling pathways. Such signaling nodes can integrate the deleterious effects of genetic aberrations, environmental stressors, or both in complex diseases, leading to cell death by various mechanisms. Here we present a systems (network) pharmacology approach that, together with transcriptomics analyses, was used to identify different G protein-coupled receptors that experimentally protected against cellular stress and death caused by linked signaling mechanisms. We describe the application of this concept to degenerative and diabetic retinopathies in appropriate mouse models as an example. Systems pharmacology also provides an attractive framework for devising strategies to combat complex diseases by using (repurposing) US Food and Drug Administration-approved pharmacological agents.

Systems biology approaches for advancing the discovery of effective drug combinations

Complex diseases like cancer are regulated by large, interconnected networks with many pathways affecting cell proliferation, invasion, and drug resistance. However, current cancer therapy predominantly relies on the reductionist approach of one gene-one disease. Combinations of drugs may overcome drug resistance by limiting mutations and induction of escape pathways, but given the enormous number of possible drug combinations, strategies to reduce the search space and prioritize experiments are needed. In this review, we focus on the use of computational modeling, bioinformatics and high-throughput experimental methods for discovery of drug combinations. We highlight cutting-edge systems approaches, including large-scale modeling of cell signaling networks, network motif analysis, statistical association-based models, identifying correlations in gene signatures, functional genomics, and high-throughput combination screens. We also present a list of publicly available data and resources to aid in discovery of drug combinations. Integration of these systems approaches will enable faster discovery and translation of clinically relevant drug combinations.

Graphical abstract

Systems pharmacology dissection of the anti-inflammatory mechanism for the medicinal herb Folium eriobotryae

Inflammation is a hallmark of many diseases like diabetes, cancers, atherosclerosis and arthritis. Thus, lots of concerns have been raised toward developing novel anti-inflammatory agents. Many alternative herbal medicines possess excellent anti-inflammatory properties, yet their precise mechanisms of action are yet to be elucidated. Here, a novel systems pharmacology approach based on a large number of chemical, biological and pharmacological data was developed and exemplified by a probe herb Folium Eriobotryae, a widely used clinical anti-inflammatory botanic drug. The results show that 11 ingredients of this herb with favorable pharmacokinetic properties are predicted as active compounds for anti-inflammatory treatment. In addition, via systematic network analyses, their targets are identified to be 43 inflammation-associated proteins including especially COX2, ALOX5, PPARG, TNF and RELA that are mainly involved in the mitogen-activated protein kinase (MAPK) signaling pathway, the rheumatoid arthritis pathway and NF-κB signaling pathway. All these demonstrate that the integrated systems pharmacology method provides not only an effective tool to illustrate the anti-inflammatory mechanisms of herbs, but also a new systems-based approach for drug discovery from, but not limited to, herbs, especially when combined with further experimental validations.

Systems pharmacology augments drug safety surveillance

Small molecule drugs are the foundation of modern medical practice yet their use is limited by the onset of unexpected and severe adverse events (AEs). Regulatory agencies rely on post-marketing surveillance to monitor safety once drugs are approved for clinical use. Despite advances in pharmacovigilance methods that address issues of confounding bias, clinical data of AEs are inherently noisy. Systems pharmacology– the integration of systems biology and chemical genomics – can illuminate drug mechanisms of action. We hypothesize that these data can improve drug safety surveillance by highlighting drugs with a mechanistic connection to the target phenotype (enriching true positives) and filtering those that do not (depleting false positives). We present an algorithm, the modular assembly of drug safety subnetworks (MADSS), to combine systems pharmacology and pharmacovigilance data and significantly improve drug safety monitoring for four clinically relevant adverse drug reactions.

Harnessing pharmacological knowledge for personalized medicine and pharmacotyping: Challenges and lessons learned

The contribution of the genetic make-up to an individual’s capacity has long been recognized in modern pharmacology as a crucial factor leading to therapy inefficiency and toxicity, negatively impacting the economic burden of healthcare and restricting the monitoring of diseases. In practical terms, and in order for drug prescription to be improved toward meeting the personalized medicine concept in drug delivery, the maximum clinical outcome for most, if not all, patients must be achieved, i.e., pharmacotyping. Such a direction although promising and of high expectation from the society, it is however hardly to be afforded for healthcare worldwide. To overcome any existed hurdles, this means that practical clinical utility of personalized medicine decisions have to be documented and validated in the clinical setting. The latter implies for drug delivery the efficient implementation of previously gained in vivo pharmacology experience with pharmacogenomics knowledge. As an approach to work faster and in a more productive way, the elaboration of advanced physiologically based pharmacokinetics models is discussed. And in better clarifying this topic, the example of tamoxifen is thoroughly presented. Overall, pharmacotyping represents a major challenge in modern therapeutics for which pharmacologists need to work in successfully fulfilling this task.

Improving postapproval drug safety surveillance: getting better information sooner

Adverse drug events (ADEs) are an important public health concern, accounting for 5% of all hospital admissions and two-thirds of all complications occurring shortly after hospital discharge. There are often long delays between when a drug is approved and when serious ADEs are identified. Recent and ongoing advances in drug safety surveillance include the establishment of government-sponsored networks of population databases, the use of data mining approaches, and the formal integration of diverse sources of drug safety information. These advances promise to reduce delays in identifying drug-related risks and in providing reassurance about the absence of such risks.

Systems biology, complexity, and the impact on antiepileptic drug discovery

The number of available anticonvulsant drugs increased in the period spanning over more than a century, amounting to the current panoply of nearly two dozen so-called antiepileptic drugs (AEDs). However, none of them actually prevents/reduces the post-brain insult development of epilepsy in man, and in no less than a third of patients with epilepsy, the seizures are not drug-controlled. Plausibly, the enduring limitation of AEDs' efficacy derives from the insufficient understanding of epileptic pathology. This review pinpoints the unbalanced reductionism of the analytic approaches that overlook the intrinsic complexity of epilepsy and of the drug resistance in epilepsy as the core conceptual flaw hampering the discovery of truly antiepileptogenic drugs. A rising awareness of the complexity of epileptic pathology is, however, brought about by the emergence of nonreductionist systems biology (SB) that considers the networks of interactions underlying the normal organismic functions and of SB-based systems (network) pharmacology that aims to restore pathological networks. By now, the systems pharmacology approaches of AED discovery are fairly meager, but their forthcoming development is both a necessity and a realistic prospect, explored in this review.

Lean Big Data integration in systems biology and systems pharmacology

Data sets from recent large-scale projects can be integrated into one unified puzzle that can provide new insights into how drugs and genetic perturbations applied to human cells are linked to whole-organism phenotypes. Data that report how drugs affect the phenotype of human cell lines and how drugs induce changes in gene and protein expression in human cell lines can be combined with knowledge about human disease, side effects induced by drugs, and mouse phenotypes. Such data integration efforts can be achieved through the conversion of data from the various resources into single-node-type networks, gene-set libraries, or multipartite graphs. This approach can lead us to the identification of more relationships between genes, drugs, and phenotypes as well as benchmark computational and experimental methods. Overall, this lean 'Big Data' integration strategy will bring us closer toward the goal of realizing personalized medicine.

Comparison of ultra-fast 2D and 3D ligand and target descriptors for side effect prediction and network analysis in polypharmacology

BACKGROUND AND PURPOSE

Some existing computational methods are used to infer protein targets of small molecules and can therefore be used to find new targets for existing drugs, with the goals of re-directing the molecule towards a different therapeutic purpose or explaining off-target effects due to multiple targeting. Inherent limitations, however, arise from the fact that chemical analogy is calculated on the basis of common frameworks or scaffolds and also because target information is neglected. The method we present addresses these issues by taking into account 3D information from both the ligand and the target.

EXPERIMENTAL APPROACH

ElectroShape is an established method for ultra-fast comparison of the shapes and charge distributions of ligands that is validated here for prediction of on-target activities, off-target profiles and adverse effects of drugs and drug-like molecules taken from the DrugBank database.

KEY RESULTS

The method is shown to predict polypharmacology profiles and relate targets from two complementary viewpoints (ligand- and target-based networks).

CONCLUSIONS AND IMPLICATIONS

The open-access web tool presented here (http://ub.cbm.uam.es/chemogenomics/) allows interactive navigation in a unified 'pharmacological space' from the viewpoints of both ligands and targets. It also enables prediction of pharmacological profiles, including likely side effects, for new compounds. We hope this web interface will help many pharmacologists to become aware of this new paradigm (up to now mostly used in the realm of the so-called 'chemical biology') and encourage its use with a view to revealing 'hidden' relationships between new and existing compounds and pharmacologically relevant targets.

Systems pharmacology of adverse event mitigation by drug combinations

Drugs are designed for therapy, but medication-related adverse events are common, and risk/benefit analysis is critical for determining clinical use. Rosiglitazone, an efficacious antidiabetic drug, is associated with increased myocardial infarctions (MIs), thus limiting its usage. Because diabetic patients are often prescribed multiple drugs, we searched for usage of a second drug ("drug B") in the Food and Drug Administration's Adverse Event Reporting System (FAERS) that could mitigate the risk of rosiglitazone ("drug A")-associated MI. In FAERS, rosiglitazone usage is associated with increased occurrence of MI, but its combination with exenatide significantly reduces rosiglitazone-associated MI. Clinical data from the Mount Sinai Data Warehouse support the observations from FAERS. Analysis for confounding factors using logistic regression showed that they were not responsible for the observed effect. Using cell biological networks, we predicted that the mitigating effect of exenatide on rosiglitazone-associated MI could occur through clotting regulation. Data we obtained from the db/db mouse model agreed with the network prediction. To determine whether polypharmacology could generally be a basis for adverse event mitigation, we analyzed the FAERS database for other drug combinations wherein drug B reduced serious adverse events reported with drug A usage such as anaphylactic shock and suicidality. This analysis revealed 19,133 combinations that could be further studied. We conclude that this type of crowdsourced approach of using databases like FAERS can help to identify drugs that could potentially be repurposed for mitigation of serious adverse events.

A scalable computational framework for establishing long-term behavior of stochastic reaction networks

Reaction networks are systems in which the populations of a finite number of species evolve through predefined interactions. Such networks are found as modeling tools in many biological disciplines such as biochemistry, ecology, epidemiology, immunology, systems biology and synthetic biology. It is now well-established that, for small population sizes, stochastic models for biochemical reaction networks are necessary to capture randomness in the interactions. The tools for analyzing such models, however, still lag far behind their deterministic counterparts. In this paper, we bridge this gap by developing a constructive framework for examining the long-term behavior and stability properties of the reaction dynamics in a stochastic setting. In particular, we address the problems of determining ergodicity of the reaction dynamics, which is analogous to having a globally attracting fixed point for deterministic dynamics. We also examine when the statistical moments of the underlying process remain bounded with time and when they converge to their steady state values. The framework we develop relies on a blend of ideas from probability theory, linear algebra and optimization theory. We demonstrate that the stability properties of a wide class of biological networks can be assessed from our sufficient theoretical conditions that can be recast as efficient and scalable linear programs, well-known for their tractability. It is notably shown that the computational complexity is often linear in the number of species. We illustrate the validity, the efficiency and the wide applicability of our results on several reaction networks arising in biochemistry, systems biology, epidemiology and ecology. The biological implications of the results as well as an example of a non-ergodic biological network are also discussed.

In many biological disciplines, computational modeling of interaction networks is the key for understanding biological phenomena. Such networks are traditionally studied using deterministic models. However, it has been recently recognized that when the populations are small in size, the inherent random effects become significant and to incorporate them, a stochastic modeling paradigm is necessary. Hence, stochastic models of reaction networks have been broadly adopted and extensively used. Such models, for instance, form a cornerstone for studying heterogeneity in clonal cell populations. In biological applications, one is often interested in knowing the long-term behavior and stability properties of reaction networks even with incomplete knowledge of the model parameters. However for stochastic models, no analytical tools are known for this purpose, forcing many researchers to use a simulation-based approach, which is highly unsatisfactory. To address this issue, we develop a theoretical and computational framework for determining the long-term behavior and stability properties for stochastic reaction networks. Our approach is based on a mixture of ideas from probability theory, linear algebra and optimization theory. We illustrate the broad applicability of our results by considering examples from various biological areas. The biological implications of our results are discussed as well.

Towards structural systems pharmacology to study complex diseases and personalized medicine

Genome-Wide Association Studies (GWAS), whole genome sequencing, and high-throughput omics techniques have generated vast amounts of genotypic and molecular phenotypic data. However, these data have not yet been fully explored to improve the effectiveness and efficiency of drug discovery, which continues along a one-drug-one-target-one-disease paradigm. As a partial consequence, both the cost to launch a new drug and the attrition rate are increasing. Systems pharmacology and pharmacogenomics are emerging to exploit the available data and potentially reverse this trend, but, as we argue here, more is needed. To understand the impact of genetic, epigenetic, and environmental factors on drug action, we must study the structural energetics and dynamics of molecular interactions in the context of the whole human genome and interactome. Such an approach requires an integrative modeling framework for drug action that leverages advances in data-driven statistical modeling and mechanism-based multiscale modeling and transforms heterogeneous data from GWAS, high-throughput sequencing, structural genomics, functional genomics, and chemical genomics into unified knowledge. This is not a small task, but, as reviewed here, progress is being made towards the final goal of personalized medicines for the treatment of complex diseases.

A network pharmacology study of Chinese medicine QiShenYiQi to reveal its underlying multi-compound, multi-target, multi-pathway mode of action

Chinese medicine is a complex system guided by traditional Chinese medicine (TCM) theories, which has proven to be especially effective in treating chronic and complex diseases. However, the underlying modes of action (MOA) are not always systematically investigated. Herein, a systematic study was designed to elucidate the multi-compound, multi-target and multi-pathway MOA of a Chinese medicine, QiShenYiQi (QSYQ), on myocardial infarction. QSYQ is composed of Astragalus membranaceus (Huangqi), Salvia miltiorrhiza (Danshen), Panax notoginseng (Sanqi), and Dalbergia odorifera (Jiangxiang). Male Sprague Dawley rat model of myocardial infarction were administered QSYQ intragastrically for 7 days while the control group was not treated. The differentially expressed genes (DEGs) were identified from myocardial infarction rat model treated with QSYQ, followed by constructing a cardiovascular disease (CVD)-related multilevel compound-target-pathway network connecting main compounds to those DEGs supported by literature evidences and the pathways that are functionally enriched in ArrayTrack. 55 potential targets of QSYQ were identified, of which 14 were confirmed in CVD-related literatures with experimental supporting evidences. Furthermore, three sesquiterpene components of QSYQ, Trans-nerolidol, (3S,6S,7R)-3,7,11-trimethyl-3,6-epoxy-1,10-dodecadien-7-ol and (3S,6R,7R)-3,7,11-trimethyl-3,6-epoxy-1,10-dodecadien-7-ol from Dalbergia odorifera T. Chen, were validated experimentally in this study. Their anti-inflammatory effects and potential targets including extracellular signal-regulated kinase-1/2, peroxisome proliferator-activated receptor-gamma and heme oxygenase-1 were identified. Finally, through a three-level compound-target-pathway network with experimental analysis, our study depicts a complex MOA of QSYQ on myocardial infarction.

Rectifying cancer drug discovery through network pharmacology

In spite of the expensive preclinical testing, the consistent failure to translate many promising targeted drugs from the laboratory bench to the clinic raises the question of whether the single-pathway drug-discovery strategies offer the correct perspective. As revealed by network biology, cancers harbor robust biological networks that are inherently resistant to changes, such as those induced by drugs with very narrow mechanisms of action. Therefore, network pharmacology strategies, the treatment of cancer by modulating more than one target, are needed. Different promiscuous approaches targeting multiple avenues within cancer-associated networks, such as the pleiotropic natural products, are emerging. Nevertheless, there is a long way before such 'proof-of-concept strategies' can be successfully applied in the clinical setting. This article provides a perspective on the current challenges in drug discovery, the reasons for high failure rates and how network pharmacology can aid the successful design of agents against cancer.

Roles of raft-anchored adaptor Cbp/PAG1 in spatial regulation of c-Src kinase

The tyrosine kinase c-Src is upregulated in numerous human cancers, implying a role for c-Src in cancer progression. Previously, we have shown that sequestration of activated c-Src into lipid rafts via a transmembrane adaptor, Cbp/PAG1, efficiently suppresses c-Src-induced cell transformation in Csk-deficient cells, suggesting that the transforming activity of c-Src is spatially regulated via Cbp in lipid rafts. To dissect the molecular mechanisms of the Cbp-mediated regulation of c-Src, a combined analysis was performed that included mathematical modeling and in vitro experiments in a c-Src- or Cbp-inducible system. c-Src activity was first determined as a function of c-Src or Cbp levels, using focal adhesion kinase (FAK) as a crucial c-Src substrate. Based on these experimental data, two mathematical models were constructed, the sequestration model and the ternary model. The computational analysis showed that both models supported our proposal that raft localization of Cbp is crucial for the suppression of c-Src function, but the ternary model, which includes a ternary complex consisting of Cbp, c-Src, and FAK, also predicted that c-Src function is dependent on the lipid-raft volume. Experimental analysis revealed that c-Src activity is elevated when lipid rafts are disrupted and the ternary complex forms in non-raft membranes, indicating that the ternary model accurately represents the system. Moreover, the ternary model predicted that, if Cbp enhances the interaction between c-Src and FAK, Cbp could promote c-Src function when lipid rafts are disrupted. These findings underscore the crucial role of lipid rafts in the Cbp-mediated negative regulation of c-Src-transforming activity, and explain the positive role of Cbp in c-Src regulation under particular conditions where lipid rafts are perturbed.

A network pharmacology approach to understanding the mechanisms of action of traditional medicine: Bushenhuoxue formula for treatment of chronic kidney disease

Traditional Chinese medicine (TCM) has unique therapeutic effects for complex chronic diseases. However, for the lack of an effective systematic approach, the research progress on the effective substances and pharmacological mechanism of action has been very slow. In this paper, by incorporating network biology, bioinformatics and chemoinformatics methods, an integrated approach was proposed to systematically investigate and explain the pharmacological mechanism of action and effective substances of TCM. This approach includes the following main steps: First, based on the known drug targets, network biology was used to screen out putative drug targets; Second, the molecular docking method was used to calculate whether the molecules from TCM and drug targets related to chronic kidney diseases (CKD) interact or not; Third, according to the result of molecular docking, natural product-target network, main component-target network and compound-target network were constructed; Finally, through analysis of network characteristics and literature mining, potential effective multi-components and their synergistic mechanism were putatively identified and uncovered. Bu-shen-Huo-xue formula (BSHX) which was frequently used for treating CKD, was used as the case to demonstrate reliability of our proposed approach. The results show that BSHX has the therapeutic effect by using multi-channel network regulation, such as regulating the coagulation and fibrinolytic balance, and the expression of inflammatory factors, inhibiting abnormal ECM accumulation. Tanshinone IIA, rhein, curcumin, calycosin and quercetin may be potential effective ingredients of BSHX. This research shows that the integration approach can be an effective means for discovering active substances and revealing their pharmacological mechanisms of TCM.

Metabolomics and systems pharmacology: why and how to model the human metabolic network for drug discovery

Metabolism represents the 'sharp end' of systems biology, because changes in metabolite concentrations are necessarily amplified relative to changes in the transcriptome, proteome and enzyme activities, which can be modulated by drugs. To understand such behaviour, we therefore need (and increasingly have) reliable consensus (community) models of the human metabolic network that include the important transporters. Small molecule 'drug' transporters are in fact metabolite transporters, because drugs bear structural similarities to metabolites known from the network reconstructions and from measurements of the metabolome. Recon2 represents the present state-of-the-art human metabolic network reconstruction; it can predict inter alia: (i) the effects of inborn errors of metabolism; (ii) which metabolites are exometabolites, and (iii) how metabolism varies between tissues and cellular compartments. However, even these qualitative network models are not yet complete. As our understanding improves so do we recognise more clearly the need for a systems (poly)pharmacology.

Functional Module Connectivity Map (FMCM): a framework for searching repurposed drug compounds for systems treatment of cancer and an application to colorectal adenocarcinoma

Drug repurposing has become an increasingly attractive approach to drug development owing to the ever-growing cost of new drug discovery and frequent withdrawal of successful drugs caused by side effect issues. Here, we devised Functional Module Connectivity Map (FMCM) for the discovery of repurposed drug compounds for systems treatment of complex diseases, and applied it to colorectal adenocarcinoma. FMCM used multiple functional gene modules to query the Connectivity Map (CMap). The functional modules were built around hub genes identified, through a gene selection by trend-of-disease-progression (GSToP) procedure, from condition-specific gene-gene interaction networks constructed from sets of cohort gene expression microarrays. The candidate drug compounds were restricted to drugs exhibiting predicted minimal intracellular harmful side effects. We tested FMCM against the common practice of selecting drugs using a genomic signature represented by a single set of individual genes to query CMap (IGCM), and found FMCM to have higher robustness, accuracy, specificity, and reproducibility in identifying known anti-cancer agents. Among the 46 drug candidates selected by FMCM for colorectal adenocarcinoma treatment, 65% had literature support for association with anti-cancer activities, and 60% of the drugs predicted to have harmful effects on cancer had been reported to be associated with carcinogens/immune suppressors. Compounds were formed from the selected drug candidates where in each compound the component drugs collectively were beneficial to all the functional modules while no single component drug was harmful to any of the modules. In cell viability tests, we identified four candidate drugs: GW-8510, etacrynic acid, ginkgolide A, and 6-azathymine, as having high inhibitory activities against cancer cells. Through microarray experiments we confirmed the novel functional links predicted for three candidate drugs: phenoxybenzamine (broad effects), GW-8510 (cell cycle), and imipenem (immune system). We believe FMCM can be usefully applied to repurposed drug discovery for systems treatment of other types of cancer and other complex diseases.

Integrative approaches for predicting in vivo effects of chemicals from their structural descriptors and the results of short-term biological assays

Cheminformatics approaches such as Quantitative Structure Activity Relationship (QSAR) modeling have been used traditionally for predicting chemical toxicity. In recent years, high throughput biological assays have been increasingly employed to elucidate mechanisms of chemical toxicity and predict toxic effects of chemicals in vivo. The data generated in such assays can be considered as biological descriptors of chemicals that can be combined with molecular descriptors and employed in QSAR modeling to improve the accuracy of toxicity prediction. In this review, we discuss several approaches for integrating chemical and biological data for predicting biological effects of chemicals in vivo and compare their performance across several data sets. We conclude that while no method consistently shows superior performance, the integrative approaches rank consistently among the best yet offer enriched interpretation of models over those built with either chemical or biological data alone. We discuss the outlook for such interdisciplinary methods and offer recommendations to further improve the accuracy and interpretability of computational models that predict chemical toxicity.

Network pharmacology strategies toward multi-target anticancer therapies: from computational models to experimental design principles

Polypharmacology has emerged as novel means in drug discovery for improving treatment response in clinical use. However, to really capitalize on the polypharmacological effects of drugs, there is a critical need to better model and understand how the complex interactions between drugs and their cellular targets contribute to drug efficacy and possible side effects. Network graphs provide a convenient modeling framework for dealing with the fact that most drugs act on cellular systems through targeting multiple proteins both through on-target and off-target binding. Network pharmacology models aim at addressing questions such as how and where in the disease network should one target to inhibit disease phenotypes, such as cancer growth, ideally leading to therapies that are less vulnerable to drug resistance and side effects by means of attacking the disease network at the systems level through synergistic and synthetic lethal interactions. Since the exponentially increasing number of potential drug target combinations makes pure experimental approach quickly unfeasible, this review depicts a number of computational models and algorithms that can effectively reduce the search space for determining the most promising combinations for experimental evaluation. Such computational-experimental strategies are geared toward realizing the full potential of multi-target treatments in different disease phenotypes. Our specific focus is on system-level network approaches to polypharmacology designs in anticancer drug discovery, where we give representative examples of how network-centric modeling may offer systematic strategies toward better understanding and even predicting the phenotypic responses to multi-target therapies.

CVDHD: a cardiovascular disease herbal database for drug discovery and network pharmacology

Background

Cardiovascular disease (CVD) is the leading cause of death and associates with multiple risk factors. Herb medicines have been used to treat CVD long ago in china and several natural products or derivatives (e.g., aspirin and reserpine) are most common drugs all over the world. The objective of this work was to construct a systematic database for drug discovery based on natural products separated from CVD-related medicinal herbs and to research on action mechanism of herb medicines.

Description

The cardiovascular disease herbal database (CVDHD) was designed to be a comprehensive resource for virtual screening and drug discovery from natural products isolated from medicinal herbs for cardiovascular-related diseases. CVDHD comprises 35230 distinct molecules and their identification information (chemical name, CAS registry number, molecular formula, molecular weight, international chemical identifier (InChI) and SMILES), calculated molecular properties (AlogP, number of hydrogen bond acceptor and donors, etc.), docking results between all molecules and 2395 target proteins, cardiovascular-related diseases, pathways and clinical biomarkers. All 3D structures were optimized in the MMFF94 force field and can be freely accessed.

Conclusions

CVDHD integrated medicinal herbs, natural products, CVD-related target proteins, docking results, diseases and clinical biomarkers. By using the methods of virtual screening and network pharmacology, CVDHD will provide a platform to streamline drug/lead discovery from natural products and explore the action mechanism of medicinal herbs. CVDHD is freely available at http://pkuxxj.pku.edu.cn/CVDHD.

Predicting drug-target interactions using probabilistic matrix factorization

Quantitative analysis of known drug-target interactions emerged in recent years as a useful approach for drug repurposing and assessing side effects. In the present study, we present a method that uses probabilistic matrix factorization (PMF) for this purpose, which is particularly useful for analyzing large interaction networks. DrugBank drugs clustered based on PMF latent variables show phenotypic similarity even in the absence of 3D shape similarity. Benchmarking computations show that the method outperforms those recently introduced provided that the input data set of known interactions is sufficiently large--which is the case for enzymes and ion channels, but not for G-protein coupled receptors (GPCRs) and nuclear receptors. Runs performed on DrugBank after hiding 70% of known interactions show that, on average, 88 of the top 100 predictions hit the hidden interactions. De novo predictions permit us to identify new potential interactions. Drug-target pairs implicated in neurobiological disorders are overrepresented among de novo predictions.

A network inference method for large-scale unsupervised identification of novel drug-drug interactions

Characterizing interactions between drugs is important to avoid potentially harmful combinations, to reduce off-target effects of treatments and to fight antibiotic resistant pathogens, among others. Here we present a network inference algorithm to predict uncharacterized drug-drug interactions. Our algorithm takes, as its only input, sets of previously reported interactions, and does not require any pharmacological or biochemical information about the drugs, their targets or their mechanisms of action. Because the models we use are abstract, our approach can deal with adverse interactions, synergistic/antagonistic/suppressing interactions, or any other type of drug interaction. We show that our method is able to accurately predict interactions, both in exhaustive pairwise interaction data between small sets of drugs, and in large-scale databases. We also demonstrate that our algorithm can be used efficiently to discover interactions of new drugs as part of the drug discovery process.

Tamoxifen integromics and personalized medicine: dynamic modular transformations underpinning response to tamoxifen in breast cancer treatment

Recent advances in pharmacogenomics technologies allow bold steps to be taken towards personalized medicine, more accurate health planning, and personalized drug development. In this framework, systems pharmacology network-based approaches offer an appealing way for integrating multi-omics data and set the basis for defining systems-level drug response biomarkers. On the road to individualized tamoxifen treatment in estrogen receptor-positive breast cancer patients, we examine the dynamics of the attendant pharmacological response mechanisms. By means of an "integromics" network approach, we assessed the tamoxifen effect through the way the high-order organization of interactome (i.e., the modules) is perturbed. To accomplish that, first we integrated the time series transcriptome data with the human protein interaction data, and second, an efficient module-detecting algorithm was applied onto the composite graphs. Our findings show that tamoxifen induces severe modular transformations on specific areas of the interactome. Our modular biomarkers in response to tamoxifen attest to the immunomodulatory role of tamoxifen, and further reveal that it deregulates cell cycle and apoptosis pathways, while coordinating the proteasome and basal transcription factors. To the best of our knowledge, this is the first report that informs the fields of personalized medicine and clinical pharmacology about the actual dynamic interactome response to tamoxifen administration.

Deciphering the combination principles of Traditional Chinese Medicine from a systems pharmacology perspective based on Ma-huang Decoction

ETHNOPHARMACOLOGY RELEVANCE

The main therapeutic concept in Traditional Chinese Medicine (TCM) is herb formula, which treats various diseases via potential herb interactions to maximize the efficacy and minimize the adverse effects. However, the combination principle of herb formula still remains a mystery due to the lack of appropriate methods.

METHODS

A systems pharmacology method integrating the pharmacokinetic analysis, drug targeting, and drug-target-disease network is developed to dissect this rule embedded in the herbal formula. All these are exemplified by a representative TCM formula, Ma-huang decoction, made up of four botanic herbs.

RESULTS

Based on the deep investigation of the function and compatibility of each herb, in a molecular/systems level, we demonstrate the different pharmacological roles that each herb might play in the prescription. By the way of enhancing the bioavailability and/or making the pharmacological synergy among different herbs, the four herbs effectively combine together to be suitable for treating diseases.

CONCLUSIONS

The present work lays foundations for a more comprehensive understanding of the combination rule of TCM, which might also be beneficial to drug development and applications.