Pervasive Protein Thermal Stability Variation during the Cell Cycle.

Quantitative mass spectrometry has established proteome-wide regulation of protein abundance and post-translational modifications in various biological processes. Here, we used quantitative mass spectrometry to systematically analyze the thermal stability and solubility of proteins on a proteome-wide scale during the eukaryotic cell cycle. We demonstrate pervasive variation of these biophysical parameters with most changes occurring in mitosis and G1. Various cellular pathways and components vary in thermal stability, such as cell-cycle factors, polymerases, and chromatin remodelers. We demonstrate that protein thermal stability serves as a proxy for enzyme activity, DNA binding, and complex formation in situ. Strikingly, a large cohort of intrinsically disordered and mitotically phosphorylated proteins is stabilized and solubilized in mitosis, suggesting a fundamental remodeling of the biophysical environment of the mitotic cell. Our data represent a rich resource for cell, structural, and systems biologists interested in proteome regulation during biological transitions.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

A comprehensive comparison of tools for differential ChIP-seq analysis.

ChIP-seq has become a widely adopted genomic assay in recent years to determine binding sites for transcription factors or enrichments for specific histone modifications. Beside detection of enriched or bound regions, an important question is to determine differences between conditions. While this is a common analysis for gene expression, for which a large number of computational approaches have been validated, the same question for ChIP-seq is particularly challenging owing to the complexity of ChIP-seq data in terms of noisiness and variability. Many different tools have been developed and published in recent years. However, a comprehensive comparison and review of these tools is still missing. Here, we have reviewed 14 tools, which have been developed to determine differential enrichment between two conditions. They differ in their algorithmic setups, and also in the range of applicability. Hence, we have benchmarked these tools on real data sets for transcription factors and histone modifications, as well as on simulated data sets to quantitatively evaluate their performance. Overall, there is a great variety in the type of signal detected by these tools with a surprisingly low level of agreement. Depending on the type of analysis performed, the choice of method will crucially impact the outcome.

Neural Circuits Trained with Standard Reinforcement Learning Can Accumulate Probabilistic Information during Decision Making.

Much experimental evidence suggests that during decision making, neural circuits accumulate evidence supporting alternative options. A computational model well describing this accumulation for choices between two options assumes that the brain integrates the log ratios of the likelihoods of the sensory inputs given the two options. Several models have been proposed for how neural circuits can learn these log-likelihood ratios from experience, but all of these models introduced novel and specially dedicated synaptic plasticity rules. Here we show that for a certain wide class of tasks, the log-likelihood ratios are approximately linearly proportional to the expected rewards for selecting actions. Therefore, a simple model based on standard reinforcement learning rules is able to estimate the log-likelihood ratios from experience and on each trial accumulate the log-likelihood ratios associated with presented stimuli while selecting an action. The simulations of the model replicate experimental data on both behavior and neural activity in tasks requiring accumulation of probabilistic cues. Our results suggest that there is no need for the brain to support dedicated plasticity rules, as the standard mechanisms proposed to describe reinforcement learning can enable the neural circuits to perform efficient probabilistic inference.

MapMyFlu: visualizing spatio-temporal relationships between related influenza sequences.

Understanding the molecular dynamics of viral spreading is crucial for anticipating the epidemiological implications of disease outbreaks. In the case of influenza, reassortments or point mutations affect the adaption to new hosts or resistance to anti-viral drugs and can determine whether a new strain will result in a pandemic infection or a less severe progression. To this end, tools integrating molecular information with epidemiological parameters are important to understand how molecular characteristics reflect in the infection dynamics. We present a new web tool, MapMyFlu, which allows to spatially and temporally display influenza viruses related to a query sequence on a Google Map based on BLAST results against the NCBI Influenza Database. Temporal and geographical trends appear clearly and may help in reconstructing the evolutionary history of a particular sequence. The tool is accessible through a web server, hence without the need for local installation. The website has an intuitive design and provides an easy-to-use service, and is available at http://mapmyflu.ipmb.uni-heidelberg.de.

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.

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.

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.