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

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

Bassoon contributes to tau-seed propagation and neurotoxicity

Tau aggregation is a defining histopathological feature of Alzheimer's disease and other tauopathies. However, the cellular mechanisms involved in tau propagation remain unclear. Here, we performed an unbiased quantitative proteomic study to identify proteins that specifically interact with this tau seed. We identified Bassoon (BSN), a presynaptic scaffolding protein, as an interactor of the tau seed isolated from a mouse model of tauopathy, and from Alzheimer's disease and progressive supranuclear palsy postmortem samples. We show that BSN exacerbates tau seeding and toxicity in both mouse and Drosophila models for tauopathy, and that BSN downregulation decreases tau spreading and overall disease pathology, rescuing synaptic and behavioral impairments and reducing brain atrophy. Our findings improve the understanding of how tau seeds can be stabilized by interactors such as BSN. Inhibiting tau-seed interactions is a potential new therapeutic approach for neurodegenerative tauopathies.

Tumor collection/processing under physioxia uncovers highly relevant signaling networks and drug sensitivity

Preclinical studies of primary cancer cells are typically done after tumors are removed from patients or animals at ambient atmospheric oxygen (O2, ~21%). However, O2 concentrations in organs are in the ~3 to 10% range, with most tumors in a hypoxic or 1 to 2% O2 environment in vivo. Although effects of O2 tension on tumor cell characteristics in vitro have been studied, these studies are done only after tumors are first collected and processed in ambient air. Similarly, sensitivity of primary cancer cells to anticancer agents is routinely examined at ambient O2. Here, we demonstrate that tumors collected, processed, and propagated at physiologic O2 compared to ambient air display distinct differences in key signaling networks including LGR5/WNT, YAP, and NRF2/KEAP1, nuclear reactive oxygen species, alternative splicing, and sensitivity to targeted therapies. Therefore, evaluating cancer cells under physioxia could more closely recapitulate their physiopathologic status in the in vivo microenvironment.

Ubiquitination and degradation of SUMO1 by small-molecule degraders extends survival of mice with patient-derived tumors

[Figure: see text].

Boosting Detection of Low-Abundance Proteins in Thermal Proteome Profiling Experiments by Addition of an Isobaric Trigger Channel to TMT Multiplexes

The study of low-abundance proteins is a challenge to discovery-based proteomics. Mass spectrometry (MS) applications, such as thermal proteome profiling (TPP), face specific challenges in the detection of the whole proteome as a consequence of the use of nondenaturing extraction buffers. TPP is a powerful method for the study of protein thermal stability, but quantitative accuracy is highly dependent on consistent detection. Therefore, TPP can be limited in its amenability to study low-abundance proteins that tend to have stochastic or poor detection by MS. To address this challenge, we incorporated an affinity-purified protein complex sample at submolar concentrations as an isobaric trigger channel into a mutant TPP (mTPP) workflow to provide reproducible detection and quantitation of the low-abundance subunits of the cleavage and polyadenylation factor (CPF) complex. The inclusion of an isobaric protein complex trigger channel increased detection an average of 40× for previously detected subunits and facilitated detection of CPF subunits that were previously below the limit of detection. Importantly, these gains in CPF detection did not cause large changes in melt temperature (T_m) calculations for other unrelated proteins in the samples, with a high positive correlation between T_m estimates in samples with and without isobaric trigger channel addition. Overall, the incorporation of an affinity-purified protein complex as an isobaric trigger channel within a tandem mass tag (TMT) multiplex for mTPP experiments is an effective and reproducible way to gather thermal profiling data on proteins that are not readily detected using the original TPP or mTPP protocols.

Atractylenolide I enhances responsiveness to immune checkpoint blockade therapy by activating tumor antigen presentation

One of the primary mechanisms of tumor cell immune evasion is the loss of antigenicity, which arises due to lack of immunogenic tumor antigens as well as dysregulation of the antigen processing machinery. In a screen for small-molecule compounds from herbal medicine that potentiate T cell-mediated cytotoxicity, we identified atractylenolide I (ATT-I), which substantially promotes tumor antigen presentation of both human and mouse colorectal cancer (CRC) cells and thereby enhances the cytotoxic response of CD8+ T cells. Cellular thermal shift assay (CETSA) with multiplexed quantitative mass spectrometry identified the proteasome 26S subunit non-ATPase 4 (PSMD4), an essential component of the immunoproteasome complex, as a primary target protein of ATT-I. Binding of ATT-I with PSMD4 augments the antigen-processing activity of immunoproteasome, leading to enhanced MHC-I-mediated antigen presentation on cancer cells. In syngeneic mouse CRC models and human patient-derived CRC organoid models, ATT-I treatment promotes the cytotoxicity of CD8+ T cells and thus profoundly enhances the efficacy of immune checkpoint blockade therapy. Collectively, we show here that targeting the function of immunoproteasome with ATT-I promotes tumor antigen presentation and empowers T cell cytotoxicity, thus elevating the tumor response to immunotherapy.

Inflect: Optimizing Computational Workflows for Thermal Proteome Profiling Data Analysis

The CETSA and Thermal Proteome Profiling (TPP) analytical methods are invaluable for the study of protein-ligand interactions and protein stability in a cellular context. These tools have increasingly been leveraged in work ranging from understanding signaling paradigms to drug discovery. Consequently, there is an important need to optimize the data analysis pipeline that is used to calculate protein melt temperatures (Tm) and relative melt shifts from proteomics abundance data. Here, we report a user-friendly analysis of the melt shift calculation workflow where we describe the impact of each individual calculation step on the final output list of stabilized and destabilized proteins. This report also includes a description of how key steps in the analysis workflow quantitatively impact the list of stabilized/destabilized proteins from an experiment. We applied our findings to develop a more optimized analysis workflow that illustrates the dramatic sensitivity of chosen calculation steps on the final list of reported proteins of interest in a study and have made the R based program Inflect available for research community use through the CRAN repository [McCracken, N. Inflect: Melt Curve Fitting and Melt Shift Analysis. R package version 1.0.3, 2021]. The Inflect outputs include melt curves for each protein which passes filtering criteria in addition to a data matrix which is directly compatible with downstream packages such as UpsetR for replicate comparisons and identification of biologically relevant changes. Overall, this work provides an essential resource for scientists as they analyze data from TPP and CETSA experiments and implement their own analysis pipelines geared toward specific applications.

Mutant thermal proteome profiling for characterization of missense protein variants and their associated phenotypes within the proteome

Temperature-sensitive (TS) missense mutants have been foundational for characterization of essential gene function. However, an unbiased approach for analysis of biochemical and biophysical changes in TS missense mutants within the context of their functional proteomes is lacking. We applied MS-based thermal proteome profiling (TPP) to investigate the proteome-wide effects of missense mutations in an application that we refer to as mutant thermal proteome profiling (mTPP). This study characterized global impacts of temperature sensitivity-inducing missense mutations in two different subunits of the 26S proteasome. The majority of alterations identified by RNA-Seq and global proteomics were similar between the mutants, which could suggest that a similar functional disruption is occurring in both missense variants. Results from mTPP, however, provide unique insights into the mechanisms that contribute to the TS phenotype in each mutant, revealing distinct changes that were not obtained using only steady-state transcriptome and proteome analyses. Computationally, multisite λ-dynamics simulations add clear support for mTPP experimental findings. This work shows that mTPP is a precise approach to measure changes in missense mutant-containing proteomes without the requirement for large amounts of starting material, specific antibodies against proteins of interest, and/or genetic manipulation of the biological system. Although experiments were performed under permissive conditions, mTPP provided insights into the underlying protein stability changes that cause dramatic cellular phenotypes observed at nonpermissive temperatures. Overall, mTPP provides unique mechanistic insights into missense mutation dysfunction and connection of genotype to phenotype in a rapid, nonbiased fashion.

Phosphoproteomic analysis identifies phospho-Threonine-17 site of phospholamban important in low molecular weight isoform of fibroblast growth factor 2-induced protection against post-ischemic cardiac dysfunction

RATIONALE

Among its many biological roles, fibroblast growth factor 2 (FGF2) protects the heart from dysfunction and damage associated with an ischemic attack. Our laboratory demonstrated that its protection against myocardial dysfunction occurs by the low molecular weight (LMW) isoform of FGF2, while the high molecular weight (HMW) isoforms are associated with a worsening in post-ischemic recovery of cardiac function. LMW FGF2-mediated cardioprotection is facilitated by activation of multiple kinases, including PKCalpha, PKCepsilon, and ERK, and inhibition of p38 and JNK.

OBJECTIVE

Yet, the substrates of those kinases associated with LMW FGF2-induced cardioprotection against myocardial dysfunction remain to be elucidated.

METHODS AND RESULTS

To identify substrates in LMW FGF2 improvement of post-ischemic cardiac function, mouse hearts expressing only LMW FGF2 were subjected to ischemia-reperfusion (I/R) injury and analyzed by a mass spectrometry (MS)-based quantitative phosphoproteomic strategy. MS analysis identified 50 phosphorylation sites from 7 sarcoendoplasmic reticulum (SR) proteins that were significantly altered in I/R-treated hearts only expressing LMW FGF2 compared to those hearts lacking FGF2. One of those phosphorylated SR proteins identified was phospholamban (PLB), which exhibited rapid, increased phosphorylation at Threonine-17 (Thr17) after I/R in hearts expressing only LMW FGF2; this was further validated using Selected Reaction Monitoring-based MS workflow. To demonstrate a mechanistic role of phospho-Thr17 PLB in LMW FGF2-mediated cardioprotection, hearts only expressing LMW FGF2 and those expressing only LMW FGF2 with a mutant PLB lacking phosphorylatable Thr17 (Thr17Ala PLB) were subjected to I/R. Hearts only expressing LMW FGF2 showed significantly improved recovery of cardiac function following I/R (p < 0.05), and this functional improvement was significantly abrogated in hearts expressing LMW FGF2 and Thr17Ala PLB (p < 0.05).

CONCLUSION

The findings indicate that LMW FGF2 modulates intracellular calcium handling/cycling via regulatory changes in SR proteins essential for recovery from I/R injury, and thereby protects the heart from post-ischemic cardiac dysfunction.

Endoplasmic reticulum stress alters ryanodine receptor function in the murine pancreatic β cell

Alterations in endoplasmic reticulum (ER) calcium (Ca2+) levels diminish insulin secretion and reduce β-cell survival in both major forms of diabetes. The mechanisms responsible for ER Ca2+ loss in β cells remain incompletely understood. Moreover, a specific role for either ryanodine receptor (RyR) or inositol 1,4,5-triphosphate receptor (IP3R) dysfunction in the pathophysiology of diabetes remains largely untested. To this end, here we applied intracellular and ER Ca2+ imaging techniques in INS-1 β cells and isolated islets to determine whether diabetogenic stressors alter RyR or IP3R function. Our results revealed that the RyR is sensitive mainly to ER stress-induced dysfunction, whereas cytokine stress specifically alters IP3R activity. Consistent with this observation, pharmacological inhibition of the RyR with ryanodine and inhibition of the IP3R with xestospongin C prevented ER Ca2+ loss under ER and cytokine stress conditions, respectively. However, RyR blockade distinctly prevented β-cell death, propagation of the unfolded protein response (UPR), and dysfunctional glucose-induced Ca2+ oscillations in tunicamycin-treated INS-1 β cells and mouse islets and Akita islets. Monitoring at the single-cell level revealed that ER stress acutely increases the frequency of intracellular Ca2+ transients that depend on both ER Ca2+ leakage from the RyR and plasma membrane depolarization. Collectively, these findings indicate that RyR dysfunction shapes ER Ca2+ dynamics in β cells and regulates both UPR activation and cell death, suggesting that RyR-mediated loss of ER Ca2+ may be an early pathogenic event in diabetes.

Phosphopeptide separation using radially aligned titania nanotubes on titanium wire

Phosphoproteomic analysis offers a unique view of cellular function and regulation in biological systems by providing global measures of a key cellular regulator in the form of protein phosphorylation. Understanding the phosphorylation changes between normal and diseased cells or tissues offers a window into the mechanism of disease and thus potential targets for therapeutic intervention. A key step in these studies is the enrichment of phosphorylated peptides that are typically separated and analyzed by using liquid chromatography mass spectrometry. The mesoporous titania beads/particles (e.g., Titansphere TiO2 beads from GL Sciences Inc., Japan) that are widely used for phosphopeptide enrichment are expensive and offer very limited opportunities for further performance improvement. Titiania nanotube arrays have shown promising characteristics for phosphopeptide separation. Here we report a proof-of-concept study to evaluate the efficacy of nanotubes on Ti-wire for phosphoproteomics research. We used titania nanotubes radially grown on titanium wires as well as the commercial beads to separate phosphopeptides generated from mouse liver complex tissue extracts. Our studies revealed that the nanotubes on metal wire provide comparable efficacy for enrichment of phophopeptides and offer an ease of use advantage versus mesoporous beads, thus having the potential to become a low cost and more practical material/methodology for phosphopeptide enrichment in biological studies.

Reversed phase liquid chromatography hyphenated to continuous flow-extractive desorption electrospray ionization-mass spectrometry for analysis and charge state manipulation of undigested proteins

The application of continuous flow-extractive desorption electrospray ionization (CF-EDESI), an ambient ionization source demonstrated previously for use with intact protein analysis, is expanded here for the coupling of reversed phase protein separations to mass spectrometry. This configuration allows the introduction of charging additives to enhance detection without affecting the chromatographic separation mechanism. Two demonstrations of the advantages of CF-EDESI are presented in this work. First, a proof-of- principle is presented to demonstrate the applicability of hyphenation of liquid chromatography (LC) to CF- EDESI. LC-CF-EDESI-MS has good sensitivity compared to LC-electrospray ionization (ESI)-mass spectrometry. Second, the supercharging mechanism investigated in CF-EDESI provides an insight into a highly debated supercharging process in ESI. The results indicate that the mechanism of protein charging seen in HPLC-CF-EDESI is different from supercharging phenomena in conventional ESI. The surface tension mechanism and binding mechanism may both contribute to protein supercharging in ESI.

Quantitative phosphoproteomics using acetone-based peptide labeling: method evaluation and application to a cardiac ischemia/reperfusion model

Mass spectrometry (MS) techniques to globally profile protein phosphorylation in cellular systems that are relevant to physiological or pathological changes have been of significant interest in biological research. An MS-based strategy utilizing an inexpensive acetone-based peptide-labeling technique known as reductive alkylation by acetone (RABA) for quantitative phosphoproteomics was explored to evaluate its capacity. Because the chemistry for RABA labeling for phosphorylation profiling had not been previously reported, it was first validated using a standard phosphoprotein and identical phosphoproteomes from cardiac tissue extracts. A workflow was then utilized to compare cardiac tissue phosphoproteomes from mouse hearts not expressing FGF2 versus hearts expressing low-molecular-weight fibroblast growth factor-2 (LMW FGF2) to relate low-molecular-weight fibroblast growth factor-2 (LMW FGF2)-mediated cardioprotective phenomena induced by ischemia/reperfusion injury of hearts, with downstream phosphorylation changes in LMW FGF2 signaling cascades. Statistically significant phosphorylation changes were identified at 14 different sites on 10 distinct proteins, including some with mechanisms already established for LMW FGF2-mediated cardioprotective signaling (e.g., connexin-43), some with new details linking LMW FGF2 to the cardioprotective mechanisms (e.g., cardiac myosin binding protein C or cMyBPC), and also several new downstream effectors not previously recognized for cardio-protective signaling by LMW FGF2. Additionally, one of the phosphopeptides, cMyBPC/pSer-282, identified was further verified with site-specific quantification using an SRM (selected reaction monitoring)-based approach that also relies on isotope labeling of a synthetic phosphopeptide with deuterated acetone as an internal standard. Overall, this study confirms that the inexpensive acetone-based peptide labeling can be used in both exploratory and targeted quantification phosphoproteomic studies to identify and verify biologically relevant phosphorylation changes in whole tissues.

Affinity mesh screen materials for selective extraction and analysis of antibiotics using transmission mode desorption electrospray ionization mass spectrometry

The extraction of active compounds from natural sources has shown to be an effective approach to drug discovery. However, the isolation and identification of natural products from complex extracts can be an arduous task. A novel approach to drug discovery is presented through the use of polymer screens functionalized with an l-lysine-d-alanine-d-alanine (Kaa) peptide to create new affinity capture mesh screen materials. The Kaa sequence is a well-characterized specific binding site for antibiotics that inhibit cell wall synthesis in Gram-positive bacteria. The detailed synthesis and characterization of these novel screen materials are presented in this work. Polypropylene mesh screens were first coated with a poly(acrylic acid) film by pulsed plasma polymerization. The synthesized Kaa peptide was then covalently attached to carboxylic acid groups through a condensation reaction. An analysis of captured compounds was performed in a rapid fashion with transmission-mode desorption electrospray ionization (TM-DESI) mass spectrometry. A proof of principle was demonstrated to show the ability of the novel affinity capture materials to select for a macrocyclic antibiotic, vancomycin, over a negative control compound, spectinomycin. With further development, this method may provide a rapid screening technique for new antibacterial compounds, for example, those extracted from natural product sources having a limited supply. Here, we show that the screen can capture vancomycin preferentially over spectinomycin in a spiked extract of tea leaves.

Antimony(III)-D, L-tartrates exhibit proton-assisted enantioselective binding in solution and in the gas phase

The negative ion mode ESI mass spectral analysis of antimony(III)-D- and -L-tartrate ("tartar emetic"), in association with leucine enantiomeric isotopomers, revealed remarkable proton-assisted enantioselective molecular recognition phenomena. The current study infers that recognition of amino acids by antimony(III)-D,L-tartrate complexes requires that the chiral selector associate a proton to become enantioselective. The dianionic selector itself failed to show enantiomeric discrimination capacity. This observation was shown to be consistent both in solution-phase targeting full scan and gas-phase targeting collision threshold dissociation (CTD) experiments. Importantly, this disparity in enantioselective binding capacity between the dianionic and the protonated monoanionic representatives of antimony(III)-D- and -L-tartrates could only be clearly revealed by ESI-MS and tandem mass spectrometry experiments as described herein. This finding urges a more in-depth study of mechanisms associated with exhibited enantiomeric resolving capacity of antimony tartrates in HPLC and CE applications, as well as in former ESI-MS association studies.

Dynamic solvation in phosphonium ionic liquids: comparison of bulk and micellar systems and considerations for the construction of the solvation correlation function, C(t)

Dynamic solvation of the dye coumarin 153 is studied in a phosphonium ionic liquid: hexadecyltributylphosphonium bromide, [(C4)3C16P+][Br-]. It forms micelles in water, and the bulk also exists as a liquid under our experimental conditions. This system permits a comparison with an imidazolium ionic liquid studied earlier, which also formed micelles in water (J. Phys. Chem. A 2006, 110, 10725-10730). We conclude that our analysis of the comparable situation in a phosphonium liquid is not as definitive as we had proposed earlier, i.e., that the majority of the early-time solvation arises from the organic cation. Part of the difficulty in performing this analysis is most likely due to the amount of water that is associated with the micelle. In the course of this work, we have focused on the calculation of the solvation correlation function, C(t), and investigated how it depends upon the methods with which the "zero-time" spectrum is constructed.