Deep thermal profiling for detection of functional proteoform groups.

The complexity of the functional proteome extends considerably beyond the coding genome, resulting in millions of proteoforms. Investigation of proteoforms and their functional roles is important to understand cellular physiology and its deregulation in diseases but challenging to perform systematically. Here we applied thermal proteome profiling with deep peptide coverage to detect functional proteoform groups in acute lymphoblastic leukemia cell lines with different cytogenetic aberrations. We detected 15,846 proteoforms, capturing differently spliced, cleaved and post-translationally modified proteins expressed from 9,290 genes. We identified differential co-aggregation of proteoform pairs and established links to disease biology. Moreover, we systematically made use of measured biophysical proteoform states to find specific biomarkers of drug sensitivity. Our approach, thus, provides a powerful and unique tool for systematic detection and functional annotation of proteoform groups.

Systematic discovery of biomolecular condensate-specific protein phosphorylation.

Reversible protein phosphorylation is an important mechanism for regulating (dis)assembly of biomolecular condensates. However, condensate-specific phosphosites remain largely unknown, thereby limiting our understanding of the underlying mechanisms. Here, we combine solubility proteome profiling with phosphoproteomics to quantitatively map several hundred phosphosites enriched in either soluble or condensate-bound protein subpopulations, including a subset of phosphosites modulating protein-RNA interactions. We show that multi-phosphorylation of the C-terminal disordered segment of heteronuclear ribonucleoprotein A1 (HNRNPA1), a key RNA-splicing factor, reduces its ability to locate to nuclear clusters. For nucleophosmin 1 (NPM1), an essential nucleolar protein, we show that phosphorylation of S254 and S260 is crucial for lowering its partitioning to the nucleolus and additional phosphorylation of distal sites enhances its retention in the nucleoplasm. These phosphorylation events decrease RNA and protein interactions of NPM1 to regulate its condensation. Our dataset is a rich resource for systematically uncovering the phosphoregulation of biomolecular condensates.

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.

Transcriptional and Post-Transcriptional Polar Effects in Bacterial Gene Deletion Libraries.

Single-gene deletions can affect the expression levels of other genes in the same operon in bacterial genomes. Here, we used proteomics for 133 Escherichia coli gene deletion mutants and transcriptome sequencing (RNA-seq) data from 71 mutants to probe the extent of transcriptional and post-transcriptional effects of gene deletions in operons. Transcriptional effects were common on genes located downstream of the deletion and were consistent across all operon members, with nearly 40% of operons showing greater than 2-fold up- or downregulation. Surprisingly, we observed an additional post-transcriptional effect that leads to the downregulation of the gene located directly downstream of the targeted gene. This effect was correlated with their intergenic distance, despite the ribosome binding site of the gene downstream remaining intact during library construction. Overall, the data presented can guide future library construction and highlight the importance of follow-up experiments for assessing direct effects of single-gene deletions in operons. IMPORTANCE Single-gene deletion libraries have allowed genome-wide characterization of gene function and interactions. While each mutant intends to disrupt the function of a single gene, it can unintentionally target other genes, such as those located in the same operon as the deletion. The extent to which such polar effects occur in deletion libraries has not been assessed. In this work, we use proteomics and transcriptomics data to show that transcript level changes lead to nearly 40% of deletions in operons affecting the protein levels of genes located downstream by at least 2-fold. Furthermore, we observed a post-transcriptional effect on the gene located directly downstream of the deletion. These results can guide the design of future gene deletion libraries and emphasizes the importance of follow-up work when linking genotypes to phenotypes.

SARS-CoV-2 infection remodels the host protein thermal stability landscape.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a global threat to human health and has compromised economic stability. In addition to the development of an effective vaccine, it is imperative to understand how SARS-CoV-2 hijacks host cellular machineries on a system-wide scale so that potential host-directed therapies can be developed. In situ proteome-wide abundance and thermal stability measurements using thermal proteome profiling (TPP) can inform on global changes in protein activity. Here we adapted TPP to high biosafety conditions amenable to SARS-CoV-2 handling. We discovered pronounced temporal alterations in host protein thermostability during infection, which converged on cellular processes including cell cycle, microtubule and RNA splicing regulation. Pharmacological inhibition of host proteins displaying altered thermal stability or abundance during infection suppressed SARS-CoV-2 replication. Overall, this work serves as a framework for expanding TPP workflows to globally important human pathogens that require high biosafety containment and provides deeper resolution into the molecular changes induced by SARS-CoV-2 infection.

Cell surface thermal proteome profiling tracks perturbations and drug targets on the plasma membrane.

Numerous drugs and endogenous ligands bind to cell surface receptors leading to modulation of downstream signaling cascades and frequently to adaptation of the plasma membrane proteome. In-depth analysis of dynamic processes at the cell surface is challenging due to biochemical properties and low abundances of plasma membrane proteins. Here we introduce cell surface thermal proteome profiling for the comprehensive characterization of ligand-induced changes in protein abundances and thermal stabilities at the plasma membrane. We demonstrate drug binding to extracellular receptors and transporters, discover stimulation-dependent remodeling of T cell receptor complexes and describe a competition-based approach to measure target engagement of G-protein-coupled receptor antagonists. Remodeling of the plasma membrane proteome in response to treatment with the TGFB receptor inhibitor SB431542 leads to partial internalization of the monocarboxylate transporters MCT1/3 explaining the antimetastatic effects of the drug.

Rtpca: an R package for differential thermal proximity coaggregation analysis.

SUMMARY: Rtpca is an R package implementing methods for inferring protein-protein interactions (PPIs) based on thermal proteome profiling experiments of a single condition or in a differential setting via an approach called thermal proximity coaggregation. It offers user-friendly tools to explore datasets for their PPI predictive performance and easily integrates with available R packages. AVAILABILITY AND IMPLEMENTATION: Rtpca is available from Bioconductor (https://bioconductor.org/packages/Rtpca). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ShinyButchR: Interactive NMF-based decomposition workflow of genome-scale datasets.

Non-negative matrix factorization (NMF) has been widely used for the analysis of genomic data to perform feature extraction and signature identification due to the interpretability of the decomposed signatures. However, running a basic NMF analysis requires the installation of multiple tools and dependencies, along with a steep learning curve and computing time. To mitigate such obstacles, we developed ShinyButchR, a novel R/Shiny application that provides a complete NMF-based analysis workflow, allowing the user to perform matrix decomposition using NMF, feature extraction, interactive visualization, relevant signature identification, and association to biological and clinical variables. ShinyButchR builds upon the also novel R package ButchR, which provides new TensorFlow solvers for algorithms of the NMF family, functions for downstream analysis, a rational method to determine the optimal factorization rank and a novel feature selection strategy.

The functional proteome landscape of Escherichia coli.

Recent developments in high-throughput reverse genetics1,2 have revolutionized our ability to map gene function and interactions3-6. The power of these approaches depends on their ability to identify functionally associated genes, which elicit similar phenotypic changes across several perturbations (chemical, environmental or genetic) when knocked out7-9. However, owing to the large number of perturbations, these approaches have been limited to growth or morphological readouts10. Here we use a high-content biochemical readout, thermal proteome profiling11, to measure the proteome-wide protein abundance and thermal stability in response to 121 genetic perturbations in Escherichia coli. We show that thermal stability, and therefore the state and interactions of essential proteins, is commonly modulated, raising the possibility of studying a protein group that is particularly inaccessible to genetics. We find that functionally associated proteins have coordinated changes in abundance and thermal stability across perturbations, owing to their co-regulation and physical interactions (with proteins, metabolites or cofactors). Finally, we provide mechanistic insights into previously determined growth phenotypes12 that go beyond the deleted gene. These data represent a rich resource for inferring protein functions and interactions.

A computational method for detection of ligand-binding proteins from dose range thermal proteome profiles.

Detecting ligand-protein interactions in living cells is a fundamental challenge in molecular biology and drug research. Proteome-wide profiling of thermal stability as a function of ligand concentration promises to tackle this challenge. However, current data analysis strategies use preset thresholds that can lead to suboptimal sensitivity/specificity tradeoffs and limited comparability across datasets. Here, we present a method based on statistical hypothesis testing on curves, which provides control of the false discovery rate. We apply it to several datasets probing epigenetic drugs and a metabolite. This leads us to detect off-target drug engagement, including the finding that the HDAC8 inhibitor PCI-34051 and its analog BRD-3811 bind to and inhibit leucine aminopeptidase 3. An implementation is available as an R package from Bioconductor ( https://bioconductor.org/packages/TPP2D ). We hope that our method will facilitate prioritizing targets from thermal profiling experiments.

Aggregation and disaggregation features of the human proteome.

Protein aggregates have negative implications in disease. While reductionist experiments have increased our understanding of aggregation processes, the systemic view in biological context is still limited. To extend this understanding, we used mass spectrometry-based proteomics to characterize aggregation and disaggregation in human cells after non-lethal heat shock. Aggregation-prone proteins were enriched in nuclear proteins, high proportion of intrinsically disordered regions, high molecular mass, high isoelectric point, and hydrophilic amino acids. During recovery, most aggregating proteins disaggregated with a rate proportional to the aggregation propensity: larger loss in solubility was counteracted by faster disaggregation. High amount of intrinsically disordered regions were associated with faster disaggregation. However, other characteristics enriched in aggregating proteins did not correlate with the disaggregation rates. In addition, we analyzed changes in protein thermal stability after heat shock. Soluble remnants of aggregated proteins were more thermally stable compared with control condition. Therefore, our results provide a rich resource of heat stress-related protein solubility data and can foster further studies related to protein aggregation diseases.

Dissecting the sequence determinants for dephosphorylation by the catalytic subunits of phosphatases PP1 and PP2A.

The phosphatases PP1 and PP2A are responsible for the majority of dephosphorylation reactions on phosphoserine (pSer) and phosphothreonine (pThr), and are involved in virtually all cellular processes and numerous diseases. The catalytic subunits exist in cells in form of holoenzymes, which impart substrate specificity. The contribution of the catalytic subunits to the recognition of substrates is unclear. By developing a phosphopeptide library approach and a phosphoproteomic assay, we demonstrate that the specificity of PP1 and PP2A holoenzymes towards pThr and of PP1 for basic motifs adjacent to the phosphorylation site are due to intrinsic properties of the catalytic subunits. Thus, we dissect this amino acid specificity of the catalytic subunits from the contribution of regulatory proteins. Furthermore, our approach enables discovering a role for PP1 as regulator of the GRB-associated-binding protein 2 (GAB2)/14-3-3 complex. Beyond this, we expect that this approach is broadly applicable to detect enzyme-substrate recognition preferences.

Meltome atlas-thermal proteome stability across the tree of life.

We have used a mass spectrometry-based proteomic approach to compile an atlas of the thermal stability of 48,000 proteins across 13 species ranging from archaea to humans and covering melting temperatures of 30-90 °C. Protein sequence, composition and size affect thermal stability in prokaryotes and eukaryotic proteins show a nonlinear relationship between the degree of disordered protein structure and thermal stability. The data indicate that evolutionary conservation of protein complexes is reflected by similar thermal stability of their proteins, and we show examples in which genomic alterations can affect thermal stability. Proteins of the respiratory chain were found to be very stable in many organisms, and human mitochondria showed close to normal respiration at 46 °C. We also noted cell-type-specific effects that can affect protein stability or the efficacy of drugs. This meltome atlas broadly defines the proteome amenable to thermal profiling in biology and drug discovery and can be explored online at http://meltomeatlas.proteomics.wzw.tum.de:5003/ and http://www.proteomicsdb.org.

Thermal proteome profiling for interrogating protein interactions.

Thermal proteome profiling (TPP) is based on the principle that, when subjected to heat, proteins denature and become insoluble. Proteins can change their thermal stability upon interactions with small molecules (such as drugs or metabolites), nucleic acids or other proteins, or upon post-translational modifications. TPP uses multiplexed quantitative mass spectrometry-based proteomics to monitor the melting profile of thousands of expressed proteins. Importantly, this approach can be performed in vitro, in situ, or in vivo. It has been successfully applied to identify targets and off-targets of drugs, or to study protein-metabolite and protein-protein interactions. Therefore, TPP provides a unique insight into protein state and interactions in their native context and at a proteome-wide level, allowing to study basic biological processes and their underlying mechanisms.

Identifying drug targets in tissues and whole blood with thermal-shift profiling.

Monitoring drug-target interactions with methods such as the cellular thermal-shift assay (CETSA) is well established for simple cell systems but remains challenging in vivo. Here we introduce tissue thermal proteome profiling (tissue-TPP), which measures binding of small-molecule drugs to proteins in tissue samples from drug-treated animals by detecting changes in protein thermal stability using quantitative mass spectrometry. We report organ-specific, proteome-wide thermal stability maps and derive target profiles of the non-covalent histone deacetylase inhibitor panobinostat in rat liver, lung, kidney and spleen and of the B-Raf inhibitor vemurafenib in mouse testis. In addition, we devised blood-CETSA and blood-TPP and applied it to measure target and off-target engagement of panobinostat and the BET family inhibitor JQ1 directly in whole blood. Blood-TPP analysis of panobinostat confirmed its binding to known targets and also revealed thermal stabilization of the zinc-finger transcription factor ZNF512. These methods will help to elucidate the mechanisms of drug action in vivo.

Nonparametric Analysis of Thermal Proteome Profiles Reveals Novel Drug-binding Proteins.

Detecting the targets of drugs and other molecules in intact cellular contexts is a major objective in drug discovery and in biology more broadly. Thermal proteome profiling (TPP) pursues this aim at proteome-wide scale by inferring target engagement from its effects on temperature-dependent protein denaturation. However, a key challenge of TPP is the statistical analysis of the measured melting curves with controlled false discovery rates at high proteome coverage and detection power. We present nonparametric analysis of response curves (NPARC), a statistical method for TPP based on functional data analysis and nonlinear regression. We evaluate NPARC on five independent TPP data sets and observe that it is able to detect subtle changes in any region of the melting curves, reliably detects the known targets, and outperforms a melting point-centric, single-parameter fitting approach in terms of specificity and sensitivity. NPARC can be combined with established analysis of variance (ANOVA) statistics and enables flexible, factorial experimental designs and replication levels. An open source software implementation of NPARC is provided.

Proteome-wide solubility and thermal stability profiling reveals distinct regulatory roles for ATP.

Adenosine triphosphate (ATP) plays fundamental roles in cellular biochemistry and was recently discovered to function as a biological hydrotrope. Here, we use mass spectrometry to interrogate ATP-mediated regulation of protein thermal stability and protein solubility on a proteome-wide scale. Thermal proteome profiling reveals high affinity interactions of ATP as a substrate and as an allosteric modulator that has widespread influence on protein complexes and their stability. Further, we develop a strategy for proteome-wide solubility profiling, and discover ATP-dependent solubilization of at least 25% of the insoluble proteome. ATP increases the solubility of positively charged, intrinsically disordered proteins, and their susceptibility for solubilization varies depending on their localization to different membrane-less organelles. Moreover, a few proteins, exhibit an ATP-dependent decrease in solubility, likely reflecting polymer formation. Our data provides a proteome-wide, quantitative insight into how ATP influences protein structure and solubility across the spectrum of physiologically relevant concentrations.

Thermal proteome profiling in bacteria: probing protein state in vivo.

Increasing antibiotic resistance urges for new technologies for studying microbes and antimicrobial mechanism of action. We adapted thermal proteome profiling (TPP) to probe the thermostability of Escherichia coli proteins in vivoE. coli had a more thermostable proteome than human cells, with protein thermostability depending on subcellular location-forming a high-to-low gradient from the cell surface to the cytoplasm. While subunits of protein complexes residing in one compartment melted similarly, protein complexes spanning compartments often had their subunits melting in a location-wise manner. Monitoring the E. coli meltome and proteome at different growth phases captured changes in metabolism. Cells lacking TolC, a component of multiple efflux pumps, exhibited major physiological changes, including differential thermostability and levels of its interaction partners, signaling cascades, and periplasmic quality control. Finally, we combined in vitro and in vivo TPP to identify targets of known antimicrobial drugs and to map their downstream effects. In conclusion, we demonstrate that TPP can be used in bacteria to probe protein complex architecture, metabolic pathways, and intracellular drug target engagement.