BackgroundCells process signals using complex and dynamic networks. Studying how this is performed in a context and cell type specific way is essential to understand signaling both in physiological and diseased situations. Context-specific medium/high throughput proteomic data measured upon perturbation is now relatively easy to obtain but formalisms that can take advantage of these features to build models of signaling are still comparatively scarce.ResultsHere we present CellNOptR, an open-source R software package for building predictive logic models of signaling networks by training networks derived from prior knowledge to signaling (typically phosphoproteomic) data. CellNOptR features different logic formalisms, from Boolean models to differential equations, in a common framework. These different logic model representations accommodate state and time values with increasing levels of detail. We provide in addition an interface via Cytoscape (CytoCopteR) to facilitate use and integration with Cytoscape network-based capabilities.ConclusionsModels generated with this pipeline have two key features. First, they are constrained by prior knowledge about the network but trained to data. They are therefore context and cell line specific, which results in enhanced predictive and mechanistic insights. Second, they can be built using different logic formalisms depending on the richness of the available data. Models built with CellNOptR are useful tools to understand how signals are processed by cells and how this is altered in disease. They can be used to predict the effect of perturbations (individual or in combinations), and potentially to engineer therapies that have differential effects/side effects depending on the cell type or context.

Our understanding of physiology and disease is hampered by the difficulty of measuring the circuitry and plasticity of signaling networks that regulate cell biology, and how these relate to phenotypes. Here, using mass spectrometry-based phosphoproteomics, we systematically characterized the topology of a network comprising the PI3K/Akt/mTOR and MEK/ERK signaling axes and confirmed its biological relevance by assessing its dynamics upon EGF and IGF1 stimulation. Measuring the activity of this network in models of acquired drug resistance revealed that cells chronically treated with PI3K or mTORC1/2 inhibitors differed in the way their networks were remodeled. Unexpectedly, we also observed a degree of heterogeneity in the network state between cells resistant to the same inhibitor, indicating that even identical and carefully controlled experimental conditions can give rise to the evolution of distinct kinase network statuses. These data suggest that the initial conditions of the system do not necessarily determine the mechanism by which cancer cells become resistant to PI3K/mTOR targeted therapies. The patterns of signaling network activity observed in the resistant cells mirrored the patterns of response to several drug combination treatments, suggesting that the activity of the defined signaling network truly reflected the evolved phenotypic diversity.signaling network | phosphoproteomics | systems biology | PI3K | MEK

Motivation: High-throughput screens (HTS) by RNAi or small molecules are among the most promising tools in functional genomics. They enable researchers to observe detailed reactions to experimental perturbations on a genome-wide scale. While there is a core set of computational approaches used in many publications to analyze these data, a specialized software combining them and making them easily accessible has so far been missing.Results: Here we describe , a flexible software to build integrated analysis pipelines for HTS data that contains over-representation analysis, gene set enrichment analysis, comparative gene set analysis and rich sub-network identification. interfaces with commonly used pre-processing packages for HTS data and presents its results as HTML pages and network plots.Availability: Our software is written in the R language and freely available via the Bioconductor project at http://www.bioconductor.org.Contact: florian.markowetz@cancer.org.uk

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.

Despite the current wealth of high-throughput data, our understanding of signal transduction is still incomplete. Mathematical modeling can be a tool to gain an insight into such processes. Detailed biochemical modeling provides deep understanding, but does not scale well above relatively a few proteins. In contrast, logic modeling can be used where the biochemical knowledge of the system is sparse and, because it is parameter free (or, at most, uses relatively a few parameters), it scales well to large networks that can be derived by manual curation or retrieved from public databases. Here, we present an overview of logic modeling formalisms in the context of training logic models to data, and specifically the different approaches to modeling qualitative to quantitative data (state) and dynamics (time) of signal transduction. We use a toy model of signal transduction to illustrate how different logic formalisms (Boolean, fuzzy logic and differential equations) treat state and time. Different formalisms allow for different features of the data to be captured, at the cost of extra requirements in terms of computational power and data quality and quantity. Through this demonstration, the assumptions behind each formalism are discussed, as well as their advantages and disadvantages and possible future developments.

Cellular communication and information processing is performed by complex, dynamic, and context specific signaling networks. Mathematical modeling is a very useful tool to make sense of this complexity. Building a model relies on two main ingredients: data and an adequate model formalism. In the case of signaling networks, we build mainly upon data at the proteome level, in particular about the phosphorylation of proteins. In this chapter we review recent developments in both data acquisition and computational analysis. We describe two approaches, antibody based technologies and mass spectrometry (MS), along with their main features and limitations. We then go on to describe some model formalisms that have been applied to such high-throughput phospho-proteomics data sets. We consider a variety of formalisms from clustering and data mining approaches to differential equation-based mechanistic models, rule-based, and logic based models, and on through Bayesian network inference and linear regressions.

Advances in mass spectrometry have made the quantitative measurement of proteins across multiple samples a reality, allowing for the study of complex biological systems such as the metabolic syndrome. Although the deregulation of lipid metabolism and increased hepatic storage of triacylglycerides are known to play a part in the onset of the metabolic syndrome, its molecular basis and dependency on dietary and genotypic factors are poorly characterized. Here, we used an experimental design with two different mouse strains and dietary and metabolic perturbations to generate a compendium of quantitative proteome data using three mass spectrometric techniques. The data reproduce known properties of the metabolic system and indicate differential molecular adaptation of the two mouse strains to perturbations, contributing to a better understanding of the metabolic syndrome. We show that high-quality, high-throughput proteomic data sets provide an unbiased broad overview of the behavior of complex systems after perturbation.

SUMMARYAdvances in mass spectrometry have made the quantitative measurement of proteins across multiple samples a reality, allowing for the study of complex biological systems such as the metabolic syndrome. Although deregulation of lipid metabolism and increased hepatic storage of triacylglycerides are known to play a part in the onset of the metabolic syndrome, its molecular basis and dependency on dietary and genotypic factors are poorly characterized. Here, we used a rich experimental design with two different mouse strains, dietary and metabolic perturbations to generate a compendium of quantitative proteome data, using three mass spectrometry strategies.The data recapitulates known properties of the metabolic system and indicate differential molecular adaptation of the two mouse strains to perturbations, contributing to a better understanding of the metabolic syndrome. We show that high-quality, high-throughput proteomic datasets provide an unbiased broad overview of complex systems upon perturbation.