PeptideDepot: flexible relational database for visual analysis of quantitative proteomic data and integration of existing protein information

SILAC Phosphoproteomics Reveals Unique Signaling Circuits in CAR-T Cells and the Inhibition of B Cell-Activating Phosphorylation in Target Cells

Chimeric antigen receptor (CAR) is a single-pass transmembrane receptor designed to specifically target and eliminate cancers. While CARs prove highly efficacious against B cell malignancies, the intracellular signaling events which promote CAR T cell activity remain elusive. To gain further insight into both CAR T cell signaling and the potential signaling response of cells targeted by CAR, we analyzed phosphopeptides captured by two separate phosphoenrichment strategies from third generation CD19-CAR T cells cocultured with SILAC labeled Raji B cells by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Here, we report that CD19-CAR T cells upregulated several key phosphorylation events also observed in canonical T cell receptor (TCR) signaling, while Raji B cells exhibited a significant decrease in B cell receptor-signaling related phosphorylation events in response to coculture. Our data suggest that CD19-CAR stimulation activates a mixture of unique CD19-CAR-specific signaling pathways and canonical TCR signaling, while global phosphorylation in Raji B cells is reduced after association with the CD19-CAR T cells.

Cyclin-Dependent Kinase-Mediated Phosphorylation of FANCD2 Promotes Mitotic Fidelity

Fanconi anemia (FA) is a rare genetic disease characterized by increased risk for bone marrow failure and cancer. The FA proteins function together to repair damaged DNA. A central step in the activation of the FA pathway is the monoubiquitination of the FANCD2 and FANCI proteins, which occurs upon exposure to DNA-damaging agents and during the S phase of the cell cycle. The regulatory mechanisms governing S-phase monoubiquitination, in particular, are poorly understood. In this study, we have identified a cyclin-dependent kinase (CDK) regulatory phosphosite (S592) proximal to the site of FANCD2 monoubiquitination. FANCD2 S592 phosphorylation was detected by liquid chromatography-tandem mass spectrometry (LC-MS/MS) and by immunoblotting with an S592 phospho-specific antibody. Mutation of S592 leads to abrogated monoubiquitination of FANCD2 during the S phase. Furthermore, FA-D2 (FANCD2^-/-) patient cells expressing S592 mutants display reduced proliferation under conditions of replication stress and increased mitotic aberrations, including micronuclei and multinucleated cells. Our findings describe a novel cell cycle-specific regulatory mechanism for the FANCD2 protein that promotes mitotic fidelity.

Extracellular Vesicles Promote Arteriogenesis in Chronically Ischemic Myocardium in the Setting of Metabolic Syndrome

Background

Ischemic heart disease continues to be a leading cause of mortality in patients. Extracellular vesicles (EVs) provide a potential for treatment that may induce collateral vessel growth to increase myocardial perfusion.

Methods and Results

Nineteen male Yorkshire pigs were given a high‐fat diet for 4 weeks, then underwent placement of an ameroid constrictor on the left circumflex artery to induce chronic myocardial ischemia. Two weeks later, the pigs received either intramyocardial vehicle (n=6), EVs (high‐fat diet with myocardial EV injection [HVM]; n=8), or HVM and calpain inhibition (n=5). Five weeks later, myocardial function, perfusion, coronary vascular density, and cell signaling were examined. Perfusion in the collateral‐dependent myocardium was increased during rapid ventricular pacing in the HVM group in both nonischemic (P=0.04) and ischemic areas of the ventricle (P=0.05). Cardiac output and stroke volume were significantly improved in the HVM group compared with the control group during ventricular pacing (P=0.006). Increased arteriolar density was seen in the HVM group in both nonischemic and ischemic myocardium (P=0.003 for both). However, no significant changes in the capillary density were observed between the control, HVM, and HVM and calpain inhibition groups (P=0.07). The group that received EVs with oral calpain inhibition had neither increased vessel density (P>0.99) nor improvement in blood flow or cardiac function (P=0.48) when compared with the control group.

Conclusions

These findings suggest that EVs promote angiogenesis in areas of chronic myocardial ischemia and improve cardiac function under conditions of diet‐induced metabolic syndrome.

Septin-2 is overexpressed in epithelial ovarian cancer and mediates proliferation via regulation of cellular metabolic proteins

Epithelial Ovarian Cancer (EOC) is associated with dismal survival rates due to the fact that patients are frequently diagnosed at an advanced stage and eventually become resistant to traditional chemotherapeutics. Hence, there is a crucial need for new and innovative therapies. Septin-2, a member of the septin family of GTP binding proteins, has been characterized in EOC for the first time and represents a potential future target. Septin-2 was found to be overexpressed in serous and clear cell human patient tissue compared to benign disease. Stable septin-2 knockdown clones developed in an ovarian cancer cell line exhibited a significant decrease in proliferation rates. Comparative label-free proteomic analysis of septin-2 knockdown cells revealed differential protein expression of pathways associated with the TCA cycle, acetyl CoA, proteasome and spliceosome. Further validation of target proteins indicated that septin-2 plays a predominant role in post-transcriptional and translational modifications as well as cellular metabolism, and suggested the potential novel role of septin-2 in promoting EOC tumorigenesis through these mechanisms.

Lck promotes Zap70-dependent LAT phosphorylation by bridging Zap70 to LAT

T cell-antigen receptor (TCR) signaling requires the sequential activities of the kinases Lck and Zap70. Upon TCR stimulation, Lck phosphorylates the TCR, thus leading to the recruitment, phosphorylation, and activation of Zap70. Lck binds and stabilizes phosho-Zap70 by using its SH2 domain, and Zap70 phosphorylates the critical adaptors LAT and SLP76, which coordinate downstream signaling. It is unclear whether phosphorylation of these adaptors occurs through passive diffusion or active recruitment. We report the discovery of a conserved proline-rich motif in LAT that mediates efficient LAT phosphorylation. Lck associates with this motif via its SH3 domain, and with phospho-Zap70 via its SH2 domain, thereby acting as a molecular bridge that facilitates the colocalization of Zap70 and LAT. Elimination of this proline-rich motif compromises TCR signaling and T cell development. These results demonstrate the remarkable multifunctionality of Lck, wherein each of its domains has evolved to orchestrate a distinct step in TCR signaling.

Extracellular Vesicle Injection Improves Myocardial Function and Increases Angiogenesis in a Swine Model of Chronic Ischemia

BACKGROUND

Mesenchymal stem cell-derived extracellular vesicles (EVs) are believed to be cardioprotective in myocardial infarct. The objective of this study was to examine the effects of human mesenchymal cell-derived EV injection on cardiac function, myocardial blood flow, and vessel density in the setting of chronic myocardial ischemia.

METHODS AND RESULTS

Twenty-three Yorkshire swine underwent placement of an ameroid constrictor on their left circumflex artery. Two weeks later, the animals were split into 2 groups: the control group (CON; n=7) and the EV myocardial injection group (MVM; n=10). The MVM group underwent myocardial injection of 50 μg of EVs in 2 mL 0.9% saline into the ischemic myocardium. Five weeks later, the pigs underwent a harvest procedure, and the left ventricular myocardium was analyzed. Absolute blood flow and the ischemic/nonischemic myocardial perfusion ratio were increased in the ischemic myocardium in the MVM group compared with the CON group. Pigs in the MVM group had increased capillary and arteriolar density in the ischemic myocardial tissue compared with CON pigs. There was an increase in expression of the phospho-mitogen-activated protein kinase/mitogen-activated protein kinase ratio, the phospho-endothelial nitric oxide synthase/endothelial nitric oxide synthase ratio, and total protein kinase B in the MVM group compared with CON. There was an increase in cardiac output and stroke volume in the MVM group compared with CON.

CONCLUSIONS

In the setting of chronic myocardial ischemia, myocardial injection of human mesenchymal cell-derived EVs increases blood flow to ischemic myocardial tissue by induction of capillary and arteriolar growth via activation of the protein kinase B/endothelial nitric oxide synthase and mitogen-activated protein kinase signaling pathways resulting in increased cardiac output and stroke volume.

Proteomic analysis of laser capture microdissected focal lesions in a rat model of progenitor marker-positive hepatocellular carcinoma

We have shown previously that rapamycin, the canonical inhibitor of the mechanistic target of rapamycin (mTOR) complex 1, markedly inhibits the growth of focal lesions in the resistant hepatocyte (Solt-Farber) model of hepatocellular carcinoma (HCC) in the rat. In the present study, we characterized the proteome of persistent, pre-neoplastic focal lesions in this model. One group was administered rapamycin by subcutaneous pellet for 3 weeks following partial hepatectomy and euthanized 4 weeks after the cessation of rapamycin. A second group received placebo pellets. Results were compared to unmanipulated control animals and to animals that underwent an incomplete Solt-Farber protocol to activate hepatic progenitor cells. Regions of formalin-fixed, paraffin-embedded tissue were obtained by laser capture microdissection (LCM). Proteomic analysis yielded 11,070 unique peptides representing 2,227 proteins. Quantitation of the peptides showed increased abundance of known HCC markers (e.g., glutathione S-transferase-P, epoxide hydrolase, 6 others) and potential markers (e.g., aflatoxin aldehyde reductase, glucose 6-phosphate dehydrogenase, 10 others) in foci from placebo-treated and rapamycin-treated rats. Peptides derived from cytochrome P450 enzymes were generally reduced. Comparisons of the rapamycin samples to normal liver and to the progenitor cell model indicated that rapamycin attenuated a loss of differentiation relative to placebo. We conclude that early administration of rapamycin in the Solt-Farber model not only inhibits the growth of pre-neoplastic foci but also attenuates the loss of differentiated function. In addition, we have demonstrated that the combination of LCM and mass spectrometry-based proteomics is an effective approach to characterize focal liver lesions.

Drosophila Insulin-Like Peptides DILP2 and DILP5 Differentially Stimulate Cell Signaling and Glycogen Phosphorylase to Regulate Longevity

Insulin and IGF signaling (IIS) is a complex system that controls diverse processes including growth, development, metabolism, stress responses, and aging. Drosophila melanogaster IIS is propagated by eight Drosophila insulin-like peptides (DILPs), homologs of both mammalian insulin and IGFs, with various spatiotemporal expression patterns and functions. DILPs 1-7 are thought to act through a single Drosophila insulin/IGF receptor, InR, but it is unclear how the DILPs thereby mediate a range of physiological phenotypes. We determined the distinct cell signaling effects of DILP2 and DILP5 stimulation upon Drosophila S2 cells. DILP2 and DILP5 induced similar transcriptional patterns but differed in signal transduction kinetics. DILP5 induced sustained phosphorylation of Akt, while DILP2 produced acute, transient Akt phosphorylation. Accordingly, we used phosphoproteomic analysis to identify distinct patterns of non-genomic signaling induced by DILP2 and DILP5. Across all treatments and replicates, 5,250 unique phosphopeptides were identified, representing 1,575 proteins. Among these peptides, DILP2, but not DILP5, dephosphorylated Ser15 on glycogen phosphorylase (GlyP), and DILP2, but not DILP5, was subsequently shown to repress enzymatic GlyP activity in S2 cells. The functional consequences of this difference were evaluated in adult Drosophila dilp mutants: dilp2 null adults have elevated GlyP enzymatic activity relative to wild type, while dilp5 mutants have reduced GlyP activity. In flies with intact insulin genes, GlyP overexpression extended lifespan in a Ser15 phosphorylation-dependent manner. In dilp2 mutants, that are otherwise long-lived, longevity was repressed by expression of phosphonull GlyP that is enzymatically inactive. Overall, DILP2, unlike DILP5, signals to affect longevity in part through its control of phosphorylation to deactivate glycogen phosphorylase, a central modulator of glycogen storage and gluconeogenesis.

mzStudio: A Dynamic Digital Canvas for User-Driven Interrogation of Mass Spectrometry Data

Although not yet truly ‘comprehensive’, modern mass spectrometry-based experiments can generate quantitative data for a meaningful fraction of the human proteome. Importantly for large-scale protein expression analysis, robust data pipelines are in place for identification of un-modified peptide sequences and aggregation of these data to protein-level quantification. However, interoperable software tools that enable scientists to computationally explore and document novel hypotheses for peptide sequence, modification status, or fragmentation behavior are not well-developed. Here, we introduce mzStudio, an open-source Python module built on our multiplierz project. This desktop application provides a highly-interactive graphical user interface (GUI) through which scientists can examine and annotate spectral features, re-search existing PSMs to test different modifications or new spectral matching algorithms, share results with colleagues, integrate other domain-specific software tools, and finally create publication-quality graphics. mzStudio leverages our common application programming interface (mzAPI) for access to native data files from multiple instrument platforms, including ion trap, quadrupole time-of-flight, Orbitrap, matrix-assisted laser desorption ionization, and triple quadrupole mass spectrometers and is compatible with several popular search engines including Mascot, Proteome Discoverer, X!Tandem, and Comet. The mzStudio toolkit enables researchers to create a digital provenance of data analytics and other evidence that support specific peptide sequence assignments.

Proteomic responses to elevated ocean temperature in ovaries of the ascidian Ciona intestinalis

Ciona intestinalis, a common sea squirt, exhibits lower reproductive success at the upper extreme of the water temperatures it experiences in coastal New England. In order to understand the changes in protein expression associated with elevated temperatures, and possible response to global temperature change, we reared C. intestinalis from embryos to adults at 18°C (a temperature at which they reproduce normally at our collection site in Rhode Island) and 22°C (the upper end of the local temperature range). We then dissected ovaries from animals at each temperature, extracted protein, and measured proteomic levels using shotgun mass spectrometry (LC-MS/MS). 1532 proteins were detected at a 1% false discovery rate present in both temperature groups by our LC-MS/MS method. 62 of those proteins are considered up- or down-regulated according to our statistical criteria. Principal component analysis shows a clear distinction in protein expression pattern between the control (18°C) group and high temperature (22°C) group. Similar to previous studies, cytoskeletal and chaperone proteins are upregulated in the high temperature group. Unexpectedly, we find evidence that proteolysis is downregulated at the higher temperature. We propose a working model for the high temperature response in C. intestinalis ovaries whereby increased temperature induces upregulation of signal transduction pathways involving PTPN11 and CrkL, and activating coordinated changes in the proteome especially in large lipid transport proteins, cellular stress responses, cytoskeleton, and downregulation of energy metabolism.

A PLC-γ1 Feedback Pathway Regulates Lck Substrate Phosphorylation at the T-Cell Receptor and SLP-76 Complex

Phospholipase C gamma 1 (PLC-γ1) occupies a critically important position in the T-cell signaling pathway. While its functions as a regulator of both Ca²⁺ signaling and PKC-family kinases are well characterized, PLC-γ1's role in the regulation of early T-cell receptor signaling events is incompletely understood. Activation of the T-cell receptor leads to the formation of a signalosome complex between SLP-76, LAT, PLC-γ1, Itk, and Vav1. Recent studies have revealed the existence of both positive and negative feedback pathways from SLP-76 to the apical kinase in the pathway, Lck. To determine if PLC-γ1 contributes to the regulation of these feedback networks, we performed a quantitative phosphoproteomic analysis of PLC-γ1-deficient T cells. These data revealed a previously unappreciated role for PLC-γ1 in the positive regulation of Zap-70 and T-cell receptor tyrosine phosphorylation. Conversely, PLC-γ1 negatively regulated the phosphorylation of SLP-76-associated proteins, including previously established Lck substrate phosphorylation sites within this complex. While the positive and negative regulatory phosphorylation sites on Lck were largely unchanged, Tyr¹⁹² phosphorylation was elevated in Jgamma1. The data supports a model wherein Lck's targeting, but not its kinase activity, is altered by PLC-γ1, possibly through Lck Tyr¹⁹² phosphorylation and increased association of the kinase with protein scaffolds SLP-76 and TSAd.

Highly reproducible improved label-free quantitative analysis of cellular phosphoproteome by optimization of LC-MS/MS gradient and analytical column construction

Expanding the sequencing depth of the peptides with a statistically significant quantitative change derived from a biological stimulation is critical. Here we demonstrate that optimization of LC gradient and analytical column construction can reveal over 30,000 unique peptides and 23,000 phosphopeptides at high confidence. The quantitative reproducibility of different analytical workflows was evaluated by comparing the phosphoproteome of CD3/4 stimulated and unstimulated T-cells as a model system. A fritless, 50cm-long column packed with 1.9μm particles operated with a standard pressure HPLC significantly improved the sequencing depth 51% and decreased the selected ion chromatogram peak spreading. Most importantly, under the optimal workflow we observed an improvement of over 300% in detection of significantly changed phosphopeptides in the stimulated cells compared with the other workflows. The discovery power of the optimized column configuration was illustrated by identification of significantly altered phosphopeptides harboring novel sites from proteins previously established as important in T cell signaling including A-Raf, B-Raf, c-Myc, CARMA1, Fyn, ITK, LAT, NFAT1/2/3, PKCα, PLCγ1/2, RAF1, and SOS1. Taken together, our results reveal the analytical power of optimized chromatography using sub 2μm particles for the analysis of the T cell phosphoproteome to reveal a vast landscape of significantly altered phosphorylation changes in response to T cell receptor stimulation.

Quantitative Phosphoproteomic Analysis of T-Cell Receptor Signaling

TCR signaling critically depends on protein phosphorylation across many proteins. Localization of each phosphorylation event relative to the T-cell receptor (TCR) and canonical T-cell signaling proteins will provide clues about the structure of TCR signaling networks. Quantitative phosphoproteomic analysis by mass spectrometry provides a wide-scale view of cellular phosphorylation networks. However, analysis of phosphorylation by mass spectrometry is still challenging due to the relative low abundance of phosphorylated proteins relative to all proteins and the extraordinary diversity of phosphorylation sites across the proteome. Highly selective enrichment of phosphorylated peptides is essential to provide the most comprehensive view of the phosphoproteome. Optimization of phosphopeptide enrichment methods coupled with highly sensitive mass spectrometry workflows significantly improves the sequencing depth of the phosphoproteome to over 10,000 unique phosphorylation sites from complex cell lysates. Here we describe a step-by-step method for phosphoproteomic analysis that has achieved widespread success for identification of serine, threonine, and tyrosine phosphorylation. Reproducible quantification of relative phosphopeptide abundance is provided by intensity-based label-free quantitation. An ideal set of mass spectrometry analysis parameters is also provided that optimize the yield of identified sites. We also provide guidelines for the bioinformatic analysis of this type of data to assess the quality of the data and to comply with proteomic data reporting requirements.

Vav1 Regulates T-Cell Activation through a Feedback Mechanism and Crosstalk between the T-Cell Receptor and CD28

Vav1, a Rac/Rho guanine nucleotide exchange factor and a critical component of the T-cell receptor (TCR) signaling cascade is tyrosine phosphorylated rapidly in response to T-cell activation. Vav1 has established roles in proliferation, cytokine secretion, Ca(2+) responses, and actin cytoskeleton regulation; however, its function in the regulation of phosphorylation of TCR components, including the ζ chain, the CD3 δ, ε, γ chains, and the associated kinases Lck and ZAP-70, is not well established. To obtain a more comprehensive picture of the role of Vav1 in receptor proximal signaling, we performed a wide-scale characterization of Vav1-dependent tyrosine phosphorylation events using quantitative phosphoproteomic analysis of Vav1-deficient T cells across a time course of TCR stimulation. Importantly, this study revealed a new function for Vav1 in the negative feedback regulation of the phosphorylation of immunoreceptor tyrosine-based activation motifs within the ζ chains, CD3 δ, ε, γ chains, as well as activation sites on the critical T cell tyrosine kinases Itk, Lck, and ZAP-70. Our study also uncovered a previously unappreciated role for Vav1 in crosstalk between the CD28 and TCR signaling pathways.

The MARK-AGE extended database: data integration and pre-processing

MARK-AGE is a recently completed European population study, where bioanalytical and anthropometric data were collected from human subjects at a large scale. To facilitate data analysis and mathematical modelling, an extended database had to be constructed, integrating the data sources that were part of the project. This step involved checking, transformation and documentation of data. The success of downstream analysis mainly depends on the preparation and quality of the integrated data. Here, we present the pre-processing steps applied to the MARK-AGE data to ensure high quality and reliability in the MARK-AGE Extended Database. Various kinds of obstacles that arose during the project are highlighted and solutions are presented.

MARK-AGE data management: Cleaning, exploration and visualization of data

Databases are an organized collection of data and necessary to investigate a wide spectrum of research questions. For data evaluation analyzers should be aware of possible data quality problems that can compromise results validity. Therefore data cleaning is an essential part of the data management process, which deals with the identification and correction of errors in order to improve data quality. In our cross-sectional study, biomarkers of ageing, analytical, anthropometric and demographic data from about 3000 volunteers have been collected in the MARK-AGE database. Although several preventive strategies were applied before data entry, errors like miscoding, missing values, batch problems etc., could not be avoided completely. Such errors can result in misleading information and affect the validity of the performed data analysis. Here we present an overview of the methods we applied for dealing with errors in the MARK-AGE database. We especially describe our strategies for the detection of missing values, outliers and batch effects and explain how they can be handled to improve data quality. Finally we report about the tools used for data exploration and data sharing between MARK-AGE collaborators.

The catalytic activity of the kinase ZAP-70 mediates basal signaling and negative feedback of the T cell receptor pathway

T cell activation by antigens binding to the T cell receptor (TCR) must be properly regulated to ensure normal T cell development and effective immune responses to pathogens and transformed cells while avoiding autoimmunity. The Src family kinase Lck and the Syk family kinase ZAP-70 (ζ chain-associated protein kinase of 70 kD) are sequentially activated in response to TCR engagement and serve as critical components of the TCR signaling machinery that leads to T cell activation. We performed a mass spectrometry-based phosphoproteomic study comparing the quantitative differences in the temporal dynamics of phosphorylation in stimulated and unstimulated T cells with or without inhibition of ZAP-70 catalytic activity. The data indicated that the kinase activity of ZAP-70 stimulates negative feedback pathways that target Lck and thereby modulate the phosphorylation patterns of the immunoreceptor tyrosine-based activation motifs (ITAMs) of the CD3 and ζ chain components of the TCR and of signaling molecules downstream of Lck, including ZAP-70. We developed a computational model that provides a mechanistic explanation for the experimental findings on ITAM phosphorylation in wild-type cells, ZAP-70-deficient cells, and cells with inhibited ZAP-70 catalytic activity. This model incorporated negative feedback regulation of Lck activity by the kinase activity of ZAP-70 and predicted the order in which tyrosines in the ITAMs of TCR ζ chains must be phosphorylated to be consistent with the experimental data.

Wide-scale quantitative phosphoproteomic analysis reveals that cold treatment of T cells closely mimics soluble antibody stimulation

The activation of T lymphocytes through antigen-mediated T cell receptor (TCR) clustering is vital in regulating the adaptive immune response. Although T cell receptor signaling has been extensively studied, the fundamental mechanisms for signal initiation are not fully understood. Reduced temperatures have initiated some of the hallmarks of TCR signaling, such as increased phosphorylation and activation on ERK and calcium release from the endoplasmic reticulum, as well as coalesced the T cell membrane microdomains. The precise mechanism of the TCR signaling initiation due to temperature change remains obscure. One critical question is whether the signaling initiated by the cold treatment of T cells differs from the signaling initiated by the cross-linking of the T cell receptor. To address this uncertainty, we performed a wide-scale, quantitative mass-spectrometry-based phosphoproteomic analysis on T cells stimulated either by temperature shifts or through the cross-linking of the TCR. Careful statistical comparisons between the two stimulations revealed a striking level of identity among the subset of 339 sites that changed significantly with both stimulations. This study demonstrates for the first time, in unprecedented detail, that T cell cold treatment was sufficient to initiate signaling patterns that were nearly identical to those of soluble antibody stimulation, shedding new light on the mechanism of activation of these critically important immune cells.

SRC homology 2 domain-containing leukocyte phosphoprotein of 76 kDa (SLP-76) N-terminal tyrosine residues regulate a dynamic signaling equilibrium involving feedback of proximal T-cell receptor (TCR) signaling

SRC homology 2 domain-containing leukocyte phosphoprotein of 76 kDa (SLP-76) is a cytosolic adaptor protein that plays an important role in the T-cell receptor-mediated T-cell signaling pathway. SLP-76 links proximal receptor stimulation to downstream effectors through interaction with many signaling proteins. Previous studies showed that mutation of three tyrosine residues, Tyr(112), Tyr(128), and Tyr(145), in the N terminus of SLP-76 results in severely impaired phosphorylation and activation of Itk and PLCγ1, which leads to defective calcium mobilization, Erk activation, and NFAT activation. To expand our knowledge of the role of N-terminal phosphorylation of SLP-76 from these three tyrosine sites, we characterized nearly 1000 tyrosine phosphorylation sites via mass spectrometry in SLP-76 reconstituted wild-type cells and SLP-76 mutant cells in which three tyrosine residues were replaced with phenylalanines (Y3F mutant). Mutation of the three N-terminal tyrosine residues of SLP-76 phenocopied SLP-76-deficient cells for the majority of tyrosine phosphorylation sites observed, including feedback on proximal T-cell receptor signaling proteins. Meanwhile, reversed phosphorylation changes were observed on Tyr(192) of Lck when we compared mutants to the complete removal of SLP-76. In addition, N-terminal tyrosine sites of SLP-76 also perturbed phosphorylation of Tyr(440) of Fyn, Tyr(702) of PLCγ1, Tyr(204), Tyr(397), and Tyr(69) of ZAP-70, revealing new modes of regulation on these sites. All these findings confirmed the central role of N-terminal tyrosine sites of SLP-76 in the pathway and also shed light on novel signaling events that are uniquely regulated by SLP-76 N-terminal tyrosine residues.

PhosFox: a bioinformatics tool for peptide-level processing of LC-MS/MS-based phosphoproteomic data

Background

It is possible to identify thousands of phosphopeptides and –proteins in a single experiment with mass spectrometry-based phosphoproteomics. However, a current bottleneck is the downstream data analysis which is often laborious and requires a number of manual steps.

Results

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

Conclusions

PhosFox facilitates efficient and more in-depth comparisons between phosphoproteins in case–control settings. The open source implementation is easily extendable to accommodate additional features for widespread application use cases.

Quantitative phosphoproteomic analysis identifies activation of the RET and IGF-1R/IR signaling pathways in neuroblastoma

Neuroblastoma is an embryonal tumor of childhood with a heterogenous clinical presentation that reflects differences in activation of complex biological signaling pathways. Protein phosphorylation is a key component of cellular signal transduction and plays a critical role in processes that control cancer cell growth and survival. We used shotgun LC/MS to compare phosphorylation between a human MYCN amplified neuroblastoma cell line (NB10), modeling a resistant tumor, and a human neural precursor cell line (NPC), modeling a normal baseline neural crest cell. 2181 unique phosphorylation sites representing 1171 proteins and 2598 phosphopeptides were found. Protein kinases accounted for 6% of the proteome, with a predominance of tyrosine kinases, supporting their prominent role in oncogenic signaling pathways. Highly abundant receptor tyrosine kinase (RTK) phosphopeptides in the NB10 cell line relative to the NPC cell line included RET, insulin-like growth factor 1 receptor/insulin receptor (IGF-1R/IR), and fibroblast growth factor receptor 1 (FGFR1). Multiple phosphorylated peptides from downstream mediators of the PI3K/AKT/mTOR and RAS pathways were also highly abundant in NB10 relative to NPC. Our analysis highlights the importance of RET, IGF-1R/IR and FGFR1 as RTKs in neuroblastoma and suggests a methodology that can be used to identify potential novel biological therapeutic targets. Furthermore, application of this previously unexploited technology in the clinic opens the possibility of providing a new wide-scale molecular signature to assess disease progression and prognosis.

ERK positive feedback regulates a widespread network of tyrosine phosphorylation sites across canonical T cell signaling and actin cytoskeletal proteins in Jurkat T cells

Competing positive and negative signaling feedback pathways play a critical role in tuning the sensitivity of T cell receptor activation by creating an ultrasensitive, bistable switch to selectively enhance responses to foreign ligands while suppressing signals from self peptides. In response to T cell receptor agonist engagement, ERK is activated to positively regulate T cell receptor signaling through phosphorylation of Ser⁵⁹ Lck. To obtain a wide-scale view of the role of ERK in propagating T cell receptor signaling, a quantitative phosphoproteomic analysis of 322 tyrosine phosphorylation sites by mass spectrometry was performed on the human Jurkat T cell line in the presence of U0126, an inhibitor of ERK activation. Relative to controls, U0126-treated cells showed constitutive decreases in phosphorylation through a T cell receptor stimulation time course on tyrosine residues found on upstream signaling proteins (CD3 chains, Lck, ZAP-70), as well as downstream signaling proteins (VAV1, PLCγ1, Itk, NCK1). Additional constitutive decreases in phosphorylation were found on the majority of identified proteins implicated in the regulation of actin cytoskeleton pathway. Although the majority of identified sites on T cell receptor signaling proteins showed decreases in phosphorylation, Tyr⁵⁹⁸ of ZAP-70 showed elevated phosphorylation in response to U0126 treatment, suggesting differential regulation of this site via ERK feedback. These findings shed new light on ERK’s role in positive feedback in T cell receptor signaling and reveal novel signaling events that are regulated by this kinase, which may fine tune T cell receptor activation.

Phosphoproteomic analysis of liver homogenates

Regulation of protein function via reversible phosphorylation is an essential component of cell signaling. Our ability to understand complex phosphorylation networks in the physiological context of a whole organism or tissue remains limited. This is largely due to the technical challenge of isolating serine/threonine phosphorylated peptides from a tissue sample. In the present study, we developed a phosphoproteomic strategy to purify and identify phosphopeptides from a tissue sample by employing protein gel filtration, protein strong anion exchange and strong cation exchange (SCX) chromatography, peptide SCX chromatography, and TiO(2) affinity purification. By applying this strategy to the mass spectrometry-based analysis of rat liver homogenates, we were able to identify with high confidence and quantify over 4,000 unique phosphopeptides. Finally, the reproducibility of our methodology was demonstrated by its application to analysis of the mammalian Target of Rapamycin (mTOR) signaling pathways in liver samples obtained from rats in which hepatic mTOR was activated by refeeding following a period of fasting.

VESPA: software to facilitate genomic annotation of prokaryotic organisms through integration of proteomic and transcriptomic data

Background

The procedural aspects of genome sequencing and assembly have become relatively inexpensive, yet the full, accurate structural annotation of these genomes remains a challenge. Next-generation sequencing transcriptomics (RNA-Seq), global microarrays, and tandem mass spectrometry (MS/MS)-based proteomics have demonstrated immense value to genome curators as individual sources of information, however, integrating these data types to validate and improve structural annotation remains a major challenge. Current visual and statistical analytic tools are focused on a single data type, or existing software tools are retrofitted to analyze new data forms. We present Visual Exploration and Statistics to Promote Annotation (VESPA) is a new interactive visual analysis software tool focused on assisting scientists with the annotation of prokaryotic genomes though the integration of proteomics and transcriptomics data with current genome location coordinates.

Results

VESPA is a desktop Java™ application that integrates high-throughput proteomics data (peptide-centric) and transcriptomics (probe or RNA-Seq) data into a genomic context, all of which can be visualized at three levels of genomic resolution. Data is interrogated via searches linked to the genome visualizations to find regions with high likelihood of mis-annotation. Search results are linked to exports for further validation outside of VESPA or potential coding-regions can be analyzed concurrently with the software through interaction with BLAST. VESPA is demonstrated on two use cases (Yersinia pestis Pestoides F and Synechococcus sp. PCC 7002) to demonstrate the rapid manner in which mis-annotations can be found and explored in VESPA using either proteomics data alone, or in combination with transcriptomic data.

Conclusions

VESPA is an interactive visual analytics tool that integrates high-throughput data into a genomic context to facilitate the discovery of structural mis-annotations in prokaryotic genomes. Data is evaluated via visual analysis across multiple levels of genomic resolution, linked searches and interaction with existing bioinformatics tools. We highlight the novel functionality of VESPA and core programming requirements for visualization of these large heterogeneous datasets for a client-side application. The software is freely available at https://www.biopilot.org/docs/Software/Vespa.php.

Phosphoproteomic profiling of in vivo signaling in liver by the mammalian target of rapamycin complex 1 (mTORC1)

BACKGROUND

Our understanding of signal transduction networks in the physiological context of an organism remains limited, partly due to the technical challenge of identifying serine/threonine phosphorylated peptides from complex tissue samples. In the present study, we focused on signaling through the mammalian target of rapamycin (mTOR) complex 1 (mTORC1), which is at the center of a nutrient- and growth factor-responsive cell signaling network. Though studied extensively, the mechanisms involved in many mTORC1 biological functions remain poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

We developed a phosphoproteomic strategy to purify, enrich and identify phosphopeptides from rat liver homogenates. Using the anticancer drug rapamycin, the only known target of which is mTORC1, we characterized signaling in liver from rats in which the complex was maximally activated by refeeding following 48 hr of starvation. Using protein and peptide fractionation methods, TiO(2) affinity purification of phosphopeptides and mass spectrometry, we reproducibly identified and quantified over four thousand phosphopeptides. Along with 5 known rapamycin-sensitive phosphorylation events, we identified 62 new rapamycin-responsive candidate phosphorylation sites. Among these were PRAS40, gephyrin, and AMP kinase 2. We observed similar proportions of increased and reduced phosphorylation in response to rapamycin. Gene ontology analysis revealed over-representation of mTOR pathway components among rapamycin-sensitive phosphopeptide candidates.

CONCLUSIONS/SIGNIFICANCE

In addition to identifying potential new mTORC1-mediated phosphorylation events, and providing information relevant to the biology of this signaling network, our experimental and analytical approaches indicate the feasibility of large-scale phosphoproteomic profiling of tissue samples to study physiological signaling events in vivo.

mzResults: an interactive viewer for interrogation and distribution of proteomics results

The growing use of mass spectrometry in the context of biomedical research has been accompanied by an increased demand for distribution of results in a format that facilitates rapid and efficient validation of claims by reviewers and other interested parties. However, the continued evolution of mass spectrometry hardware, sample preparation methods, and peptide identification algorithms complicates standardization and creates hurdles related to compliance with journal submission requirements. Moreover, the recently announced Philadelphia Guidelines (1, 2) suggest that authors provide native mass spectrometry data files in support of their peer-reviewed research articles. These trends highlight the need for data viewers and other tools that work independently of manufacturers' proprietary data systems and seamlessly connect proteomics results with original data files to support user-driven data validation and review. Based upon our recently described API(1)-based framework for mass spectrometry data analysis (3, 4), we created an interactive viewer (mzResults) that is built on established database standards and enables efficient distribution and interrogation of results associated with proteomics experiments, while also providing a convenient mechanism for authors to comply with data submission standards as described in the Philadelphia Guidelines. In addition, the architecture of mzResults supports in-depth queries of the native mass spectrometry files through our multiplierz software environment. We use phosphoproteomics data to illustrate the features and capabilities of mzResults.

HTAPP: high-throughput autonomous proteomic pipeline

Recent advances in the speed and sensitivity of mass spectrometers and in analytical methods, the exponential acceleration of computer processing speeds, and the availability of genomic databases from an array of species and protein information databases have led to a deluge of proteomic data. The development of a lab-based automated proteomic software platform for the automated collection, processing, storage, and visualization of expansive proteomic data sets is critically important. The high-throughput autonomous proteomic pipeline described here is designed from the ground up to provide critically important flexibility for diverse proteomic workflows and to streamline the total analysis of a complex proteomic sample. This tool is composed of a software that controls the acquisition of mass spectral data along with automation of post-acquisition tasks such as peptide quantification, clustered MS/MS spectral database searching, statistical validation, and data exploration within a user-configurable lab-based relational database. The software design of high-throughput autonomous proteomic pipeline focuses on accommodating diverse workflows and providing missing software functionality to a wide range of proteomic researchers to accelerate the extraction of biological meaning from immense proteomic data sets. Although individual software modules in our integrated technology platform may have some similarities to existing tools, the true novelty of the approach described here is in the synergistic and flexible combination of these tools to provide an integrated and efficient analysis of proteomic samples.