Proteomic profiling across breast cancer cell lines and models.

We performed quantitative proteomics on 60 human-derived breast cancer cell line models to a depth of ~13,000 proteins. The resulting high-throughput datasets were assessed for quality and reproducibility. We used the datasets to identify and characterize the subtypes of breast cancer and showed that they conform to known transcriptional subtypes, revealing that molecular subtypes are preserved even in under-sampled protein feature sets. All datasets are freely available as public resources on the LINCS portal. We anticipate that these datasets, either in isolation or in combination with complimentary measurements such as genomics, transcriptomics and phosphoproteomics, can be mined for the purpose of predicting drug response, informing cell line specific context in models of signalling pathways, and identifying markers of sensitivity or resistance to therapeutics.

Tandem mass tag-based quantitative proteomic profiling identifies candidate serum biomarkers of drug-induced liver injury in humans.

Diagnosis of drug-induced liver injury (DILI) and its distinction from other liver diseases are significant challenges in drug development and clinical practice. Here, we identify, confirm, and replicate the biomarker performance characteristics of candidate proteins in patients with DILI at onset (DO; n = 133) and follow-up (n = 120), acute non-DILI at onset (NDO; n = 63) and follow-up (n = 42), and healthy volunteers (HV; n = 104). Area under the receiver operating characteristic curve (AUC) for cytoplasmic aconitate hydratase, argininosuccinate synthase, carbamoylphosphate synthase, fumarylacetoacetase, fructose-1,6-bisphosphatase 1 (FBP1) across cohorts achieved near complete separation (range: 0.94-0.99) of DO and HV. In addition, we show that FBP1, alone or in combination with glutathione S-transferase A1 and leukocyte cell-derived chemotaxin 2, could potentially assist in clinical diagnosis by distinguishing NDO from DO (AUC range: 0.65-0.78), but further technical and clinical validation of these candidate biomarkers is needed.

A Comparison of Two Stability Proteomics Methods for Drug Target Identification in OnePot 2D Format.

Stability proteomics techniques that do not require drug modifications have emerged as an attractive alternative to affinity purification methods in drug target engagement studies. Two representative techniques include the chemical-denaturation-based SPROX (Stability of Proteins from Rates of Oxidation), which utilizes peptide-level quantification and thermal-denaturation-based TPP (Thermal Proteome Profiling), which utilizes protein-level quantification. Recently, the "OnePot" strategy was adapted for both SPROX and TPP to increase the throughput. When combined with the 2D setup which measures both the denaturation and the drug dose dimensions, the OnePot 2D format offers improved analysis specificity with higher resource efficiency. However, a systematic evaluation of the OnePot 2D format and a comparison between SPROX and TPP are still lacking. Here, we performed SPROX and TPP to identify protein targets of a well-studied pan-kinase inhibitor staurosporine with K562 lysate, in curve-fitting and OnePot 2D formats. We found that the OnePot 2D format provided ∼10× throughput, achieved ∼1.6× protein coverage and involves more straightforward data analysis. We also compared SPROX with the current "gold-standard" stability proteomics technique TPP in the OnePot 2D format. The protein coverage of TPP is ∼1.5 fold of SPROX; however, SPROX offers protein domain-level information, identifies comparable numbers of kinase hits, has higher signal (R value), and requires ∼3× less MS time. Unique SPROX hits encompass higher-molecular-weight proteins, compared to the unique TPP hits, and include atypical kinases. We also discuss hit stratification and prioritization strategies to promote the efficiency of hit followup.

Multiomics Profiling Establishes the Polypharmacology of FDA-Approved CDK4/6 Inhibitors and the Potential for Differential Clinical Activity.

The target profiles of many drugs are established early in their development and are not systematically revisited at the time of FDA approval. Thus, it is often unclear whether therapeutics with the same nominal targets but different chemical structures are functionally equivalent. In this paper we use five different phenotypic and biochemical assays to compare approved inhibitors of cyclin-dependent kinases 4/6-collectively regarded as breakthroughs in the treatment of hormone receptor-positive breast cancer. We find that transcriptional, proteomic, and phenotypic changes induced by palbociclib, ribociclib, and abemaciclib differ significantly; abemaciclib in particular has advantageous activities partially overlapping those of alvocidib, an older polyselective CDK inhibitor. In cells and mice, abemaciclib inhibits kinases other than CDK4/6 including CDK2/cyclin A/E-implicated in resistance to CDK4/6 inhibition-and CDK1/cyclin B. The multifaceted experimental and computational approaches described here therefore uncover underappreciated differences in CDK4/6 inhibitor activities with potential importance in treating human patients.

Multi omics analysis of fibrotic kidneys in two mouse models.

Kidney fibrosis represents an urgent unmet clinical need due to the lack of effective therapies and an inadequate understanding of the molecular pathogenesis. We have generated a comprehensive and combined multi-omics dataset (proteomics, mRNA and small RNA transcriptomics) of fibrotic kidneys that is searchable through a user-friendly web application: http://hbcreports.med.harvard.edu/fmm/ . Two commonly used mouse models were utilized: a reversible chemical-induced injury model (folic acid (FA) induced nephropathy) and an irreversible surgically-induced fibrosis model (unilateral ureteral obstruction (UUO)). mRNA and small RNA sequencing, as well as 10-plex tandem mass tag (TMT) proteomics were performed with kidney samples from different time points over the course of fibrosis development. The bioinformatics workflow used to process, technically validate, and combine the single omics data will be described. In summary, we present temporal multi-omics data from fibrotic mouse kidneys that are accessible through an interrogation tool (Mouse Kidney Fibromics browser) to provide a searchable transcriptome and proteome for kidney fibrosis researchers.

Dual Sensing of Physiologic pH and Calcium by EFCAB9 Regulates Sperm Motility.

Varying pH of luminal fluid along the female reproductive tract is a physiological cue that modulates sperm motility. CatSper is a sperm-specific, pH-sensitive calcium channel essential for hyperactivated motility and male fertility. Multi-subunit CatSper channel complexes organize linear Ca2+ signaling nanodomains along the sperm tail. Here, we identify EF-hand calcium-binding domain-containing protein 9 (EFCAB9) as a bifunctional, cytoplasmic machine modulating the channel activity and the domain organization of CatSper. Knockout mice studies demonstrate that EFCAB9, in complex with the CatSper subunit, CATSPERζ, is essential for pH-dependent and Ca2+-sensitive activation of the CatSper channel. In the absence of EFCAB9, sperm motility and fertility is compromised, and the linear arrangement of the Ca2+ signaling domains is disrupted. EFCAB9 interacts directly with CATSPERζ in a Ca2+-dependent manner and dissociates at elevated pH. These observations suggest that EFCAB9 is a long-sought, intracellular, pH-dependent Ca2+ sensor that triggers changes in sperm motility.

PF-06651600, a Dual JAK3/TEC Family Kinase Inhibitor.

PF-06651600 was developed as an irreversible inhibitor of JAK3 with selectivity over the other three JAK isoforms. A high level of selectivity toward JAK3 is achieved by the covalent interaction of PF-06651600 with a unique cysteine residue (Cys-909) in the catalytic domain of JAK3, which is replaced by a serine residue in the other JAK isoforms. Importantly, 10 other kinases in the kinome have a cysteine at the equivalent position of Cys-909 in JAK3. Five of those kinases belong to the TEC kinase family including BTK, BMX, ITK, RLK, and TEC and are also inhibited by PF-06651600. Preclinical data demonstrate that inhibition of the cytolytic function of CD8+ T cells and NK cells by PF-06651600 is driven by the inhibition of TEC kinases. On the basis of the underlying pathophysiology of inflammatory diseases such as rheumatoid arthritis, inflammatory bowel disease, alopecia areata, and vitiligo, the dual activity of PF-06651600 toward JAK3 and the TEC kinase family may provide a beneficial inhibitory profile for therapeutic intervention.

Adaptation of Human iPSC-Derived Cardiomyocytes to Tyrosine Kinase Inhibitors Reduces Acute Cardiotoxicity via Metabolic Reprogramming.

Tyrosine kinase inhibitors (TKIs) are widely used to treat solid tumors but can be cardiotoxic. The molecular basis for this toxicity and its relationship to therapeutic mechanisms remain unclear; we therefore undertook a systems-level analysis of human cardiomyocytes (CMs) exposed to four TKIs. CMs differentiated from human induced pluripotent stem cells (hiPSCs) were exposed to sunitinib, sorafenib, lapatinib, or erlotinib, and responses were assessed by functional assays, microscopy, RNA sequencing, and mass spectrometry (GEO: GSE114686; PRIDE: PXD012043). TKIs have diverse effects on hiPSC-CMs distinct from inhibition of tyrosine-kinase-mediated signal transduction; cardiac metabolism is particularly sensitive. Following sorafenib treatment, oxidative phosphorylation is downregulated, resulting in a profound defect in mitochondrial energetics. Cells adapt by upregulating aerobic glycolysis. Adaptation makes cells less acutely sensitive to sorafenib but may have long-term negative consequences. Thus, CMs exhibit adaptive responses to anti-cancer drugs conceptually similar to those previously shown in tumors to mediate drug resistance.

Leveraging Compound Promiscuity to Identify Targetable Cysteines within the Kinome.

Covalent kinase inhibitors, which typically target cysteine residues, represent an important class of clinically relevant compounds. Approximately 215 kinases are known to have potentially targetable cysteines distributed across 18 spatially distinct locations proximal to the ATP-binding pocket. However, only 40 kinases have been covalently targeted, with certain cysteine sites being the primary focus. To address this disparity, we have developed a strategy that combines the use of a multi-targeted acrylamide-modified inhibitor, SM1-71, with a suite of complementary chemoproteomic and cellular approaches to identify additional targetable cysteines. Using this single multi-targeted compound, we successfully identified 23 kinases that are amenable to covalent inhibition including MKNK2, MAP2K1/2/3/4/6/7, GAK, AAK1, BMP2K, MAP3K7, MAPKAPK5, GSK3A/B, MAPK1/3, SRC, YES1, FGFR1, ZAK (MLTK), MAP3K1, LIMK1, and RSK2. The identification of nine of these kinases previously not targeted by a covalent inhibitor increases the number of targetable kinases and highlights opportunities for covalent kinase inhibitor development.

The NAD+ Salvage Pathway Supports PHGDH-Driven Serine Biosynthesis.

NAD+ is a key metabolic redox cofactor that is regenerated from nicotinamide through the NAD+ salvage pathway. Here, we find that inhibiting the NAD+ salvage pathway depletes serine biosynthesis from glucose by impeding the NAD+-dependent protein, 3-phosphoglycerate dehydrogenase (PHGDH). Importantly, we find that PHGDHhigh breast cancer cell lines are exquisitely sensitive to inhibition of the NAD+ salvage pathway. Further, we find that PHGDH protein levels and those of the rate-limiting enzyme of NAD+ salvage, NAMPT, correlate in ER-negative, basal-like breast cancers. Although NAD+ salvage pathway inhibitors are actively being pursued in cancer treatment, their efficacy has been poor, and our findings suggest that they may be effective for PHGDH-dependent cancers.

Role of Selenof as a Gatekeeper of Secreted Disulfide-Rich Glycoproteins.

Selenof (15-kDa selenoprotein; Sep15) is an endoplasmic reticulum (ER)-resident thioredoxin-like oxidoreductase that occurs in a complex with UDP-glucose:glycoprotein glucosyltransferase. We found that Selenof deficiency in mice leads to elevated levels of non-functional circulating plasma immunoglobulins and increased secretion of IgM during in vitro splenic B cell differentiation. However, Selenof knockout animals show neither enhanced bacterial killing capacity nor antigen-induced systemic IgM activity, suggesting that excess immunoglobulins are not functional. In addition, ER-to-Golgi transport of a target glycoprotein was delayed in Selenof knockout embryonic fibroblasts, and proteomic analyses revealed that Selenof deficiency is primarily associated with antigen presentation and ER-to-Golgi transport. Together, the data suggest that Selenof functions as a gatekeeper of immunoglobulins and, likely, other client proteins that exit the ER, thereby supporting redox quality control of these proteins.

MS3-IDQ: Utilizing MS3 Spectra beyond Quantification Yields Increased Coverage of the Phosphoproteome in Isobaric Tag Experiments.

Protein phosphorylation is critically important for many cellular processes, including progression through the cell cycle, cellular metabolism, and differentiation. Isobaric labeling, for example, tandem mass tags (TMT), in phosphoproteomics workflows enables both relative and absolute quantitation of these phosphorylation events. Traditional TMT workflows identify peptides using fragment ions at the MS2 level and quantify reporter ions at the MS3 level. However, in addition to the TMT reporter ions, MS3 spectra also include fragment ions that can be used to identify peptides. Here we describe using MS3 spectra for both phosphopeptide identification and quantification, a process that we term MS3-IDQ. To maximize quantified phosphopeptides, we optimize several instrument parameters, including the modality of mass analyzer (i.e., ion trap or Orbitrap), MS2 automatic gain control (AGC), and MS3 normalized collision energy (NCE), to achieve the best balance of identified and quantified peptides. Our optimized MS3-IDQ method included the following parameters for the MS3 scan: NCE = 37.5 and AGC target = 1.5 × 105, and scan range = 100-2000. Data from the MS3 scan were complementary to those of the MS2 scan, and the combination of these scans can increase phosphoproteome coverage by >50%, thereby yielding a greater number of quantified and accurately localized phosphopeptides.

Aminoglycoside-driven biosynthesis of selenium-deficient Selenoprotein P.

Selenoprotein biosynthesis relies on the co-translational insertion of selenocysteine in response to UGA codons. Aminoglycoside antibiotics interfere with ribosomal function and may cause codon misreading. We hypothesized that biosynthesis of the selenium (Se) transporter selenoprotein P (SELENOP) is particularly sensitive to antibiotics due to its ten in frame UGA codons. As liver regulates Se metabolism, we tested the aminoglycosides G418 and gentamicin in hepatoma cell lines (HepG2, Hep3B and Hepa1-6) and in experimental mice. In vitro, SELENOP levels increased strongly in response to G418, whereas expression of the glutathione peroxidases GPX1 and GPX2 was marginally affected. Se content of G418-induced SELENOP was dependent on Se availability, and was completely suppressed by G418 under Se-poor conditions. Selenocysteine residues were replaced mainly by cysteine, tryptophan and arginine in a codon-specific manner. Interestingly, in young healthy mice, antibiotic treatment failed to affect Selenop biosynthesis to a detectable degree. These findings suggest that the interfering activity of aminoglycosides on selenoprotein biosynthesis can be severe, but depend on the Se status, and other parameters likely including age and general health. Focused analyses with aminoglycoside-treated patients are needed next to evaluate a possible interference of selenoprotein biosynthesis by the antibiotics and elucidate potential side effects.

A Strategy to Combine Sample Multiplexing with Targeted Proteomics Assays for High-Throughput Protein Signature Characterization.

Targeted mass spectrometry assays for protein quantitation monitor peptide surrogates, which are easily multiplexed to target many peptides in a single assay. However, these assays have generally not taken advantage of sample multiplexing, which allows up to ten analyses to occur in parallel. We present a two-dimensional multiplexing workflow that utilizes synthetic peptides for each protein to prompt the simultaneous quantification of >100 peptides from up to ten mixed sample conditions. We demonstrate that targeted analysis of unfractionated lysates (2 hr) accurately reproduces the quantification of fractionated lysates (72 hr analysis) while obviating the need for peptide detection prior to quantification. We targeted 131 peptides corresponding to 69 proteins across all 60 National Cancer Institute cell lines in biological triplicate, analyzing 180 samples in only 48 hr (the equivalent of 16 min/sample). These data further elucidated a correlation between the expression of key proteins and their cellular response to drug treatment.

Neutral Loss Is a Very Common Occurrence in Phosphotyrosine-Containing Peptides Labeled with Isobaric Tags.

While developing a multiplexed phosphotyrosine peptide quantification assay, an unexpected observation was made: significant neutral loss from phosphotyrosine (pY) containing peptides. Using a 2000-member peptide library, we sought to systematically investigate this observation by comparing unlabeled peptides with the two highest-plex isobaric tags (iTRAQ8 and TMT10) across CID, HCD, and ETD fragmentation using high resolution high mass accuracy Orbitrap instrumentation. We found pY peptide neutral loss behavior was consistent with reduced proton mobility, and does not occur during ETD. The site of protonation at the peptide N-terminus changes from a primary to a tertiary amine as a result of TMT labeling which would increase the gas phase basicity and reduce proton mobility at this site. This change in fragmentation behavior has implications during instrument method development and interpretation of MS/MS spectra, and therefore ensuing follow-up studies. We show how sites not localized to tyrosine by search and site localization algorithms can be confidently reassigned to tyrosine using neutral loss and phosphotyrosine immonium ions. We believe these findings will be of general interest to those studying pY signal transduction using isobaric tags.

Quantitative mass spectrometry-based multiplexing compares the abundance of 5000 S. cerevisiae proteins across 10 carbon sources.

UNLABELLED: The budding yeast Saccharomyces cerevisiae is a model system for investigating biological processes. Cellular events are known to be dysregulated due to shifts in carbon sources. However, the comprehensive proteomic alterations thereof have not been fully investigated. Here we examined proteomic alterations in S. cerevisiae due to the adaptation of yeast from glucose to nine different carbon sources - maltose, trehalose, fructose, sucrose, glycerol, acetate, pyruvate, lactic acid, and oleate. Isobaric tag-based mass spectrometry techniques are at the forefront of global proteomic investigations. As such, we used a TMT10-plex strategy to study multiple growth conditions in a single experiment. The SPS-MS3 method on an Orbitrap Fusion Lumos mass spectrometer enabled the quantification of over 5000 yeast proteins across ten carbon sources at a 1% protein-level FDR. On average, the proteomes of yeast cultured in fructose and sucrose deviated the least from those cultured in glucose. As expected, gene ontology classification revealed the major alteration in protein abundances occurred in metabolic pathways and mitochondrial proteins. Our protocol lays the groundwork for further investigation of carbon source-induced protein alterations. Additionally, these data offer a hypothesis-generating resource for future studies aiming to investigate both characterized and uncharacterized genes. BIOLOGICAL SIGNIFICANCE: We investigate the proteomic alterations in S. cerevisiae resulting from adaptation of yeast from glucose to nine different carbon sources - maltose, trehalose, fructose, sucrose, glycerol, acetate, pyruvate, lactic acid, and oleate. SPS-MS3 TMT10plex analysis is used for quantitative proteomic analysis. We showcase a technique that allows the quantification of over 5000 yeast proteins, the highest number to date in S. cerevisiae, across 10 growth conditions in a single experiment. As expected, gene ontology classification of proteins with the major alterations in abundances occurred in metabolic pathways and mitochondrial proteins, reflecting the degree of metabolic stress when cellular machinery shifts from growth on glucose to an alternative carbon source. Our protocol lays the groundwork for further investigation of carbon source-induced protein alterations. Improving depth of coverage - measuring abundance changes of over 5000 proteins - increases our understanding of difficult-to-study genes in the model system S. cerevisiae and by homology human cell biology. We submit this highly comprehensive dataset as a hypothesis generating resource for targeted studies on uncharacterized genes.

Regulation of Selenocysteine Content of Human Selenoprotein P by Dietary Selenium and Insertion of Cysteine in Place of Selenocysteine.

Selenoproteins are a unique group of proteins that contain selenium in the form of selenocysteine (Sec) co-translationally inserted in response to a UGA codon with the help of cis- and trans-acting factors. Mammalian selenoproteins contain single Sec residues, with the exception of selenoprotein P (SelP) that has 7-15 Sec residues depending on species. Assessing an individual's selenium status is important under various pathological conditions, which requires a reliable selenium biomarker. Due to a key role in organismal selenium homeostasis, high Sec content, regulation by dietary selenium, and availability of robust assays in human plasma, SelP has emerged as a major biomarker of selenium status. Here, we found that Cys is present in various Sec positions in human SelP. Treatment of cells expressing SelP with thiophosphate, an analog of the selenium donor for Sec synthesis, led to a nearly complete replacement of Sec with Cys, whereas supplementation of cells with selenium supported Sec insertion. SelP isolated directly from human plasma had up to 8% Cys inserted in place of Sec, depending on the Sec position. These findings suggest that a change in selenium status may be reflected in both SelP concentration and its Sec content, and that availability of the SelP-derived selenium for selenoprotein synthesis may be overestimated under conditions of low selenium status due to replacement of Sec with Cys.

Quantitative Proteomic Atlas of Ubiquitination and Acetylation in the DNA Damage Response.

Execution of the DNA damage response (DDR) relies upon a dynamic array of protein modifications. Using quantitative proteomics, we have globally profiled ubiquitination, acetylation, and phosphorylation in response to UV and ionizing radiation. To improve acetylation site profiling, we developed the strategy FACET-IP. Our datasets of 33,500 ubiquitination and 16,740 acetylation sites provide valuable insight into DDR remodeling of the proteome. We find that K6- and K33-linked polyubiquitination undergo bulk increases in response to DNA damage, raising the possibility that these linkages are largely dedicated to DDR function. We also show that Cullin-RING ligases mediate 10% of DNA damage-induced ubiquitination events and that EXO1 is an SCF-Cyclin F substrate in the response to UV radiation. Our extensive datasets uncover additional regulated sites on known DDR players such as PCNA and identify previously unknown DDR targets such as CENPs, underscoring the broad impact of the DDR on cellular physiology.

Comprehensive Temporal Protein Dynamics during the Diauxic Shift in Saccharomyces cerevisiae.

Yeast (Saccharomyces cerevisiae) has served as a key model system in biology and as a benchmark for "omics" technology. Although near-complete proteomes of log phase yeast have been measured, protein abundance in yeast is dynamic, particularly during the transition from log to stationary phase. Defining the dynamics of proteomic changes during this transition, termed the diauxic shift, is important to understand the basic biology of proliferative versus quiescent cells. Here, we perform temporal quantitative proteomics to fully capture protein induction and repression during the diauxic shift. Accurate and sensitive quantitation at a high temporal resolution and depth of proteome coverage was achieved using TMT10 reagents and LC-MS3 analysis on an Orbitrap Fusion tribrid mass spectrometer deploying synchronous precursor selection. Triplicate experiments were analyzed using the time-course R package and a simple template matching strategy was used to reveal groups of proteins with similar temporal patterns of protein induction and repression. Within these groups are functionally distinct types of proteins such as those of glyoxylate metabolism and many proteins of unknown function not previously associated with the diauxic shift (e.g. YNR034W-A and FMP16). We also perform a dual time-course experiment to determine Hap2-dependent proteins during the diauxic shift. These data serve as an important basic model for fermentative versus respiratory growth of yeast and other eukaryotes and are a benchmark for temporal quantitative proteomics.

Mechanism-based proteomic screening identifies targets of thioredoxin-like proteins.

Thioredoxin (Trx)-fold proteins are protagonists of numerous cellular pathways that are subject to thiol-based redox control. The best characterized regulator of thiols in proteins is Trx1 itself, which together with thioredoxin reductase 1 (TR1) and peroxiredoxins (Prxs) comprises a key redox regulatory system in mammalian cells. However, there are numerous other Trx-like proteins, whose functions and redox interactors are unknown. It is also unclear if the principles of Trx1-based redox control apply to these proteins. Here, we employed a proteomic strategy to four Trx-like proteins containing CXXC motifs, namely Trx1, Rdx12, Trx-like protein 1 (Txnl1) and nucleoredoxin 1 (Nrx1), whose cellular targets were trapped in vivo using mutant Trx-like proteins, under conditions of low endogenous expression of these proteins. Prxs were detected as key redox targets of Trx1, but this approach also supported the detection of TR1, which is the Trx1 reductant, as well as mitochondrial intermembrane proteins AIF and Mia40. In addition, glutathione peroxidase 4 was found to be a Rdx12 redox target. In contrast, no redox targets of Txnl1 and Nrx1 could be detected, suggesting that their CXXC motifs do not engage in mixed disulfides with cellular proteins. For some Trx-like proteins, the method allowed distinguishing redox and non-redox interactions. Parallel, comparative analyses of multiple thiol oxidoreductases revealed differences in the functions of their CXXC motifs, providing important insights into thiol-based redox control of cellular processes.