Simplistic pathways or complex networks?

Analyzing the Impact of Diesel Exhaust Particles on Lung Fibrosis Using Dual PCR Array and Proteomics: YWHAZ Signaling

Air pollutants are associated with exacerbations of asthma, chronic bronchitis, and airway inflammation. Diesel exhaust particles (DEPs) can induce and worsen lung diseases. However, there are insufficient data to guide polymerase chain reaction (PCR) array proteomics studies regarding the impacts of DEPs on respiratory diseases. This study was performed to identify genes and proteins expressed in normal human bronchial epithelial (NHBE) cells. MicroRNAs (miRNAs) and proteins expressed in NHBE cells exposed to DEPs at 1 μg/cm2 for 8 h and 24 h were identified using PCR array analysis and 2D PAGE/LC-MS/MS, respectively. YWHAZ gene expression was estimated using PCR, immunoblotting, and immunohistochemical analyses. Genes discovered through an overlap analysis were validated in DEP-exposed mice. Proteomics approaches showed that exposing NHBE cells to DEPs led to changes in 32 protein spots. A transcriptomics PCR array analysis showed that 6 of 84 miRNAs were downregulated in the DEP exposure groups compared to controls. The mRNA and protein expression levels of YWHAZ, β-catenin, vimentin, and TGF-β were increased in DEP-treated NHBE cells and DEP-exposed mice. Lung fibrosis was increased in mice exposed to DEPs. Our combined PCR array-omics analysis demonstrated that DEPs can induce airway inflammation and lead to lung fibrosis through changes in the expression levels of YWHAZ, β-catenin, vimentin, and TGF-β. These findings suggest that dual approaches can help to identify biomarkers and therapeutic targets involved in pollutant-related respiratory diseases.

Exploring absent protein function in yeast: assaying post translational modification and human genetic variation

Yeast is a valuable eukaryotic model organism that has evolved many processes conserved up to humans, yet many protein functions, including certain DNA and protein modifications, are absent. It is this absence of protein function that is fundamental to approaches using yeast as an in vivo test system to investigate human proteins. Functionality of the heterologous expressed proteins is connected to a quantitative, selectable phenotype, enabling the systematic analyses of mechanisms and specificity of DNA modification, post-translational protein modifications as well as the impact of annotated cancer mutations and coding variation on protein activity and interaction. Through continuous improvements of yeast screening systems, this is increasingly carried out on a global scale using deep mutational scanning approaches. Here we discuss the applicability of yeast systems to investigate absent human protein function with a specific focus on the impact of protein variation on protein-protein interaction modulation.

Pathway-focused PCR array profiling of CAL-27 cell with over-expressed ZNF750

Zinc-finger protein 750 (ZNF750) is the potential anti-cancer gene in oral squamous cell carcinoma (OSCC). The present study was to investigate the expression changes of ZNF750 in OSCC tissue and to reveal the induction of altered mRNA expression profiles caused by over-expressed ZNF750 in CAL-27 cell. The expression level of ZNF750 in tissue specimens from OSCC patients was detected by immunohistochemistry. Gene expression profiling was performed using Human Signal Transduction PathwayFinder RT² Profiler™ PCR Array. The expression changes of 84 key genes representing 10 signal transduction pathways in human following over-expressed ZNF750 in CAL-27 cell was examined. The expression of ZNF750 protein was reduced in OSCC tissues. The R² PCR Array analysis revealed that 39 of the 84 examined genes that changed at least a two-fold between control and ZNF750 groups. These genes related to oxidative stress, WNT, JAK/STAT, TGFβ, NF-kappaB (NFκB), p53, Notch, Hedgehog, PPAR and Hypoxia signaling. ZNF750 could inhibit the candidate genes ATF4, SQSTM1, HMOX1, CCND1, TNF-alpha, TNFSF10 and FOSL1 but activate CDKN1A and EMP1. Our studies suggest that ZNF750 can regulate signaling pathways that related to proliferation, cell cycle, inflammation and oxidative stress in CAL-27 cell.

Phosphoproteomics-Based Profiling of Kinase Activities in Cancer Cells

Cellular signaling, predominantly mediated by phosphorylation through protein kinases, is found to be deregulated in most cancers. Accordingly, protein kinases have been subject to intense investigations in cancer research, to understand their role in oncogenesis and to discover new therapeutic targets. Despite great advances, an understanding of kinase dysfunction in cancer is far from complete.A powerful tool to investigate phosphorylation is mass-spectrometry (MS)-based phosphoproteomics, which enables the identification of thousands of phosphorylated peptides in a single experiment. Since every phosphorylation event results from the activity of a protein kinase, high-coverage phosphoproteomics data should indirectly contain comprehensive information about the activity of protein kinases.In this chapter, we discuss the use of computational methods to predict kinase activity scores from MS-based phosphoproteomics data. We start with a short explanation of the fundamental features of the phosphoproteomics data acquisition process from the perspective of the computational analysis. Next, we briefly review the existing databases with experimentally verified kinase-substrate relationships and present a set of bioinformatic tools to discover novel kinase targets. We then introduce different methods to infer kinase activities from phosphoproteomics data and these kinase-substrate relationships. We illustrate their application with a detailed protocol of one of the methods, KSEA (Kinase Substrate Enrichment Analysis). This method is implemented in Python within the framework of the open-source Kinase Activity Toolbox (kinact), which is freely available at http://github.com/saezlab/kinact/ .

A parallel metaheuristic for large mixed-integer dynamic optimization problems, with applications in computational biology

Background

We consider a general class of global optimization problems dealing with nonlinear dynamic models. Although this class is relevant to many areas of science and engineering, here we are interested in applying this framework to the reverse engineering problem in computational systems biology, which yields very large mixed-integer dynamic optimization (MIDO) problems. In particular, we consider the framework of logic-based ordinary differential equations (ODEs).

Methods

We present saCeSS2, a parallel method for the solution of this class of problems. This method is based on an parallel cooperative scatter search metaheuristic, with new mechanisms of self-adaptation and specific extensions to handle large mixed-integer problems. We have paid special attention to the avoidance of convergence stagnation using adaptive cooperation strategies tailored to this class of problems.

Results

We illustrate its performance with a set of three very challenging case studies from the domain of dynamic modelling of cell signaling. The simpler case study considers a synthetic signaling pathway and has 84 continuous and 34 binary decision variables. A second case study considers the dynamic modeling of signaling in liver cancer using high-throughput data, and has 135 continuous and 109 binaries decision variables. The third case study is an extremely difficult problem related with breast cancer, involving 690 continuous and 138 binary decision variables. We report computational results obtained in different infrastructures, including a local cluster, a large supercomputer and a public cloud platform. Interestingly, the results show how the cooperation of individual parallel searches modifies the systemic properties of the sequential algorithm, achieving superlinear speedups compared to an individual search (e.g. speedups of 15 with 10 cores), and significantly improving (above a 60%) the performance with respect to a non-cooperative parallel scheme. The scalability of the method is also good (tests were performed using up to 300 cores).

Conclusions

These results demonstrate that saCeSS2 can be used to successfully reverse engineer large dynamic models of complex biological pathways. Further, these results open up new possibilities for other MIDO-based large-scale applications in the life sciences such as metabolic engineering, synthetic biology, drug scheduling.

Kinase activity ranking using phosphoproteomics data (KARP) quantifies the contribution of protein kinases to the regulation of cell viability

Cell survival is regulated by a signaling network driven by the activity of protein kinases; however, determining the contribution that each kinase in the network makes to such regulation remains challenging. Here, we report a computational approach that uses mass spectrometry-based phosphoproteomics data to rank protein kinases based on their contribution to cell regulation. We found that the scores returned by this algorithm, which we have termed kinase activity ranking using phosphoproteomics data (KARP), were a quantitative measure of the contribution that individual kinases make to the signaling output. Application of KARP to the analysis of eight hematological cell lines revealed that cyclin-dependent kinase (CDK) 1/2, casein kinase (CK) 2, extracellular signal-related kinase (ERK), and p21-activated kinase (PAK) were the most frequently highly ranked kinases in these cell models. The patterns of kinase activation were cell-line specific yet showed a significant association with cell viability as a function of kinase inhibitor treatment. Thus, our study exemplifies KARP as an untargeted approach to empirically and systematically identify regulatory kinases within signaling networks.

Knowledge-guided fuzzy logic modeling to infer cellular signaling networks from proteomic data

Modeling of signaling pathways is crucial for understanding and predicting cellular responses to drug treatments. However, canonical signaling pathways curated from literature are seldom context-specific and thus can hardly predict cell type-specific response to external perturbations; purely data-driven methods also have drawbacks such as limited biological interpretability. Therefore, hybrid methods that can integrate prior knowledge and real data for network inference are highly desirable. In this paper, we propose a knowledge-guided fuzzy logic network model to infer signaling pathways by exploiting both prior knowledge and time-series data. In particular, the dynamic time warping algorithm is employed to measure the goodness of fit between experimental and predicted data, so that our method can model temporally-ordered experimental observations. We evaluated the proposed method on a synthetic dataset and two real phosphoproteomic datasets. The experimental results demonstrate that our model can uncover drug-induced alterations in signaling pathways in cancer cells. Compared with existing hybrid models, our method can model feedback loops so that the dynamical mechanisms of signaling networks can be uncovered from time-series data. By calibrating generic models of signaling pathways against real data, our method supports precise predictions of context-specific anticancer drug effects, which is an important step towards precision medicine.

Large-scale models of signal propagation in human cells derived from discovery phosphoproteomic data

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.

Phosphoproteomics can offer significant insight into cell signalling and how signalling is modified in response to perturbations. Here the authors develop a new tool for the analysis of high-content phosphoproteomics in the context of kinase/phosphatase-substrate knowledge, which is used to train logic models.

Empirical inference of circuitry and plasticity in a kinase signaling network

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.

Systems biology and the integration of mechanistic explanation and mathematical explanation

The paper discusses how systems biology is working toward complex accounts that integrate explanation in terms of mechanisms and explanation by mathematical models-which some philosophers have viewed as rival models of explanation. Systems biology is an integrative approach, and it strongly relies on mathematical modeling. Philosophical accounts of mechanisms capture integrative in the sense of multilevel and multifield explanations, yet accounts of mechanistic explanation (as the analysis of a whole in terms of its structural parts and their qualitative interactions) have failed to address how a mathematical model could contribute to such explanations. I discuss how mathematical equations can be explanatorily relevant. Several cases from systems biology are discussed to illustrate the interplay between mechanistic research and mathematical modeling, and I point to questions about qualitative phenomena (rather than the explanation of quantitative details), where quantitative models are still indispensable to the explanation. Systems biology shows that a broader philosophical conception of mechanisms is needed, which takes into account functional-dynamical aspects, interaction in complex networks with feedback loops, system-wide functional properties such as distributed functionality and robustness, and a mechanism's ability to respond to perturbations (beyond its actual operation). I offer general conclusions for philosophical accounts of explanation.

Ovarian cancer evolution through stochastic genome alterations: defining the genomic role in ovarian cancer

Ovarian cancer is the fifth leading cause of death among women worldwide. Characterized by complex etiology and multi-level heterogeneity, its origins are not well understood. Intense research efforts over the last decade have furthered our knowledge by identifying multiple risk factors that are associated with the disease. However, it is still unclear how genetic heterogeneity contributes to tumor formation, and more specifically, how genome-level heterogeneity acts as the key driving force of cancer evolution. Most current genomic approaches are based on 'average molecular profiling.' While effective for data generation, they often fail to effectively address the issue of high level heterogeneity because they mask variation that exists in a cell population. In this synthesis, we hypothesize that genome-mediated cancer evolution can effectively explain diverse factors that contribute to ovarian cancer. In particular, the key contribution of genome replacement can be observed during major transitions of ovarian cancer evolution including cellular immortalization, transformation, and malignancy. First, we briefly review major updates in the literature, and illustrate how current gene-mediated research will offer limited insight into cellular heterogeneity and ovarian cancer evolution. We next explain a holistic framework for genome-based ovarian cancer evolution and apply it to understand the genomic dynamics of a syngeneic ovarian cancer mouse model. Finally, we employ single cell assays to further test our hypothesis, discuss some predictions, and report some recent findings.

Machine learning prediction of cancer cell sensitivity to drugs based on genomic and chemical properties

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.

Mapping the human phosphatome on growth pathways

Phosphatases control cell growth by a variety of mechanisms. A novel strategy is presented that combines multiparametric analysis of cell perturbations with logic modeling to achieve a detailed mapping of human phosphatase function on growth pathways.

siRNA-mediated downregulation of 298 phosphatase and phosphatase-related genes coupled to automated microscopy was used to characterize their impact on key growth pathways.

In parallel, a literature-derived signed directed network was derived and optimized by training with experimental data.

The resulting logic-based growth model was used to infer the cell state upon perturbation of each signaling node and compare it with the profiles obtained upon phosphatase perturbation.

Mapping of 67% of the protein phosphatase onto the growth model shows that phosphatases are key modulators of growth pathways and affect cell-cycle progression.

This novel approach is general and enables to efficiently map proteins onto complex pathways.

Large-scale siRNA screenings allow linking the function of poorly characterized genes to phenotypic readouts. According to this strategy, genes are associated with a function of interest if the alteration of their expression perturbs the phenotypic readouts. However, given the intricacy of the cell regulatory network, the mapping procedure is low resolution and the resulting models provide little mechanistic insights. We have developed a new strategy that combines multiparametric analysis of cell perturbation with logic modeling to achieve a more detailed functional mapping of human genes onto complex pathways. A literature-derived optimized model is used to infer the cell activation state following upregulation or downregulation of the model entities. By matching this signature with the experimental profile obtained in the high-throughput siRNA screening it is possible to infer the target of each protein, thus defining its ‘entry point' in the network. By this novel approach, 41 phosphatases that affect key growth pathways were identified and mapped onto a human epithelial cell-specific growth model, thus providing insights into the mechanisms underlying their function.

Construction of cell type-specific logic models of signaling networks using CellNOpt

Mathematical models are useful tools for understanding protein signaling networks because they provide an integrated view of pharmacological and toxicological processes at the molecular level. Here we describe an approach previously introduced based on logic modeling to generate cell-specific, mechanistic and predictive models of signal transduction. Models are derived from a network encoding prior knowledge that is trained to signaling data, and can be either binary (based on Boolean logic) or quantitative (using a recently developed formalism, constrained fuzzy logic). The approach is implemented in the freely available tool CellNetOptimizer (CellNOpt). We explain the process CellNOpt uses to train a prior knowledge network to data and illustrate its application with a toy example as well as a realistic case describing signaling networks in the HepG2 liver cancer cell line.

PROTEINCHALLENGE: crowd sourcing in proteomics analysis and software development

In large-scale proteomics studies there is a temptation, after months of experimental work, to plug resulting data into a convenient-if poorly implemented-set of tools, which may neither do the data justice nor help answer the scientific question. In this paper we have captured key concerns, including arguments for community-wide open source software development and "big data" compatible solutions for the future. For the meantime, we have laid out ten top tips for data processing. With these at hand, a first large-scale proteomics analysis hopefully becomes less daunting to navigate. However there is clearly a real need for robust tools, standard operating procedures and general acceptance of best practises. Thus we submit to the proteomics community a call for a community-wide open set of proteomics analysis challenges--PROTEINCHALLENGE--that directly target and compare data analysis workflows, with the aim of setting a community-driven gold standard for data handling, reporting and sharing.

CellNOptR: a flexible toolkit to train protein signaling networks to data using multiple logic formalisms

Background

Cells 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.

Results

Here 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.

Conclusions

Models 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.

State-time spectrum of signal transduction logic models

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.

A directed protein interaction network for investigating intracellular signal transduction

Cellular signal transduction is a complex process involving protein-protein interactions (PPIs) that transmit information. For example, signals from the plasma membrane may be transduced to transcription factors to regulate gene expression. To obtain a global view of cellular signaling and to predict potential signal modulators, we searched for protein interaction partners of more than 450 signaling-related proteins by means of automated yeast two-hybrid interaction mating. The resulting PPI network connected 1126 proteins through 2626 PPIs. After expansion of this interaction map with publicly available PPI data, we generated a directed network resembling the signal transduction flow between proteins with a naïve Bayesian classifier. We exploited information on the shortest PPI paths from membrane receptors to transcription factors to predict input and output relationships between interacting proteins. Integration of directed PPI with time-resolved protein phosphorylation data revealed network structures that dynamically conveyed information from the activated epidermal growth factor and extracellular signal-regulated kinase (EGF/ERK) signaling cascade to directly associated proteins and more distant proteins in the network. From the model network, we predicted 18 previously unknown modulators of EGF/ERK signaling, which we validated in mammalian cell-based assays. This generic experimental and computational approach provides a framework for elucidating causal connections between signaling proteins and facilitates the identification of proteins that modulate the flow of information in signaling networks.

NMR analysis of a stress response metabolic signaling network

We previously hypothesized that Staphylococcus epidermidis senses a diverse set of environmental and nutritional factors associated with biofilm formation through a modulation in the activity of the tricarboxylic acid (TCA) cycle. Herein, we report our further investigation of the impact of additional environmental stress factors on TCA cycle activity and provide a detailed description of our NMR methodology. S. epidermidis wild-type strain 1457 was treated with stressors that are associated with biofilm formation, a sublethal dose of tetracycline, 5% NaCl, 2% glucose, and autoinducer-2 (AI-2). As controls and to integrate our current data with our previous study, 4% ethanol stress and iron-limitation were also used. Consistent with our prior observations, the effect of many environmental stress factors on the S. epidermidis metabolome was essentially identical to the effect of TCA cycle inactivation in the aconitase mutant strain 1457-acnA::tetM. A detailed quantitative analysis of metabolite concentration changes using 2D (1)H-(13)C HSQC and (1)H-(1)H TOCSY spectra identified a network of 37 metabolites uniformly affected by the stressors and TCA cycle inactivation. We postulate that the TCA cycle acts as the central pathway in a metabolic signaling network.

Defining pluripotent stem cells through quantitative proteomic analysis

Embryonic stem cells (ESCs) are at the center stage of intense research, inspired by their potential to give rise to all cell types of the adult individual. This property makes ESCs suitable candidates for generating specialized cells to replace damaged tissue lost after injury or disease. However, such clinical applications require a detailed insight of the molecular mechanisms underlying the self-renewal, expansion and differentiation of stem cells. This has gained further relevance since the introduction of induced pluripotent stem cells (iPSCs), which are functionally very similar to ESCs. The key property that iPSCs can be derived from somatic cells lifts some of the major ethical issues related to the need for embryos to generate ESCs. Yet, this has only increased the need to define the similarity of iPSCs and ESCs at the molecular level, both before and after they are induced to differentiate. In this article, we describe the proteomic approaches that have been used to characterize ESCs with regard to self-renewal and differentiation, with an emphasis on signaling cascades and histone modifications. We take this as a lead to discuss how quantitative proteomics can be deployed to study reprogramming and iPSC identity. In addition, we discuss how emerging proteomic technologies can become a useful tool to monitor the (de)differentiation status of ESCs and iPSCs.

Spatial phosphoprotein profiling reveals a compartmentalized extracellular signal-regulated kinase switch governing neurite growth and retraction

Brain development and spinal cord regeneration require neurite sprouting and growth cone navigation in response to extension and collapsing factors present in the extracellular environment. These external guidance cues control neurite growth cone extension and retraction processes through intracellular protein phosphorylation of numerous cytoskeletal, adhesion, and polarity complex signaling proteins. However, the complex kinase/substrate signaling networks that mediate neuritogenesis have not been investigated. Here, we compare the neurite phosphoproteome under growth and retraction conditions using neurite purification methodology combined with mass spectrometry. More than 4000 non-redundant phosphorylation sites from 1883 proteins have been annotated and mapped to signaling pathways that control kinase/phosphatase networks, cytoskeleton remodeling, and axon/dendrite specification. Comprehensive informatics and functional studies revealed a compartmentalized ERK activation/deactivation cytoskeletal switch that governs neurite growth and retraction, respectively. Our findings provide the first system-wide analysis of the phosphoprotein signaling networks that enable neurite growth and retraction and reveal an important molecular switch that governs neuritogenesis.

Training signaling pathway maps to biochemical data with constrained fuzzy logic: quantitative analysis of liver cell responses to inflammatory stimuli

Predictive understanding of cell signaling network operation based on general prior knowledge but consistent with empirical data in a specific environmental context is a current challenge in computational biology. Recent work has demonstrated that Boolean logic can be used to create context-specific network models by training proteomic pathway maps to dedicated biochemical data; however, the Boolean formalism is restricted to characterizing protein species as either fully active or inactive. To advance beyond this limitation, we propose a novel form of fuzzy logic sufficiently flexible to model quantitative data but also sufficiently simple to efficiently construct models by training pathway maps on dedicated experimental measurements. Our new approach, termed constrained fuzzy logic (cFL), converts a prior knowledge network (obtained from literature or interactome databases) into a computable model that describes graded values of protein activation across multiple pathways. We train a cFL-converted network to experimental data describing hepatocytic protein activation by inflammatory cytokines and demonstrate the application of the resultant trained models for three important purposes: (a) generating experimentally testable biological hypotheses concerning pathway crosstalk, (b) establishing capability for quantitative prediction of protein activity, and (c) prediction and understanding of the cytokine release phenotypic response. Our methodology systematically and quantitatively trains a protein pathway map summarizing curated literature to context-specific biochemical data. This process generates a computable model yielding successful prediction of new test data and offering biological insight into complex datasets that are difficult to fully analyze by intuition alone.

Biological signalling activity measurements using mass spectrometry

MS (mass spectrometry) techniques are rapidly evolving to high levels of performance and robustness. This is allowing the application of these methods to the interrogation of signalling networks with unprecedented depth and accuracy. In the present review we discuss how MS-based multiplex quantification of kinase activities and phosphoproteomics provide complementary means to assess biological signalling activity. In addition, we discuss how a wider application of these analytical concepts to quantify kinase signalling will result in a more comprehensive understanding of normal and disease biology at the system level.

Phospho.ELM: a database of phosphorylation sites--update 2011

The Phospho.ELM resource (http://phospho.elm.eu.org) is a relational database designed to store in vivo and in vitro phosphorylation data extracted from the scientific literature and phosphoproteomic analyses. The resource has been actively developed for more than 7 years and currently comprises 42 574 serine, threonine and tyrosine non-redundant phosphorylation sites. Several new features have been implemented, such as structural disorder/order and accessibility information and a conservation score. Additionally, the conservation of the phosphosites can now be visualized directly on the multiple sequence alignment used for the score calculation. Finally, special emphasis has been put on linking to external resources such as interaction networks and other databases.

Rapid cytosolic delivery of luminescent nanocrystals in live cells with endosome-disrupting polymer colloids

Luminescent nanocrystals hold great potential for bioimaging because of their exceptional optical properties, but their use in live cells has been limited. When nanocrystals enter live cells, they are taken up in vesicles. This vesicular sequestration is persistent and precludes nanocrystals from reaching intracellular targets. Here, we describe a unique, cationic core-shell polymer colloid that translocates nanocrystals to the cytosol by disrupting endosomal membranes via a low-pH triggered mechanism. Confocal fluorescence microscopy and flow cytometry indicate that picomolar concentrations of quantum dots are sufficient for cytosolic labeling, with the process occurring within a few hours of incubation. We anticipate a host of advanced applications arising from efficient cytosolic delivery of nanocrystal imaging probes: from single particle tracking experiments to monitoring protein-protein interactions in live cells for extended periods.