High-throughput mass spectrometry maps the sepsis plasma proteome and differences in patient response.

Sepsis, the dysregulated host response to infection causing life-threatening organ dysfunction, is a global health challenge requiring better understanding of pathophysiology and new therapeutic approaches. Here, we applied high-throughput tandem mass spectrometry to delineate the plasma proteome for sepsis and comparator groups (noninfected critical illness, postoperative inflammation, and healthy volunteers) involving 2612 samples (from 1611 patients) and 4553 liquid chromatography-mass spectrometry analyses acquired through a single batch of continuous measurements, with a throughput of 100 samples per day. We show how this scale of data can delineate proteins, pathways, and coexpression modules in sepsis and be integrated with paired leukocyte transcriptomic data (837 samples from n = 649 patients). We mapped the plasma proteomic landscape of the host response in sepsis, including changes over time, and identified features relating to etiology, clinical phenotypes (including organ failures), and severity. This work reveals subphenotypes informative for sepsis response state, disease processes, and outcome; identifies potential biomarkers; and advances opportunities for a precision medicine approach to sepsis.

Quantitative proteome-wide O-glycoproteomics analysis with FragPipe.

Identification of O-glycopeptides from tandem mass spectrometry data is complicated by the near complete dissociation of O-glycans from the peptide during collisional activation and by the combinatorial explosion of possible glycoforms when glycans are retained intact in electron-based activation. The recent O-Pair search method provides an elegant solution to these problems, using a collisional activation scan to identify the peptide sequence and total glycan mass, and a follow-up electron-based activation scan to localize the glycosite(s) using a graph-based algorithm in a reduced search space. Our previous O-glycoproteomics methods with MSFragger-Glyco allowed for extremely fast and sensitive identification of O-glycopeptides from collisional activation data but had limited support for site localization of glycans and quantification of glycopeptides. Here, we report an improved pipeline for O-glycoproteomics analysis that provides proteome-wide, site-specific, quantitative results by incorporating the O-Pair method as a module within FragPipe. In addition to improved search speed and sensitivity, we add flexible options for oxonium ion-based filtering of glycans and support for a variety of MS acquisition methods and provide a comparison between all software tools currently capable of O-glycosite localization in proteome-wide searches.

diaTracer enables spectrum-centric analysis of diaPASEF proteomics data.

Data-independent acquisition (DIA) has become a widely used strategy for peptide and protein quantification in mass spectrometry-based proteomics studies. The integration of ion mobility separation into DIA analysis, such as the diaPASEF technology available on Bruker's timsTOF platform, further improves the quantification accuracy and protein depth achievable using DIA. We introduce diaTracer, a new spectrum-centric computational tool optimized for diaPASEF data. diaTracer performs three-dimensional (m/z, retention time, ion mobility) peak tracing and feature detection to generate precursor-resolved "pseudo-MS/MS" spectra, facilitating direct ("spectral-library free") peptide identification and quantification from diaPASEF data. diaTracer is available as a stand-alone tool and is fully integrated into the widely used FragPipe computational platform. We demonstrate the performance of diaTracer and FragPipe using diaPASEF data from cerebrospinal fluid (CSF) and plasma samples, data from phosphoproteomics and HLA immunopeptidomics experiments, and low-input data from a spatial proteomics study. We also show that diaTracer enables unrestricted identification of post-translational modifications from diaPASEF data using open/mass offset searches.

An MSstats workflow for detecting differentially abundant proteins in large-scale data-independent acquisition mass spectrometry experiments with FragPipe processing.

Technological advances in mass spectrometry and proteomics have made it possible to perform larger-scale and more-complex experiments. The volume and complexity of the resulting data create major challenges for downstream analysis. In particular, next-generation data-independent acquisition (DIA) experiments enable wider proteome coverage than more traditional targeted approaches but require computational workflows that can manage much larger datasets and identify peptide sequences from complex and overlapping spectral features. Data-processing tools such as FragPipe, DIA-NN and Spectronaut have undergone substantial improvements to process spectral features in a reasonable time. Statistical analysis tools are needed to draw meaningful comparisons between experimental samples, but these tools were also originally designed with smaller datasets in mind. This protocol describes an updated version of MSstats that has been adapted to be compatible with large-scale DIA experiments. A very large DIA experiment, processed with FragPipe, is used as an example to demonstrate different MSstats workflows. The choice of workflow depends on the user's computational resources. For datasets that are too large to fit into a standard computer's memory, we demonstrate the use of MSstatsBig, a companion R package to MSstats. The protocol also highlights key decisions that have a major effect on both the results and the processing time of the analysis. The MSstats processing can be expected to take 1-3 h depending on the usage of MSstatsBig. The protocol can be run in the point-and-click graphical user interface MSstatsShiny or implemented with minimal coding expertise in R.

Integrated stress response plasticity governs normal cell adaptation to chronic stress via the PP2A-TFE3-ATF4 pathway.

The integrated stress response (ISR) regulates cell fate during conditions of stress by leveraging the cell's capacity to endure sustainable and efficient adaptive stress responses. Protein phosphatase 2A (PP2A) activity modulation has been shown to be successful in achieving both therapeutic efficacy and safety across various cancer models; however, the molecular mechanisms driving its selective antitumor effects remain unclear. Here, we show for the first time that ISR plasticity relies on PP2A activation to regulate drug response and dictate cellular fate under conditions of chronic stress. We demonstrate that genetic and chemical modulation of the PP2A leads to chronic proteolytic stress and triggers an ISR to dictate cell fate. More specifically, we uncovered that the PP2A-TFE3-ATF4 pathway governs ISR cell plasticity during endoplasmic reticular and cellular stress independent of the unfolded protein response. We further show that normal cells reprogram their genetic signatures to undergo ISR-mediated adaptation and homeostatic recovery thereby successfully avoiding toxicity following PP2A-mediated stress. Conversely, oncogenic specific cytotoxicity induced by chemical modulation of PP2A is achieved by activating chronic and irreversible ISR in cancer cells. Our findings propose that a differential response to chemical modulation of PP2A is determined by intrinsic ISR plasticity, providing a novel biological vulnerability to selectively induce cancer cell death and improve targeted therapeutic efficacy.

Immunopeptidomics-based identification of naturally presented non-canonical circRNA-derived peptides.

Circular RNAs (circRNAs) are covalently closed non-coding RNAs lacking the 5' cap and the poly-A tail. Nevertheless, it has been demonstrated that certain circRNAs can undergo active translation. Therefore, aberrantly expressed circRNAs in human cancers could be an unexplored source of tumor-specific antigens, potentially mediating anti-tumor T cell responses. This study presents an immunopeptidomics workflow with a specific focus on generating a circRNA-specific protein fasta reference. The main goal of this workflow is to streamline the process of identifying and validating human leukocyte antigen (HLA) bound peptides potentially originating from circRNAs. We increase the analytical stringency of our workflow by retaining peptides identified independently by two mass spectrometry search engines and/or by applying a group-specific FDR for canonical-derived and circRNA-derived peptides. A subset of circRNA-derived peptides specifically encoded by the region spanning the back-splice junction (BSJ) are validated with targeted MS, and with direct Sanger sequencing of the respective source transcripts. Our workflow identifies 54 unique BSJ-spanning circRNA-derived peptides in the immunopeptidome of melanoma and lung cancer samples. Our approach enlarges the catalog of source proteins that can be explored for immunotherapy.

Analysis and visualization of quantitative proteomics data using FragPipe-Analyst.

The FragPipe computational proteomics platform is gaining widespread popularity among the proteomics research community because of its fast processing speed and user-friendly graphical interface. Although FragPipe produces well-formatted output tables that are ready for analysis, there is still a need for an easy-to-use and user-friendly downstream statistical analysis and visualization tool. FragPipe-Analyst addresses this need by providing an R shiny web server to assist FragPipe users in conducting downstream analyses of the resulting quantitative proteomics data. It supports major quantification workflows including label-free quantification, tandem mass tags, and data-independent acquisition. FragPipe-Analyst offers a range of useful functionalities, such as various missing value imputation options, data quality control, unsupervised clustering, differential expression (DE) analysis using Limma, and gene ontology and pathway enrichment analysis using Enrichr. To support advanced analysis and customized visualizations, we also developed FragPipeAnalystR, an R package encompassing all FragPipe-Analyst functionalities that is extended to support site-specific analysis of post-translational modifications (PTMs). FragPipe-Analyst and FragPipeAnalystR are both open-source and freely available.

Proteome-scale tissue mapping using mass spectrometry based on label-free and multiplexed workflows.

Multiplexed bimolecular profiling of tissue microenvironment, or spatial omics, can provide deep insight into cellular compositions and interactions in both normal and diseased tissues. Proteome-scale tissue mapping, which aims to unbiasedly visualize all the proteins in whole tissue section or region of interest, has attracted significant interest because it holds great potential to directly reveal diagnostic biomarkers and therapeutic targets. While many approaches are available, however, proteome mapping still exhibits significant technical challenges in both protein coverage and analytical throughput. Since many of these existing challenges are associated with mass spectrometry-based protein identification and quantification, we performed a detailed benchmarking study of three protein quantification methods for spatial proteome mapping, including label-free, TMT-MS2, and TMT-MS3. Our study indicates label-free method provided the deepest coverages of ~3500 proteins at a spatial resolution of 50 µm and the largest quantification dynamic range, while TMT-MS2 method holds great benefit in mapping throughput at >125 pixels per day. The evaluation also indicates both label-free and TMT-MS2 provide robust protein quantifications in terms of identifying differentially abundant proteins and spatially co-variable clusters. In the study of pancreatic islet microenvironment, we demonstrated deep proteome mapping not only enables to identify protein markers specific to different cell types, but more importantly, it also reveals unknown or hidden protein patterns by spatial co-expression analysis.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Integrated proteogenomic characterization of glioblastoma evolution.

The evolutionary trajectory of glioblastoma (GBM) is a multifaceted biological process that extends beyond genetic alterations alone. Here, we perform an integrative proteogenomic analysis of 123 longitudinal glioblastoma pairs and identify a highly proliferative cellular state at diagnosis and replacement by activation of neuronal transition and synaptogenic pathways in recurrent tumors. Proteomic and phosphoproteomic analyses reveal that the molecular transition to neuronal state at recurrence is marked by post-translational activation of the wingless-related integration site (WNT)/ planar cell polarity (PCP) signaling pathway and BRAF protein kinase. Consistently, multi-omic analysis of patient-derived xenograft (PDX) models mirror similar patterns of evolutionary trajectory. Inhibition of B-raf proto-oncogene (BRAF) kinase impairs both neuronal transition and migration capability of recurrent tumor cells, phenotypic hallmarks of post-therapy progression. Combinatorial treatment of temozolomide (TMZ) with BRAF inhibitor, vemurafenib, significantly extends the survival of PDX models. This study provides comprehensive insights into the biological mechanisms of glioblastoma evolution and treatment resistance, highlighting promising therapeutic strategies for clinical intervention.

Sulfoproteomics Workflow with Precursor Ion Accurate Mass Shift Analysis Reveals Novel Tyrosine Sulfoproteins in the Golgi.

Tyrosine sulfation in the Golgi of secreted and membrane proteins is an important post-translational modification (PTM). However, its labile nature has limited analysis by mass spectrometry (MS), a major reason why no sulfoproteome studies have been previously reported. Here, we show that a phosphoproteomics experimental workflow, which includes serial enrichment followed by high resolution, high mass accuracy MS, and tandem MS (MS/MS) analysis, enables sulfopeptide coenrichment and identification via accurate precursor ion mass shift open MSFragger database search. This approach, supported by manual validation, allows the confident identification of sulfotyrosine-containing peptides in the presence of high levels of phosphorylated peptides, thus enabling these two sterically and ionically similar isobaric PTMs to be distinguished and annotated in a single proteomic analysis. We applied this approach to isolated interphase and mitotic rat liver Golgi membranes and identified 67 tyrosine sulfopeptides, corresponding to 26 different proteins. This work discovered 23 new sulfoproteins with functions related to, for example, Ca2+-binding, glycan biosynthesis, and exocytosis. In addition, we report the first preliminary evidence for crosstalk between sulfation and phosphorylation in the Golgi, with implications for functional control.

Autoinhibited kinesin-1 adopts a hierarchical folding pattern.

Conventional kinesin-1 is the primary anterograde motor in cells for transporting cellular cargo. While there is a consensus that the C-terminal tail of kinesin-1 inhibits motility, the molecular architecture of a full-length autoinhibited kinesin-1 remains unknown. Here, we combine crosslinking mass spectrometry (XL-MS), electron microscopy (EM), and AlphaFold structure prediction to determine the architecture of the full-length autoinhibited kinesin-1 homodimer (kinesin-1 heavy chain [KHC]) and kinesin-1 heterotetramer (KHC bound to kinesin light chain 1 [KLC1]). Our integrative analysis shows that kinesin-1 forms a compact, bent conformation through a break in coiled-coil 3. Moreover, our XL-MS analysis demonstrates that kinesin light chains stabilize the folded inhibited state rather than inducing a new structural state. Using our structural model, we show that disruption of multiple interactions between the motor, stalk, and tail domains is required to activate the full-length kinesin-1. Our work offers a conceptual framework for understanding how cargo adaptors and microtubule-associated proteins relieve autoinhibition to promote activation.

Efficient Analysis of Proteome-Wide FPOP Data by FragPipe.

Monitoring protein structure before and after environmental alterations (e.g., different cell states) can give insights into the role and function of proteins. Fast photochemical oxidation of proteins (FPOP) coupled with mass spectrometry (MS) allows for monitoring of structural rearrangements by exposing proteins to OH radicals that oxidize solvent-accessible residues, indicating protein regions undergoing movement. Some of the benefits of FPOP include high throughput and a lack of scrambling due to label irreversibility. However, the challenges of processing FPOP data have thus far limited its proteome-scale uses. Here, we present a computational workflow for fast and sensitive analysis of FPOP data sets. Our workflow, implemented as part of the FragPipe computational platform, combines the speed of the MSFragger search with a unique hybrid search method to restrict the large search space of FPOP modifications. Together, these features enable more than 10-fold faster FPOP searches that identify 150% more modified peptide spectra than previous methods. We hope this new workflow will increase the accessibility of FPOP to enable more protein structure and function relationships to be explored.

Identification of LMAN1- and SURF4-Dependent Secretory Cargoes.

Most proteins secreted into the extracellular space are first recruited from the endoplasmic reticulum into coat protein complex II (COPII)-coated vesicles or tubules that facilitate their transport to the Golgi apparatus. Although several secreted proteins have been shown to be actively recruited into COPII vesicles and tubules by the cargo receptors LMAN1 and SURF4, the full cargo repertoire of these receptors is unknown. We now report mass spectrometry analysis of conditioned media and cell lysates from HuH7 cells CRISPR targeted to inactivate the LMAN1 or SURF4 gene. We found that LMAN1 has limited clients in HuH7 cells, whereas SURF4 traffics a broad range of cargoes. Analysis of putative SURF4 cargoes suggests that cargo recognition is governed by complex mechanisms rather than interaction with a universal binding motif..

The Intrinsically Disordered N Terminus in Atg12 from Yeast Is Necessary for the Functional Structure of the Protein.

The Atg12 protein in yeast is an indispensable polypeptide in the highly conserved ubiquitin-like conjugation system operating in the macroautophagy/autophagy pathway. Atg12 is covalently conjugated to Atg5 through the action of Atg7 and Atg10; the Atg12-Atg5 conjugate binds Atg16 to form an E3 ligase that functions in a separate conjugation pathway involving Atg8. Atg12 is comprised of a ubiquitin-like (UBL) domain preceded at the N terminus by an intrinsically disordered protein region (IDPR), a domain that comprises a major portion of the protein but remains elusive in its conformation and function. Here, we show that the IDPR in unconjugated Atg12 is positioned in proximity to the UBL domain, a configuration that is important for the functional structure of the protein. A major deletion in the IDPR disrupts intactness of the UBL domain at the unconjugated C terminus, and a mutation in the predicted α0 helix in the IDPR prevents Atg12 from binding to Atg7 and Atg10, which ultimately affects the protein function in the ubiquitin-like conjugation cascade. These findings provide evidence that the IDPR is an indispensable part of the Atg12 protein from yeast.

Solid-Phase Compatible Silane-Based Cleavable Linker Enables Custom Isobaric Quantitative Chemoproteomics.

Mass spectrometry-based chemoproteomics has emerged as an enabling technology for functional biology and drug discovery. To address limitations of established chemoproteomics workflows, including cumbersome reagent synthesis and low throughput sample preparation, here, we established the silane-based cleavable isotopically labeled proteomics (sCIP) method. The sCIP method is enabled by a high yielding and scalable route to dialkoxydiphenylsilane fluorenylmethyloxycarbonyl (DADPS-Fmoc)-protected amino acid building blocks, which enable the facile synthesis of customizable, isotopically labeled, and chemically cleavable biotin capture reagents. sCIP is compatible with both MS1- and MS2-based quantitation, and the sCIP-MS2 method is distinguished by its click-assembled isobaric tags in which the reporter group is encoded in the sCIP capture reagent and balancer in the pan cysteine-reactive probe. The sCIP-MS2 workflow streamlines sample preparation with early stage isobaric labeling and sample pooling, allowing for high coverage and increased sample throughput via customized low cost six-plex sample multiplexing. When paired with a custom FragPipe data analysis workflow and applied to cysteine-reactive fragment screens, sCIP proteomics revealed established and unprecedented cysteine-ligand pairs, including the discovery that mitochondrial uncoupling agent FCCP acts as a covalent-reversible cysteine-reactive electrophile.

Mass Spectrometric Profiling of HLA-B44 Peptidomes Provides Evidence for Tapasin-Mediated Tryptophan Editing.

The extreme polymorphisms of HLA class I proteins result in structural variations in their peptide binding sites to achieve diversity in Ag presentation. External factors could independently constrict or alter HLA class I peptide repertoires. Such effects of the assembly factor tapasin were assessed for HLA-B*44:05 (Y116) and a close variant, HLA-B*44:02 (D116), which have low and high tapasin dependence, respectively, for their cell surface expression. Analyses of the HLA-B*44:05 peptidomes in the presence and absence of tapasin reveal that peptides with C-terminal tryptophans and higher predicted affinities are preferentially selected by tapasin, coincident with reduced frequencies of peptides with other C-terminal amino acids, including leucine. Comparisons of the HLA-B*44:05 and HLA-B*44:02 peptidomes indicate the expected structure-based alterations near the peptide C termini, but also C-terminal amino acid frequency and predicted affinity changes among the unique and shared peptide groups for B*44:02 and B*44:05. Overall, these findings indicate that the presence of tapasin and the tapasin dependence of assembly alter HLA class I peptide-binding preferences at the peptide C terminus. The particular C-terminal amino acid preferences that are altered by tapasin are expected to be determined by the intrinsic peptide-binding specificities of HLA class I allotypes. Additionally, the findings suggest that tapasin deficiency and reduced tapasin dependence expand the permissive affinities of HLA class I-bound peptides, consistent with prior findings that HLA class I allotypes with low tapasin dependence have increased breadth of CD8+ T cell epitope presentation and are more protective in HIV infections.

Multi-omic stratification of the missense variant cysteinome.

Cancer genomes are rife with genetic variants; one key outcome of this variation is gain-ofcysteine, which is the most frequently acquired amino acid due to missense variants in COSMIC. Acquired cysteines are both driver mutations and sites targeted by precision therapies. However, despite their ubiquity, nearly all acquired cysteines remain uncharacterized. Here, we pair cysteine chemoproteomics-a technique that enables proteome-wide pinpointing of functional, redox sensitive, and potentially druggable residues-with genomics to reveal the hidden landscape of cysteine acquisition. For both cancer and healthy genomes, we find that cysteine acquisition is a ubiquitous consequence of genetic variation that is further elevated in the context of decreased DNA repair. Our chemoproteogenomics platform integrates chemoproteomic, whole exome, and RNA-seq data, with a customized 2-stage false discovery rate (FDR) error controlled proteomic search, further enhanced with a user-friendly FragPipe interface. Integration of CADD predictions of deleteriousness revealed marked enrichment for likely damaging variants that result in acquisition of cysteine. By deploying chemoproteogenomics across eleven cell lines, we identify 116 gain-of-cysteines, of which 10 were liganded by electrophilic druglike molecules. Reference cysteines proximal to missense variants were also found to be pervasive, 791 in total, supporting heretofore untapped opportunities for proteoform-specific chemical probe development campaigns. As chemoproteogenomics is further distinguished by sample-matched combinatorial variant databases and compatible with redox proteomics and small molecule screening, we expect widespread utility in guiding proteoform-specific biology and therapeutic discovery.