Abstract
In response to an ever-increasing demand of new small molecules therapeutics, numerous chemical and genetic tools have been developed to interrogate compound mechanism of action. Owing to its ability to characterize compound-dependent changes in thermal stability, the proteome-wide thermal shift assay has emerged as a powerful tool in this arsenal. The most recent iterations have drastically improved the overall efficiency of these assays, providing an opportunity to screen compounds at a previously unprecedented rate. Taking advantage of this advance, we quantified 1.498 million thermal stability measurements in response to multiple classes of therapeutic and tool compounds (96 compounds in living cells and 70 compounds in lysates). When interrogating the dataset as a whole, approximately 80% of compounds (with quantifiable targets) caused a significant change in the thermal stability of an annotated target. There was also a wealth of evidence portending off-target engagement despite the extensive use of the compounds in the laboratory and/or clinic. Finally, the combined application of cell- and lysate-based assays, aided in the classification of primary (direct ligand binding) and secondary (indirect) changes in thermal stability. Overall, this study highlights the value of these assays in the drug development process by affording an unbiased and reliable assessment of compound mechanism of action.
eLife assessment
This study is a useful showcase of a workflow to perform large-scale characterization of drug mechanisms of action using proteomics. The work is backed by solid evidence, however, more statistical analyses and a user-friendly interface to enhance data mining by the readers are recommended. The strengths of this study include the large number of compounds tested within a common workflow and well-described experimental protocols. This will be of broad interest to medicinal chemists, toxicologists, and biochemists.
Introduction
A critical aspect of any drug development effort is the ability to define the mechanism of action. Most small molecule therapeutics function by engaging a protein (target) through a direct physical interaction, thereby inhibiting or activating a specific target in order to modulate enzymatic activity, cell signaling pathways, or gene expression. Because of its unique ability to probe cell function on a proteome-wide scale, mass spectrometry (MS) has emerged as a powerful tool for interrogating compound mechanism of action by aiding in the validation of known targets and the identification of new targets (Moellering and Cravatt, 2012; Daub, 2015; Schürmann et al., 2016; Frantzi et al., 2019; Robers et al., 2020; Kuljanin et al., 2021; Meissner et al., 2022). In addition to traditional proteomic approaches that profile protein expression and post-translational modifications (PTMs), the mass spectrometer has also become indispensable for unbiased assessments of compound target engagement and, by extension, mechanism of action. This has led to the development of numerous MS-based methods to screen the proteome for evidence of ligand binding (Molina et al., 2013; Jafari et al., 2014; Savitski et al., 2014; Liu et al., 2017; Gaetani et al., 2019; Ball et al., 2020; Van Vranken et al., 2021; Beusch et al., 2022).
All proteins have a characteristic melting temperature, and this physical property can be influenced by numerous factors (ligand-binding, protein-protein interactions, and PTMs), which can make a given protein more or less resistant to thermal denaturation (Koshland, 1958). Classically, the thermal shift assay was used to study a single protein of interest, however, by coupling the traditional thermal shift assay with modern quantitative proteomics and sample multiplexing, it is possible to assess changes in protein thermal stability on a proteome-wide scale (Molina et al., 2013; Jafari et al., 2014; Savitski et al., 2014). As such, MS-coupled thermal shift assays have emerged as a powerful tool for studying compound target engagement and mechanism of action.
The cellular thermal shift assay (CETSA) and thermal proteome profiling (TPP) use tandem mass tag (TMT)-based quantitative proteomics to quantify complete melting curves and assign melting temperatures to thousands of proteins, simultaneously, making it possible to screen the proteome for ligand-induced changes in protein thermal stability (Jafari et al., 2014; Molina et al., 2013; Savitski et al., 2014). Recently, the proteome integral solubility alteration (PISA) assay was proposed to improve the efficiency of the proteome-wide thermal shift assay (Gaetani et al., 2019). Rather than building complete melting curves, PISA uses TMT-based quantitative proteomics to estimate, or integrate, the area under a protein melting curve (Figure 1 – figure supplement 1A). Instead of TMT labelling the soluble fraction from each individual temperature—as in CETSA and TPP—PISA pools the soluble fractions of multiple samples heated across a thermal gradient such that a single TMT reporter represents an entire integrated melting curve. In the end, any change in protein thermal stability (or a thermal shift) can be quantified as a fold change in the abundance of soluble protein in a compound-treated sample vs. a vehicle-treated control after thermal denaturation and high-speed centrifugation. Therefore, an increase in melting temperature, which one could determine using CETSA or TPP, will lead to an increase in the area under the curve and an increase in the soluble protein abundance relative to controls (positive log2 fold change). Conversely, a decrease in melting temperature will result in a decrease in the area under the curve and a decrease in the soluble protein abundance relative to controls (negative log2 fold change). In the end, the ability to compress entire melting curves into a single TMT channel enables the synchronized interrogation of multiple compounds at multiple concentrations with multiple replicates in one experiment.
Owing to improved scalability, PISA results in an 8-fold theoretical improvement in throughput compared to TPP and CETSA. Therefore, PISA makes it feasible to readily approach large-scale chemical perturbation experiments, such as chemical library screening. With this in mind, we sought to establish an efficient platform for screening chemical libraries for drug-induced changes in protein thermal stability. We selected 96 compounds with known mechanisms of action for PISA screening in live K562 cells and further assayed 70 of the compounds in native lysates. Using these data, we identified evidence of on-target binding, as well as putative examples of off-target engagement. In addition to direct drug binding, we were able to identify secondary changes in protein thermal stability that did not result from direct engagement and, instead, likely stemmed from other factors such as PTMs or macromolecular interactions. Overall, this study serves as a roadmap for the scalable implementation of protein thermal stability quantitation to interrogate compound mechanism of action.
Results
Establishing a workflow for large-scale chemical perturbation screening
To enable large-scale chemical perturbation screening, we first sought to establish a robust workflow for assessing protein thermal stability changes in live K562 cells. In a PISA experiment differential thermal stability (i.e., a change in melting temperature or a thermal shift) is defined as a significant deviation in protein abundance following thermal melting and high-speed centrifugation. Most commonly this is manifested as a log2 fold change comparing a compound treated sample to a vehicle-treated control (Figure 1 – figure supplement 1A). The ranges of temperatures used in a PISA experiment can impact the magnitude of the fold-change measurements, with narrower temperature ranges often resulting in larger fold changes (Li et al., 2020). After testing a number of ranges experimentally, we ultimately settled on a temperature range of 48°C-58°C, which encompasses the back half of most protein melting curves in K562 cells (Figure 1 – figure supplement 1B; Jarzab et al., 2020).
To highlight the advantage of this window, K562 cells were treated with each of the clinically-available CDK4/6 inhibitors—ribociclib, abemaciclib, and palbociclib—or DMSO (10 µM, 30 minutes; Figure 1 – source data 1). Notably, all three compounds and a DMSO control were able to be assayed with four biological replicates in a single TMTPro 16plex experiment, requiring just 24 hours of instrument time. The cells were either heated across a range of 37°C-62°C (to approximate a full melting curve) or 48°C-58°C (thermal window) and any changes in protein thermal stability were assessed by PISA. In cells heated across the full melting curve, CDK4 and CDK6 experienced minimal log2 fold changes, whereas utilizing the narrower range yielded a statistically-significant change in CDK4 and CDK6 abundance following thermal denaturation and centrifugation, which, in a PISA experiment, is consistent with a ligand-induced thermal shift or, more specifically, a thermal stabilization (Figure 1 – figure supplement 1C). In addition to CDK4/6, several other kinases experienced a significant change in thermal stability in response to treatment with each compound (Figure 1 – figure supplement 1D-G). In fact, each of the three CDK4/6 inhibitors appear to engage a unique set of kinases, which might contribute to their disparate clinical and molecular phenotypes (Figure 1 – figure supplement 1G; Hafner et al., 2019). Finally, there was a thermal destabilization (negative log2 fold change) of RB1, a well-established phosphorylation target of CDK4/6 (Figure 1 – figure supplement 1C-F; Narasimha et al., 2014). The simplest explanation of this shift is that inhibition of CDK4/6 prevents phosphorylation of RB1 thereby inducing a change in its thermal stability.
A large-scale chemical perturbation screen in K562 cells
We next sought to execute a large-scale chemical screen in live K562 cells. To that end, we curated a custom chemical library comprised of 96 commonly used cancer drugs and tool compounds. Each compound selected for the library had at least one well-annotated target and a known mechanism of action (Table 1). Of these 96 compounds, 70 targeted protein kinases, while the remaining 26 compounds targeted other classes of proteins including histone deacetylases (HDACs), lysine demethylases, poly (ADP-ribose) polymerases (PARPs), and others (Figure 1A). The library also contained multiple compounds targeting certain proteins— eight HDAC inhibitors, four aurora kinase inhibitors, and six BRAF inhibitors, etc. To initiate the screen, K562 cells were treated with each compound in biological duplicate at 10 µM for 30 minutes (Figure 1A). The live cells were heated across a thermal gradient (48°C-58°C) and any changes in protein thermal stability were assessed using PISA (Figure 1 – figure supplement 2). This resulted in a total of 256 samples, which were arranged into 16 TMTpro 16-plex experiments, which required a minimum of 384 hours of instrument time (Figure 1 – figure supplement 3A).
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig1.tif/full/max/0/default.jpg)
Establishing a robust and efficient workflow for assessing changes in protein thermal stability in living cells.
A. Schematic of the cell-based PISA workflow utilized in the large-scale chemical perturbation screen. B. Histogram plotting the log2 fold change values for all thermal stability measurements (main panel). Box plot depicting the %CV of curated targets vs. all other proteins (inset). C. Plots highlighting the significant changes in protein thermal stability following treatment with OTS167 (top) or GSK-LSD1 (bottom). The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. Blue boxes contain proteins that have been stabilized (increase in melting temperature) and orange boxes contain proteins that have been destabilized (decrease in melting temperature). Orange points represent the known target of each compound. D. The total of significant thermal stability changes quantified upon treatment with each compound. The number of stabilizing events is plotted on the y-axis and destabilizing events are plotted on the x-axis. E. Plot comparing the log2 fold changes measurements for each protein in BI 2536 (y-axis)- and volaserib (x-axis)-treated K562 cells.
In total, we quantified approximately 6,800 proteins per treatment for a total of 871,120 drug-protein thermal stability measurements (Figure 1 – figure supplement 3B,C and Figure 1 – source data 2). Compared to the full dataset, annotated targets of the library compounds had significantly different thermal stability compared to non-targets (Wilcoxon rank sum test, p=3.56×10-17). To define high-confidence thermal stability changes, we developed an empirically derived framework (see Methods; Figure 1 – figure supplement 3D). Briefly, we determined the log2 fold change of each compound treatment in reference to the vehicle-treated control and then quantified the trimmed standard deviation of all stability measurements for a given protein across all treatments. In order to be considered a hit, the log2 fold change needed to be 3.5 standard deviations from the mean and have an absolute value greater than or equal to 0.2 (|log2Cmpd/DMSO|>0.2, |nSD|>3.5; Figure 1 – figure supplement 3D-G). The fold change filtering removed background proteins while the per protein standard deviation filtering removed proteins that were observed to engage in large numbers of non-specific interactions, such as the kinase GAK (Figure 1 – source data 2). From the 871,120 total protein thermal stability measurements, filtering the data resulted in 3,156 putative hits (|log2Cmpd/DMSO|>0.2, |nSD|>3.5; Figure 1B).
Upon assembling compound-centric thermal stability profiles, known compound targets could be easily resolved. For example, in cells treated with OTS167, the primary target, MELK, had one of the largest negative log2 fold changes (log2Cmpd/DMSO: −0.840; nSD: −5.688), which is consistent with a thermal destabilization (Figure 1C). In addition to its intended target, OTS167 also induced consistent protein thermal stability changes in other kinases (WEE1, PRKD3, CDK11B, AURKA, and others). These findings are consistent with previous reports demonstrating that the mechanism of action for OTS167 is in part driven by off-target protein engagement (Giuliano et al., 2018). Similarly, treatment with GSK-LSD1 caused a marked thermal stabilization of its primary target—the lysine demethylase, KDM1A (log2Cmpd/DMSO: 1.726; nSD: 26.825; Figure 1C). In addition to the primary target, GSK-LSD1 also stabilized several proteins known to form a complex with KDM1A, including GSE1 and RCOR1/3 (Figure 1C). Finally, we observed disparities in how different molecules targeting the same protein impact the target. The library contained two compounds known to target AKT1—MK-2206, an allosteric inhibitor and CCT128930, an ATP- competitive inhibitor (Amiran et al., 2023; Hirai et al., 2010; Yap et al., 2011). Interestingly, engagement of the allosteric ligand, MK-2206, induced a large change in AKT1 thermal stability, while the ATP-competitive inhibitor appeared to have no effect (Figure 1 – figure supplement 4A,B).
Beyond the known targets, these data define a set of high-confidence ligand-induced thermal stability changes for each of the 96 compounds assayed in this screen (Figure 1D). OTS167 caused a significant change in the thermal stability of 154 proteins, including other protein kinases, which likely stem from off-target engagement (Figure 1C,D). GSK-LSD1 treatment, on the other hand, resulted in only ten significant changes (Figure 1C,D). While OTS167 was among the most active compounds in the screen (as defined by total filtered changes), treatment with the PLK1 inhibitor volasertib resulted in the most thermal stability changes (435, Figure 1D). While there were several other compounds that caused >100 changes in thermal stability (ER-27319, BS-181, SGI-1776, etc.), the median number of thermally stabilized or destabilization proteins per compound was 14.5 (Figure 1D).
Our screening library contained two dual PLK1/BRD4 inhibitors—volasertib and BI 2536—which have similar chemical structures (Figure 1E). Appropriately, the thermal stability profiles that resulted from treatment with each compound correlated well (rspearman = 0.49, p < 1.0×10-15), highlighted by consistent changes in the thermal stability of PLK1, BRD4, and the proteome in general (Figure 1E). In addition to these dual PLK1/BRD4 inhibitors, the library also contained two compounds—JQ-1 and CPI-203—which target BRD4, but not PLK1. While the thermal stability profiles of the dual PLK1/BRD4 inhibitors were correlated, this was not the case for volasertib and JQ-1. Indeed, there was clear evidence that JQ-1 could engage BRD4 but not PLK1 however, there was little overall correlation, suggesting that the observed volasertib-dependent changes are largely PLK1-depepndent (Figure 1 – figure supplement 4C).
Finally, we wanted to determine if the set of high-confidence thermal stability alterations could provide insights related to compound mechanism of action. PF-3758309 is designated as a PAK4 inhibitor. Including PAK4, treatment with this inhibitor caused a total of 102 proteins to undergo a significant change in protein thermal stability (Figure 1D). Expectedly, this list included a number of protein kinases including MELK, TBK1, WEE1, and others (Figure 1 – figure supplement 4D). Unexpectedly, there was also a significant enrichment of spliceosome subunits among these proteins (Figure 1 – figure supplement 4D). These changes could stem from upstream inhibition of a kinase that regulates splicing and that the mechanism of action of this compound might, at least, in part, works through an alteration of this process. This result is consistent with a recent study that identified PF-3758309 as a potent modulator of pre-mRNA splicing (Shi et al., 2020). Overall, these data provide further evidence that PF-3758309 is a promiscuous kinase inhibitor capable of impacting the assembly and/or activity of the spliceosome.
A protein-centric view of compound-dependent changes in thermal stability
Having established criteria for identifying significant compound-dependent changes in protein thermal stability, we focused on the known targets of the library compounds in order to identify key factors affecting the interpretation of PISA data. In total, we quantified at least one known target for 79 of the 96 compounds assayed in the primary screen (84%) and found that 56 (71%) of these compounds induced a thermal shift of an annotated target (Figure 2A). While compound treatment more commonly resulted in a thermal stabilization of known targets (positive log2 fold change), we also observed a number of examples in which inhibitor binding stimulated a thermal destabilization (negative log2 fold change). Consistent thermal stabilization was observed for proteins including BRD4, HDAC1, and AURKA. Conversely, compound engagement by PLK1 resulted in a consistent thermal destabilization (Figure 2A). Therefore, both positive and negative log2 fold changes represent evidence of target engagement and specific proteins are consistently stabilized or destabilized in response to inhibitor binding.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig2.tif/full/max/0/default.jpg)
Using cell-based PISA data to assess compound target engagement.
A. A plot of the mean log2 fold change (x-axis) and mean nSD (y-axis) of each known compound-target pair that were quantified in the screen. The blue box contains proteins that have been stabilized (increase in melting temperature) and the orange box contains proteins that have been destabilized (decrease in melting temperature). B and C. Plots depicting the mean log2 fold change of p38α (B) and PIK3CB (C) in response to each of the compound and DMSO treatments. Orange points indicate compounds known to target each protein. Blue points indicate other compounds that induce a significant change in thermal stability for each protein. Green dashed lines mark a SD of 3.5 for each protein. D. A protein-centric view of PLK1 thermal stability in response to all treatments. Log2 fold change is plotted on the y-axis. The points represent each of the 256 treatments that were performed. Orange points indicate compounds known to target PLK1. Blue points indicate other compounds that result in a significant thermal destabilization of PLK1. Green dashed lines mark a SD of 3.5 cutoff for PLK1. E and F. K562 cells were treated with the indicated concentrations of palbociclib or BI 2536 for 15 minutes. Changes in protein thermal stability are represented as a log2 fold change in protein abundance for each treatment in reference to a DMSO-treated control (E). PLK1 activity was determined using a western blot for p-TCTP (S46) and total TCTP (F). G. Chemical structures of BI 2536, palbociclib and NVP- TAE-226. Common structural features with BI 2536 are highlighted in blue for Palbociclib and NVP-TAE-226. H. Co-crystal structure of CDK6 bound to palbociclib (yellow) with BI 2536 (blue) modeled into the active site. I. Co-crystal structure of PLK1 bound to BI 2536 (blue) with palbociclib (yellow) modeled into the active site.
Next, we characterized the magnitude of thermal stability changes that we measured for known targets. Having treated cells with each compound at 10 µM, we assumed that each target was fully saturated with a given compound. Therefore, the measurements likely reflect the maximal ligand-induced thermal stability changes for each target. We observed a large range of thermal stability measurements for known compound-target pairs, from a four-fold reduction in protein stability to a four-fold increase in protein stability upon compound engagement (Figure 2A). We also observed many proteins with small (15%) but consistent stabilization following compound treatment. p38α-targeting compounds ralimetinib and doramapimod caused a large apparent increase in p38α thermal stability (log2Cmpd/DMSO = 1.52 and 1.94, respectively), while compounds targeting PI3K (pictilisib, duvesilib, and IC-87114), on the other hand, generated smaller thermal stability changes (log2Cmpd/DMSO = 0.258, 0.319, and 0.213, respectively; Figure 2B,C and Figure 2 – figure supplement 1A,B). In all cases, each of these measurements surpassed the nSD cutoff of 3.5 for each particular target. Taken at face value, one might conclude that the greater thermal stability effect that doramapimod exerts on p38α is more meaningful than the much smaller change that pictilisib exerts on PIK3CB, however, this is not necessarily true as both compounds can engage their targets at nanomolar concentrations. Instead, we posit that maximum ligand-induced change in thermal stability is target-specific. Therefore, a small thermal stability change for one protein (PIK3CB, log2Cmpd/DMSO ≥ 0.2), might be just as meaningful as a large change for another (p38α, log2Cmpd/DMSO ≥ 2).
Using the cell-based screening data to identify off-target compound engagement
Having defined a set of high-confidence thermal stability changes for each compound, we next looked for potential examples of off-target engagement. In addition to volasertib and BI 2536 (PLK1/BRD inhibitors), several other compounds from the screen also thermally destabilized PLK1 to a similar extent (Figure 2D). These included two promiscuous kinase inhibitors—OTS167 and PF-3758309. In addition to these compounds, two other ATP-competitive kinase inhibitors also impacted the thermal stability of PLK1—NVP-TAE-226, a FAK inhibitor and palbociclib, one of the CDK4/6 inhibitors used to establish our workflow. Because these two compounds induced a change in thermal stability to a similar magnitude as BI 2536 and volasertib, we hypothesized that these molecules might bind the active site and inhibit PLK1.
To further explore the impact of NVP-TAE-226 and palbociclib on PLK1 thermal stability, we treated K562 cells with increasing concentrations of BI 2536, NVP-TAE-226, and palbociclib, ranging from 100 nM to 10 µM and performed cell-based PISA profiling to quantify the changes in thermal stability that occurred at each inhibitor concentration (Figure 2 – source data 1 and source data 2). Ideally, this would provide a relative assessment of the ability of each molecule to engage (and inhibit) PLK1. Appropriately, BI 2536 caused a dose-dependent decrease in PLK1 thermal stability down to the lowest dose tested (100 nM; Figure 2E and Figure 2 – figure supplement 1C). Treatment with palbociclib and NVP-TAE-226 also caused a dose-dependent decrease in PLK1 thermal stability, however, each of these compounds appeared to impact PLK1 thermal stability to a lesser extent than BI 2536. In fact, 100 nM BI 2536 treatment impacted PLK1 thermal stability to the same degree as 10 µM treatments with palbociclib and NVP-TAE-226 (Figure 2E and Figure 2 – figure supplement 1C). Next, we synchronized HCT116 cells in G2/M phase to test the ability of each compound to directly inhibit PLK1 in a cell-based assay. During this phase of the cell cycle, PLK1 is active and known to phosphorylate TCTP (at serine 46), a protein involved in microtubule stabilization (Ulise Cucchi et al., 2010). Following synchronization, these cells were treated with increasing concentrations of each compound and phosphorylation of TCTP at S46 was assayed to read out PLK1 activity. BI 2536 completely abolished PLK1- dependent phosphorylation of TCTP at all tested doses. Palbociclib and NVP-TAE-226 also inhibited TCTP phosphorylation however, these compounds were much less potent than BI 2536 and were only able to inhibit PLK1 at doses between 1 µM and 10 µM (Figure 2F and Figure 2 – figure supplement 1D). Importantly, these results closely mirror the PISA data, providing further evidence that proteomic thermal stability measurements accurately recapitulate target engagement. Therefore, palbociclib and NVP-TAE-226 appear to be weak inhibitors of PLK1, which is consistent with previous studies that suggested these compounds as off-target inhibitors of this kinase (Hafner et al., 2019).
Based on the PISA data, palbociclib was able to engage both CDK4/6 and PLK1, while BI 2536 was only able to engage PLK1 (Figure 2E). This is despite the fact that both compounds have a similar chemical structure (Figure 2G). In order to further interrogate this disparity, we modeled palbociclib and BI 2536 into co-crystal structures of PLK1 (PDB: 2rku) and CDK6 (PDB: 5l2i) with their ‘specific’ inhibitors (Figure 2H,I; Chen et al., 2016; Kothe et al., 2007). When we overlaid BI 2536 in the active site of CDK6, we observed three points of steric clashing with the protein density of the CDK6 active site. These findings were consistent with the lack of thermal stabilization of CDK4/6 when cells were treated with BI 2536. From overlaid projections of palbociclib and BI 2536 in PLK1, we did not observe steric interference in binding suggesting that both of these inhibitors could readily bind to the active site of PLK1. These data are again consistent with our observed thermal destabilization of PLK1 when cells were treated with either BI 2536 or palbociclib.
PISA screening in crude cell extracts
Thus far we have focused exclusively on protein thermal stability in living cells, but PISA can also be used to screen for changes in cell lysates or crude extracts. Having already assembled one of the largest drug-based thermal denaturation studies in cells, we wanted to further profile a subset of these compounds in native extracts. We selected 70 compounds from the cell-based screen for further interrogation in native extracts. For this screen, we prepared crude extracts by dounce homogenizing K562 cells in PBS (Sridharan et al., 2019). The crude extracts were treated with each of the 70 compounds in biological duplicate at 10 µM for 30 minutes and any changes in thermal stability were determined using PISA. In total, the 70 compounds and DMSO controls generated a total of 160 samples, which were arrayed across 10 TMT 16-plex experiments and quantified using approximately 240 hours of instrument time (Figure 3A and Figure 3 – figure supplement 1A).
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig3.tif/full/max/0/default.jpg)
A chemical perturbation screen of 70 compounds in K562 native extracts.
A. Schematic of the lysate-based PISA workflow utilized in the large-scale chemical perturbation screen. B. Histogram plotting the log2 fold change values for all thermal stability measurements in lysate-based screen. C. Plots highlighting the significant changes in protein thermal stability following treatment with MK-2206 (top) or GSK-LSD1 (bottom). The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. Blue boxes contain proteins that have been stabilized (increase in melting temperature) and orange boxes contain proteins that have been destabilized (decrease in melting temperature). Orange points represent the known target of each compound. D. Plot depicting the total thermal stability changes quantified upon treatment with each compound. The number of stabilizing events is plotted on the x-axis and destabilizing events are plotted on the y-axis. E. Plot depicting the number of compounds for which a target was quantified (gray), the number of compounds that caused a significant change in a known target (orange), and the number of compounds that failed to cause a significant change in a known target (blue). F. Plot depicting the total percentage of compounds that caused a significant change in thermal stability of a known target using each experimental approach.
We quantified an average of 7,840 proteins per treatment (Figure 3B and Figure 3 – figure supplement 1B), which resulted in a total of 627,176 thermal stability measurements from all K562 lysate-based experiments (Figure 1 – source data 2). After applying the same filters used in the cell-based screen, our dataset consisted of 2,176 protein thermal stability changes (Figure 3B and Figure 3 – figure supplement 1C). Similar to the cell-based studies, the known targets of library compounds frequently experienced a significant change in thermal stability (Figure 3 – figure supplement 1C-G). For example, treatment of K562 lysates with MK-2206 and GSK-LSD1 stabilized AKT1/2 and KDM1A, respectively (Figure 3C). In addition to AKT1/2, MK-2206 also stabilized additional kinases, including subunits of RPS6K and PI3K, which highlights the ability of this approach to find putative off-target engagements (Figure 3C). Finally, we observed a GSK- LSD1-dependent stabilization of several KDM1 binding partners, which echoes the similar observation made in cell-based experiments (Figure 3C).
Since changes in thermal stability in lysate-based experiments stem almost exclusively from direct ligand binding events, we expected to observe fewer overall stability changes in this setting compared to cell-based experiments. As expected, the lysate-based screen (2,176) yielded fewer significant changes than the cell-based screen (3,156). It is also important to note that three compounds (vorinostat, Y-33075, and GSK-1070916) contribute 1,200 (55%) of all significant changes in the lysate-based screen. It is highly unlikely that these three molecules actively engage so many proteins and, therefore, the 2,176 hits in the lysate-based screen were likely affected in part by artefactual effects of lysate-based analyses. Overall, the median number of protein thermal stability changes in lysates was 66% of those observed in cells, with 9.5 proteins in lysate-based experiments compared to 14.5 in cell-based experiments (Figure 3D).
We directly compared the overall ability of the cell- and lysate-based approaches to identify the known target(s) of each compound. In the cell-based screen, we assayed 96 molecules and quantified known targets of 79 of these inhibitors. Of the 70 molecules assayed in K562 lysates, we quantified protein targets for 64 (91%) of these molecules and found that 45 (70%) of these compounds induced a change in an annotated target (Figure 3E,F). In each case we observed a thermal stability change of a known target for about 70% of the inhibitors we screened (Figure 3F). However, when we combined the data from the two approaches, we were able to quantify targets for 59 of the 69 compounds that were assayed in both settings and were able to observe a significant change in thermal stability for 82% of these compounds (Figure 3E,F). Therefore, combining the two approaches allows for greater coverage of the cellular proteome and provides a better chance of observing the protein target for a compound of interest.
Cell- and lysate-based PISA are complimentary approaches to determine compound engagement
Having re-screened dozens of compounds in lysate-based experiments, we wanted to assess the complementarity of these approaches. We focused on the known compound targets quantified in both cell- and lysate-based experiments and compared the thermal stability changes measured using each approach (Figure 4A). Encouragingly, we observed compounds that impacted the thermal stability of their known targets in both cell- and lysate-based approaches. These include GSK-1070916 and MK-2206, which stabilized AURKA and AKT1, respectively, in both cell- and lysate-based experiments (Figure 4A,B). Conversely, we observed compounds that failed to impact the thermal stability of their target proteins in either the cell-based or lysate-based experiments (Figure 4A and Figure 4 – figure supplement 1A). These include compounds such as IRAK4-Inhibitor-1 and entospletinib which are known to target IRAK4 and SYK kinase, respectively.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig4.tif/full/max/0/default.jpg)
Cell- and lysate-based PISA are complimentary approaches for assessing mechanism of action.
A. A plot depicting the mean log2 fold change for all compound-target pairs quantified in cell-(x-axis) and lysate-based (y-axis) PISA. Dashed lines indicate a log2 fold change of +/-0.2. B. The log2 fold change of a selected number of compound-target pair in cell- and lysate-based PISA. C-E. Plots depicting the log2 fold change values for BRD4 (C), BRAF (D), and RIPK1 (E) in response to treatment with all compounds assayed using cell- and lysate-based PISA. Dashed lines indicate a log2 fold change of +/-0.2. F. HCT116 cell lysates (N=3) were treated with each compound or DMSO for 15 minutes and assay using PISA. The plot indicates the log2 fold change of several proteins in response to treatment with each compound. G. HCT116 cells were treated with the indicated concentration of each compound for 1 hour following the initiation of necroptosis using TSZ treatment. RIPK1 autophosphorylation was assayed by western blot using antibodies targeting p-RIPK1 (S116) or total RIPK1.
We note that the magnitudes of the thermal stability changes varied between the cell- and lysate-based approaches. Yet, compounds that stabilized their annotated protein targets in cells also tended to stabilize their targets in lysates with one notable exception. Sotrastaurin induces a thermal stabilization of its target, PRKCA, in lysates (log2Cmpd/DMSO = 0.924), but a thermal destabilization in cells (log2Cmpd/DMSO = −0.502; Figure 4A,B). These data suggest that the thermal destabilization observed in cell-based experiments might stem from a complex biophysical rearrangement.
While many compounds significantly altered their known targets using both approaches, there were some compounds that preferentially altered the thermal stability of known targets in one experimental setting but not the other (Figure 4A). For example, there were four compounds known to target BRD4 (JQ-1, CPI-203, BI 2536, and volasertib) that were assayed in both cell- and lysate-based experiments. Treatment of cells with any of these four inhibitors resulted in a log2Cmpd/DMSO ∼ 2, while the lysate-based assays resulted in smaller thermal stability changes (Figure 4C and Figure 4 – figure supplement 1A). Conversely, multiple inhibitors caused a significant shift in BRAF thermal stability in lysates (AZ 628, PLX-4720, L-779450) but did not impact BRAF thermal stability in cell-based assays (Figure 4D and Figure 4 – figure supplement 1A). Thus, certain protein targets are more prone to thermal stability changes in one experimental setting compared to the other.
Combining cell- and lysate-based data to discover off-target engagement
Having demonstrated the complementarity of the cell- and lysate-based approaches, we sought to explore the corroborative value of integrating these data. We reasoned that any unexpected ligand-induced changes in thermal stability that were shared between the two approaches would provide strong evidence of compound engagement. Despite an absence in the library of compounds designated as RIPK1 kinase inhibitors, there were four compounds that significantly stabilized RIPK1 in cells and lysates (Figure 4E and Figure 4 – figure supplement 1B). Although three of the compounds—GSK-2656167, GSK-2606414, and tozasertib—had previously been shown to be off-target inhibitors of RIPK, AZD-5438, a cyclin-dependent kinase (CDK) inhibitor, had never been attributed such activity (Byth et al., 2009; Martens et al., 2018; Rojas-Rivera et al., 2017). Nonetheless, this compound was capable of stabilizing RIPK1 to a similar extent as the other three known RIPK1 off-target inhibitors (Figure 4E and Figure 4 – figure supplement 1B).
In order to gain a sense of the potency of AZD-5438 (and the other molecules), we assayed these compounds alongside necrostatin-2, a bona fide RIPK1 inhibitor, in HCT116 lysates (Figure 4 – source data 1). Treatment with necrostatin-2 or GSK-2606414 resulted in a similar thermal stability change for RIPK1, that was far greater than both tozasertib and AZD-5438 (Figure 4F and Figure 4 – figure supplement 1C-F). Treatment of HCT116 cells with TNFα, a Smad mimetic, and a pan-caspase inhibitor (zVAD) will lead to an initiation of necroptosis and, importantly, autophosphorylation of S116 by RIPK1 (Laurien et al., 2020). While each of our putative RIPK1 inhibitors was capable of inhibiting this autophosphorylation, AZD-5438 and tozasertib were less potent than necrostatin-2 and GSK-2606414 (Figure 4G). Therefore, the thermal stability measurements made in our lysate-based PISA experiment correlated with compound potency. This demonstrates that the combination of cell- and lysate-based data can be employed to identify instances of off-target engagement, even for weak inhibitors.
In addition to compound potency, the combined cell and lysate data provided valuable information regarding compound specificity. While it is difficult to separate necrostatin-2 and GSK-2606414 based on potency, necrostatin-2 had a greater apparent specificity with respect to RIPK1 engagement. Indeed, necrostatin-2 only impacted the thermal stability of RIPK1 in the lysate-based experiments (Figure 4 – figure supplement 1C). GSK-2606414, on the other hand, engaged other kinases including CSK, LIMK1, and TBK1 (Figure 4 – figure supplement 1D). Finally, AZD- 5438 and tozasertib were less potent and less specific for RIPK1 than necrostatin-2 (Figure 4 – figure supplement 1E,F).
Disparities in cell- and lysate-based data pinpoint secondary changes in protein thermal stability
While the data suggests that cell- and lysate-based PISA are equally valuable in screening the proteome for evidence of target engagement (Figure 3F), cell-based experiments have added potential to identify secondary changes in protein thermal stability that occur independent of direct ligand binding. These changes could occur as a result of changes in interactions with other proteins, ligands, nucleic acids, or the effects of PTMs. Importantly, all of these factors can help to define compound mechanisms of action. This type of change was previously highlighted in the context of the CDK4/6 inhibitors. Indeed, treatment of cells with CDK4/6 inhibitors induced a thermal destabilization of RB1, which we attributed to a change in phosphorylation due to upstream target inhibition (Figure 1 – Figure supplement 1C). Other CDK inhibitors (AZD-5438 and flavopiridol) also lead to a significant destabilization of RB1. In lysates, on the other hand, when cellular signaling networks are disrupted, RB1 thermal stability was unaffected by these compounds (Figure 5A). This concept is also apparent when comparing the PF-3758309- dependent thermal stability profiles in each experimental setting (Figure 5 – figure supplement 1A). Indeed, there is a striking difference in the number of significant changes in cells (102) vs. lysates (21). This disparity is driven by splicesome subunits, which are significantly impacted in cell-based experiments but unaffected in lysates. This supports the interpretation that the apparent thermal destabilization of the spliceosome in cells is the result of secondary changes rather than direct ligand binding.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig5.tif/full/max/0/default.jpg)
Disparities in cell- and lysate-based PISA can point toward secondary changes in protein thermal stability.
A and B. Plots depicting the mean log2 fold change values for RB1 (A) and CRKL (B) in response to treatment with all compounds assayed using cell (x)- and lysate (y)-based PISA. Dashed lines indicate a log2 fold change of +/-0.2. C. K562 cells were treated with each indicated compound at 10 µM for 15 minutes. Western blots were used to assess levels of p-CRKL (Y207) and total CRKL. D. A plot of the log2 fold change of CRKL in response to each indicated compound. E. Plots highlighting the significant changes in protein thermal stability following treatment with AZD-7762. The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. Blue boxes contain proteins that have been stabilized (increase in melting temperature) and orange boxes contain proteins that have been destabilized (decrease in melting temperature). Orange points represent the known target of AZD-7762 and the blue points represents tyrosine kinases that experience a significant change in thermal stability. F. A plot of the log2 fold change of each protein in response to AZD- 7762 treatment.
Exploiting the disparities between cell- and lysate-based data can be a powerful tool in determining compound mechanism of action in the absence of reliable evidence of target engagement. These studies were performed in K562 cells/lysates, a chronic myelogenous leukemia cell line that expresses the BCR-ABL fusion. Consistent with previous reports, treatment of K562 cells with BCR-ABL-targeted compounds did not induce any apparent change in the thermal stability their primary target (Molina et al., 2013). We were, however, able to detect a change in thermal stability for well-known BCR-ABL substrates—most notably CRKL and CRK (Figure 5B and Figure 5 – figure supplement 1B-D). Importantly, the change in CRKL thermal stability only occurred in cell-based experiments and not in lysates (Figure 5 – figure supplement 1B,C). These data are consistent with a distal effector of thermal stability. A number of other compounds that largely target tyrosine kinases also caused a significant thermal destabilization of CRKL in K562 cells but not lysates (Figure 5B). We treated K562 cells with a small panel of kinase inhibitors and found that compounds that induce a thermal destabilization of CRKL also inhibit phosphorylation at Y207 (Figure 5C). Therefore, the change in CRKL thermal stability correlates with a change in phosphorylation and is likely dependent on the inhibition of a primary compound target—such as BCR-ABL—that then inhibits CRKL phosphorylation. These data highlight that protein thermal stability measurements can reveal critical insights into primary and secondary effectors of compound mechanism of action, including kinase-substrate relationships.
Using secondary changes to identify putative off-target effects
Previously, we demonstrated that the combined cell and lysate data could identify off-target inhibition through direct target engagement. Next, we wanted to determine if secondary effects could be used to determine putative off-target engagement. In addition to tyrosine kinase inhibitors, we also observed a strong thermal destabilization of CRKL in response to treatment with the CHEK1 inhibitor AZD-7762 (Figure 5B,C and Figure 5 – figure supplement 5B). Importantly, our screen contained two additional inhibitors known to engage CHEK1 (CCT241533 and BML-227), neither of which had any impact on CRKL thermal stability (Figure 5D). Therefore, the AZD-7762-dependent destabilization of CRKL is unlikely to be related to on-target inhibition of CHEK1 and, instead, likely stems from off-target engagement of other kinases. Consistent with our hypothesis, AZD-7762 induced a change in the thermal stability of many tyrosine kinases in lysate-based experiment including FER, LYN, CSK, and YES1 (Figure 5E,F and Figure 5 – figure supplement 5E). These putative interactors were not limited to tyrosine kinases but also included dozens of serine/threonine kinases. Overall, these data illustrate how secondary changes in protein thermal stability can be exploited to identify putative examples of off-target engagement.
Library-scale assays connect consistent thermal stability changes
Binary comparisons of compounds with shared targets such as volasertib and BI 2536 (Figure 1E) revealed consistent thermal stability responses for known targets PLK1 and bromodomain proteins. Next, we set out to determine if small but consistent effects on the proteome due to drug treatment could reveal new information concerning compound engagement of proteins, protein classes, or protein complexes. To address this, we implemented a network-based approach using the correlation of thermal stability changes for each protein. To this end, we mapped the all-by-all correlation of proteins in the cell-based dataset. We filtered the 70,392,100 binary comparisons to include only the correlations in the top ∼5% of absolute magnitude (rspearman > 0.35). Because 70 of the 96 compounds used were kinase inhibitors, we further focused on correlated thermal stability changes within the 396 human kinases we quantified in our dataset (Figure 6A).
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig6.tif/full/max/0/default.jpg)
Concordance of protein engagement for protein complex members and structurally similar kinases.
A. Correlation network graph for all observed kinases in the cell-based PISA assays. Human kinases (nodes) are connected based on correlation coefficients (RSpearman > 0.35). Nodes are colored based on kinase families. B. Example of significant correlation between RPS6KA1, RPS6KA2, and RPS6KA3 with consistent but sometimes small changes in thermal stability. C. Sub-graph of the network in A showing CMGC-family p38α/MAPK14 (isoforms 1 and 2) correlation with CAMK-family member kinases MAPKAPK2 and MAPKAPK3. Edges are colored based on RSpearman. D. Correlation plots for the subnetwork proteins in C. E. Protein interaction network derived from the BioPlex interactome highlighting that MAPKAPK proteins directly interact with p38α/MAPK14 in multiple cell lines. F. Sub-graph of the network in A highlighting the AGC protein kinases. Edges are colored based on RSpearman. G. AGC kinases AKT1 and AKT2 have a significant correlation driven in part by thermal stability effects of MK-2206. H. Structural comparison of the CMGC kinases GSK3A (Alphafold) and CDK16 complexed with the TKI rebastinib (PDB: 5g6v). Alignment generated an RMSD of 0.89 Å. I. Plot of the significant correlation between CDK16 and GSK3A derived from the network in A.
We first sought to determine if the network-based approach could reveal consistent thermal stability responses for kinases with conserved sequences or general kinase families (Eid et al., 2017). Strikingly, for the highly conserved 90 kDa ribosomal S6 kinases, we observed highly correlated thermal stability for three members, RPS6KA1, RPS6KA2, and RPS6KA3 (Figure 6B). We reasoned that this may be a function of the sequence conservation between these proteins or direct binding interactions (Huttlin et al., 2021, 2017). To this end we compared the cell and lysate data. In cells, RPS6KA1, RPS6KA2, and RPS6KA3 were generally destabilized by compound treatment. However, in lysates, these proteins were generally stabilized, even by the same compound that destabilized the proteins in cells, such as NVP-TAE-226. These data were in line with what we observed for sotrastaurin engagement of PRKCA (Figure 4B). Interestingly, within our thermal stability network, we observed PRKCA thermal stability was significantly correlated with PRKCB and PRKCQ (Figure 6 – figure supplement 1A,B).
Both RPS6K and PRKC subnetworks reveal correlation between kinases of the same kinome group. RPS6K proteins belong to the CAMK kinases, and PRKC proteins belong to the AGC kinases. To explore this within group correlation, we annotated the network of quantified human protein kinases based on their established groupings (AGC, Atypical, CAMK, CK1, CMGC, STE, TK, TKL, Other) along with a set of non-protein kinases (Figure 6A). Within the kinase network, we observed correlated thermal stability changes within and between families of kinases. For example, we observed correlated profiles between isoforms of p38α/MAPK14 (a CMGC kinase) and the MAPKAPK2/3 (CAMK kinases, Figure 6C,D). MAPKAPK2 and MAPKAPK3 were both stabilized with treatment of the p38 kinase inhibitor doramapimod. Indeed, these proteins are known to both act as substrates of p38 and to form complexes in divergent human cancer cells (Figure 6E). We further observed interconnected correlation networks within larger kinase groups. The AGC kinases, including AKT1/2 and ROCK1/2 were particularly well correlated (Figure 6F). AKT1 and AKT2 was driven in part by their thermal stability change due to AKT inhibitor MK-2206 (Figure 6G).
We reasoned that because protein kinases and non-protein kinases share similar structures and cofactor binding pockets, that non-protein kinases may represent an important class of off-target effectors for our library compounds. Within the correlated thermal stability profile network, we observed strong correlation due to compound-based thermal stability changes within the non-protein kinase group. PIP4K2A, PIP4K2B, and PIP4K2C thermal stabilities were highly correlated in both cells and lysates (Figure 6 – figure supplement 1C-E). For example, the PLK1 inhibitor BI 2536 stabilized all three PIP4K proteins in both cells and lysates. It is worth noting that the magnitude of thermal stability changes for these PIP4K proteins was higher in cells than in lysates consistent with loss of membranes and cellular metabolite concentrations (Sridharan et al., 2019). Yet the thermal stability responses in both assays were highly consistent for all three proteins using either approach. We explored the observation of small but consistent changes further and observed significant correlation of the SRPK1, SRPK2, and RIOK2 kinases (Figure 6 – figure supplement 1F,G). In our final datasets, we did not call any of these proteins as hits due to either the large numbers of thermal stability effects (SRPK1, SRPK2) or the small magnitude of their thermal stability changes (RIOK2). However, with correlations of 0.95 and 0.77 for SPRK1-SPRK2 and SPRK1-RIOK2, it is clear that these proteins respond consistently to chemical perturbation suggesting that our initial filters for full PISA datasets may mask some real and important protein-compound engagements.
Finally, we observed a thermal stability correlation between the CMGC kinases CDK16, GSK3A, and GSK3B (Figure 6A). CDK inhibitors have been reported as active for GSK3 kinases owing to phylogenetic and structural conservation (Figure 6H; Fischer, 2003; Meijer et al., 2003). In our cell-based approach the inhibitor AZD5438 has dual efficacy for GSK3A and CDK16 and we observed that thermal stability responses for CDK16, GSK3A, and GSK3B were correlated (Figure 6H). This effect was driven in part by annotated inhibitors of CDKs (AZD5438, flavopiridol) or GSKs (AZD5438, LY2090314). Yet, known inhibitors of several kinases stabilized both GSKs and CDK16, including MELK (OTSSP167), JAK (LY2784544), FAK (defactinib), and PAK4 (PF- 3738309), as well as the PARP inhibitor rucaparib (Figure 6I). Interestingly, compounds such as VE821 (ATR), A205804 (ICAM-1), CC-401 (JNK), and L-779450 (Raf) exhibited GSK3 specific thermal stabilization in cells. The consistent and differential responses of CDK16 and the GSKs offers clues into compound specificity or the utilization of combined effects in cells.
Discussion
The determination of protein-compound engagement and mechanism of action is an essential aspect of the drug development process. Herein we present approachable methods for the implementation and interpretation of proteome-wide thermal shifts assays on at 96-well plate scale. Taking advantage of the 8-fold increase in efficiency afforded by PISA, we catalogued compound-dependent changes in protein thermal stability in response to 96 compounds in living K562 cells and 70 compounds in native extracts. This study resulted in a total of 871,120 thermal stability measurements in living cells and an additional 627,176 measurements in lysates. Exploiting this wealth of data, we were able to find clear evidence of on-target engagement for the majority of compounds tested in both cell- and lysate-based experiments. Furthermore, there were numerous potential examples of off-target engagement. In addition to direct target engagement, there were many secondary changes in protein thermal stability that occurred independent of direct binding. Finally, using an all-by-all analysis of compounds and proteins, we demonstrated how structurally related proteins and protein complexes engage consistently with our library of compounds. Our study serves as an unbiased assessment of the current technology and serves as a comprehensive guide for ongoing work to enable library-scale study of protein-compound and protein-ligand interactions.
One of the most valuable features of the proteome-wide thermal shift assay, including PISA, is the ability to screen compounds in both living cells and native lysates (Molina et al., 2013; Savitski et al., 2014). In this study, we demonstrate that these two approaches are not redundant, but, rather, provide complimentary data that can be employed to help establish compound mechanism of action. First and foremost, each approach was equally proficient in identifying direct ligand (inhibitor) binding events (Figures 1 and 3). In both cell- and lysate-based PISA experiments, about 70% of all compounds (whose target was quantified) induced a significant change in the thermal stability of its well-established target(s). Furthermore, when the two datasets were combined, that figure increased to over 80%. This is because certain targets, like BRAF, for example were more prone to undergo a change in thermal stability in one experimental setting rather than the other (Figure 4D).
In general, a compound-dependent thermal shift that occurs in a lysate-based experiment is almost certain to stem from direct target engagement. This is because cell signaling pathways and cellular structures are disrupted and diluted. Cell-based studies, on the other hand, have the added potential to identify secondary changes that occur independent of direct engagement. As a result of this dichotomy, we found that, on average, compounds caused a greater number of thermal shifts in cells than lysates. This is further evidenced by the PF-3758309-dependent thermal destabilization of the spliceosome. Finally, we documented a number of examples of phosphorylation-dependent thermal shifts. For example, compounds targeting CDKs tended to impact the thermal stability of RB1. Likewise, compounds targeting tyrosine kinases induced a thermal destabilization of downstream phosphorylation targets like CRKL and CRK. In these examples, an inhibitor engages and inhibits a direct target, which, in turn, prevents the phosphorylation of a downstream target thereby inducing a change in thermal stability.
The idea of area under the curve analysis (PISA) remains a relatively recent modification of the traditional proteome-wide thermal shift assay (Gaetani et al., 2019). Using this dataset, we were able to benchmark many aspects of this advance and highlight many important considerations when designing and interpreting these experiments. First, we highlight the importance of selecting an appropriate thermal window (Li et al., 2020). In this study we selected a temperature range of 48°C-58°C, however, it is possible that deploying multiple thermal windows could be a viable strategy for focused interrogation of a small number of high-value compounds. In utilizing thermal windows, the goal is to maximize the potential magnitudes of the final fold change measurements. Still, we observed a large range in the magnitudes of the log2 fold changes that were quantified for known compound targets. This is despite the fact, that we utilized a high treatment dose (10 µM) to maximize the magnitude of the fold change of the primary target. Ultimately, the data seemed to indicate that the maximal PISA fold change is protein-specific. Therefore, a log2 fold change of 2 for one compound-protein pair could be just as meaningful as a log2 fold change of 0.2 for another. In either case though we find that it is essential to optimize the thermal denaturation to ensure maximally useful effect sizes when comparing compounds.
While the primary screen was carried out at fixed dose, the increased throughput of PISA allowed for certain compounds to be assayed at multiple doses in a single experiment. In these instances, there was a clear does-dependent change in thermal stability of primary targets, off-targets, and secondary targets. This not only helped corroborate observations from the primary screen, but also could provide a qualitative assessment of relative compound potency. Specifically, the compounds that most strongly impacted the thermal stability of targets, also acted as the most potent inhibitors. In order to be a candidate for this type of study, a target must have a large maximal thermal shift (magnitude of log2 fold change) because there must be a large enough dynamic range to clearly resolve different doses. From this study, PIK3CB (maximal log2 fold change ∼0.2) would make a bad candidate, while p38α (maximal log2 fold change ∼2.0) would be a good candidate. Importantly, however, it should be possible to tune the thermal window in order to specifically increase the maximal log2 fold change of PI3KCB for dose-dependency studies.
Each of the compounds used in this study have a well-defined target and known mechanism of action. Furthermore, many of the compounds have been used extensively as tool compounds in basic research and as therapeutics in the treatment of multiple cancers. Despite the frequency with which many of these compounds have been used, the data generated in this study suggests that the current mechanistic understanding of many (if not all) of these compounds remains incomplete. For example, we found that the clinically-deployed CDK4/6 inhibitor palbociclib is capable of directly engaging and inhibiting PLK1. Despite being a weak inhibitor of PLK1, this activity likely still contributes to the mechanism of action of palbociclib and further supports the conclusions of a previous study that implicated palbociclib as a potential inhibitor of PLK1 (Hafner et al., 2019). Similarly, PISA helped reveal that the CHEK1 inhibitor AZD-7762 is, in fact, a highly promiscuous kinase inhibitor that, amongst many other kinases, appears to engage a number of tyrosine kinases and inhibits downstream signaling. While this compound has not been approved for clinical use, it has been used frequently as a tool compound to study the biological impact of CHEK1 inhibition. Additionally, the PLK1 inhibitor BI 2536 engages non-protein kinases such as PIP4K2A/B/C. As inhibition of PIP4K2s has been shown un hematological cells, this potentially suggests the mechanism of dose-limiting hematologic events upon patient treatment with BI 2536 (Frost et al., 2012; Lima et al., 2022). Overall, these examples, and the dataset as a whole, suggest that proteome-wide thermal shift assays and PISA, specifically, are valuable tools to help resolve the full range of compound activities. Despite this utility, however, it is important to note that these assays, alone, are unlikely to paint a full picture of mechanism of action. Instead, they represent a valuable tool that can be coupled with other MS-based and genetic approaches.
In the end, this study highlights the immense power of thermal stability assays in screening large collections of molecules to help define target engagement and mechanism of action on a proteome-wide scale. While this type of study is ideal for finding the targets of compounds that have been identified in high-throughput screening, they can provide valuable information in almost any context. This includes compounds that seemingly already have a well-defined target. Despite finding many ways to utilize this data, we expect that this dataset contains many interesting observations that can not only be used to better understand specific drugs, but also has the potential to reveal new biology. Moving forward, we hope that this data will serve as a valuable resource for those interested in drug development and the research community at large.
Acknowledgements
We would like to acknowledge the Gygi and Schweppe labs for advice and technical assistance concerning methods development, experimental implementation, and data interpretation. We would also like to acknowledge our funding sources: R35GM150919-01 (DKS), R01GM067945- 20 (SPG), Andy Hill CARE Distinguished Researcher Award (DKS), Cancer Consortium New Investigator Award (DKS). JGV is the Mark Foundation for Cancer Research Fellow of the Damon Runyon Cancer Research Institute (DRG-2359-19).
Methods
Cell culture
K562 cells were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum. For cell-based PISA experiments, cells were grown to approximately 1×106 cells/mL and immediately used for the assay. For lysate-based experiments, cells were grown to approximately 1×106 cells/mL, washed with phosphate-buffered saline, and flash frozen in liquid nitrogen. Cell pellets were stored at −80°C until ready for use.
HCT116 cell were cultured in Dulbecco’s Modified Eagle’s Medium supplemented with 10% fetal bovine serum and 1X penicillin-streptomycin. Cells were grown until ∼80% confluency, harvested by scraping, washed with phosphate-buffered saline, and flash frozen in liquid nitrogen. Cell pellets were stored at −80°C until ready for use.
SDS-PAGE and western blotting
K562 and HCT116 cell lysates were combined with Laemmli buffer and resolved on Novex WedgeWell 4-20% Tris-Glycine gels (Invitrogen). Gels were transferred to an Immuno-Blot PVDF membrane (Bio-Rad). Membranes were immunoblotted with antibodies against CRKL (Cell Signaling Technologies (CST), 38710), p-CRKL Y207 (CST, 3490), RIPK1 (CST, 3493), p-RIPK1 S116 (CST, 65746), TCTP (CST, 5128), and p-TCTP S46 (CST, 5251). Following primary antibody, membranes were incubated in goat anti-rabbit IgG-HRP secondary antibody (Santa Cruz Biotechnology, sc-2004).
Cell cycle synchronization
HCT116 cells were grown to approximately 50% confluency. The media was removed and replaced with media containing 100 ng/mL nocodazole (Cayman, 13857). The cells were allowed to incubate for 20 hours. After 20 hours, each compound was added at the desired concentration. The cells were allowed to further incubate for 1 hour. After 1 hour the cells were washed three times with PBS and lysed in RIPA buffer. 20 ug of each lysate was separated by SDS-PAGE and analyzed with immunoblot.
Initiation of necroptosis
HCT116 cells were grown to approximately 50% confluency. The media was removed and replaced with media containing 20 μM Z-VAD(Ome)-FMK (Cayman, 14463), 100 nM Birinapant (Cayman, 19699), and 10 ng/mL soluble recombinant human TNF-α (Cayman, 32020) (Laurien et al., 2020). Cells were incubated for 6 hours in the presence of various inhibitors or vehicle (DMSO). In the end, cells were washed three times with PBS and lysed in RIPA buffer. 20 ug of each lysate was separated by SDS-PAGE and analyzed with immunoblot.
Cell-based PISA
K562 cells were grown to a concentration of 1×106 cells/mL. The cells were pelleted and resuspended on a 1:1 mixture of conditioned media and fresh media for a final concentration of 6×106 cells/mL. Each compound was added to fresh media at a 3X concentration (30 μM). In order to initiate the experiment, 500 μL of cell suspension was mixed with 1 mL of 3X treatment media and plated in a 24 well untreated tissue culture plate to achieve a final cell concentration of 2×106 cells/mL and a compound concentration of 10 μM. The cells were allowed to incubate for 30 minutes. After incubation, an equal volume of each culture was transferred to 10 PCR tubes. The PCR tubes were heated across a thermal gradient ranging from 48°C-58°C (or 37°C-62°C) for 3 minutes to induce thermal denaturation. The samples were allowed to cool to room temperature and then an equal volume from each PCR tube was pooled. The cells in each pooled sample were washed once with PBS and then an equal volume of extraction buffer (1X PBS pH 7.4, 0.5% NP-40, protease inhibitors) was added to added to each pellet. Samples were incubated for 10 minutes at 4°C on a roller. Extracted samples were spun at 21,000xg for 90 minutes to separate insoluble aggregates from soluble protein. An equal volume from each soluble fraction was collected and prepared for LC-MS/MS analysis.
Crude extract PISA (Lysate-based PISA)
Frozen K562 or HCT116 pellets were thawed on ice and resuspended in lysis buffer (1X PBS pH 7.4, 1 mM MgCl2, protease inhibitor). The proteomes were extracted using a dounce homogenizer (20 strokes). The extracts were spun at 300xg for 3 minutes to remove any unbroken cells. The resulting crude extract was diluted to 2 mg/mL in lysis buffer. Each compound was added to lysis buffer at a 2X concentration (usually 20 μM). In order to initiate the experiment, an equal volume of crude extract and treatment buffer were combined, to achieve a final protein concentration of 1 mg/mL and compound concentration of 10 μM, and incubated for 30 minutes. After incubation, an equal volume of each sample was transferred to 10 PCR tubes. The PCR tubes were heated across a thermal gradient ranging from 48°C-58°C for 3 minutes to induce thermal denaturation. An equal volume from each PCR tube was pooled. An equal volume of extraction buffer (1X PBS pH 7.4, 1% NP-40, protease inhibitors) was added to added to each pooled sample to achieve a final NP-40 concentration of 0.5%. Samples were incubated for 10 minutes at 4°C on a roller. Extracted samples were spun at 21,000xg for 90 minutes to separate insoluble aggregates from soluble protein. An equal volume from each soluble fraction was collected and prepared for LC- MS/MS analysis.
LC-MS sample preparation
Samples (20 µg protein) were diluted in prep buffer (400 mM EPPS pH 8.5, 1% SDS, 10 mM tris(2-carboxyethyl)phosphine hydrochloride) and incubated at room temperature for 10 minutes. Iodoacetimide was added to a final concentration of 10 mM to each sample and incubated for 25 minutes in the dark. Finally, DTT was added to each sample to a final concentration of 10 mM. A buffer exchange was carried out using a modified SP3 protocol (Hughes et al., 2019, 2014). Briefly, ∼250 µg of Cytiva SpeedBead Magnetic Carboxylate Modified Particles (65152105050250 and 4515210505250), mixed at a 1:1 ratio, were added to each sample. 100% ethanol was added to each sample to achieve a final ethanol concentration of at least 50%. Samples were incubated with gentle shaking for 15 mins. Samples were washed three times with 80% ethanol. Protein was eluted from SP3 beads using 200 mM EPPS pH 8.5 containing trypsin (ThermoFisher Scientific, 90305R20) and Lys-C (Wako, 129-02541). Samples were digested overnight at 37°C with vigorous shaking. Acetonitrile was added to each sample to achieve a final concentration of ∼33%. Each sample was labelled, in the presence of SP3 beads, with ∼60 µg of TMTPro 16plex reagents (ThermoFisher Scientific). Following confirmation of satisfactory labelling (>97%), excess TMT was quenched by addition of hydroxylamine to a final concentration of 0.3%. The full volume from each sample was pooled and acetonitrile was removed by vacuum centrifugation for 1 hour. The pooled sample was acidified using formic acid and peptides were de-salted using a Sep-Pak 50mg tC18 cartridge (Waters). Peptides were eluted in 70% acetonitrile, 1% formic acid and dried by vacuum centrifugation. The peptides were resuspended in 10mM ammonium bicarbonate pH 8, 5% acetonitrile and fractionated by basic pH reverse phase HPLC. In total 24 fractions were collected. The fractions were dried in a vacuum centrifuge, resuspended in 5% acetonitrile, 1% formic acid and desalted by stage-tip. Finally, peptides were eluted in, 70% acetonitrile, 1% formic acid, dried, resuspended in 5% actetonitrile, 5% formic acid, and analyzed by LC-MS/MS.
Offline basic reversed phase fractionation
TMT labeled peptides were solubilized in 5% acetonitrile/10 mM ammonium bicarbonate, pH 8.0 and ∼300 µg of TMT labeled peptides were separated by an Agilent 300 Extend C18 column (3.5 μm particles, 4.6 mm ID and 250 mm in length). An Agilent 1260 binary pump coupled with a photodiode array (PDA) detector (ThermoFisher Scientific) was used to separate the peptides. A 45-minute linear gradient from 10% to 40% acetonitrile in 10 mM ammonium bicarbonate pH 8.0 (flow rate of 0.25 mL/min) separated the peptide mixtures into a total of 96 fractions (36 seconds). A total of 96 Fractions were consolidated into 24 samples in a checkerboard fashion and vacuum dried to completion.
Mass spectrometry data acquisition
Total proteome data were collected on Orbitrap Eclipse mass spectrometer (ThermoFisher Scientific) coupled to a Proxeon EASY-nLC 1000 (or 1200) LC pump (ThermoFisher Scientific). Peptides were separated using a 90-120 minute gradient at 500-550 nL/min on a 30 cm column (i.d. 100 μm, Accucore, 2.6 μm, 150 Å) packed inhouse. High-field asymmetric-waveform ion mobility spectroscopy (FAIMS) was enabled during data acquisition with compensation voltages (CVs) set as −40 V, −60 V, and −80 V (Schweppe et al., 2019). MS1 data were collected using the Orbitrap (Resolution – 60,000; Scan range – 400-1600 Th; Automatic gain control (AGC) −4×105; Normalized AGC target – 100%; maximum ion injection time – 50 ms). Determined charge states between 2 and 6 were required for sequencing, and a 90 s dynamic exclusion window was used. Data dependent mode was set as cycle time (1 s). MS2 scans were collected in the orbitrap after high-energy collision dissociation (HCD) fragmentation (Resolution – 50,000; AGC target – 1×105; Normalized AGC target – 200%; Normalized collision energy – 36; Isolation window – 0.5 Th; Maximum ion injection time – 86 ms).
Mass spectrometry data analysis
Raw files were first converted to mzXML, and monoisotopic peaks were assigned using Monocle (Rad et al., 2021). Database searching included all human entries from Uniprot (downloaded on February 25th, 2020). The database was concatenated with one composed of all protein sequences in the reversed order (Elias and Gygi, 2007). Sequences of common contaminant proteins (e.g., trypsin, keratins, etc.) were appended as well. Searches were performed with Comet (Eng et al., 2013) using a 50 ppm precursor ion tolerance and 0.02 Da product ion tolerance. TMT on lysine residues and peptide N termini (+304.207 Da) and carbamidomethylation of cysteine residues (+57.0215 Da) were set as static modifications, while oxidation of methionine residues (+15.9949 Da) was set as a variable modification. Peptide-spectrum matches (PSMs) were adjusted to a 1% false discovery rate (FDR; Elias and Gygi, 2007) PSM filtering was performed using linear discriminant analysis (LDA) as described previously (Huttlin et al., 2010) while considering the following parameters: XCorr, ΔCn, missed cleavages, peptide length, charge state, and precursor mass accuracy. Each run was filtered separately. Protein-level FDR was subsequently estimated at a data set level. For each protein across all samples, the posterior probabilities reported by the LDA model for each peptide were multiplied to give a protein-level probability estimate. Using the Picked FDR method (Savitski et al., 2015) proteins were filtered to the target 1% FDR level. TMT reporter ion intensities were measured using a 0.003 Da window around the theoretical m/z of each reporter ion. Proteins were quantified by summing reporter ion counts across all matching PSMs. Reporter ion intensities were adjusted to correct for the isotopic impurities of the different TMT reagents according to manufacturer specifications. Peptides were filtered to exclude those with a summed signal-to-noise (SN) > 160 across all TMT channels. To control for different total protein loading within a TMT experiment, the summed protein quantities of each channel were adjusted to be equal within the experiment.
Protein engagement hit calling
Compound engagement was determined based on relative thermal stability to DMSO controls. Owing to the fact that the duplicate analyses used in the initial screen would result in common statistical tests (Welch’s t-test) being underpowered, we used a combination of fold changes compared to DMSO and individual protein variance to call hits. First, relative thermal stability to the DMSO controls was determined. Second, for each protein across all cells or lysate assays, the number of standard deviations away from the mean thermal stability measurement (z-score) for a given protein was quantified. Cutoffs for fold change and z-score were determined to limit the number of hits derived from DMSO-treated samples. We maintained the same cutoffs across both datasets. We considered proteins to engage a compound if both replicates of the compound treatments resulted when the thermal stability fold change compared to DMSO greater than an absolute value of log20.2 and an absolute z-score greater than 3.5. This resulted in a false hit rate (nDMSO-Hits) of 1% across all lysate-based assays and 3% for all cell-based assays.
Statistical analyses
Follow-up PISA data were analyzed using Perseus (Tyanova et al., 2016). Significant changes were determined using a permutation-based FDR with the following settings – FDR – 0.05, S0 – 0.1, and number of randomizations – 250. Individual fold change values were calculated in reference to the mean of the vehicle-treated samples.
Correlation analyses were run in R version 4.3.1 using rank-based Spearman’s rho (rspearman) to minimize the effects of thermal stability-based outliers. The list of binary comparisons was filtered to include the top 5% of highly correlated and anticorrelate protein pairs. Significant pairwise correlations for individual comparisons were determined using the psych package.
Structural modeling
Protein structures and associated ligands were modeled using the ICMBrowser v3.9 (Molsoft). PDB files were downloaded from RSCB PDB for the following identifiers: 2rku, 5l2i, 5g6v.
Data availability
The mass spectrometry data have been deposited to the ProteomeXchange Consortium with the data set identifier PDX048009. Information regarding TMT channel assignments and plex layout can be found in Table 2.
Table 1. Library compounds. A list of all compounds utilized in the cell- and lysate-based PISA screens.
Table 2. Sample multiplexing layout. This table contains the TMTPro channel assignments for all samples generated in the cell- and lysate-based screens.
Figure 1 – source data 1. This table contains cell-based PISA data for K562 cells treated with ribociclib, palbociclib, and abemaciclib.
Figure 1 – source data 2. This table contains the PISA data for cell- and lysate-based screens.
Figure 2 – source data 1. This table contains cell-based PISA data for K562 cells treated with increasing concentrations of BI 2536 and palbocilicb.
Figure 2 – source data 2. This table contains cell-based PISA data for K562 cells treated with increasing concentrations of BI 2536 and NVP-TAE-226.
Figure 2 – source data 3. This document contains the unedited western scans for the panels displayed in Figure 2F and Figure 2 – figure supplement 1D.
Figure 4 – source data 1. This table contains lysate-based PISA data for HCT116 lysates treated with necrostatin-2, GSK2606414, AZD-5438, and tozasertib.
Figure 4 – source data 2. This document contains the unedited western scans for the panel displayed in Figure 4G.
Figure 5 – source data 1. This document contains the unedited western scans for the panels displayed in Figure 5C.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig1-1.tif/full/max/0/default.jpg)
Establishing a robust and efficient workflow for assessing changes in protein thermal stability in living cells.
A. Schematic depicting how thermal shifts are measured in a PISA experiment. An increase in melting temperature (thermal stabilization) will result in an increase in the AUC and a positive log2 fold change. A decrease in melting temperature (thermal destabilization) will result in a decrease in the AUC and a negative log2 fold change. B. PISA schematic depicting a thermal window vs. a full melting curve. C-G. K562 cells (N=4) were treated with 10 µM palbociclib (D), ribociclib (E), or abemaciclib (F) for 15 minutes and assayed using PISA. (C) The log2 fold changes for specific proteins are plotted for each treatment using a full melting curve (left) or thermal window (48°C-58°C; right). (D-F). Data is presented as a volcano plot to highlight significant changes in abundance utilizing the thermal window (48°C-58°C). Significant changes were determined using a permutation-based FDR (FDR – 0.05, S0 – 0.1). (G). Log2 fold change values for selected proteins following treatment with each compound in D-F.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig1-2.tif/full/max/0/default.jpg)
Establishing a robust and efficient workflow for assessing changes in protein thermal stability in living cells.
K562 cells were treated with each of the 96 compounds at 10 µM for 30 minutes. Three compounds and a DMSO control were assayed in duplicate (batch). All 96 compounds and 32 DMSO controls were assayed using thirty-two treatment batches. Following treatment, an equal number of cells were transferred to 10 PCR tubes. The cells were placed in a thermal cycler and heated across a thermal gradient from 48°C- 58°C for 3 minutes. The cells were allowed to cool to room temperature for 5 minutes. An equal volume from each PCR tube was pooled and spun at 300 x g for 3 minutes to pellet cells. Cells were washed one with PBS and lysed in a buffer containing 0.5% NP-40, which will dissolve membranes without disrupting heat-induced protein aggregates. The lysates were centrifuged for 90 minutes at 21,000 x g to separate soluble protein from aggregates. An equal volume of each soluble fraction (∼20 µg) was prepared for LC-MS/MS analysis. Two treatment batches were combined for each TMTpro 16-plex. Each soluble fraction was reduced and alkylated. Each sample was precipitated onto SP3 carboxylate-coated beads to facilitate a buffer exchange. Proteomes were eluted off the SP3 beads into digestion buffer and digested with a combination of Lys-C and trypsin. Peptides from each sample were labeled with a unique TMTpro reagent. TMT-labelled peptides were pooled into a single sample, which was desalted using a sep-pak. Dried peptides were resuspended in HPLC buffer A and fractionated by basic reverse-phase HPLC. Twelve to twenty-four fractions were stage-tipped and analyzed on an Orbitrap Eclipse with a FAIMS device enabled (Thermo Fisher). Changes in thermal stability were determined by comparing compound-treated samples to vehicle-treated controls.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig1-3.tif/full/max/0/default.jpg)
Establishing a robust and efficient workflow for assessing changes in protein thermal stability in living cells.
A. Minimal instrument time (assuming 12 fractions per plex) required to assay 96 compounds and 32 DMSO controls in duplicate using CETSA (left) or PISA (right). B. Total proteins quantified per TMTpro 16-plex. C. Total thermal stability measurements. D. Schematic depiction of the filtering scheme used to define significant changes. E. Histogram depicting the log2 fold change measurements that were greater than 3.5 standard deviations from the mean. F. Total number of proteins passing each filter. G. Pie charts depicting the fraction of all measurements that pass each filter.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig1-4.tif/full/max/0/default.jpg)
Establishing a robust and efficient workflow for assessing changes in protein thermal stability in living cells.
A-B. K562 cells (N=2) were treated with MK-2206 (top) or CCT128930 (bottom) at 10 µM for 30 minutes and assayed using PISA. The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. Blue boxes contain proteins that have been stabilized (increase in melting temperature) and orange boxes contain proteins that have been destabilized (decrease in melting temperature). Orange points represent the known target of each compound. B. The plot displays the log2 fold change values for selected proteins following treatment with each compound in A. C. Plot highlighting the significant changes in protein thermal stability following treatment with PF- 3758309. The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. The blue box contains proteins that have been stabilized (increase in melting temperature) and the orange box contains proteins that have been destabilized (decrease in melting temperature). The blue point represents the known target. D. Plot comparing the log2 fold changes measurements for each protein in JQ-1 (x-axis)- and volaserib (y-axis)-treated K562 cells. Orange points represent the known targets of the compounds.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig2-1.tif/full/max/0/default.jpg)
Using cell-based PISA data to assess compound target engagement.
A-B. A protein-centric view of p38α (A) and PIK3CB (B) thermal stability in response to all treatments. Log2 fold change is plotted on the y-axis. The points represent each of the 256 treatments that were performed. Orange points indicate compounds known to target each respective protein. Blue points indicate other compounds that result in a significant thermal stabilization. Green dashed lines mark a SD of 3.5 cutoff for each treatment. C-D. K562 cells were treated with the indicated concentrations of NVP-TAE-226 or BI 2536 for 15 minutes. Changes in protein thermal stability are represented as a log2 fold change in protein abundance for each treatment in reference to a DMSO-treated control (C). PLK1 activity was determined using a western blot for p-TCTP (S46) and total TCTP (D).
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig3-1.tif/full/max/0/default.jpg)
A chemical perturbation screen of 70 compounds in K562 native extracts.
A. Minimal instrument time required to assay 70 compounds and 10 DMSO controls in duplicate using TPP (left) or PISA (right). B. Total number of proteins quantified in each plex in the lysate-based screen. C. A plot of the mean log2 fold change (x-axis) and mean nSD (y-axis) of each known compound-target pair that were quantified in the lysate-based screen. The blue box contains proteins that have been stabilized (increase in melting temperature) and the orange box contains proteins that have been destabilized (decrease in melting temperature). D-G. Protein-centric view of HDAC1 (D), p38α (E), CDK2 (F), and CDK7 (G) thermal stability in response to all treatments. Orange points indicate compounds known to target each respective protein. Blue points indicate other compounds that result in a significant thermal destabilization of each respective protein. Green dashed lines mark a SD of 3.5 cutoff for each compound.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig4-1.tif/full/max/0/default.jpg)
Cell- and lysate-based PISA are complimentary approaches for assessing mechanism of action.
A. Changes in protein thermal stability for the indicated proteins in cell- and lysate-based PISA are represented as a log2 fold change in protein abundance for each treatment in reference to a DMSO-treated control. Dashed lines indicate a log2 fold change of +/-0.2. B. A protein-centric view of RIPK1 thermal stability in response to all treatments. Log2 fold change is plotted on the y-axis. The points represent each of the 256 treatments that were performed. Blue points highlight compounds that result in a significant thermal stabilization. Green dashed lines mark a SD of 3.5 cutoff for RIPK1. C-F. HCT116 lysates were treated with necrostatin-2 (C), GSK2606414 (D), AZD-5438 (E), or tozasertib (F) (N=3) or DMSO (N=4) for 15 minutes and any changes in thermal stability were determined using PISA. Data is presented as a volcano plot to highlight significant changes in abundance. Significant changes were determined using a permutation-based FDR (FDR – 0.05, S0 – 0.1).
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig5-1.tif/full/max/0/default.jpg)
Disparities in cell- and lysate-based PISA can point toward secondary changes in protein thermal stability.
A. Plot highlighting the significant changes in protein thermal stability following treatment with PF-3758309 in cells (no fill) and lysates (black points). The mean log2 fold change of duplicate measurements is plotted on the x-axis and the mean nSD is plotted on the y-axis. The blue box contains proteins that have been stabilized (increase in melting temperature) and the orange box contains proteins that have been destabilized (decrease in melting temperature). B and C. A protein-centric view of CRKL thermal stability in K562 cells (B) and lysates (C) in response to all treatments. Log2 fold change is plotted on the y-axis. The points represent each treatment that was performed in each screen. Green dashed lines mark a SD of 3.5 cutoff for each treatment. D. Plots depicting the log2 fold change values for CRK in response to treatment with all compounds assayed using cell- and lysate-based PISA. E. Changes in protein thermal stability for the indicated proteins in cell- and lysate-based PISA are represented as a log2 fold change in protein abundance for each treatment in reference to a DMSO-treated control.
![](https://prod--epp.elifesciences.org/iiif/2/95595%2Fv1%2Fcontent%2F577428v1_fig6-1.tif/full/max/0/default.jpg)
Concordance of protein engagement for protein complex members and structurally similar kinases.
A. Subgraph from Figure 7A highlighting the correlation between PRKC kinases. Edges are colored based on RSpearman. B. Correlation of PRKCB and either PRKCA or PRKCQ. These correlations are driven by consistent but relatively small changes in thermal stability. C. Non-protein kinases of the correlation network graph. Library compounds did not specifically target any non-protein kinases though these compounds elicited consistent thermal stability changes in these kinases. Comparison of compound engagement for PIP4K2 non-protein kinases (D) PIP4K2A and PIP4K2B as well as (E) PIP4K2B and PIP4K2C. Compound engagement generates statistically significant correlations in both cells and lysates, with a larger thermal stability effect size in cell-based assays. F-G. Significant correlation with small thermal stability fold changes for the non-protein kinases SRPK1, SRPK2, and RIOK2.
References
- Molecular insights into the behavior of the allosteric and ATP-competitive inhibitors in interaction with AKT1 protein: A molecular dynamics studyInternational Journal of Biological Macromolecules 242https://doi.org/10.1016/j.ijbiomac.2023.124853
- An isothermal shift assay for proteome scale drug-target identificationCommunications Biology 3https://doi.org/10.1038/s42003-020-0795-6
- Ion-Based Proteome-Integrated Solubility Alteration Assays for Systemwide Profiling of Protein–Molecule InteractionsAnal. Chem 94:7066–7074https://doi.org/10.1021/acs.analchem.2c00391
- AZD5438, a potent oral inhibitor of cyclin-dependent kinases 1, 2, and 9, leads to pharmacodynamic changes and potent antitumor effects in human tumor xenograftsMolecular Cancer Therapeutics 8:1856–1866https://doi.org/10.1158/1535-7163.MCT-08-0836
- Spectrum and Degree of CDK Drug Interactions Predicts Clinical PerformanceMolecular Cancer Therapeutics 15:2273–2281https://doi.org/10.1158/1535-7163.MCT-16-0300
- Quantitative Proteomics of Kinase Inhibitor Targets and MechanismsACS Chem. Biol 10:201–212https://doi.org/10.1021/cb5008794
- KinMap: a web-based tool for interactive navigation through human kinome dataBMC Bioinformatics 18https://doi.org/10.1186/s12859-016-1433-7
- Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometryNature Methods 4:207–214https://doi.org/10.1038/nmeth1019
- Comet: An open-source MS/MS sequence database search toolPROTEOMICS 13:22–24https://doi.org/10.1002/pmic.201200439
- CDK versus GSK-3 Inhibition: a Purple Haze No Longer?Chemistry & Biology 10:1144–1146https://doi.org/10.1016/j.chembiol.2003.12.009
- Proteomics in Drug Development: The Dawn of a New Era?PROTEOMICS – Clinical Applications 13https://doi.org/10.1002/prca.201800087
- Phase I Study of the Plk1 Inhibitor BI 2536 Administered Intravenously on Three Consecutive Days in Advanced Solid TumoursCurrent Oncology 19:28–35https://doi.org/10.3747/co.19.866
- Proteome Integral Solubility Alteration: A High-Throughput Proteomics Assay for Target DeconvolutionJ. Proteome Res 18:4027–4037https://doi.org/10.1021/acs.jproteome.9b00500
- MELK expression correlates with tumor mitotic activity but is not required for cancer growtheLife 7https://doi.org/10.7554/eLife.32838
- Multiomics Profiling Establishes the Polypharmacology of FDA-Approved CDK4/6 Inhibitors and the Potential for Differential Clinical ActivityCell Chemical Biology 26:1067–1080https://doi.org/10.1016/j.chembiol.2019.05.005
- MK-2206, an Allosteric Akt Inhibitor, Enhances Antitumor Efficacy by Standard Chemotherapeutic Agents or Molecular Targeted Drugs In vitro and In vivoMol Cancer Ther 9https://doi.org/10.1158/1535-7163.MCT-09-1012
- Ultrasensitive proteome analysis using paramagnetic bead technologyMolecular Systems Biology 10https://doi.org/10.15252/msb.20145625
- Single-pot, solid-phase-enhanced sample preparation for proteomics experimentsNature Protocols 14:68–85https://doi.org/10.1038/s41596-018-0082-x
- Dual proteome-scale networks reveal cell-specific remodeling of the human interactomeCell 184:3022–3040https://doi.org/10.1016/j.cell.2021.04.011
- Architecture of the human interactome defines protein communities and disease networksNature 545:505–509https://doi.org/10.1038/nature22366
- A Tissue-Specific Atlas of Mouse Protein Phosphorylation and ExpressionCell 143:1174–1189https://doi.org/10.1016/j.cell.2010.12.001
- The cellular thermal shift assay for evaluating drug target interactions in cellsNature Protocols 9:2100–2122https://doi.org/10.1038/nprot.2014.138
- Meltome atlas—thermal proteome stability across the tree of lifeNature Methods 17:495–503https://doi.org/10.1038/s41592-020-0801-4
- Application of a Theory of Enzyme Specificity to Protein Synthesis*Proceedings of the National Academy of Sciences 44:98–104https://doi.org/10.1073/pnas.44.2.98
- Research Article: Selectivity-determining Residues in Plk1Chemical Biology & Drug Design 70:540–546https://doi.org/10.1111/j.1747-0285.2007.00594.x
- Reimagining high-throughput profiling of reactive cysteines for cell-based screening of large electrophile librariesNature Biotechnology 39:630–641https://doi.org/10.1038/s41587-020-00778-3
- Autophosphorylation at serine 166 regulates RIP kinase 1-mediated cell death and inflammationNature Communications 11https://doi.org/10.1038/s41467-020-15466-8
- Selection of Heating Temperatures Improves the Sensitivity of the Proteome Integral Solubility Alteration AssayJ. Proteome Res 19:2159–2166https://doi.org/10.1021/acs.jproteome.0c00063
- The PIP4K2 inhibitor THZ-P1-2 exhibits antileukemia activity by disruption of mitochondrial homeostasis and autophagyBlood Cancer Journal 12https://doi.org/10.1038/s41408-022-00747-w
- Large-Scale Analysis of Breast Cancer-Related Conformational Changes in Proteins Using SILAC-SPROXJ. Proteome Res 16:3277–3286https://doi.org/10.1021/acs.jproteome.7b00283
- RIPK1-dependent cell death: a novel target of the Aurora kinase inhibitor Tozasertib (VX-680)Cell Death & Disease 9https://doi.org/10.1038/s41419-017-0245-7
- GSK-3-Selective Inhibitors Derived from Tyrian Purple IndirubinsChemistry & Biology 10:1255–1266https://doi.org/10.1016/j.chembiol.2003.11.010
- The emerging role of mass spectrometry-based proteomics in drug discoveryNature Reviews Drug Discovery 21:637–654https://doi.org/10.1038/s41573-022-00409-3
- How Chemoproteomics Can Enable Drug Discovery and DevelopmentChemistry & Biology 19:11–22https://doi.org/10.1016/j.chembiol.2012.01.001
- Monitoring Drug Target Engagement in Cells and Tissues Using the Cellular Thermal Shift AssayScience 341https://doi.org/10.1126/science.1233606
- Cyclin D activates the Rb tumor suppressor by mono-phosphorylationeLife 3https://doi.org/10.7554/eLife.02872
- Improved Monoisotopic Mass Estimation for Deeper Proteome CoverageJ. Proteome Res 20:591–598https://doi.org/10.1021/acs.jproteome.0c00563
- Quantifying Target Occupancy of Small Molecules Within Living CellsAnnu. Rev. Biochem 89:557–581https://doi.org/10.1146/annurev-biochem-011420-092302
- When PERK inhibitors turn out to be new potent RIPK1 inhibitors: critical issues on the specificity and use of GSK2606414 and GSK2656157Cell Death & Differentiation 24:1100–1110https://doi.org/10.1038/cdd.2017.58
- Tracking cancer drugs in living cells by thermal profiling of the proteomeScience 346https://doi.org/10.1126/science.1255784
- A Scalable Approach for Protein False Discovery Rate Estimation in Large Proteomic Data SetsMol Cell Proteomics 14:2394–2404https://doi.org/10.1074/mcp.M114.046995
- Small-Molecule Target Engagement in CellsCell Chemical Biology 23:435–441https://doi.org/10.1016/j.chembiol.2016.03.008
- Characterization and Optimization of Multiplexed Quantitative Analyses Using High-Field Asymmetric-Waveform Ion Mobility Mass SpectrometryAnal. Chem 91:4010–4016https://doi.org/10.1021/acs.analchem.8b05399
- An exon skipping screen identifies antitumor drugs that are potent modulators of pre-mRNA splicing, suggesting new therapeutic applicationsPLOS ONE 15https://doi.org/10.1371/journal.pone.0233672
- Proteome-wide solubility and thermal stability profiling reveals distinct regulatory roles for ATPNature Communications 10https://doi.org/10.1038/s41467-019-09107-y
- The Perseus computational platform for comprehensive analysis of (prote)omics dataNature Methods 13:731–740https://doi.org/10.1038/nmeth.3901
- Phosphorylation of TCTP as a Marker for Polo-like Kinase-1 ActivityVivo. Anticancer Res 30
- Assessing target engagement using proteome-wide solvent shift assayseLife 10https://doi.org/10.7554/eLife.70784
- Preclinical Pharmacology, Antitumor Activity, and Development of Pharmacodynamic Markers for the Novel, Potent AKT Inhibitor CCT128930Molecular Cancer Therapeutics 10:360–371https://doi.org/10.1158/1535-7163.MCT-10-0760
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Van Vranken et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 402
- downloads
- 32
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.