Generating synthetic signaling networks for in silico modeling studies.

Predictive models of signaling pathways have proven to be difficult to develop. Traditional approaches to developing mechanistic models rely on collecting experimental data and fitting a single model to that data. This approach works for simple systems but has proven unreliable for complex systems such as biological signaling networks. Thus, there is a need to develop new approaches to create predictive mechanistic models of complex systems. To meet this need, we developed a method for generating artificial signaling networks that were reasonably realistic and thus could be treated as ground truth models. These synthetic models could then be used to generate synthetic data for developing and testing algorithms designed to recover the underlying network topology and associated parameters. We defined the reaction degree and reaction distance to measure the topology of reaction networks, especially to consider enzymes. To determine whether our generated signaling networks displayed meaningful behavior, we compared them with signaling networks from the BioModels Database. This comparison indicated that our generated signaling networks had high topological similarities with BioModels signaling networks with respect to the reaction degree and distance distributions. In addition, our synthetic signaling networks had similar behavioral dynamics with respect to both steady states and oscillations, suggesting that our method generated synthetic signaling networks comparable with BioModels and thus could be useful for building network evaluation tools.

Signal integration and integral feedback control with biochemical reaction networks.

Biochemical reaction networks perform a variety of signal processing functions, one of which is computing the integrals of signal values. This is often used in integral feedback control, where it enables a system's output to respond to changing inputs, but to then return exactly back to some pre-determined setpoint value afterward. To gain a deeper understanding of how biochemical networks are able to both integrate signals and perform integral feedback control, we investigated these abilities for several simple reaction networks. We found imperfect overlap between these categories, with some networks able to perform both tasks, some able to perform integration but not integral feedback control, and some the other way around. Nevertheless, networks that could either integrate or perform integral feedback control shared key elements. In particular, they included a chemical species that was neutrally stable in the open loop system (no feedback), meaning that this species does not have a unique stable steady-state concentration. Neutral stability could arise from zeroth order decay reactions, binding to a partner that was produced at a constant rate (which occurs in antithetic control), or through a long chain of covalent cycles. Mathematically, it arose from rate equations for the reaction network that were underdetermined when evaluated at steady-state.

Design patterns of biological cells.

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

Design Patterns of Biological Cells.

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology in order to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

A Phosphoproteomics Data Resource for Systems-level Modeling of Kinase Signaling Networks.

Building mechanistic models of kinase-driven signaling pathways requires quantitative measurements of protein phosphorylation across physiologically relevant conditions, but this is rarely done because of the insensitivity of traditional technologies. By using a multiplexed deep phosphoproteome profiling workflow, we were able to generate a deep phosphoproteomics dataset of the EGFR-MAPK pathway in non-transformed MCF10A cells across physiological ligand concentrations with a time resolution of <12 min and in the presence and absence of multiple kinase inhibitors. An improved phosphosite mapping technique allowed us to reliably identify >46,000 phosphorylation sites on >6600 proteins, of which >4500 sites from 2110 proteins displayed a >2-fold increase in phosphorylation in response to EGF. This data was then placed into a cellular context by linking it to 15 previously published protein databases. We found that our results were consistent with much, but not all previously reported data regarding the activation and negative feedback phosphorylation of core EGFR-ERK pathway proteins. We also found that EGFR signaling is biphasic with substrates downstream of RAS/MAPK activation showing a maximum response at <3ng/ml EGF while direct substrates, such as HGS and STAT5B, showing no saturation. We found that RAS activation is mediated by at least 3 parallel pathways, two of which depend on PTPN11. There appears to be an approximately 4-minute delay in pathway activation at the step between RAS and RAF, but subsequent pathway phosphorylation was extremely rapid. Approximately 80 proteins showed a >2-fold increase in phosphorylation across all experiments and these proteins had a significantly higher median number of phosphorylation sites (~18) relative to total cellular phosphoproteins (~4). Over 60% of EGF-stimulated phosphoproteins were downstream of MAPK and included mediators of cellular processes such as gene transcription, transport, signal transduction and cytoskeletal arrangement. Their phosphorylation was either linear with respect to MAPK activation or biphasic, corresponding to the biphasic signaling seen at the level of the EGFR. This deep, integrated phosphoproteomics data resource should be useful in building mechanistic models of EGFR and MAPK signaling and for understanding how downstream responses are regulated.

Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways.

Post-translational modifications (PTMs) are key regulatory mechanisms that can control protein function. Of these, phosphorylation is the most common and widely studied. Because of its importance in regulating cell signaling, precise and accurate measurements of protein phosphorylation across wide dynamic ranges are crucial to understanding how signaling pathways function. Although immunological assays are commonly used to detect phosphoproteins, their lack of sensitivity, specificity, and selectivity often make them unreliable for quantitative measurements of complex biological samples. Recent advances in Mass Spectrometry (MS)-based targeted proteomics have made it a more useful approach than immunoassays for studying the dynamics of protein phosphorylation. Selected reaction monitoring (SRM)-also known as multiple reaction monitoring (MRM)-and parallel reaction monitoring (PRM) can quantify relative and absolute abundances of protein phosphorylation in multiplexed fashions targeting specific pathways. In addition, the refinement of these tools by enrichment and fractionation strategies has improved measurement of phosphorylation of low-abundance proteins. The quantitative data generated are particularly useful for building and parameterizing mathematical models of complex phospho-signaling pathways. Potentially, these models can provide a framework for linking analytical measurements of clinical samples to better diagnosis and treatment of disease.

Addressing the genetic/nongenetic duality in cancer with systems biology.

The dogma that cancer is a genetic disease is being questioned. Recent findings suggest that genetic/nongenetic duality is necessary for cancer progression. A think tank organized by the Shraman Foundation's Institute for Theoretical Biology compiled key challenges and opportunities that theoreticians, experimentalists, and clinicians can explore from a systems biology perspective to provide a better understanding of the disease as well as help discover new treatment options and therapeutic strategies.

SBbadger: biochemical reaction networks with definable degree distributions.

MOTIVATION: An essential step in developing computational tools for the inference, optimization and simulation of biochemical reaction networks is gauging tool performance against earlier efforts using an appropriate set of benchmarks. General strategies for the assembly of benchmark models include collection from the literature, creation via subnetwork extraction and de novo generation. However, with respect to biochemical reaction networks, these approaches and their associated tools are either poorly suited to generate models that reflect the wide range of properties found in natural biochemical networks or to do so in numbers that enable rigorous statistical analysis. RESULTS: In this work, we present SBbadger, a python-based software tool for the generation of synthetic biochemical reaction or metabolic networks with user-defined degree distributions, multiple available kinetic formalisms and a host of other definable properties. SBbadger thus enables the creation of benchmark model sets that reflect properties of biological systems and generate the kinetics and model structures typically targeted by computational analysis and inference software. Here, we detail the computational and algorithmic workflow of SBbadger, demonstrate its performance under various settings, provide sample outputs and compare it to currently available biochemical reaction network generation software. AVAILABILITY AND IMPLEMENTATION: SBbadger is implemented in Python and is freely available at https://github.com/sys-bio/SBbadger and via PyPI at https://pypi.org/project/SBbadger/. Documentation can be found at https://SBbadger.readthedocs.io. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamics and Sensitivity of Signaling Pathways.

Purpose of Review: Signaling pathways serve to communicate information about extracellular conditions into the cell, to both the nucleus and cytoplasmic processes to control cell responses. Genetic mutations in signaling network components are frequently associated with cancer and can result in cells acquiring an ability to divide and grow uncontrollably. Because signaling pathways play such a significant role in cancer initiation and advancement, their constituent proteins are attractive therapeutic targets. In this review, we discuss how signaling pathway modeling can assist with identifying effective drugs for treating diseases, such as cancer. An achievement that would facilitate the use of such models is their ability to identify controlling biochemical parameters in signaling pathways, such as molecular abundances and chemical reaction rates, because this would help determine effective points of attack by therapeutics. Recent Findings: We summarize the current state of understanding the sensitivity of phosphorylation cycles with and without sequestration. We also describe some basic properties of regulatory motifs including feedback and feedforward regulation. Summary: Although much recent work has focused on understanding the dynamics and particularly the sensitivity of signaling networks in eukaryotic systems, there is still an urgent need to build more scalable models of signaling networks that can appropriately represent their complexity across different cell types and tumors.

Facile One-Pot Nanoproteomics for Label-Free Proteome Profiling of 50-1000 Mammalian Cells.

Recent advances in sample preparation enable label-free mass spectrometry (MS)-based proteome profiling of small numbers of mammalian cells. However, specific devices are often required to downscale sample processing volume from the standard 50-200 µL to sub-µL for effective nanoproteomics, which greatly impedes the implementation of current nanoproteomics methods by the proteomics research community. Herein, we report a facile one-pot nanoproteomics method termed SOPs-MS (surfactant-assisted one-pot sample processing at the standard volume coupled with MS) for convenient robust proteome profiling of 50-1000 mammalian cells. Building upon our recent development of SOPs-MS for label-free single-cell proteomics at a low µL volume, we have systematically evaluated its processing volume at 10-200 µL using 100 human cells. The processing volume of 50 µL that is in the range of volume for standard proteomics sample preparation has been selected for easy sample handling with a benchtop micropipette. SOPs-MS allows for reliable label-free quantification of ∼1200-2700 protein groups from 50 to 1000 MCF10A cells. When applied to small subpopulations of mouse colon crypt cells, SOPs-MS has revealed protein signatures between distinct subpopulation cells with identification of ∼1500-2500 protein groups for each subpopulation. SOPs-MS may pave the way for routine deep proteome profiling of small numbers of cells and low-input samples.

Receptor-Driven ERK Pulses Reconfigure MAPK Signaling and Enable Persistence of Drug-Adapted BRAF-Mutant Melanoma Cells.

Targeted inhibition of oncogenic pathways can be highly effective in halting the rapid growth of tumors but often leads to the emergence of slowly dividing persister cells, which constitute a reservoir for the selection of drug-resistant clones. In BRAFV600E melanomas, RAF and MEK inhibitors efficiently block oncogenic signaling, but persister cells emerge. Here, we show that persister cells escape drug-induced cell-cycle arrest via brief, sporadic ERK pulses generated by transmembrane receptors and growth factors operating in an autocrine/paracrine manner. Quantitative proteomics and computational modeling show that ERK pulsing is enabled by rewiring of mitogen-activated protein kinase (MAPK) signaling: from an oncogenic BRAFV600E monomer-driven configuration that is drug sensitive to a receptor-driven configuration that involves Ras-GTP and RAF dimers and is highly resistant to RAF and MEK inhibitors. Altogether, this work shows that pulsatile MAPK activation by factors in the microenvironment generates a persistent population of melanoma cells that rewires MAPK signaling to sustain non-genetic drug resistance.

Publisher Correction: Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt.

Intestinal stem cells are non-quiescent, dividing epithelial cells that rapidly differentiate into progenitor cells of the absorptive and secretory cell lineages. The kinetics of this process is rapid such that the epithelium is replaced weekly. To determine how the transcriptome and proteome keep pace with rapid differentiation, we developed a new cell sorting method to purify mouse colon epithelial cells. Here we show that alternative mRNA splicing and polyadenylation dominate changes in the transcriptome as stem cells differentiate into progenitors. In contrast, as progenitors differentiate into mature cell types, changes in mRNA levels dominate the transcriptome. RNA processing targets regulators of cell cycle, RNA, cell adhesion, SUMOylation, and Wnt and Notch signaling. Additionally, global proteome profiling detected >2,800 proteins and revealed RNA:protein patterns of abundance and correlation. Paired together, these data highlight new potentials for autocrine and feedback regulation and provide new insights into cell state transitions in the crypt.

Carrier-Assisted Single-Tube Processing Approach for Targeted Proteomics Analysis of Low Numbers of Mammalian Cells.

Heterogeneity in composition is inherent in all cell populations, even those containing a single cell type. Single-cell proteomics characterization of cell heterogeneity is currently achieved by antibody-based technologies, which are limited by the availability of high-quality antibodies. Herein we report a simple, easily implemented, mass spectrometry (MS)-based targeted proteomics approach, termed cLC-SRM (carrier-assisted liquid chromatography coupled to selected reaction monitoring), for reliable multiplexed quantification of proteins in low numbers of mammalian cells. We combine a new single-tube digestion protocol to process low numbers of cells with minimal loss together with sensitive LC-SRM for protein quantification. This single-tube protocol builds upon trifluoroethanol digestion and further minimizes sample losses by tube pretreatment and the addition of carrier proteins. We also optimized the denaturing temperature and trypsin concentration to significantly improve digestion efficiency. cLC-SRM was demonstrated to have sufficient sensitivity for reproducible detection of most epidermal growth factor receptor (EGFR) pathway proteins expressed at levels ≥30 000 and ≥3000 copies per cell for 10 and 100 mammalian cells, respectively. Thus, cLC-SRM enables reliable quantification of low to moderately abundant proteins in less than 100 cells and could be broadly useful for multiplexed quantification of important proteins in small subpopulations of cells or in size-limited clinical samples. Further improvements of this method could eventually enable targeted single-cell proteomics when combined with either SRM or other emerging ultrasensitive MS detection.

Facile carrier-assisted targeted mass spectrometric approach for proteomic analysis of low numbers of mammalian cells.

There is an unmet technical challenge for mass spectrometry (MS)-based proteomic analysis of single mammalian cells. Quantitative proteomic analysis of single cells has been previously achieved by antibody-based immunoassays but is limited by the availability of high-quality antibodies. Herein we report a facile targeted MS-based proteomics method, termed cPRISM-SRM (carrier-assisted high-pressure, high-resolution separations with intelligent selection and multiplexing coupled to selected reaction monitoring), for reliable analysis of low numbers of mammalian cells. The method capitalizes on using "carrier protein" to assist processing of low numbers of cells with minimal loss, high-resolution PRISM separation for target peptide enrichment, and sensitive SRM for protein quantification. We have demonstrated that cPRISM-SRM has sufficient sensitivity to quantify proteins expressed at ≥200,000 copies per cell at the single-cell level and ≥3000 copies per cell in 100 mammalian cells. We envision that with further improvement cPRISM-SRM has the potential to move toward targeted MS-based single-cell proteomics.

A systems perspective of heterocellular signaling.

Signal exchange between different cell types is essential for development and function of multicellular organisms, and its dysregulation is causal in many diseases. Unfortunately, most cell-signaling work has employed single cell types grown under conditions unrelated to their native context. Recent technical developments have started to provide the tools needed to follow signaling between multiple cell types, but gaps in the information they provide have limited their usefulness in building realistic models of heterocellular signaling. Currently, only targeted assays have the necessary sensitivity, selectivity, and spatial resolution to usefully probe heterocellular signaling processes, but these are best used to test specific, mechanistic models. Decades of systems biology research with monocultures has provided a solid foundation for building models of heterocellular signaling, but current models lack a realistic description of regulated proteolysis and the feedback processes triggered within and between cells. Identification and understanding of key regulatory processes in the extracellular environment and of recursive signaling patterns between cells will be essential to building predictive models of heterocellular systems.

Blazing Signature Filter: a library for fast pairwise similarity comparisons.

BACKGROUND: Identifying similarities between datasets is a fundamental task in data mining and has become an integral part of modern scientific investigation. Whether the task is to identify co-expressed genes in large-scale expression surveys or to predict combinations of gene knockouts which would elicit a similar phenotype, the underlying computational task is often a multi-dimensional similarity test. As datasets continue to grow, improvements to the efficiency, sensitivity or specificity of such computation will have broad impacts as it allows scientists to more completely explore the wealth of scientific data. RESULTS: The Blazing Signature Filter (BSF) is a highly efficient pairwise similarity algorithm which enables extensive data mining within a reasonable amount of time. The algorithm transforms datasets into binary metrics, allowing it to utilize the computationally efficient bit operators and provide a coarse measure of similarity. We demonstrate the utility of our algorithm using two common bioinformatics tasks: identifying data sets with similar gene expression profiles, and comparing annotated genomes. CONCLUSIONS: The BSF is a highly efficient pairwise similarity algorithm that can scale to billions of comparisons without the need for specialized hardware.

Targeted Quantification of Phosphorylation Dynamics in the Context of EGFR-MAPK Pathway.

Large-scale phosphoproteomics with coverage of over 10,000 sites of phosphorylation have now been routinely achieved with advanced mass spectrometry (MS)-based workflows. However, accurate targeted MS-based quantification of phosphorylation dynamics, an important direction for gaining quantitative understanding of signaling pathways or networks, has been much less investigated. Herein, we report an assessment of the targeted workflow in the context of signal transduction pathways, using the epidermal growth factor receptor (EGFR)-mitogen-activated protein kinase (MAPK) pathway as our model. A total of 43 phosphopeptides from the EGFR-MAPK pathway were selected for the study. The recovery and sensitivity of two commonly used enrichment methods, immobilized metal affinity chromatography (IMAC) and titanium oxide (TiO2), combined with selected reaction monitoring (SRM)-MS were evaluated. The recovery of phosphopeptides by IMAC and TiO2 enrichment was quantified to be 38 ± 5% and 58 ± 20%, respectively, based on internal standards. Moreover, both enrichment methods provided comparable sensitivity from 1 to 100 µg starting peptides. Robust quantification was consistently achieved for most targeted phosphopeptides when starting with 25-100 µg peptides. However, the numbers of quantified targets significantly dropped when peptide samples were in the 1-25 µg range. Finally, IMAC-SRM was applied to quantify signaling dynamics of EGFR-MAPK pathway in Hs578T cells following 10 ng/mL EGF treatment. The kinetics of phosphorylation clearly revealed early and late phases of phosphorylation, even for very low abundance proteins. These results demonstrate the feasibility of robust targeted quantification of phosphorylation dynamics for specific pathways, even starting with relatively small amounts of protein.

How low can you go?

Extremely low numbers of active epidermal growth factor receptors are sufficient to drive tumor growth.

Decoding Signal Processing at the Single-Cell Level.

The feedforward circuitry regulating ERK-dependent early response genes acts as a signal integrator rather than a signal persistence detector.