Improved determination of protein turnover rate with heavy water labeling by mass isotopomer ratio selection.

The synthesis and degradation rates of proteins form an essential component of gene expression control. Heavy water labeling has been used in conjunction with mass spectrometry to measure protein turnover rates, but the optimal analytical approaches to derive turnover rates from the isotopomer patterns of deuterium labeled peptides continue to be a subject of research. Here we describe a method, which comprises a reverse lookup of numerically approximated peptide isotope envelopes, coupled to the selection of optimal isotopomer pairs based on peptide sequence, to calculate the molar fraction of new peptide synthesis in heavy water labeling mass spectrometry experiments. We validated this approach using an experimental calibration curve comprising mixtures of fully unlabeled and fully labeled proteomes. We then re-analyzed 17 proteome-wide turnover experiments from four mouse organs, and showed that the method increases the coverage of well-fitted peptides in protein turnover experiments by 25-82%. The method is implemented in the Riana software tool for protein turnover analysis, and may avail ongoing efforts to study the synthesis and degradation kinetics of proteins in animals on a proteome-wide scale. What's new: We describe a reverse lookup method to calculate the molar fraction of new synthesis from numerically approximated peptide isotopomer profiles in heavy water labeling mass spectrometry experiments. Using an experimental calibration curve comprising mixtures of fully unlabeled and fully labeled proteomes at various proportions, we show that this method provides a straightforward way to calculate the proportion of new proteins in a protein pool from arbitrarily chosen isotopomer ratios. We next analyzed which of the isotopomer pairs within the peptide isotope envelope yielded isotopomer time courses that fit most closely to kinetic models, and found that the identity of the isotopomer pair depends partially on the number of deuterium accessible labeling sites of the peptide. We next derived a strategy to automatically select the isotopomer pairs to calculate turnover rates based on peptide sequence, and showed that this increases the coverage of existing proteome-wide turnover experiments in multiple data sets of the mouse heart, liver, kidney, and skeletal muscle by up to 25-82%.

Simultaneous proteome localization and turnover analysis reveals spatiotemporal features of protein homeostasis disruptions.

The spatial and temporal distributions of proteins are critical to protein function, but cannot be directly assessed by measuring protein bundance. Here we describe a mass spectrometry-based proteomics strategy, Simultaneous Proteome Localization and Turnover (SPLAT), to measure concurrently protein turnover rates and subcellular localization in the same experiment. Applying the method, we find that unfolded protein response (UPR) has different effects on protein turnover dependent on their subcellular location in human AC16 cells, with proteome-wide slowdown but acceleration among stress response proteins in the ER and Golgi. In parallel, UPR triggers broad differential localization of proteins including RNA-binding proteins and amino acid transporters. Moreover, we observe newly synthesized proteins including EGFR that show a differential localization under stress than the existing protein pools, reminiscent of protein trafficking disruptions. We next applied SPLAT to an induced pluripotent stem cell derived cardiomyocyte (iPSC-CM) model of cancer drug cardiotoxicity upon treatment with the proteasome inhibitor carfilzomib. Paradoxically, carfilzomib has little effect on global average protein half-life, but may instead selectively disrupt sarcomere protein homeostasis. This study provides a view into the interactions of protein spatial and temporal dynamics and demonstrates a method to examine protein homeostasis regulations in stress and drug response.

Simultaneous proteome localization and turnover analysis reveals spatiotemporal features of protein homeostasis disruptions.

The functions of proteins depend on their spatial and temporal distributions, which are not directly measured by static protein abundance. Under endoplasmic reticulum (ER) stress, the unfolded protein response (UPR) pathway remediates proteostasis in part by altering the turnover kinetics and spatial distribution of proteins. A global view of these spatiotemporal changes has yet to emerge and it is unknown how they affect different cellular compartments and pathways. Here we describe a mass spectrometry-based proteomics strategy and data analysis pipeline, termed Simultaneous Proteome Localization and Turnover (SPLAT), to measure concurrently the changes in protein turnover and subcellular distribution in the same experiment. Investigating two common UPR models of thapsigargin and tunicamycin challenge in human AC16 cells, we find that the changes in protein turnover kinetics during UPR varies across subcellular localizations, with overall slowdown but an acceleration in endoplasmic reticulum and Golgi proteins involved in stress response. In parallel, the spatial proteomics component of the experiment revealed an externalization of amino acid transporters and ion channels under UPR, as well as the migration of RNA-binding proteins toward an endosome co-sedimenting compartment. The SPLAT experimental design classifies heavy and light SILAC labeled proteins separately, allowing the observation of differential localization of new and old protein pools and capturing a partition of newly synthesized EGFR and ITGAV to the ER under stress that suggests protein trafficking disruptions. Finally, application of SPLAT toward human induced pluripotent stem cell derived cardiomyocytes (iPSC-CM) exposed to the cancer drug carfilzomib, identified a selective disruption of proteostasis in sarcomeric proteins as a potential mechanism of carfilzomib-mediated cardiotoxicity. Taken together, this study provides a global view into the spatiotemporal dynamics of human cardiac cells and demonstrates a method for inferring the coordinations between spatial and temporal proteome regulations in stress and drug response.

Co-culture with normal and senescent AC16 cardiomyocytes modulates fibrotic response in primary human ventricular fibroblasts.

Transwell co-culture with human AC16 cardiomyocyte-like cells modifies the response of primary human ventricular fibroblasts to TGF-ß stimulation. Fibrotic response markers including collagen I (COL1A1) and ɑ-smooth muscle actin (ACTA2) are amplified in the presence of AC16 cells, whereas others including periostin (POSTN) and fibronectin (FN1) are suppressed. Similar modulation is observed when the ventricular fibroblasts are co-cultured with AC16 cells under baseline and induced senescence conditions. Given that the response to TGF-ß stimulation is commonly measured to study fibrotic signaling and drug treatments in vitro, the results here suggest that the effect of cellular crosstalk should be more broadly considered.

Protein prediction models support widespread post-transcriptional regulation of protein abundance by interacting partners.

Protein and mRNA levels correlate only moderately. The availability of proteogenomics data sets with protein and transcript measurements from matching samples is providing new opportunities to assess the degree to which protein levels in a system can be predicted from mRNA information. Here we examined the contributions of input features in protein abundance prediction models. Using large proteogenomics data from 8 cancer types within the Clinical Proteomic Tumor Analysis Consortium (CPTAC) data set, we trained models to predict the abundance of over 13,000 proteins using matching transcriptome data from up to 958 tumor or normal adjacent tissue samples each, and compared predictive performances across algorithms, data set sizes, and input features. Over one-third of proteins (4,648) showed relatively poor predictability (elastic net r ≤ 0.3) from their cognate transcripts. Moreover, we found widespread occurrences where the abundance of a protein is considerably less well explained by its own cognate transcript level than that of one or more trans locus transcripts. The incorporation of additional trans-locus transcript abundance data as input features increasingly improved the ability to predict sample protein abundance. Transcripts that contribute to non-cognate protein abundance primarily involve those encoding known or predicted interaction partners of the protein of interest, including not only large multi-protein complexes as previously shown, but also small stable complexes in the proteome with only one or few stable interacting partners. Network analysis further shows a complex proteome-wide interdependency of protein abundance on the transcript levels of multiple interacting partners. The predictive model analysis here therefore supports that protein-protein interaction including in small protein complexes exert post-transcriptional influence on proteome compositions more broadly than previously recognized. Moreover, the results suggest mRNA and protein co-expression analysis may have utility for finding gene interactions and predicting expression changes in biological systems.