Identification of phosphosites that alter protein thermal stability

Profiling Proteins and Phosphorylation Sites During T Cell Activation Using an Integrated Thermal Shift Assay

T cell activation is a complex biological process of naive cells maturing into effector cells. Proteomic and phospho-proteomic approaches have provided critical insights into this process, yet it is not always clear how changes in individual proteins or phosphorylation sites have functional significance. Here, we developed the Phosphorylation Integrated Thermal Shift Assay (PITSA) that combines the measurement of protein or phosphorylation site abundance and thermal stability into a single tandem mass tags experiment and apply this method to study T cell activation. We quantified the abundance and thermal stability of over 7500 proteins and 5000 phosphorylation sites and identified significant differences in chromatin-related, TCR signaling, DNA repair, and proliferative phosphoproteins. PITSA may be applied to a wide range of biological contexts to generate hypotheses as to which proteins or phosphorylation sites are functionally regulated in a given system as well as the mechanisms by which this regulation may occur.

Discovery of metal-binding proteins by thermal proteome profiling

Metal-binding proteins (MBPs) have various and important biological roles in all living species and many human diseases are intricately linked to dysfunctional MBPs. Here, we report a chemoproteomic method named 'metal extraction-triggered agitation logged by thermal proteome profiling' (METAL-TPP) to globally profile MBPs in proteomes. The method involves the extraction of metals from MBPs using chelators and monitoring the resulting protein stability changes through thermal proteome profiling. Applying METAL-TPP to the human proteome with a broad-spectrum chelator, EDTA, revealed a group of proteins with reduced thermal stability that contained both previously known MBPs and currently unannotated MBP candidates. Biochemical characterization of one potential target, glutamine-fructose-6-phosphate transaminase 2 (GFPT2), showed that zinc bound the protein, inhibited its enzymatic activity and modulated the hexosamine biosynthesis pathway. METAL-TPP profiling with another chelator, TPEN, uncovered additional MBPs in proteomes. Collectively, this study developed a robust tool for proteomic discovery of MBPs and provides a rich resource for functional studies of metals in cell biology.

Poly-pharmacology of existing drugs: How to crack the code?

Drug development in oncology is highly challenging, with less than 5% success rate in clinical trials. This alarming figure points out the need to study in more details the multiple biological effects of drugs in specific contexts. Indeed, the comprehensive assessment of drug poly-pharmacology can provide insights into their therapeutic and adverse effects, to optimize their utilization and maximize the success rate of clinical trials. Recent technological advances have made possible in-depth investigation of drug poly-pharmacology. This review first highlights high-throughput methodologies that have been used to unveil new mechanisms of action of existing drugs. Then, we discuss how emerging chemo-proteomics strategies allow effectively dissecting the poly-pharmacology of drugs in an unsupervised manner.

Stability-based approaches in chemoproteomics

Target deconvolution can help understand how compounds exert therapeutic effects and can accelerate drug discovery by helping optimise safety and efficacy, revealing mechanisms of action, anticipate off-target effects and identifying opportunities for therapeutic expansion. Chemoproteomics, a combination of chemical biology with mass spectrometry has transformed target deconvolution. This review discusses modification-free chemoproteomic approaches that leverage the change in protein thermodynamics induced by small molecule ligand binding. Unlike modification-based methods relying on enriching specific protein targets, these approaches offer proteome-wide evaluations, driven by advancements in mass spectrometry sensitivity, increasing proteome coverage and quantitation methods. Advances in methods based on denaturation/precipitation by thermal or chemical denaturation, or by protease degradation are evaluated, emphasising the evolving landscape of chemoproteomics and its potential impact on future drug-development strategies.

Experimental and data analysis advances in thermal proteome profiling

Method development for mass spectrometry (MS)-based thermal shift proteomic assays have advanced to probe small molecules with known and unknown protein-ligand interaction mechanisms and specificity, which is predominantly used in characterization of drug-protein interactions. In the discovery of target and off-target protein-ligand interactions, a thorough investigation of method development and their impact on the sensitivity and accuracy of protein-small molecule and protein-protein interactions is warranted. In this review, we discuss areas of improvement at each stage of thermal proteome profiling data analysis that includes processing of MS-based data, method development, and their effect on the overall quality of thermal proteome profiles. We also overview the optimization of experimental strategies and prioritization of an increased number of independent biological replicates over the number of evaluated temperatures.

Scaled-Down Thermal Profiling and Coaggregation Analysis of the Proteome for Drug Target and Protein Interaction Analysis

Thermal proteome profiling (TPP), an experimental technique combining the cellular thermal shift assay (CETSA) with quantitative protein mass spectrometry (MS), identifies interactions of drugs and chemicals with endogenous proteins. Thermal proximity coaggregation (TPCA) profiling extended TPP to study the intracellular dynamics of protein complexes. In TPP and TPCA, samples are subjected to multiple denaturing temperatures, each requiring over 100 μg of proteins, which restricts their applications for rare cells and precious clinical samples. We developed a workflow termed STASIS (scaled-down thermal profiling and coaggregation analysis with SISPROT) that scales down the required protein to as low as 1 μg per temperature. This is achieved by heating and centrifugation using the same PCR tube, processing samples with the SISPROT technology (simple and integrated spintip-based proteomics technology), and tip-based manual fractionation of TMT-labeled peptides. We evaluate the STASIS workflow with starting protein quantities of 10, 5, and 1 μg per temperature prior to heating, identifying between 4000 and 5000 proteins with 6 h of acquisition time. Importantly, we observed a high correlation in the Tm of proteins with minimal difference in TPCA performance for predicting protein complexes. Moreover, STASIS could identify the targets of methotrexate and panobinostat with high precision with 1 μg of proteins per temperature. In conclusion, STASIS is a robust cost-effective technique for target deconvolution and extended TPCA to rare primary cells and precious clinical samples for the analysis of protein complexes.

Obtaining Functional Proteomics Insights From Thermal Proteome Profiling Through Optimized Melt Shift Calculation and Statistical Analysis With InflectSSP

Thermal proteome profiling (TPP) is an invaluable tool for functional proteomics studies that has been shown to discover changes associated with protein-ligand, protein-protein, and protein-RNA interaction dynamics along with changes in protein stability resulting from cellular signaling. The increasing number of reports employing this assay has not been met concomitantly with new approaches leading to advancements in the quality and sensitivity of the corresponding data analysis. The gap between data acquisition and data analysis tools is important to fill as TPP findings have reported subtle melt shift changes related to signaling events such as protein posttranslational modifications. In this study, we have improved the Inflect data analysis pipeline (now referred to as InflectSSP, available at https://CRAN.R-project.org/package=InflectSSP) to increase the sensitivity of detection for both large and subtle changes in the proteome as measured by TPP. Specifically, InflectSSP now has integrated statistical and bioinformatic functions to improve objective functional proteomics findings from the quantitative results obtained from TPP studies through increasing both the sensitivity and specificity of the data analysis pipeline. InflectSSP incorporates calculation of a "melt coefficient" into the pipeline with production of average melt curves for biological replicate studies to aid in identification of proteins with significant melts. To benchmark InflectSSP, we have reanalyzed two previously reported datasets to demonstrate the performance of our publicly available R-based program for TPP data analysis. We report new findings following temporal treatment of human cells with the small molecule thapsigargin that induces the unfolded protein response as a consequence of inhibition of sarcoplasmic/endoplasmic reticulum calcium ATPase 2A. InflectSSP analysis of our unfolded protein response study revealed highly reproducible and statistically significant target engagement over a time course of treatment while simultaneously providing new insights into the possible mechanisms of action of the small molecule thapsigargin.

Toward the analysis of functional proteoforms using mass spectrometry-based stability proteomics

Functional proteomics aims to elucidate biological functions, mechanisms, and pathways of proteins and proteoforms at the molecular level to examine complex cellular systems and disease states. A series of stability proteomics methods have been developed to examine protein functionality by measuring the resistance of a protein to chemical or thermal denaturation or proteolysis. These methods can be applied to measure the thermal stability of thousands of proteins in complex biological samples such as cell lysate, intact cells, tissues, and other biological fluids to measure proteome stability. Stability proteomics methods have been popularly applied to observe stability shifts upon ligand binding for drug target identification. More recently, these methods have been applied to characterize the effect of structural changes in proteins such as those caused by post-translational modifications (PTMs) and mutations, which can affect protein structures or interactions and diversify protein functions. Here, we discussed the current application of a suite of stability proteomics methods, including thermal proteome profiling (TPP), stability of proteomics from rates of oxidation (SPROX), and limited proteolysis (LiP) methods, to observe PTM-induced structural changes on protein stability. We also discuss future perspectives highlighting the integration of top-down mass spectrometry and stability proteomics methods to characterize intact proteoform stability and understand the function of variable protein modifications.

Decoding protein methylation function with thermal stability analysis

Protein methylation is an important modification beyond epigenetics. However, systems analyses of protein methylation lag behind compared to other modifications. Recently, thermal stability analyses have been developed which provide a proxy of a protein functional status. Here, we show that molecular and functional events closely linked to protein methylation can be revealed by the analysis of thermal stability. Using mouse embryonic stem cells as a model, we show that Prmt5 regulates mRNA binding proteins that are enriched in intrinsically disordered regions and involved in liquid-liquid phase separation mechanisms, including the formation of stress granules. Moreover, we reveal a non-canonical function of Ezh2 in mitotic chromosomes and the perichromosomal layer, and identify Mki67 as a putative Ezh2 substrate. Our approach provides an opportunity to systematically explore protein methylation function and represents a rich resource for understanding its role in pluripotency.

Recent advances in kinase signaling network profiling by mass spectrometry

Mass spectrometry-based phosphoproteomics is currently the leading methodology for the study of global kinase signaling. The scientific community is continuously releasing technological improvements for sensitive and fast identification of phosphopeptides, and their accurate quantification. To interpret large-scale phosphoproteomics data, numerous bioinformatic resources are available that help understanding kinase network functional role in biological systems upon perturbation. Some of these resources are databases of phosphorylation sites, protein kinases and phosphatases; others are bioinformatic algorithms to infer kinase activity, predict phosphosite functional relevance and visualize kinase signaling networks. In this review, we present the latest experimental and bioinformatic tools to profile protein kinase signaling networks and provide examples of their application in biomedicine.

Protein-Peptide Turnover Profiling reveals the order of PTM addition and removal during protein maturation

Post-translational modifications (PTMs) regulate various aspects of protein function, including degradation. Mass spectrometric methods relying on pulsed metabolic labeling are popular to quantify turnover rates on a proteome-wide scale. Such data have traditionally been interpreted in the context of protein proteolytic stability. Here, we combine theoretical kinetic modeling with experimental pulsed stable isotope labeling of amino acids in cell culture (pSILAC) for the study of protein phosphorylation. We demonstrate that metabolic labeling combined with PTM-specific enrichment does not measure effects of PTMs on protein stability. Rather, it reveals the relative order of PTM addition and removal along a protein's lifetime-a fundamentally different metric. This is due to interconversion of the measured proteoform species. Using this framework, we identify temporal phosphorylation sites on cell cycle-specific factors and protein complex assembly intermediates. Our results thus allow tying PTMs to the age of the modified proteins.

Thermal proteome profiling: Insights into protein modifications, associations, and functions

Tracking proteins' biophysical characteristics on a proteome-wide scale can provide valuable information on their functions and interactions. Thermal proteome profiling (TPP) is a multiplexed quantitative proteomics approach that measures changes in protein thermal stability-a key biophysical property-across different cellular states. Developed in 2014, as a target-deconvolution assay for drugs and other small molecules, TPP has since evolved to a system-level biochemical omics technique providing insights into context-dependent changes in protein states. In this review, we summarise key advances in the experimental and data analysis pipeline that have aided this transformation and discuss the recent developments and applications of TPP.

Rewiring of the protein-protein-metabolite interactome during the diauxic shift in yeast

In budding yeast Saccharomyces cerevisiae, the switch from aerobic fermentation to respiratory growth is separated by a period of growth arrest, known as the diauxic shift, accompanied by a significant metabolic rewiring, including the derepression of gluconeogenesis and the establishment of mitochondrial respiration. Previous studies reported hundreds of proteins and tens of metabolites accumulating differentially across the diauxic shift transition. To assess the differences in the protein-protein (PPIs) and protein-metabolite interactions (PMIs) yeast samples harvested in the glucose-utilizing, fermentative phase, ethanol-utilizing and early stationary respiratory phases were analysed using isothermal shift assay (iTSA) and a co-fractionation mass spectrometry approach, PROMIS. Whereas iTSA monitors changes in protein stability and is informative towards protein interaction status, PROMIS uses co-elution to delineate putative PPIs and PMIs. The resulting dataset comprises 1627 proteins and 247 metabolites, hundreds of proteins and tens of metabolites characterized by differential thermal stability and/or fractionation profile, constituting a novel resource to be mined for the regulatory PPIs and PMIs. The examples discussed here include (i) dissociation of the core and regulatory particle of the proteasome in the early stationary phase, (ii) the differential binding of a co-factor pyridoxal phosphate to the enzymes of amino acid metabolism and (iii) the putative, phase-specific interactions between proline-containing dipeptides and enzymes of central carbon metabolism.

IsobaricQuant enables cross-platform quantification, visualization, and filtering of isobarically-labeled peptides

In mass spectrometry (MS)-based quantitative proteomics, labeling with isobaric mass tags such as iTRAQ and TMT can substantially improve sample throughput and reduce peptide missing values. Nonetheless, the quantification of labeled peptides tends to suffer from reduced accuracy due to the co-isolation of co-eluting precursors of similar mass-to-charge. Acquisition approaches such as multistage MS3 or ion mobility separation address this problem, yet are difficult to audit and limited to expensive instrumentation. Here we introduce IsobaricQuant, an open-source software tool for quantification, visualization, and filtering of peptides labeled with isobaric mass tags, with specific focus on precursor interference. IsobaricQuant is compatible with MS2 and MS3 acquisition strategies, has a viewer that allows assessing interference, and provides several scores to aid the filtering of scans with compression. We demonstrate that IsobaricQuant quantifications are accurate by comparing it with commonly used software. We further show that its QC scores can successfully filter out scans with reduced quantitative accuracy at MS2 and MS3 levels, removing inaccurate peptide quantifications and decreasing protein CVs. Finally, we apply IsobaricQuant to a PISA dataset and show that QC scores improve the sensitivity of the identification of protein targets of a kinase inhibitor. IsobaricQuant is available at https://github.com/Villen-Lab/isobaricquant.

The evolution of post-translational modifications

Post-translational modifications (PTMs) are chemical modifications that can regulate the activity and function of proteins. From an evolutionary perspective, they also represent a fast mechanism for the generation of phenotypic diversity and divergence. Advances in mass spectrometry have now enabled the identification of over 600 distinct PTM classes collectively spanning an order of 10⁶ unique sites. However, the chemical detection of PTMs has lagged far behind their functional characterisation, and relatively little is still known about the selective constraints that govern PTM evolution. In particular, the true fraction of PTM sites that are functional - and thus subject to selection - remains an open question. Here, I review advances made in the past two years towards understanding the evolution of PTMs and their associated enzymes.

Small Molecule Arranged Thermal Proximity Coaggregation (smarTPCA)-A Novel Approach to Characterize Protein-Protein Interactions in Living Cells by Similar Isothermal Dose-Responses

Chemical biology and the application of small molecules has proven to be a potent perturbation strategy, especially for the functional elucidation of proteins, their networks, and regulators. In recent years, the cellular thermal shift assay (CETSA) and its proteome-wide extension, thermal proteome profiling (TPP), have proven to be effective tools for identifying interactions of small molecules with their target proteins, as well as off-targets in living cells. Here, we asked the question whether isothermal dose-response (ITDR) CETSA can be exploited to characterize secondary effects downstream of the primary binding event, such as changes in post-translational modifications or protein-protein interactions (PPI). By applying ITDR-CETSA to MAPK14 kinase inhibitor treatment of living HL-60 cells, we found similar dose-responses for the direct inhibitor target and its known interaction partners MAPKAPK2 and MAPKAPK3. Extension of the dose-response similarity comparison to the proteome wide level using TPP with compound concentration range (TPP-CCR) revealed not only the known MAPK14 interaction partners MAPKAPK2 and MAPKAPK3, but also the potentially new intracellular interaction partner MYLK. We are confident that dose-dependent small molecule treatment in combination with ITDR-CETSA or TPP-CCR similarity assessment will not only allow discrimination between primary and secondary effects, but will also provide a novel method to study PPI in living cells without perturbation by protein modification, which we named "small molecule arranged thermal proximity coaggregation" (smarTPCA).

Thermal Proteome Profiling Reveals the O-GlcNAc-Dependent Meltome

Posttranslational modifications alter the biophysical properties of proteins and thereby influence cellular physiology. One emerging manner by which such modifications regulate protein functions is through their ability to perturb protein stability. Despite the increasing interest in this phenomenon, there are few methods that enable global interrogation of the biophysical effects of posttranslational modifications on the proteome. Here, we describe an unbiased proteome-wide approach to explore the influence of protein modifications on the thermodynamic stability of thousands of proteins in parallel. We apply this profiling strategy to study the effects of O-linked N-acetylglucosamine (O-GlcNAc), an abundant modification found on hundreds of proteins in mammals that has been shown in select cases to stabilize proteins. Using this thermal proteomic profiling strategy, we identify a set of 72 proteins displaying O-GlcNAc-dependent thermostability and validate this approach using orthogonal methods targeting specific proteins. These collective observations reveal that the majority of proteins influenced by O-GlcNAc are, surprisingly, destabilized by O-GlcNAc and cluster into distinct macromolecular complexes. These results establish O-GlcNAc as a bidirectional regulator of protein stability and provide a blueprint for exploring the impact of any protein modification on the meltome of, in principle, any organism.

Proteome-wide cellular thermal shift assay reveals unexpected cross-talk between brassinosteroid and auxin signaling

Chemical genetics, which investigates biological processes using small molecules, is gaining interest in plant research. However, a major challenge is to uncover the mode of action of the small molecules. Here, we applied the cellular thermal shift assay coupled with mass spectrometry (CETSA MS) to intact Arabidopsis cells and showed that bikinin, the plant-specific glycogen synthase kinase 3 (GSK3) inhibitor, changed the thermal stability of some of its direct targets and putative GSK3-interacting proteins. In combination with phosphoproteomics, we also revealed that GSK3s phosphorylated the auxin carrier PIN-FORMED1 and regulated its polarity that is required for the vascular patterning in the leaf.

Despite the growing interest in using chemical genetics in plant research, small molecule target identification remains a major challenge. The cellular thermal shift assay coupled with high-resolution mass spectrometry (CETSA MS) that monitors changes in the thermal stability of proteins caused by their interactions with small molecules, other proteins, or posttranslational modifications, allows the discovery of drug targets or the study of protein–metabolite and protein–protein interactions mainly in mammalian cells. To showcase the applicability of this method in plants, we applied CETSA MS to intact Arabidopsis thaliana cells and identified the thermal proteome of the plant-specific glycogen synthase kinase 3 (GSK3) inhibitor, bikinin. A comparison between the thermal and the phosphoproteomes of bikinin revealed the auxin efflux carrier PIN-FORMED1 (PIN1) as a substrate of the Arabidopsis GSK3s that negatively regulate the brassinosteroid signaling. We established that PIN1 phosphorylation by the GSK3s is essential for maintaining its intracellular polarity that is required for auxin-mediated regulation of vascular patterning in the leaf, thus revealing cross-talk between brassinosteroid and auxin signaling.

ProSAP: a GUI software tool for statistical analysis and assessment of thermal stability data

The Cellular Thermal Shift Assay (CETSA) plays an important role in drug-target identification, and statistical analysis is a crucial step significantly affecting conclusion. We put forward ProSAP (Protein Stability Analysis Pod), an open-source, cross-platform and user-friendly software tool, which provides multiple methods for thermal proteome profiling (TPP) analysis, nonparametric analysis (NPA), proteome integral solubility alteration and isothermal shift assay (iTSA). For testing the performance of ProSAP, we processed several datasets and compare the performance of different algorithms. Overall, TPP analysis is more accurate with fewer false positive targets, but NPA methods are flexible and free from parameters. For iTSA, edgeR and DESeq2 identify more true targets than t-test and Limma, but when it comes to ranking, the four methods show not much difference. ProSAP software is available at https://github.com/hcji/ProSAP and https://zenodo.org/record/5763315.

A peptidoform based proteomic strategy for studying functions of post-translational modifications

Protein post-translational modifications (PTMs) generate an enormous, but as yet undetermined, expansion of the produced proteoforms. In this Viewpoint, we firstly reviewed the concepts of proteoform and peptidoform. We show that many of the current PTM biological investigation and annotation studies largely follow a PTM site-specific rather than proteoform-specific approach. We further illustrate a potentially useful matching strategy in which a particular "modified peptidoform" is matched to the corresponding "unmodified peptidoform" as a reference for the quantitative analysis between samples and conditions. We suggest this strategy has the potential to provide more directly relevant information to learn the PTM site-specific biological functions. Accordingly, we advocate for the wider use of the nomenclature "peptidoform" in future bottom-up proteomic studies.

Strategies for mass spectrometry-based phosphoproteomics using isobaric tagging

INTRODUCTION

Protein phosphorylation is a primary mechanism of signal transduction in cellular systems. Isobaric tagging can be used to investigate alterations in phosphorylation events in sample multiplexing experiments where quantification extends across all conditions. As such, innovations in tandem mass tag methods can facilitate the expansion of the depth and breadth of phosphoproteomic analyses.

AREAS COVERED

This review discusses the current state of tandem mass tag-centric phosphoproteomics and highlights advances in reagent chemistry, instrumentation, data acquisition, and data analysis. We stress that approaches for phosphoproteomic investigations require high-specificity enrichment, sensitive detection, and accurate phosphorylation site localization.

EXPERT OPINION

Tandem mass tag-centric phosphoproteomics will continue to be an important conduit for our understanding of signal transduction in living organisms. We anticipate that progress in phosphopeptide enrichment methodologies, enhancements in instrumentation and data acquisition technologies, and further refinements in analytical strategies will be key to the discovery of biologically relevant findings from phosphoproteomics studies.

Drug Target Identification in Tissues by Thermal Proteome Profiling

Drug target deconvolution can accelerate the drug discovery process by identifying a drug's targets (facilitating medicinal chemistry efforts) and off-targets (anticipating toxicity effects or adverse drug reactions). Multiple mass spectrometry-based approaches have been developed for this purpose, but thermal proteome profiling (TPP) remains to date the only one that does not require compound modification and can be used to identify intracellular targets in living cells. TPP is based on the principle that the thermal stability of a protein can be affected by its interactions. Recent developments of this approach have expanded its applications beyond drugs and cell cultures to studying protein-drug interactions and biological phenomena in tissues. These developments open up the possibility of studying drug treatment or mechanisms of disease in a holistic fashion, which can result in the design of better drugs and lead to a better understanding of fundamental biology. Expected final online publication date for the Annual Review of Pharmacology and Toxicology, Volume 62 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A Bayesian semi-parametric model for thermal proteome profiling

The thermal stability of proteins can be altered when they interact with small molecules, other biomolecules or are subject to post-translation modifications. Thus monitoring the thermal stability of proteins under various cellular perturbations can provide insights into protein function, as well as potentially determine drug targets and off-targets. Thermal proteome profiling is a highly multiplexed mass-spectrommetry method for monitoring the melting behaviour of thousands of proteins in a single experiment. In essence, thermal proteome profiling assumes that proteins denature upon heating and hence become insoluble. Thus, by tracking the relative solubility of proteins at sequentially increasing temperatures, one can report on the thermal stability of a protein. Standard thermodynamics predicts a sigmoidal relationship between temperature and relative solubility and this is the basis of current robust statistical procedures. However, current methods do not model deviations from this behaviour and they do not quantify uncertainty in the melting profiles. To overcome these challenges, we propose the application of Bayesian functional data analysis tools which allow complex temperature-solubility behaviours. Our methods have improved sensitivity over the state-of-the art, identify new drug-protein associations and have less restrictive assumptions than current approaches. Our methods allows for comprehensive analysis of proteins that deviate from the predicted sigmoid behaviour and we uncover potentially biphasic phenomena with a series of published datasets.