Search

Recent reports have called into question the reproducibility, validity and translatability of the preclinical animal studies due to limitations in their experimental design and statistical analysis. To this end, we implemented a matching-based modelling approach for optimal intervention group allocation, randomization and power calculations, which takes full account of the complex animal characteristics at baseline prior to interventions. In prostate cancer xenograft studies, the method effectively normalized the confounding baseline variability, and resulted in animal allocations which were supported by RNA-seq profiling of the individual tumours. The matching information increased the statistical power to detect true treatment effects at smaller sample sizes in two castration-resistant prostate cancer models, thereby leading to saving of both animal lives and research costs. The novel modelling approach and its open-source and web-based software implementations enable the researchers to conduct adequately-powered and fully-blinded preclinical intervention studies, with the aim to accelerate the discovery of new therapeutic interventions.

Background: The estrogen receptor (ER) inhibitor tamoxifen reduces breast cancer mortality by 31 % and has served as the standard treatment for ER-positive breast cancers for decades. However, 50 % of advanced ER-positive cancers display de novo resistance to tamoxifen, and acquired resistance evolves in 40 % of patients who initially respond. Mechanisms underlying resistance development remain poorly understood and new therapeutic opportunities are urgently needed. Here, we report the generation and characterization of seven tamoxifen-resistant breast cancer cell lines from four parental strains., Methods: Using high throughput drug sensitivity and resistance testing (DSRT) with 279 approved and investigational oncology drugs, exome-sequencing and network analysis, we for the first time, systematically determine the drug response profiles specific to tamoxifen resistance., Results: We discovered emerging vulnerabilities towards specific drugs, such as ERK1/2-, proteasome- and BCL-family inhibitors as the cells became tamoxifen-resistant. Co-resistance to other drugs such as the survivin inhibitor YM155 and the chemotherapeutic agent paclitaxel also occurred., Conclusion: This study indicates that multiple molecular mechanisms dictate endocrine resistance, resulting in unexpected vulnerabilities to initially ineffective drugs, as well as in emerging co-resistances. Thus, combatting drug-resistant tumors will require patient-tailored strategies in order to identify new drug vulnerabilities, and to understand the associated co-resistance patterns.

Drug discovery is moving away from the single target-based approach towards harnessing the potential of polypharmacological agents that modulate the activity of multiple nodes in the complex networks of deregulations underlying disease phenotypes. Computational network pharmacology methods that use systems-level drug-response phenotypes, such as those originating from genome-wide transcriptomic profiles, have proved particularly effective for elucidating the mechanisms of action of multitargeted compounds. Here, we show, via the case study of the natural product pinosylvin, how the combination of two complementary network-based methods can provide novel, unexpected mechanistic insights. This case study also illustrates that elucidating the mechanism of action of multitargeted natural products through transcriptional response-based approaches is a challenging endeavor, often requiring multiple computational-experimental iterations., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Background: Systematic approaches to functionally identify key players in microRNA (miRNA)-target networks regulating prostate cancer (PCa) proliferation are still missing., Objective: To comprehensively map miRNA regulation of genes relevant for PCa proliferation through phenotypic screening and tumor expression data., Design, Setting, and Participants: Gain-of-function screening with 1129 miRNA molecules was performed in five PCa cell lines, measuring proliferation, viability, and apoptosis. These results were integrated with changes in miRNA expression from two cohorts of human PCa (188 tumors in total). For resulting miRNAs, the predicted targets were collected and analyzed for patterns with gene set enrichment analysis, and for their association with biochemical recurrence free survival., Outcome Measurements and Statistical Analysis: Rank product statistical analysis was used to evaluate miRNA effects in phenotypic screening and for expression differences in the prostate tumor cohorts. Expression data were analyzed using the significance analysis of microarrays (SAM) method and the patient material was subjected to Kaplan-Meier statistics., Results and Limitations: Functional screening identified 25 miRNAs increasing and 48 miRNAs decreasing cell viability. Data integration resulted in 14 miRNAs, with aberrant expression and effect on proliferation. These miRNAs are predicted to regulate >3700 genes, of which 28 were found up-regulated and 127 down-regulated in PCa compared with benign tissue. Seven genes, FLNC, MSRB3, PARVA, PCDH7, PRNP, RAB34, and SORBS1, showed an inverse association to their predicted miRNA, and were identified to significantly correlate with biochemical recurrence free survival in PCa patients., Conclusions: A systematic in vitro screening approach combined with in vivo expression and gene set enrichment analysis provide unbiased means for revealing novel miRNA-target links, possibly driving the oncogenic processes in PCa., Patient Summary: This study identified novel regulatory molecules, which impact on PCa proliferation and are aberrantly expressed in clinical tumors. Thus, our study reveals regulatory nodes with potential for therapy., (Copyright © 2015 European Association of Urology. Published by Elsevier B.V. All rights reserved.)

Background: Triple negative breast cancer (TNBC) is a highly heterogeneous and aggressive type of cancer that lacks effective targeted therapy. Despite detailed molecular profiling, no targeted therapy has been established. Hence, with the aim of gaining deeper understanding of the functional differences of TNBC subtypes and how that may relate to potential novel therapeutic strategies, we studied comprehensive anticancer-agent responses among a panel of TNBC cell lines., Method: The responses of 301 approved and investigational oncology compounds were measured in 16 TNBC cell lines applying a functional profiling approach. To go beyond the standard drug viability effect profiling, which has been used in most chemosensitivity studies, we utilized a multiplexed readout for both cell viability and cytotoxicity, allowing us to differentiate between cytostatic and cytotoxic responses., Results: Our approach revealed that most single-agent anti-cancer compounds that showed activity for the viability readout had no or little cytotoxic effects. Major compound classes that exhibited this type of response included anti-mitotics, mTOR, CDK, and metabolic inhibitors, as well as many agents selectively inhibiting oncogene-activated pathways. However, within the broad viability-acting classes of compounds, there were often subsets of cell lines that responded by cell death, suggesting that these cells are particularly vulnerable to the tested substance. In those cases we could identify differential levels of protein markers associated with cytotoxic responses. For example, PAI-1, MAPK phosphatase and Notch-3 levels associated with cytotoxic responses to mitotic and proteasome inhibitors, suggesting that these might serve as markers of response also in clinical settings. Furthermore, the cytotoxicity readout highlighted selective synergistic and synthetic lethal drug combinations that were missed by the cell viability readouts. For instance, the MEK inhibitor trametinib synergized with PARP inhibitors. Similarly, combination of two non-cytotoxic compounds, the rapamycin analog everolimus and an ATP-competitive mTOR inhibitor dactolisib, showed synthetic lethality in several mTOR-addicted cell lines., Conclusions: Taken together, by studying the combination of cytotoxic and cytostatic drug responses, we identified a deeper spectrum of cellular responses both to single agents and combinations that may be highly relevant for identifying precision medicine approaches in TNBC as well as in other types of cancers.

The oncoprotein Cancerous Inhibitor of Protein Phosphatase 2A (CIP2A) is overexpressed in most malignancies and is an obvious candidate target protein for future cancer therapies. However, the physiological importance of CIP2A-mediated PP2A inhibition is largely unknown. As PP2A regulates immune responses, we investigated the role of CIP2A in normal immune system development and during immune response in vivo. We show that CIP2A-deficient mice (CIP2AHOZ) present a normal immune system development and function in unchallenged conditions. However when challenged with Listeria monocytogenes, CIP2AHOZ mice display an impaired adaptive immune response that is combined with decreased frequency of both CD4+ T-cells and CD8+ effector T-cells. Importantly, the cell autonomous effect of CIP2A deficiency for T-cell activation was confirmed. Induction of CIP2A expression during T-cell activation was dependent on Zap70 activity. Thus, we reveal CIP2A as a hitherto unrecognized mediator of T-cell activation during adaptive immune response. These results also reveal CIP2AHOZ as a possible novel mouse model for studying the role of PP2A activity in immune regulation. On the other hand, the results also indicate that CIP2A targeting cancer therapies would not cause serious immunological side-effects.

Ovarian cancer is a heterogeneous disease of low prevalence, but poor survival. Early diagnosis is critical for survival, but it is often challenging because the symptoms of ovarian cancer are subtle and become apparent only during advanced stages of the disease. Therefore, the identification of robust biomarkers of early disease is a clinical priority. Metabolomic profiling is an emerging diagnostic tool enabling the detection of biomarkers reflecting alterations in tumor metabolism, a hallmark of cancer. In this study, we performed metabolomic profiling of serum and tumor tissue from 158 patients with high-grade serous ovarian cancer (HGSOC) and 100 control patients with benign or non-neoplastic lesions. We report metabolites of hydroxybutyric acid (HBA) as novel diagnostic and prognostic biomarkers associated with tumor burden and patient survival. The accumulation of HBA metabolites caused by HGSOC was also associated with reduced expression of succinic semialdehyde dehydrogenase (encoded by ALDH5A1), and with the presence of an epithelial-to-mesenchymal transition gene signature, implying a role for these metabolic alterations in cancer cell migration and invasion. In conclusion, our findings represent the first comprehensive metabolomics analysis in HGSOC and propose a new set of metabolites as biomarkers of disease with diagnostic and prognostic capabilities., (©2015 American Association for Cancer Research.)

Estrogens are suggested to lower the risk of developing metabolic syndrome in both sexes. In this study, we investigated how the increased circulating estrogen-to-androgen ratio (E/A) alters liver lipid metabolism in males. The cytochrome P450 aromatase (P450arom) is an enzyme converting androgens to estrogens. Male mice overexpressing human aromatase enzyme (AROM+ mice), and thus have high circulating E/A, were used as a model in this study. Proteomics and gene expression analyses indicated an increase in the peroxisomal β-oxidation in the liver of AROM+ mice as compared with their wild type littermates. Correspondingly, metabolomic analysis revealed a decrease in the amount of phosphatidylcholines with long-chain fatty acids in the plasma. With interest we noted that the expression of Cyp4a12a enzyme, which specifically metabolizes arachidonic acid (AA) to 20-hydroxy AA, was dramatically decreased in the AROM+ liver. As a consequence, increased amounts of phospholipids having AA as a fatty acid tail were detected in the plasma of the AROM+ mice. Overall, these observations demonstrate that high circulating E/A in males is linked to indicators of higher peroxisomal β-oxidation and lower AA metabolism in the liver. Furthermore, the plasma phospholipid profile reflects the changes in the liver lipid metabolism., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Introduction: System-wide identification of both on- and off-targets of chemical probes provides improved understanding of their therapeutic potential and possible adverse effects, thereby accelerating and de-risking drug discovery process. Given the high costs of experimental profiling of the complete target space of drug-like compounds, computational models offer systematic means for guiding these mapping efforts. These models suggest the most potent interactions for further experimental or pre-clinical evaluation both in cell line models and in patient-derived material., Areas Covered: The authors focus here on network-based machine learning models and their use in the prediction of novel compound-target interactions both in target-based and phenotype-based drug discovery applications. While currently being used mainly in complementing the experimentally mapped compound-target networks for drug repurposing applications, such as extending the target space of already approved drugs, these network pharmacology approaches may also suggest completely unexpected and novel investigational probes for drug development., Expert Opinion: Although the studies reviewed here have already demonstrated that network-centric modeling approaches have the potential to identify candidate compounds and selective targets in disease networks, many challenges still remain. In particular, these challenges include how to incorporate the cellular context and genetic background into the disease networks to enable more stratified and selective target predictions, as well as how to make the prediction models more realistic for the practical drug discovery and therapeutic applications.

High content protein interaction screens have revolutionized our understanding of protein complex assembly. However, one of the major challenges in translation of high content protein interaction data is identification of those interactions that are functionally relevant for a particular biological question. To address this challenge, we developed a relevance ranking platform (RRP), which consist of modular functional and bioinformatic filters to provide relevance rank among the interactome proteins. We demonstrate the versatility of RRP to enable a systematic prioritization of the most relevant interaction partners from high content data, highlighted by the analysis of cancer relevant protein interactions for oncoproteins Pin1 and PME-1. We validated the importance of selected interactions by demonstration of PTOV1 and CSKN2B as novel regulators of Pin1 target c-Jun phosphorylation and reveal previously unknown interacting proteins that may mediate PME-1 effects via PP2A-inhibition. The RRP framework is modular and can be modified to answer versatile research problems depending on the nature of the biological question under study. Based on comparison of RRP to other existing filtering tools, the presented data indicate that RRP offers added value especially for the analysis of interacting proteins for which there is no sufficient prior knowledge available. Finally, we encourage the use of RRP in combination with either SAINT or CRAPome computational tools for selecting the candidate interactors that fulfill the both important requirements, functional relevance, and high confidence interaction detection., (© 2015 by The American Society for Biochemistry and Molecular Biology, Inc.)

Motivation: Most data analysis tools for high-throughput screening (HTS) seek to uncover interesting hits for further analysis. They typically assume a low hit rate per plate. Hit rates can be dramatically higher in secondary screening, RNAi screening and in drug sensitivity testing using biologically active drugs. In particular, drug sensitivity testing on primary cells is often based on dose-response experiments, which pose a more stringent requirement for data quality and for intra- and inter-plate variation. Here, we compared common plate normalization and noise-reduction methods, including the B-score and the Loess a local polynomial fit method under high hit-rate scenarios of drug sensitivity testing. We generated simulated 384-well plate HTS datasets, each with 71 plates having a range of 20 (5%) to 160 (42%) hits per plate, with controls placed either at the edge of the plates or in a scattered configuration., Results: We identified 20% (77/384) as the critical hit-rate after which the normalizations started to perform poorly. Results from real drug testing experiments supported this estimation. In particular, the B-score resulted in incorrect normalization of high hit-rate plates, leading to poor data quality, which could be attributed to its dependency on the median polish algorithm. We conclude that a combination of a scattered layout of controls per plate and normalization using a polynomial least squares fit method, such as Loess helps to reduce column, row and edge effects in HTS experiments with high hit-rates and is optimal for generating accurate dose-response curves., Contact: john.mpindi@helsinki.fi., Availability and Implementation: Supplementary information: R code and Supplementary data are available at Bioinformatics online., (© The Author 2015. Published by Oxford University Press.)

Deconvoluting the molecular target signals behind observed drug response phenotypes is an important part of phenotype-based drug discovery and repurposing efforts. We demonstrate here how our network-based deconvolution approach, named target addiction score (TAS), provides insights into the functional importance of druggable protein targets in cell-based drug sensitivity testing experiments. Using cancer cell line profiling data sets, we constructed a functional classification across 107 cancer cell models, based on their common and unique target addiction signatures. The pan-cancer addiction correlations could not be explained by the tissue of origin, and only correlated in part with molecular and genomic signatures of the heterogeneous cancer cells. The TAS-based cancer cell classification was also shown to be robust to drug response data resampling, as well as predictive of the transcriptomic patterns in an independent set of cancer cells that shared similar addiction signatures with the 107 cancers. The critical protein targets identified by the integrated approach were also shown to have clinically relevant mutation frequencies in patients with various cancer subtypes, including not only well-established pan-cancer genes, such as PTEN tumor suppressor, but also a number of targets that are less frequently mutated in specific cancer types, including ABL1 oncoprotein in acute myeloid leukemia. An application to leukemia patient primary cell models demonstrated how the target deconvolution approach offers functional insights into patient-specific addiction patterns, such as those indicative of their receptor-type tyrosine-protein kinase FLT3 internal tandem duplication (FLT3-ITD) status and co-addiction partners, which may lead to clinically actionable, personalized drug treatment developments. To promote its application to the future drug testing studies, we have made available an open-source implementation of the TAS calculation in the form of a stand-alone R package., (© 2015. Published by The Company of Biologists Ltd.)

Rational design of multi-targeted drug combinations is a promising strategy to tackle the drug resistance problem for many complex disorders. A drug combination is usually classified as synergistic or antagonistic, depending on the deviation of the observed combination response from the expected effect calculated based on a reference model of non-interaction. The existing reference models were proposed originally for low-throughput drug combination experiments, which make the model assumptions often incompatible with the complex drug interaction patterns across various dose pairs that are typically observed in large-scale dose-response matrix experiments. To address these limitations, we proposed a novel reference model, named zero interaction potency (ZIP), which captures the drug interaction relationships by comparing the change in the potency of the dose-response curves between individual drugs and their combinations. We utilized a delta score to quantify the deviation from the expectation of zero interaction, and proved that a delta score value of zero implies both probabilistic independence and dose additivity. Using data from a large-scale anticancer drug combination experiment, we demonstrated empirically how the ZIP scoring approach captures the experimentally confirmed drug synergy while keeping the false positive rate at a low level. Further, rather than relying on a single parameter to assess drug interaction, we proposed the use of an interaction landscape over the full dose-response matrix to identify and quantify synergistic and antagonistic dose regions. The interaction landscape offers an increased power to differentiate between various classes of drug combinations, and may therefore provide an improved means for understanding their mechanisms of action toward clinical translation.

Many biological or chemical agents when combined interact with each other and produce a synergistic response that cannot be predicted based on the single agent responses alone. However, depending on the postulated null hypothesis of non-interaction, one may end up in different interpretations of synergy. Two popular reference models for null hypothesis include the Bliss independence model and the Loewe additivity model, each of which is formulated from different perspectives. During the last century, there has been an intensive debate on the suitability of these synergy models, both of which are theoretically justified and also in practice supported by different schools of scientists. More than 20 years ago, there was a community effort to make a consensus on the terminology one should use when claiming synergy. The agreement was formulated at a conference held in Saariselkä, Finland in 1992, stating that one should use the terms Bliss synergy or Loewe synergy to avoid ambiguity in the underlying models. We review the theoretical relationships between these models and argue that one should combine the advantages of both models to provide a more consistent definition of synergy and antagonism.

Chemical perturbation screens offer the possibility to identify actionable sets of cancer-specific vulnerabilities. However, most inhibitors of kinases or other cancer targets result in polypharmacological effects, which complicate the identification of target dependencies directly from the drug-response phenotypes. In this study, we developed a chemical systems biology approach that integrates comprehensive drug sensitivity and selectivity profiling to provide functional insights into both single and multi-target oncogenic signal addictions. When applied to 21 breast cancer cell lines, perturbed with 40 kinase inhibitors, the subtype-specific addiction patterns clustered in agreement with patient-derived subtypes, while showing considerable variability between the heterogeneous breast cancers. Experimental validation of the top predictions revealed a number of co-dependencies between kinase targets that led to unexpected synergistic combinations between their inhibitors, such as dasatinib and axitinib in the triple-negative basal-like HCC1937 cell line., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Hyperactivated RAS drives progression of many human malignancies. However, oncogenic activity of RAS is dependent on simultaneous inactivation of protein phosphatase 2A (PP2A) activity. Although PP2A is known to regulate some of the RAS effector pathways, it has not been systematically assessed how these proteins functionally interact. Here we have analyzed phosphoproteomes regulated by either RAS or PP2A, by phosphopeptide enrichment followed by mass-spectrometry-based label-free quantification. To allow data normalization in situations where depletion of RAS or PP2A inhibitor CIP2A causes a large uni-directional change in the phosphopeptide abundance, we developed a novel normalization strategy, named pairwise normalization. This normalization is based on adjusting phosphopeptide abundances measured before and after the enrichment. The superior performance of the pairwise normalization was verified by various independent methods. Additionally, we demonstrate how the selected normalization method influences the downstream analyses and interpretation of pathway activities. Consequently, bioinformatics analysis of RAS and CIP2A regulated phosphoproteomes revealed a significant overlap in their functional pathways. This is most likely biologically meaningful as we observed a synergistic survival effect between CIP2A and RAS expression as well as KRAS activating mutations in TCGA pan-cancer data set, and synergistic relationship between CIP2A and KRAS depletion in colony growth assays.

Hotspot mutations of Ras drive cell transformation and tumorigenesis. Less frequent mutations in Ras are poorly characterized for their oncogenic potential. Yet insight into their mechanism of action may point to novel opportunities to target Ras. Here, we show that several cancer-associated mutations in the switch III region moderately increase Ras activity in all isoforms. Mutants are biochemically inconspicuous, while their clustering into nanoscale signaling complexes on the plasma membrane, termed nanocluster, is augmented. Nanoclustering dictates downstream effector recruitment, MAPK-activity, and tumorigenic cell proliferation. Our results describe an unprecedented mechanism of signaling protein activation in cancer.

MicroRNAs (miRNAs) regulate a wide range of cellular signaling pathways and biological processes in both physiological and pathological states such as cancer. We have previously identified miR-135b as a direct regulator of androgen receptor (AR) protein level in prostate cancer (PCa). We wanted to further explore the relationship of miR-135b to hormonal receptors, particularly estrogen receptor α (ERα). Here we show that miR-135b expression is lower in ERα-positive breast tumors as compared to ERα-negative samples in two independent breast cancer (BCa) patient cohorts (101 and 1302 samples). Additionally, the miR-135b expression is higher in AR-low PCa patient samples (47 samples). We identify ERα as a novel miR-135b target by demonstrating miR-135b binding to the 3'UTR of the ERα and decreased ERα protein and mRNA level upon miR-135b overexpression in BCa cells. MiR-135b reduces proliferation of ERα-positive BCa cells MCF-7 and BT-474 as well as AR-positive PCa cells LNCaP and 22Rv1 when grown in 2D. To identify other genes regulated by miR-135b we performed gene expression studies and found a link to the hypoxia inducible factor 1α (HIF1α) pathway. We show that miR-135b influences the protein level of the inhibitor for hypoxia inducible factor 1α (HIF1AN) and is able to bind to HIF1AN 3'UTR. Our study demonstrates that miR-135b regulates ERα, AR and HIF1AN protein levels through interaction with their 3'UTR regions, and proliferation in ERα-positive BCa and AR-positive PCa cells., (Copyright © 2015 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

It is widely accepted that drug discovery often requires a systems-level polypharmacology approach to tackle problems such as lack of efficacy and emerging resistance of single-targeted compounds. Network pharmacology approaches are increasingly being developed and applied to find new therapeutic opportunities and to re-purpose approved drugs. However, these recent advances have been relatively slow to be translated into the field of natural products. Here, we argue that a network pharmacology approach would enable an effective mapping of the yet unexplored target space of natural products, hence providing a systematic means to extend the druggable space of proteins implicated in various complex diseases. We give an overview of the key network pharmacology concepts and recent experimental-computational approaches that have been successfully applied to natural product research, including unbiased elucidation of mechanisms of action as well as systematic prediction of effective therapeutic combinations. We focus specifically on anticancer applications that use in vivo and in vitro functional phenotypic measurements, such as genome-wide transcriptomic response profiles, which enable a global modelling of the multi-target activity at the level of the biological pathways and interaction networks. We also provide representative examples of other disease applications, databases and tools as well as existing and emerging resources, which may prove useful for future natural product research. Finally, we offer our personal view of the current limitations, prospective developments and open questions in this exciting field.

Unlabelled: Network pharmacology-based prediction of multi-targeted drug combinations is becoming a promising strategy to improve anticancer efficacy and safety. We developed a logic-based network algorithm, called Target Inhibition Interaction using Maximization and Minimization Averaging (TIMMA), which predicts the effects of drug combinations based on their binary drug-target interactions and single-drug sensitivity profiles in a given cancer sample. Here, we report the R implementation of the algorithm (TIMMA-R), which is much faster than the original MATLAB code. The major extensions include modeling of multiclass drug-target profiles and network visualization. We also show that the TIMMA-R predictions are robust to the intrinsic noise in the experimental data, thus making it a promising high-throughput tool to prioritize drug combinations in various cancer types for follow-up experimentation or clinical applications., Availability and Implementation: TIMMA-R source code is freely available at http://cran.r-project.org/web/packages/timma/., (© The Author 2015. Published by Oxford University Press.)

A number of supervised machine learning models have recently been introduced for the prediction of drug-target interactions based on chemical structure and genomic sequence information. Although these models could offer improved means for many network pharmacology applications, such as repositioning of drugs for new therapeutic uses, the prediction models are often being constructed and evaluated under overly simplified settings that do not reflect the real-life problem in practical applications. Using quantitative drug-target bioactivity assays for kinase inhibitors, as well as a popular benchmarking data set of binary drug-target interactions for enzyme, ion channel, nuclear receptor and G protein-coupled receptor targets, we illustrate here the effects of four factors that may lead to dramatic differences in the prediction results: (i) problem formulation (standard binary classification or more realistic regression formulation), (ii) evaluation data set (drug and target families in the application use case), (iii) evaluation procedure (simple or nested cross-validation) and (iv) experimental setting (whether training and test sets share common drugs and targets, only drugs or targets or neither). Each of these factors should be taken into consideration to avoid reporting overoptimistic drug-target interaction prediction results. We also suggest guidelines on how to make the supervised drug-target interaction prediction studies more realistic in terms of such model formulations and evaluation setups that better address the inherent complexity of the prediction task in the practical applications, as well as novel benchmarking data sets that capture the continuous nature of the drug-target interactions for kinase inhibitors., (© The Author 2014. Published by Oxford University Press.)

Predicting the best treatment strategy from genomic information is a core goal of precision medicine. Here we focus on predicting drug response based on a cohort of genomic, epigenomic and proteomic profiling data sets measured in human breast cancer cell lines. Through a collaborative effort between the National Cancer Institute (NCI) and the Dialogue on Reverse Engineering Assessment and Methods (DREAM) project, we analyzed a total of 44 drug sensitivity prediction algorithms. The top-performing approaches modeled nonlinear relationships and incorporated biological pathway information. We found that gene expression microarrays consistently provided the best predictive power of the individual profiling data sets; however, performance was increased by including multiple, independent data sets. We discuss the innovations underlying the top-performing methodology, Bayesian multitask MKL, and we provide detailed descriptions of all methods. This study establishes benchmarks for drug sensitivity prediction and identifies approaches that can be leveraged for the development of new methods.

New molecular information on potential therapeutic targets or tools for noninvasive diagnosis for endometriosis are important for patient care and treatment. However, surprisingly few efforts have described endometriosis at the protein level. In this work we enumerate the proteins in patient endometrium and ovarian endometrioma by extensive and comprehensive analysis of minute amounts of cryosectioned tissues in a three-tiered mass spectrometric approach. Quantitative comparison of the tissues revealed 214 differentially expressed proteins in ovarian endometrioma and endometrium. These proteins are reported here as a resource of SRM (selected reaction monitoring) assays that are unique, standardized, and openly available. Pathway analysis of the proteome measurements revealed a potential role for Transforming growth factor β-1 in ovarian endometriosis development. Subsequent mRNA microarray analysis further revealed clear ovarian endometrioma specificity for a subset of these proteins, which was also supported by further in silico studies. In this process two important proteins emerged, Calponin-1 and EMILIN-1, that were additionally confirmed in ovarian endometrioma tissues by immunohistochemistry and Western blotting. This study provides the most comprehensive molecular description of ovarian endometriosis to date and researchers with new molecular methods and tools for high throughput patient screening using the SRM assays.

Context: Aberrant sex steroid signaling is suggested to promote endometriosis growth by several mechanisms, and the tissue concentrations of sex steroids are key determinants of the hormone action. However, their concentrations are only superficially known in the endometrium and endometriosis lesions., Objective: This study sought to evaluate whether the tissue steroid hormone concentrations in endometriosis differ from the endometrium or serum., Main Outcome Measures: Steroid analysis of serum and tissue specimens of women with endometriosis (n = 60) and healthy controls (n=16) was measured, and supporting data from quantitative RT-PCR for steroidogenic enzymes and explant cultures of a subset of specimens is provided., Results: Endometrial tissue progesterone (P4) concentrations reflected the serum P4 levels during the menstrual cycle, whereas in endometriosis lesions, the cycle-dependent change was missing. Remarkably high tissue T concentrations were measured in endometriosis lesions independent of the cycle phase, being 5-19 times higher than the corresponding serum levels. Tissue/serum ratio of T was further increased in patients with contraceptive medication. The altered tissue steroid concentrations in endometriosis were in line with the expression of various steroidogenic enzymes in the lesions, of which HSD3B2 showed constantly high expression, whereas CYP11A1 expression was low. Furthermore, the high concentration of sex steroids detected in the ovarian lesions involves their production by the lesion and by the adjacent ovarian tissue., Conclusions: Endometriosis lesions present with progestin and androgen metabolism, which are different from that of the endometrium, and the lesions are characterized by high tissue T and a loss of cyclical changes in tissue P4 concentration.

The increasing coverage and sensitivity of LC-MS/MS-based proteomics have expanded its applications in systems medicine. In particular, label-free quantitation approaches are enabling biomarker discovery in terms of statistical comparison of proteomic profiles across large numbers of clinical samples. However, it still remains poorly understood how much protein markers can add novel insights compared to markers derived from mRNA transcriptomic profiling. Using paired label-free LC-MS/MS and gene expression microarray measurements from primary samples of patients with acute myeloid leukemia (AML), we demonstrate here that while the quantitative proteomic and transcriptomic profiles were highly correlated, in general, the marker panels showing statistically significant expression changes across the disease and healthy groups were profoundly different between protein and mRNA levels. In particular, the proteomic assay enabled unique links to known leukemic processes, which were missed when using the transcriptomic profiling alone, as well as identified additional links to metabolic regulators and chromatin remodelers, such as GPX1, fumarate hydratase, and SET oncogene, which have subsequently been evaluated in independent AML samples. Overall, these results highlighted the complementary and informative view obtained from the quantitative LC-MS/MS approach into the AML deregulated signaling networks., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Viral double-stranded RNA (dsRNA) is the most important viral structure recognized by cytosolic pattern-recognition receptors of the innate immune system, and its recognition results in the activation of signaling cascades that stimulate the production of antiviral cytokines and apoptosis of infected cells. 14-3-3 proteins are ubiquitously expressed regulatory molecules that participate in a variety of cellular processes, and 14-3-3 protein-mediated signaling pathways are activated by cytoplasmic dsRNA in human keratinocytes. However, the functional role of 14-3-3 protein-mediated interactions during viral dsRNA stimulation has remained uncharacterized. Here, we used functional proteomics to identify proteins whose phosphorylation and interaction with 14-3-3 is modulated by dsRNA and to characterize the signaling pathways activated during cytosolic dsRNA-induced innate immune response in human HaCaT keratinocytes. Phosphoproteome analysis showed that several MAPK- and immune-response-related signaling pathways were activated after dsRNA stimulation. Interactome analysis identified RelA-associated inhibitor, high-mobility group proteins, and several proteins associated with host responses to viral infection as novel 14-3-3 target proteins. Functional studies showed that RelA-associated inhibitor regulated dsRNA-induced apoptosis and TNF production. Integrated network analyses of proteomic data revealed that sirtuin1 was a central molecule regulated by 14-3-3s during dsRNA stimulation. Further experiments showed that sirtuin 1 negatively regulated dsRNA-induced NFκB transcriptional activity, suppressed expression of antiviral cytokines, and protected cells from apoptosis in dsRNA-stimulated and encephalomyocarditis-virus-infected keratinocytes. In conclusion, our data highlight the importance of 14-3-3 proteins in antiviral responses and identify RelA-associated inhibitor and sirtuin 1 as novel regulators of antiviral innate immune responses., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

Androgens are key factors involved in the development and progression of prostate cancer (PCa), and PCa growth can be suppressed by androgen deprivation therapy. In a considerable proportion of men receiving androgen deprivation therapy, however, PCa progresses to castration-resistant PCa (CRPC), making the development of efficient therapies challenging. We used an orthotopic VCaP human PCa xenograft model to study cellular and molecular changes in tumors after androgen deprivation therapy (castration). Tumor growth was monitored through weekly serum prostate-specific antigen measurements, and mice with recurrent tumors after castration were randomized to treatment groups. Serum prostate-specific antigen concentrations showed significant correlation with tumor volume. Castration-resistant tumors retained concentrations of intratumoral androgen (androstenedione, testosterone, and 5α-dihydrotestosterone) at levels similar to tumors growing in intact hosts. Accordingly, castration induced up-regulation of enzymes involved in androgen synthesis (CYP17A1, AKR1C3, and HSD17B6), as well as expression of full-length androgen receptor (AR) and AR splice variants (AR-V1 and AR-V7). Furthermore, AR target gene expression was maintained in castration-resistant xenografts. The AR antagonists enzalutamide (MDV3100) and ARN-509 suppressed PSA production of castration-resistant tumors, confirming the androgen dependency of these tumors. Taken together, the findings demonstrate that our VCaP xenograft model exhibits the key characteristics of clinical CRPC and thus provides a valuable tool for identifying druggable targets and for testing therapeutic strategies targeting AR signaling in CRPC., (Copyright © 2014 American Society for Investigative Pathology. Published by Elsevier Inc. All rights reserved.)

Background: It is possible to identify thousands of phosphopeptides and -proteins in a single experiment with mass spectrometry-based phosphoproteomics. However, a current bottleneck is the downstream data analysis which is often laborious and requires a number of manual steps., Results: Toward automating the analysis steps, we have developed and implemented a software, PhosFox, which enables peptide-level processing of phosphoproteomic data generated by multiple protein identification search algorithms, including Mascot, Sequest, and Paragon, as well as cross-comparison of their identification results. The software supports both qualitative and quantitative phosphoproteomics studies, as well as multiple between-group comparisons. Importantly, PhosFox detects uniquely phosphorylated peptides and proteins in one sample compared to another. It also distinguishes differences in phosphorylation sites between phosphorylated proteins in different samples. Using two case study examples, a qualitative phosphoproteome dataset from human keratinocytes and a quantitative phosphoproteome dataset from rat kidney inner medulla, we demonstrate here how PhosFox facilitates an efficient and in-depth phosphoproteome data analysis. PhosFox was implemented in the Perl programming language and it can be run on most common operating systems. Due to its flexible interface and open source distribution, the users can easily incorporate the program into their MS data analysis workflows and extend the program with new features. PhosFox source code, implementation and user instructions are freely available from https://bitbucket.org/phintsan/phosfox., Conclusions: PhosFox facilitates efficient and more in-depth comparisons between phosphoproteins in case-control settings. The open source implementation is easily extendable to accommodate additional features for widespread application use cases.

We developed a systematic algorithmic solution for quantitative drug sensitivity scoring (DSS), based on continuous modeling and integration of multiple dose-response relationships in high-throughput compound testing studies. Mathematical model estimation and continuous interpolation makes the scoring approach robust against sources of technical variability and widely applicable to various experimental settings, both in cancer cell line models and primary patient-derived cells. Here, we demonstrate its improved performance over other response parameters especially in a leukemia patient case study, where differential DSS between patient and control cells enabled identification of both cancer-selective drugs and drug-sensitive patient sub-groups, as well as dynamic monitoring of the response patterns and oncogenic driver signals during cancer progression and relapse in individual patient cells ex vivo. An open-source and easily extendable implementation of the DSS calculation is made freely available to support its tailored application to translating drug sensitivity testing results into clinically actionable treatment options.

Prostate cancer is the most common cancer of men in the Western world, and novel approaches for prostate cancer risk reduction are needed. Plant-derived phenolic compounds attenuate prostate cancer growth in preclinical models by several mechanisms, which is in line with epidemiological findings suggesting that consumption of plant-based diets is associated with low risk of prostate cancer. The objective of this study was to assess the effects of a novel lignan-stilbenoid mixture in PC-3M-luc2 human prostate cancer cells in vitro and in orthotopic xenografts. Lignan and stilbenoid -rich extract was obtained from Scots pine (Pinus sylvestris) knots. Pine knot extract as well as stilbenoids (methyl pinosylvin and pinosylvin), and lignans (matairesinol and nortrachelogenin) present in pine knot extract showed antiproliferative and proapoptotic efficacy at ≥ 40 μM concentration in vitro. Furthermore, pine knot extract derived stilbenoids enhanced tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) induced apoptosis already at ≥ 10 μM concentrations. In orthotopic PC-3M-luc2 xenograft bearing immunocompromized mice, three-week peroral exposure to pine knot extract (52 mg of lignans and stilbenoids per kg of body weight) was well tolerated and showed anti-tumorigenic efficacy, demonstrated by multivariate analysis combining essential markers of tumor growth (i.e. tumor volume, vascularization, and cell proliferation). Methyl pinosylvin, pinosylvin, matairesinol, nortrachelogenin, as well as resveratrol, a metabolite of pinosylvin, were detected in serum at total concentration of 7-73 μM, confirming the bioavailability of pine knot extract derived lignans and stilbenoids. In summary, our data indicates that pine knot extract is a novel and cost-effective source of resveratrol, methyl pinosylvin and other bioactive lignans and stilbenoids. Pine knot extract shows anticarcinogenic efficacy in preclinical prostate cancer model, and our in vitro data suggests that compounds derived from the extract may have potential as novel chemosensitizers to TRAIL. These findings promote further research on health-related applications of wood biochemicals.

We carried out a systematic evaluation of target selectivity profiles across three recent large-scale biochemical assays of kinase inhibitors and further compared these standardized bioactivity assays with data reported in the widely used databases ChEMBL and STITCH. Our comparative evaluation revealed relative benefits and potential limitations among the bioactivity types, as well as pinpointed biases in the database curation processes. Ignoring such issues in data heterogeneity and representation may lead to biased modeling of drugs' polypharmacological effects as well as to unrealistic evaluation of computational strategies for the prediction of drug-target interaction networks. Toward making use of the complementary information captured by the various bioactivity types, including IC50, K(i), and K(d), we also introduce a model-based integration approach, termed KIBA, and demonstrate here how it can be used to classify kinase inhibitor targets and to pinpoint potential errors in database-reported drug-target interactions. An integrated drug-target bioactivity matrix across 52,498 chemical compounds and 467 kinase targets, including a total of 246,088 KIBA scores, has been made freely available.

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.

Unlabelled: We present an individualized systems medicine (ISM) approach to optimize cancer drug therapies one patient at a time. ISM is based on (i) molecular profiling and ex vivo drug sensitivity and resistance testing (DSRT) of patients' cancer cells to 187 oncology drugs, (ii) clinical implementation of therapies predicted to be effective, and (iii) studying consecutive samples from the treated patients to understand the basis of resistance. Here, application of ISM to 28 samples from patients with acute myeloid leukemia (AML) uncovered five major taxonomic drug-response subtypes based on DSRT profiles, some with distinct genomic features (e.g., MLL gene fusions in subgroup IV and FLT3-ITD mutations in subgroup V). Therapy based on DSRT resulted in several clinical responses. After progression under DSRT-guided therapies, AML cells displayed significant clonal evolution and novel genomic changes potentially explaining resistance, whereas ex vivo DSRT data showed resistance to the clinically applied drugs and new vulnerabilities to previously ineffective drugs., Significance: Here, we demonstrate an ISM strategy to optimize safe and effective personalized cancer therapies for individual patients as well as to understand and predict disease evolution and the next line of therapy. This approach could facilitate systematic drug repositioning of approved targeted drugs as well as help to prioritize and de-risk emerging drugs for clinical testing., (©2013 AACR.)

Checkpoint kinase Chk1 is constitutively active in many cancer cell types and new generation Chk1 inhibitors show marked antitumor activity as single agents. Here we present a hitherto unrecognized mechanism that contributes to the response of cancer cells to Chk1-targeted therapy. Inhibiting chronic Chk1 activity in cancer cells induced the tumor suppressor activity of protein phosphatase protein phosphatase 2A (PP2A), which by dephosphorylating MYC serine 62, inhibited MYC activity and impaired cancer cell survival. Mechanistic investigations revealed that Chk1 inhibition activated PP2A by decreasing the transcription of cancerous inhibitor of PP2A (CIP2A), a chief inhibitor of PP2A activity. Inhibition of cancer cell clonogenicity by Chk1 inhibition could be rescued in vitro either by exogenous expression of CIP2A or by blocking the CIP2A-regulated PP2A complex. Chk1-mediated CIP2A regulation was extended in tumor models dependent on either Chk1 or CIP2A. The clinical relevance of CIP2A as a Chk1 effector protein was validated in several human cancer types, including neuroblastoma, where CIP2A was identified as an NMYC-independent prognostic factor. Because the Chk1-CIP2A-PP2A pathway is driven by DNA-PK activity, functioning regardless of p53 or ATM/ATR status, our results offer explanative power for understanding how Chk1 inhibitors mediate single-agent anticancer efficacy. Furthermore, they define CIP2A-PP2A status in cancer cells as a pharmacodynamic marker for their response to Chk1-targeted therapy., (©2013 AACR)

Chemical-genomic and genetic interaction profiling approaches are widely used to study mechanisms of drug action and resistance. However, there exist a number of scoring algorithms customized to different experimental assays, the relative performance of which remains poorly understood, especially with respect to different types of chemogenetic assays. Using yeast Saccharomyces cerevisiae as a test bed, we carried out a systematic evaluation among the main drug target analysis approaches in terms of predicting global drug-target interaction networks. We found drastic differences in their performance across different chemical-genomic assay types, such as those based on heterozygous and homozygous diploid or haploid deletion mutant libraries. Moreover, a relatively small overlap in the predicted targets was observed between those approaches that use either chemical-genomic screening alone or combined with genetic interaction profiling. A rank-based integration of the complementary scoring approaches led to improved overall performance, demonstrating that genetic interaction profiling provides added information on drug target prediction. Optimal performance was achieved when focusing specifically on the negative tail of the genetic interactions, suggesting that combining synthetic lethal interactions with chemical-genetic interactions provides highest information on drug-target interactions. A network view of rapamycin-interacting genes, pathways and complexes was used as an example to demonstrate the benefits of such integrated and optimized analysis of chemogenetic assays in yeast.

A central challenge in systems biology and medical genetics is to understand how interactions among genetic loci contribute to complex phenotypic traits and human diseases. While most studies have so far relied on statistical modeling and association testing procedures, machine learning and predictive modeling approaches are increasingly being applied to mining genotype-phenotype relationships, also among those associations that do not necessarily meet statistical significance at the level of individual variants, yet still contributing to the combined predictive power at the level of variant panels. Network-based analysis of genetic variants and their interaction partners is another emerging trend by which to explore how sub-network level features contribute to complex disease processes and related phenotypes. In this review, we describe the basic concepts and algorithms behind machine learning-based genetic feature selection approaches, their potential benefits and limitations in genome-wide setting, and how physical or genetic interaction networks could be used as a priori information for providing improved predictive power and mechanistic insights into the disease networks. These developments are geared toward explaining a part of the missing heritability, and when combined with individual genomic profiling, such systems medicine approaches may also provide a principled means for tailoring personalized treatment strategies in the future.

A recent trend in drug development is to identify drug combinations or multi-target agents that effectively modify multiple nodes of disease-associated networks. Such polypharmacological effects may reduce the risk of emerging drug resistance by means of attacking the disease networks through synergistic and synthetic lethal interactions. However, due to the exponentially increasing number of potential drug and target combinations, systematic approaches are needed for prioritizing the most potent multi-target alternatives on a global network level. We took a functional systems pharmacology approach toward the identification of selective target combinations for specific cancer cells by combining large-scale screening data on drug treatment efficacies and drug-target binding affinities. Our model-based prediction approach, named TIMMA, takes advantage of the polypharmacological effects of drugs and infers combinatorial drug efficacies through system-level target inhibition networks. Case studies in MCF-7 and MDA-MB-231 breast cancer and BxPC-3 pancreatic cancer cells demonstrated how the target inhibition modeling allows systematic exploration of functional interactions between drugs and their targets to maximally inhibit multiple survival pathways in a given cancer type. The TIMMA prediction results were experimentally validated by means of systematic siRNA-mediated silencing of the selected targets and their pairwise combinations, showing increased ability to identify not only such druggable kinase targets that are essential for cancer survival either individually or in combination, but also synergistic interactions indicative of non-additive drug efficacies. These system-level analyses were enabled by a novel model construction method utilizing maximization and minimization rules, as well as a model selection algorithm based on sequential forward floating search. Compared with an existing computational solution, TIMMA showed both enhanced prediction accuracies in cross validation as well as significant reduction in computation times. Such cost-effective computational-experimental design strategies have the potential to greatly speed-up the drug testing efforts by prioritizing those interventions and interactions warranting further study in individual cancer cases.

Influenza A viruses (IAVs) infect humans and cause significant morbidity and mortality. Different treatment options have been developed; however, these were insufficient during recent IAV outbreaks. Here, we conducted a targeted chemical screen in human nonmalignant cells to validate known and search for novel host-directed antivirals. The screen validated saliphenylhalamide (SaliPhe) and identified two novel anti-IAV agents, obatoclax and gemcitabine. Further experiments demonstrated that Mcl-1 (target of obatoclax) provides a novel host target for IAV treatment. Moreover, we showed that obatoclax and SaliPhe inhibited IAV uptake and gemcitabine suppressed viral RNA transcription and replication. These compounds possess broad spectrum antiviral activity, although their antiviral efficacies were virus-, cell type-, and species-specific. Altogether, our results suggest that phase II obatoclax, investigational SaliPhe, and FDA/EMEA-approved gemcitabine represent potent antiviral agents.

Aims: Extracellular ATP and ADP regulate diverse inflammatory, prothrombotic and vasoactive responses in the vasculature. Statins have been shown to modulate their signaling pathways in vitro. We hypothesized that altered intravascular nucleotide turnover modulates vasodilation in patients with type 1 diabetes (T1DM), and this can be partly restored with pravastatin therapy., Methods: In this randomized double blind study, plasma ATP and ADP levels and echocardiography-derived coronary flow velocity response to cold pressor test (CPT) were concurrently assessed in 42 normocholesterolemic patients with T1DM (age 30 ± 6 years, LDL cholesterol 2.5 ± 0.6 mmol/L) before and after four-month treatment with pravastatin 40 mg/day or placebo (n = 22 and n = 20, respectively), and in 41 healthy control subjects., Results: Compared to controls, T1DM patients had significantly higher concentrations of ATP (p < 0.01) and ADP (p < 0.01) and these levels were partly restored after treatment with pravastatin (p = 0.002 and p = 0.007, respectively), but not after placebo (p = 0.06 and p = 0.14, respectively). Coronary flow velocity acceleration was significantly lower in T1DM patients compared to control subjects, and it increased from pre- to post-intervention in the pravastatin (p = 0.02), but not in placebo group (p = 0.15)., Conclusions: Pravastatin treatment significantly reduces circulating ATP and ADP levels of T1DM patients, and concurrently improves coronary flow response to CPT. This study provides a novel insight in purinergic mechanisms involved in pleiotropic effects of pravastatin.

Purpose: Preclinical tumor growth experiments often result in heterogeneous datasets that include growing, regressing, or stable growth profiles in the treatment and control groups. Such confounding intertumor variability may mask the true treatment effects especially when less aggressive treatment alternatives are being evaluated., Experimental Design: We developed a statistical modeling approach in which the growing and poorly growing tumor categories were automatically detected by means of an expectation-maximization algorithm coupled within a mixed-effects modeling framework. The framework is implemented and distributed as an R package, which enables model estimation and statistical inference, as well as statistical power and precision analyses., Results: When applied to four tumor growth experiments, the modeling framework was shown to (i) improve the detection of subtle treatment effects in the presence of high within-group tumor variability; (ii) reveal hidden tumor subgroups associated with established or novel biomarkers, such as ERβ expression in a MCF-7 breast cancer model, which remained undetected with standard statistical analysis; (iii) provide guidance on the selection of sufficient sample sizes and most informative treatment periods; and (iv) offer flexibility to various cancer models, experimental designs, and treatment options. Model-based testing of treatment effect on the tumor growth rate (or slope) was shown as particularly informative in the preclinical assessment of treatment alternatives based on dietary interventions., Conclusions: In general, the modeling framework enables identification of such biologically significant differences in tumor growth profiles that would have gone undetected or had required considerably higher number of animals when using traditional statistical methods.

Background: Through the wealth of information contained within them, genome-wide association studies (GWAS) have the potential to provide researchers with a systematic means of associating genetic variants with a wide variety of disease phenotypes. Due to the limitations of approaches that have analyzed single variants one at a time, it has been proposed that the genetic basis of these disorders could be determined through detailed analysis of the genetic variants themselves and in conjunction with one another. The construction of models that account for these subsets of variants requires methodologies that generate predictions based on the total risk of a particular group of polymorphisms. However, due to the excessive number of variants, constructing these types of models has so far been computationally infeasible., Results: We have implemented an algorithm, known as greedy RLS, that we use to perform the first known wrapper-based feature selection on the genome-wide level. The running time of greedy RLS grows linearly in the number of training examples, the number of features in the original data set, and the number of selected features. This speed is achieved through computational short-cuts based on matrix calculus. Since the memory consumption in present-day computers can form an even tighter bottleneck than running time, we also developed a space efficient variation of greedy RLS which trades running time for memory. These approaches are then compared to traditional wrapper-based feature selection implementations based on support vector machines (SVM) to reveal the relative speed-up and to assess the feasibility of the new algorithm. As a proof of concept, we apply greedy RLS to the Hypertension - UK National Blood Service WTCCC dataset and select the most predictive variants using 3-fold external cross-validation in less than 26 minutes on a high-end desktop. On this dataset, we also show that greedy RLS has a better classification performance on independent test data than a classifier trained using features selected by a statistical p-value-based filter, which is currently the most popular approach for constructing predictive models in GWAS., Conclusions: Greedy RLS is the first known implementation of a machine learning based method with the capability to conduct a wrapper-based feature selection on an entire GWAS containing several thousand examples and over 400,000 variants. In our experiments, greedy RLS selected a highly predictive subset of genetic variants in a fraction of the time spent by wrapper-based selection methods used together with SVM classifiers. The proposed algorithms are freely available as part of the RLScore software library at http://users.utu.fi/aatapa/RLScore/.

We developed a computational procedure for optimizing the binding site detections in a given ChIP-seq experiment by maximizing their reproducibility under bootstrap sampling. We demonstrate how the procedure can improve the detection accuracies beyond those obtained with the default settings of popular peak calling software, or inform the user whether the peak detection results are compromised, circumventing the need for arbitrary re-iterative peak calling under varying parameter settings. The generic, open-source implementation is easily extendable to accommodate additional features and to promote its widespread application in future ChIP-seq studies. The peakROTS R-package and user guide are freely available at http://www.nic.funet.fi/pub/sci/molbio/peakROTS.

Influenza A viruses are important pathogens that cause acute respiratory diseases and annual epidemics in humans. Macrophages recognize influenza A virus infection with their pattern recognition receptors, and are involved in the activation of proper innate immune response. Here, we have used high-throughput subcellular proteomics combined with bioinformatics to provide a global view of host cellular events that are activated in response to influenza A virus infection in human primary macrophages. We show that viral infection regulates the expression and/or subcellular localization of more than one thousand host proteins at early phases of infection. Our data reveals that there are dramatic changes in mitochondrial and nuclear proteomes in response to infection. We show that a rapid cytoplasmic leakage of lysosomal proteins, including cathepsins, followed by their secretion, contributes to inflammasome activation and apoptosis seen in the infected macrophages. Also, our results demonstrate that P2X₇ receptor and src tyrosine kinase activity are essential for inflammasome activation during influenza A virus infection. Finally, we show that influenza A virus infection is associated with robust secretion of different danger-associated molecular patterns (DAMPs) suggesting an important role for DAMPs in host response to influenza A virus infection. In conclusion, our high-throughput quantitative proteomics study provides important new insight into host-response against influenza A virus infection in human primary macrophages.

Background: High-throughput genetic screening approaches have enabled systematic means to study how interactions among gene mutations contribute to quantitative fitness phenotypes, with the aim of providing insights into the functional wiring diagrams of genetic interaction networks on a global scale. However, it is poorly known how well these quantitative interaction measurements agree across the screening approaches, which hinders their integrated use toward improving the coverage and quality of the genetic interaction maps in yeast and other organisms., Results: Using large-scale data matrices from epistatic miniarray profiling (E-MAP), genetic interaction mapping (GIM), and synthetic genetic array (SGA) approaches, we carried out here a systematic comparative evaluation among these quantitative maps of genetic interactions in yeast. The relatively low association between the original interaction measurements or their customized scores could be improved using a matrix-based modelling framework, which enables the use of single- and double-mutant fitness estimates and measurements, respectively, when scoring genetic interactions. Toward an integrative analysis, we show how the detections from the different screening approaches can be combined to suggest novel positive and negative interactions which are complementary to those obtained using any single screening approach alone. The matrix approximation procedure has been made available to support the design and analysis of the future screening studies., Conclusions: We have shown here that even if the correlation between the currently available quantitative genetic interaction maps in yeast is relatively low, their comparability can be improved by means of our computational matrix approximation procedure, which will enable integrative analysis and detection of a wider spectrum of genetic interactions using data from the complementary screening approaches.

We present a versatile user-friendly software tool, PolyAlign, for the alignment of multiple LC-MS signal maps with the option of manual landmark setting or automated alignment. One of the spectral images is selected as a reference map, and after manually setting the landmarks, the program warps the images using either polynomial or Hermite transformation. The software provides an option for automated landmark finding. The software includes a very fast zoom-in function synchronized between the images, which facilitate detecting correspondences between the adjacent images. Such an interactive visual process enables the analyst to decide when the alignment is satisfactory and to correct known irregularities. We demonstrate that the software provides significant improvements in the alignment of LC-MALDI data, with 10-15 landmark pairs, and it is also applicable to correcting electrospray LC-MS data. The results with practical data show substantial improvement in peak alignment compared to MZmine, which was among the best analysis packages in a recent assessment. The PolyAlign software is freely available and easily accessible as an integrated component of the popular MZmine software, and also as a simpler stand-alone Perl implementation to preview data and apply landmark directed polynomial transformation.

Probe defects are a major source of noise in gene expression studies. While existing approaches detect noisy probes based on external information such as genomic alignments, we introduce and validate a targeted probabilistic method for analyzing probe reliability directly from expression data and independently of the noise source. This provides insights into the various sources of probe-level noise and gives tools to guide probe design.

The relative contribution of genetic risk factors to the progression of subclinical atherosclerosis is poorly understood. It is likely that multiple variants are implicated in the development of atherosclerosis, but the subtle genotypic and phenotypic differences are beyond the reach of the conventional case-control designs and the statistical significance testing procedures being used in most association studies. Our objective here was to investigate whether an alternative approach--in which common disorders are treated as quantitative phenotypes that are continuously distributed over a population--can reveal predictive insights into the early atherosclerosis, as assessed using ultrasound imaging-based quantitative measurement of carotid artery intima-media thickness (IMT). Using our population-based follow-up study of atherosclerosis precursors as a basis for sampling subjects with gradually increasing IMT levels, we searched for such subsets of genetic variants and their interactions that are the most predictive of the various risk classes, rather than using exclusively those variants meeting a stringent level of statistical significance. The area under the receiver operating characteristic curve (AUC) was used to evaluate the predictive value of the variants, and cross-validation was used to assess how well the predictive models will generalize to other subsets of subjects. By means of our predictive modeling framework with machine learning-based SNP selection, we could improve the prediction of the extreme classes of atherosclerosis risk and progression over a 6-year period (average AUC 0.844 and 0.761), compared to that of using conventional cardiovascular risk factors alone (average AUC 0.741 and 0.629), or when combined with the statistically significant variants (average AUC 0.762 and 0.651). The predictive accuracy remained relatively high in an independent validation set of subjects (average decrease of 0.043). These results demonstrate that the modeling framework can utilize the "gray zone" of genetic variation in the classification of subjects with different degrees of risk of developing atherosclerosis., Competing Interests: The authors have declared that no competing interests exist.

Rapid nondestructive screening of mutants is a common step in many research projects in plant biology. Here we report the development of a method that uses kinetic imaging of chlorophyll fluorescence to detect phenotypes that differ from wild-type plants. The method uses multiple fluorescence features simultaneously in order to catch different types of photosynthesis-related mutants with a single assay. The Mahalanobis distance was used to evaluate the degree of similarity in fluorescence features between the wild-type and test plants, and plants differing strongly from the wild-type were classified as mutants. The method was tested on a collection of photosynthesis-related mutants of Arabidopsis thaliana. The plants were evaluated from images in which the color of each pixel depended on the Mahalanobis distance of the fluorescence features. Two parameters of the color-coding procedure were used to adjust the trade-off between detection of true mutants and erratic classification of wild-type plants as mutants. We found that a large percentage of photosynthesis-related mutants can be detected with this method. Scripts for the free statistics software R are provided to facilitate the practical application of the method.