Acidic extracellular pH is characteristic of the cell microenvironment in several important physiological and pathological contexts. Although it is well established that acidic extracellular pH can have profound effects on processes such as cell adhesion and migration, the underlying molecular mechanisms are largely unknown. Integrin receptors physically connect cells to the extracellular matrix, and are thus likely to modulate cell responses to extracellular conditions. Here, we examine the role of acidic extracellular pH in regulating activation of integrin α(v)β(3). Through computational molecular dynamics simulations, we find that acidic extracellular pH promotes opening of the α(v)β(3) headpiece, indicating that acidic pH can thereby facilitate integrin activation. This prediction is consistent with our flow cytometry and atomic force microscope-mediated force spectroscopy assays of integrin α(v)β(3) on live cells, which both demonstrate that acidic pH promotes activation at the intact cell surface. Finally, quantification of cell morphology and migration measurements shows that acidic extracellular pH affects cell behavior in a manner that is consistent with increased integrin activation. Taken together, these computational and experimental results suggest a new and complementary mechanism of integrin activation regulation, with associated implications for cell adhesion and migration in regions of altered pH that are relevant to wound healing and cancer.

Recent advances in computational biology suggest that any perturbation to the transcriptional programme of the cell can be summarised by a proper 'signature': a set of genes combined with a pattern of expression. Therefore, it should be possible to generate proxies of clinicopathological phenotypes and drug effects through signatures acquired via DNA microarray technology. Gene expression signatures have recently been assembled and compared through genome-wide metrics, unveiling unexpected drug-disease and drug-drug 'connections' by matching corresponding signatures. Consequently, novel applications for existing drugs have been predicted and experimentally validated. Here, we describe related methods, case studies and resources while discussing challenges and benefits of exploiting existing repositories of microarray data that could serve as a search space for systematic drug repositioning.

Mass spectrometry is widely used to probe the proteome and its modifications in an untargeted manner, with unrivalled coverage. Applied to phosphoproteomics, it has tremendous potential to interrogate phospho-signalling and its therapeutic implications. However, this task is complicated by issues of undersampling of the phosphoproteome and challenges stemming from its high-content but low-sample-throughput nature. Hence, methods using such data to reconstruct signalling networks have been limited to restricted data sets and insights (for example, groups of kinases likely to be active in a sample). We propose a new method to handle high-content discovery phosphoproteomics data on perturbation by putting it in the context of kinase/phosphatase-substrate knowledge, from which we derive and train logic models. We show, on a data set obtained through perturbations of cancer cells with small-molecule inhibitors, that this method can study the targets and effects of kinase inhibitors, and reconcile insights obtained from multiple data sets, a common issue with these data.

Dysregulation of ErbB-family signaling underlies numerous pathologies and has been therapeutically targeted through inhibiting ErbB-receptors themselves or their cognate ligands. For the latter, "decoy" antibodies have been developed to sequester ligands including heparin-binding epidermal growth factor (HB-EGF); however, demonstrating sufficient efficacy has been difficult. Here, we hypothesized that this strategy depends on properties such as ligand-receptor binding affinity, which varies widely across the known ErbB-family ligands. Guided by computational modeling, we found that high-affinity ligands such as HB-EGF are more difficult to target with decoy antibodies compared to low-affinity ligands such as amphiregulin (AREG). To address this issue, we developed an alternative method for inhibiting HB-EGF activity by targeting its cleavage from the cell surface. In a model of the invasive disease endometriosis, we identified A Disintegrin and Metalloproteinase 12 (ADAM12) as a protease implicated in HB-EGF shedding. We designed a specific inhibitor of ADAM12 based on its recombinant prodomain (PA12), which selectively inhibits ADAM12 but not ADAM10 or ADAM17. In endometriotic cells, PA12 significantly reduced HB-EGF shedding and resultant cellular migration. Overall, specific inhibition of ligand shedding represents a possible alternative to decoy antibodies, especially for ligands such as HB-EGF that exhibit high binding affinity and localized signaling.

Cells respond to environmental stimuli via specialized signaling pathways. Concurrent stimuli trigger multiple pathways that integrate information, predominantly via protein phosphorylation. Budding yeast responds to NaCl and pheromone via two mitogen-activated protein kinase cascades, the high osmolarity, and the mating pathways, respectively. To investigate signal integration between these pathways, we quantified the time-resolved phosphorylation site dynamics after pathway co-stimulation. Using shotgun mass spectrometry, we quantified 2,536 phosphopeptides across 36 conditions. Our data indicate that NaCl and pheromone affect phosphorylation events within both pathways, which thus affect each other at more levels than anticipated, allowing for information exchange and signal integration. We observed a pheromone-induced down-regulation of Hog1 phosphorylation due to Gpd1, Ste20, Ptp2, Pbs2, and Ptc1. Distinct Ste20 and Pbs2 phosphosites responded differently to the two stimuli, suggesting these proteins as key mediators of the information exchange. A set of logic models was then used to assess the role of measured phosphopeptides in the crosstalk. Our results show that the integration of the response to different stimuli requires complex interconnections between signaling pathways.

The high failure rate of therapeutics showing promise in mouse models to translate to patients is a pressing challenge in biomedical science. Though retrospective studies have examined the fidelity of mouse models to their respective human conditions, approaches for prospective translation of insights from mouse models to patients remain relatively unexplored. Here, we develop a semi-supervised learning approach for inference of disease-associated human differentially expressed genes and pathways from mouse model experiments. We examined 36 transcriptomic case studies where comparable phenotypes were available for mouse and human inflammatory diseases and assessed multiple computational approaches for inferring human biology from mouse datasets. We found that semi-supervised training of a neural network identified significantly more true human biological associations than interpreting mouse experiments directly. Evaluating the experimental design of mouse experiments where our model was most successful revealed principles of experimental design that may improve translational performance. Our study shows that when prospectively evaluating biological associations in mouse studies, semi-supervised learning approaches, combining mouse and human data for biological inference, provide the most accurate assessment of human in vivo disease processes. Finally, we proffer a delineation of four categories of model system-to-human "Translation Problems" defined by the resolution and coverage of the datasets available for molecular insight translation and suggest that the task of translating insights from model systems to human disease contexts may be better accomplished by a combination of translation-minded experimental design and computational approaches.

Activation of brown adipose tissue-mediated thermogenesis is a strategy for tackling obesity and promoting metabolic health. BMP8b is secreted by brown/beige adipocytes and enhances energy dissipation. Here we show that adipocyte-secreted BMP8b contributes to adrenergic-induced remodeling of the neuro-vascular network in adipose tissue (AT). Overexpression of bmp8b in AT enhances browning of the subcutaneous depot and maximal thermogenic capacity. Moreover, BMP8b-induced browning, increased sympathetic innervation and vascularization of AT were maintained at 28 °C, a condition of low adrenergic output. This reinforces the local trophic effect of BMP8b. Innervation and vascular remodeling effects required BMP8b signaling through the adipocytes to 1) secrete neuregulin-4 (NRG4), which promotes sympathetic axon growth and branching in vitro, and 2) induce a pro-angiogenic transcriptional and secretory profile that promotes vascular sprouting. Thus, BMP8b and NRG4 can be considered as interconnected regulators of neuro-vascular remodeling in AT and are potential therapeutic targets in obesity.

Targeted cancer therapeutics have demonstrated more limited clinical efficacy than anticipated, due to both intrinsic and acquired drug resistance. Underlying mechanisms have been largely attributed to genetic changes, but a substantial proportion of resistance observations remain unexplained by genomic properties. Emerging evidence shows that receptor tyrosine kinase (RTK) reprogramming is a major alternative process causing targeted drug resistance, separate from genetic alterations. Hence, the contributions of mechanisms leading to this process need to be more rigorously assessed.

Despite the abundance of large-scale molecular and drug-response data, the insights gained about the mechanisms underlying treatment efficacy in cancer has been in general limited. Machine learning algorithms applied to those datasets most often are used to provide predictions without interpretation, or reveal single drug-gene association and fail to derive robust insights. We propose to use Macau, a bayesian multitask multi-relational algorithm to generalize from individual drugs and genes and explore the interactions between the drug targets and signaling pathways' activation. A typical insight would be: "Activation of pathway Y will confer sensitivity to any drug targeting protein X". We applied our methodology to the Genomics of Drug Sensitivity in Cancer (GDSC) screening, using gene expression of 990 cancer cell lines, activity scores of 11 signaling pathways derived from the tool PROGENy as cell line input and 228 nominal targets for 265 drugs as drug input. These interactions can guide a tissue-specific combination treatment strategy, for example suggesting to modulate a certain pathway to maximize the drug response for a given tissue. We confirmed in literature drug combination strategies derived from our result for brain, skin and stomach tissues. Such an analysis of interactions across tissues might help target discovery, drug repurposing and patient stratification strategies.

Patients with seemingly the same tumour can respond very differently to treatment. There are strong, well-established effects of somatic mutations on drug efficacy, but there is at-most anecdotal evidence of a germline component to drug response. Here, we report a systematic survey of how inherited germline variants affect drug susceptibility in cancer cell lines. We develop a joint analysis approach that leverages both germline and somatic variants, before applying it to screening data from 993 cell lines and 265 drugs. Surprisingly, we find that the germline contribution to variation in drug susceptibility can be as large or larger than effects due to somatic mutations. Several of the associations identified have a direct relationship to the drug target. Finally, using 17-AAG response as an example, we show how germline effects in combination with transcriptomic data can be leveraged for improved patient stratification and to identify new markers for drug sensitivity.

Drug versus Disease (DvD) provides a pipeline, available through R or Cytoscape, for the comparison of drug and disease gene expression profiles from public microarray repositories. Negatively correlated profiles can be used to generate hypotheses of drug-repurposing, whereas positively correlated profiles may be used to infer side effects of drugs. DvD allows users to compare drug and disease signatures with dynamic access to databases Array Express, Gene Expression Omnibus and data from the Connectivity Map.

During tumor progression, alternative splicing gives rise to different Mena protein isoforms. We analyzed how Mena11a, an isoform enriched in epithelia and epithelial-like cells, affects Mena-dependent regulation of actin dynamics and cell behavior. While other Mena isoforms promote actin polymerization and drive membrane protrusion, we find that Mena11a decreases actin polymerization and growth factor-stimulated membrane protrusion at lamellipodia. Ectopic Mena11a expression slows mesenchymal-like cell motility, while isoform-specific depletion of endogenous Mena11a in epithelial-like tumor cells perturbs cell:cell junctions and increases membrane protrusion and overall cell motility. Mena11a can dampen membrane protrusion and reduce actin polymerization in the absence of other Mena isoforms, indicating that it is not simply an inactive Mena isoform. We identify a phosphorylation site within 11a that is required for some Mena11a-specific functions. RNA-seq data analysis from patient cohorts demonstrates that the difference between mRNAs encoding constitutive Mena sequences and those containing the 11a exon correlates with metastasis in colorectal cancer, suggesting that 11a exon exclusion contributes to invasive phenotypes and leads to poor clinical outcomes.

Predicting the response of a specific cancer to a therapy is a major goal in modern oncology that should ultimately lead to a personalised treatment. High-throughput screenings of potentially active compounds against a panel of genomically heterogeneous cancer cell lines have unveiled multiple relationships between genomic alterations and drug responses. Various computational approaches have been proposed to predict sensitivity based on genomic features, while others have used the chemical properties of the drugs to ascertain their effect. In an effort to integrate these complementary approaches, we developed machine learning models to predict the response of cancer cell lines to drug treatment, quantified through IC₅₀ values, based on both the genomic features of the cell lines and the chemical properties of the considered drugs. Models predicted IC₅₀ values in a 8-fold cross-validation and an independent blind test with coefficient of determination R² of 0.72 and 0.64 respectively. Furthermore, models were able to predict with comparable accuracy (R² of 0.61) IC50s of cell lines from a tissue not used in the training stage. Our in silico models can be used to optimise the experimental design of drug-cell screenings by estimating a large proportion of missing IC₅₀ values rather than experimentally measuring them. The implications of our results go beyond virtual drug screening design: potentially thousands of drugs could be probed in silico to systematically test their potential efficacy as anti-tumour agents based on their structure, thus providing a computational framework to identify new drug repositioning opportunities as well as ultimately be useful for personalized medicine by linking the genomic traits of patients to drug sensitivity.

Epithelial-mesenchymal transition (EMT), whether in developmental morphogenesis or malignant transformation, prominently involves modified cell motility behavior. Although major advances have transpired in understanding the molecular pathways regulating the process of EMT induction per se by certain environmental stimuli, an important outstanding question is how the activities of signaling pathways governing motility yield the diverse movement behaviors characteristic of pre-induction versus postinduction states across a broad landscape of growth factor contexts. For the particular case of EMT induction in human mammary cells by ectopic expression of the transcription factor Twist, we found the migration responses to a panel of growth factors (EGF, HRG, IGF, HGF) dramatically disparate between confluent pre-Twist epithelial cells and sparsely distributed post-Twist mesenchymal cells-but that a computational model quantitatively integrating multiple key signaling node activities could nonetheless account for this full range of behavior. Moreover, motility in both conditions was successfully predicted a priori for an additional growth factor (PDGF) treatment. Although this signaling network state model could comprehend motility behavior globally, modulation of the network interactions underlying the altered pathway activities was identified by ascertaining differences in quantitative topological influences among the nodes between the two conditions.

The process of angiogenesis is under complex regulation in adult organisms, particularly as it often occurs in an inflammatory post-wound environment. As such, there are many impacting factors that will regulate the generation of new blood vessels which include not only pro-angiogenic growth factors such as vascular endothelial growth factor, but also angiostatic factors. During initial postwound hemostasis, a large initial bolus of platelet factor 4 is released into localized areas of damage before progression of wound healing toward tissue homeostasis. Because of its early presence and high concentration, the angiostatic chemokine platelet factor 4, which can induce endothelial anoikis, can strongly affect angiogenesis. In our work, we explored signaling crosstalk interactions between vascular endothelial growth factor and platelet factor 4 using phosphotyrosine-enriched mass spectrometry methods on human dermal microvascular endothelial cells cultured under conditions facilitating migratory sprouting into collagen gel matrices. We developed new methods to enable mass spectrometry-based phosphorylation analysis of primary cells cultured on collagen gels, and quantified signaling pathways over the first 48 h of treatment with vascular endothelial growth factor in the presence or absence of platelet factor 4. By observing early and late signaling dynamics in tandem with correlation network modeling, we found that platelet factor 4 has significant crosstalk with vascular endothelial growth factor by modulating cell migration and polarization pathways, centered around P38α MAPK, Src family kinases Fyn and Lyn, along with FAK. Interestingly, we found EphA2 correlational topology to strongly involve key migration-related signaling nodes after introduction of platelet factor 4, indicating an influence of the angiostatic factor on this ambiguous but generally angiogenic signal in this complex environment.

Understanding the mechanisms of cell function and drug action is a major endeavor in the pharmaceutical industry. Drug effects are governed by the intrinsic properties of the drug (i.e., selectivity and potency) and the specific signaling transduction network of the host (i.e., normal vs. diseased cells). Here, we describe an unbiased, phosphoproteomic-based approach to identify drug effects by monitoring drug-induced topology alterations. With our proposed method, drug effects are investigated under diverse stimulations of the signaling network. Starting with a generic pathway made of logical gates, we build a cell-type specific map by constraining it to fit 13 key phopshoprotein signals under 55 experimental conditions. Fitting is performed via an Integer Linear Program (ILP) formulation and solution by standard ILP solvers; a procedure that drastically outperforms previous fitting schemes. Then, knowing the cell's topology, we monitor the same key phosphoprotein signals under the presence of drug and we re-optimize the specific map to reveal drug-induced topology alterations. To prove our case, we make a topology for the hepatocytic cell-line HepG2 and we evaluate the effects of 4 drugs: 3 selective inhibitors for the Epidermal Growth Factor Receptor (EGFR) and a non-selective drug. We confirm effects easily predictable from the drugs' main target (i.e., EGFR inhibitors blocks the EGFR pathway) but we also uncover unanticipated effects due to either drug promiscuity or the cell's specific topology. An interesting finding is that the selective EGFR inhibitor Gefitinib inhibits signaling downstream the Interleukin-1alpha (IL1alpha) pathway; an effect that cannot be extracted from binding affinity-based approaches. Our method represents an unbiased approach to identify drug effects on small to medium size pathways which is scalable to larger topologies with any type of signaling interventions (small molecules, RNAi, etc). The method can reveal drug effects on pathways, the cornerstone for identifying mechanisms of drug's efficacy.

The spread of cancer during metastatic disease requires that tumor cells subvert normal regulatory networks governing cell motility to invade surrounding tissues and migrate toward blood and lymphatic vessels. Enabled (Ena)/vasodilator-stimulated phosphoprotein (VASP) proteins regulate cell motility by controlling the geometry of assembling actin networks. Mena, an Ena/VASP protein, is upregulated in the invasive subpopulation of breast cancer cells. In addition, Mena is alternately spliced to produce an invasion isoform, Mena(INV). Here we show that Mena and Mena(INV) promote carcinoma cell motility and invasiveness in vivo and in vitro, and increase lung metastasis. Mena and Mena(INV) potentiate epidermal growth factor (EGF)-induced membrane protrusion and increase the matrix degradation activity of tumor cells. Interestingly, Mena(INV) is significantly more effective than Mena in driving metastases and sensitizing cells to EGF-dependent invasion and protrusion. Upregulation of Mena(INV) could therefore enable tumor cells to invade in response to otherwise benign EGF stimulus levels.

When modeling cell signaling networks, a balance must be struck between mechanistic detail and ease of interpretation. In this paper we apply a fuzzy logic framework to the analysis of a large, systematic dataset describing the dynamics of cell signaling downstream of TNF, EGF, and insulin receptors in human colon carcinoma cells. Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation. We uncover a relationship between MK2 and ERK pathways that might account for the previously identified pro-survival influence of MK2. We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling. More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

Systematic study of cell signaling networks increasingly involves high throughput proteomics, transcriptional profiling, and automated literature mining with the aim of assembling large scale interaction networks. In contrast, functional analysis of cell signaling usually focuses on a much smaller sets of proteins and eschews computation but focuses directly on cellular responses to environment and perturbation. We sought to combine these two traditions by collecting cell response measures on a reasonably large scale and then attempting to infer differences in network topology between two cell types. Human hepatocytes and hepatocellular carcinoma cell lines were exposed to inducers of inflammation, innate immunity, and proliferation in the presence and absence of small molecule drugs, and multiplex biochemical measurement was then performed on intra- and extracellular signaling molecules. We uncovered major differences between primary and transformed hepatocytes with respect to the engagement of toll-like receptor and NF-kappaB-dependent secretion of chemokines and cytokines that prime and attract immune cells. Overall, our results serve as a proof of principle for an approach to network analysis that is systematic, comparative, and biochemically focused. More specifically, our data support the hypothesis that hepatocellular carcinoma cells down-regulate normal inflammatory and immune responses to avoid immune editing.

Proteome balance is safeguarded by the proteostasis network (PN), an intricately regulated network of conserved processes that evolved to maintain native function of the diverse ensemble of protein species, ensuring cellular and organismal health. Proteostasis imbalances and collapse are implicated in a spectrum of human diseases, from neurodegeneration to cancer. The characteristics of PN disease alterations however have not been assessed in a systematic way. Since the chaperome is among the central components of the PN, we focused on the chaperome in our study by utilizing a curated functional ontology of the human chaperome that we connect in a high-confidence physical protein-protein interaction network. Challenged by the lack of a systems-level understanding of proteostasis alterations in the heterogeneous spectrum of human cancers, we assessed gene expression across more than 10,000 patient biopsies covering 22 solid cancers. We derived a novel customized Meta-PCA dimension reduction approach yielding M-scores as quantitative indicators of disease expression changes to condense the complexity of cancer transcriptomics datasets into quantitative functional network topographies. We confirm upregulation of the HSP90 family and also highlight HSP60s, Prefoldins, HSP100s, ER- and mitochondria-specific chaperones as pan-cancer enriched. Our analysis also reveals a surprisingly consistent strong downregulation of small heat shock proteins (sHSPs) and we stratify two cancer groups based on the preferential upregulation of ATP-dependent chaperones. Strikingly, our analyses highlight similarities between stem cell and cancer proteostasis, and diametrically opposed chaperome deregulation between cancers and neurodegenerative diseases. We developed a web-based Proteostasis Profiler tool (Pro2) enabling intuitive analysis and visual exploration of proteostasis disease alterations using gene expression data. Our study showcases a comprehensive profiling of chaperome shifts in human cancers and sets the stage for a systematic global analysis of PN alterations across the human diseasome towards novel hypotheses for therapeutic network re-adjustment in proteostasis disorders.