Mikhail M. Savitski

Mass-spectrometry-based draft of the human proteome

Proteomes are characterized by large protein-abundance differences, cell-type- and time-dependent expression patterns and post-translational modifications, all of which carry biological information that is not accessible by genomics or transcriptomics. Here we present a mass-spectrometry-based draft of the human proteome and a public, high-performance, in-memory database for real-time analysis of terabytes of big data, called ProteomicsDB. The information assembled from human tissues, cell lines and body fluids enabled estimation of the size of the protein-coding genome, and identified organ-specific proteins and a large number of translated lincRNAs (long intergenic non-coding RNAs). Analysis of messenger RNA and protein-expression profiles of human tissues revealed conserved control of protein abundance, and integration of drug-sensitivity data enabled the identification of proteins predicting resistance or sensitivity. The proteome profiles also hold considerable promise for analysing the composition and stoichiometry of protein complexes. ProteomicsDB thus enables navigation of proteomes, provides biological insight and fosters the development of proteomic technology.

Inhibition of BET recruitment to chromatin as an effective treatment for MLL-fusion leukaemia

Recurrent chromosomal translocations involving the mixed lineage leukaemia (MLL) gene initiate aggressive forms of leukaemia, which are often refractory to conventional therapies. Many MLL-fusion partners are members of the super elongation complex (SEC), a critical regulator of transcriptional elongation, suggesting that aberrant control of this process has an important role in leukaemia induction. Here we use a global proteomic strategy to demonstrate that MLL fusions, as part of SEC and the polymerase-associated factor complex (PAFc), are associated with the BET family of acetyl-lysine recognizing, chromatin ‘adaptor’ proteins. These data provided the basis for therapeutic intervention in MLL-fusion leukaemia, via the displacement of the BET family of proteins from chromatin. We show that a novel small molecule inhibitor of the BET family, GSK1210151A (I-BET151), has profound efficacy against human and murine MLL-fusion leukaemic cell lines, through the induction of early cell cycle arrest and apoptosis. I-BET151 treatment in two human leukaemia cell lines with different MLL fusions alters the expression of a common set of genes whose function may account for these phenotypic changes. The mode of action of I-BET151 is, at least in part, due to the inhibition of transcription at key genes (BCL2, C-MYC and CDK6) through the displacement of BRD3/4, PAFc and SEC components from chromatin. In vivo studies indicate that I-BET151 has significant therapeutic value, providing survival benefit in two distinct mouse models of murine MLL–AF9 and human MLL–AF4 leukaemia. Finally, the efficacy of I-BET151 against human leukaemia stem cells is demonstrated, providing further evidence of its potent therapeutic potential. These findings establish the displacement of BET proteins from chromatin as a promising epigenetic therapy for these aggressive leukaemias.

Chemoproteomics profiling of HDAC inhibitors reveals selective targeting of HDAC complexes

The development of selective histone deacetylase (HDAC) inhibitors with anti-cancer and anti-inflammatory properties remains challenging in large part owing to the difficulty of probing the interaction of small molecules with megadalton protein complexes. A combination of affinity capture and quantitative mass spectrometry revealed the selectivity with which 16 HDAC inhibitors target multiple HDAC complexes scaffolded by ELM-SANT domain subunits, including a novel mitotic deacetylase complex (MiDAC). Inhibitors clustered according to their target profiles with stronger binding of aminobenzamides to the HDAC NCoR complex than to the HDAC Sin3 complex. We identified several non-HDAC targets for hydroxamate inhibitors. HDAC inhibitors with distinct profiles have correspondingly different effects on downstream targets. We also identified the anti-inflammatory drug bufexamac as a class IIb (HDAC6, HDAC10) HDAC inhibitor. Our approach enables the discovery of novel targets and inhibitors and suggests that the selectivity of HDAC inhibitors should be evaluated in the context of HDAC complexes and not purified catalytic subunits.

Confident Phosphorylation Site Localization Using the Mascot Delta Score

Large scale phosphorylation analysis is more and more getting into focus of proteomic research. Although it is now possible to identify thousands of phosphorylated peptides in a biological system, confident site localization remains challenging. Here we validate the Mascot Delta Score (MD-score) as a simple method that achieves similar sensitivity and specificity for phosphosite localization as the published Ascore, which is mainly used in conjunction with Sequest. The MD-score was evaluated using liquid chromatography-tandem MS data of 180 individually synthesized phosphopeptides with precisely known phosphorylation sites. We tested the MD-score for a wide range of commonly available fragmentation methods and found it to be applicable throughout with high statistical significance. However, the different fragmentation techniques differ strongly in their ability to localize phosphorylation sites. At 1% false localization rate, the highest number of correctly assigned phosphopeptides was achieved by higher energy collision induced dissociation in combination with an Orbitrap mass analyzer followed very closely by low resolution ion trap spectra obtained after electron transfer dissociation. Both these methods are significantly better than low resolution spectra acquired after collision induced dissociation and multi stage activation. Score thresholds determined from simple calibration functions for each fragmentation method were stable over replicate analyses of the phosphopeptide set. The MD-score outperforms the Ascore for tyrosine phosphorylated peptides and we further show that the ability to call sites correctly increases with increasing distance of two candidate sites within a peptide sequence. The MD-score does not require complex computational steps which makes it attractive in terms of practical utility. We provide all mass spectra and the synthetic peptides to the community so that the development of present and future localization software can be benchmarked and any laboratory can determine MD-scores and localization probabilities for their individual analytical set up.

Tracking cancer drugs in living cells by thermal profiling of the proteome

INTRODUCTION Understanding drug mechanism poses the daunting challenge of determining the affinity of the drug for all potential targets. Drug target engagement can be assessed by means of a cellular thermal shift assay (CETSA) based on ligand-induced changes in protein thermal stability. We combined the CETSA method with quantitative mass spectrometry to study the effect of drugs on the thermal profile of a cellular proteome comprising more than 7000 proteins. The approach enabled the monitoring of drug targets and downstream effectors. Tracking drugs in living cells. Drugs alter the thermal stability of proteins directly through compound binding or indirectly through changes in overall protein state. Thermal proteome profiling determines melting curves for thousands of proteins and tracks drug action in cells. RATIONALE We devised a method for the thermal profiling of cellular proteomes. Cells were cultured with or without drugs and heated to different temperatures so as to induce protein denaturation, and remaining soluble proteins were extracted with buffer. At each temperature, soluble proteins were quantified by means of high-resolution mass spectrometry, yielding denaturation curves. This allowed determination of thermal stability and the identification of ligand-induced shifts. To rank binding affinities among multiple targets, we determined stability profiles across a range of compound concentrations at a defined temperature. Comparison of the thermal profiles obtained after drug treatment of intact cells versus cell extract allowed us to distinguish effects induced by ligand binding from those induced by downstream modifications. RESULTS We performed thermal proteome profiling (TPP) on human K562 cells by heating intact cells or cell extracts and observed marked differences in melting properties between the two settings, with a trend toward increased protein stability in cell extract. Adenosine triphosphatase (ATP)–binding proteins showed a significant trend toward increased stability in intact cells, suggesting stabilization by the endogenous ligand. This was confirmed with the addition of ATP to cell extract, which resulted in increased stability for this protein group. The ability of TPP to identify target binding was validated by using the broad-specificity inhibitors staurosporine and GSK3182571, which induced shifts in the melting temperatures of many kinase targets and also affected the thermal profiles of other proteins, including regulatory subunits of kinase complexes. We identified the heme biosynthesis enzyme ferrochelatase (FECH) as an off-target of several kinase inhibitors and showed that the drug vemurafenib reaches full target occupancy of its cognate target BRAF and the off-target FECH within a narrow concentration window. FECH deficiency is genetically linked to protoporphyria, suggesting that the photosensitivity induced by vemurafenib and other drugs is mediated by FECH. Drug treatment of live cells affected not only direct target proteins but also downstream effectors. The ABL inhibitor dasatinib induced thermal shifts in several proteins downstream of BCR-ABL, including CRKL, and at concentrations in good agreement with the effect on cell growth. CONCLUSION Thermal profiling of cellular proteomes enables the differential assessment of protein ligand binding and other protein modifications, providing an unbiased measure of drug-target occupancy for multiple targets and facilitating the identification of markers for drug efficacy and toxicity. Mapping human drug targets in the cell To understand both the beneficial and the side effects of a drug, one would need to know its full binding profile to all cellular proteins. Savitski et al. take significant steps toward meeting this daunting challenge. They monitored the unfolding or “melting” of over 7000 human proteins and measured how small-molecule binding changes individual melting profiles. As a proof of principle, over 50 targets were identified for an inhibitor known to bind a broad spectrum of kinases. Two cancer drugs, vemurafib and Alectinib, are known to have a side effect of photosensitivity. The thermal profiling approach identified drug-protein interactions responsible for these side effects. Science, this issue 10.1126/science.1255784 Monitoring drug effects on the thermal profile of a cell’s proteins identifies drug targets and off-targets. The thermal stability of proteins can be used to assess ligand binding in living cells. We have generalized this concept by determining the thermal profiles of more than 7000 proteins in human cells by means of mass spectrometry. Monitoring the effects of small-molecule ligands on the profiles delineated more than 50 targets for the kinase inhibitor staurosporine. We identified the heme biosynthesis enzyme ferrochelatase as a target of kinase inhibitors and suggest that its inhibition causes the phototoxicity observed with vemurafenib and alectinib. Thermal shifts were also observed for downstream effectors of drug treatment. In live cells, dasatinib induced shifts in BCR-ABL pathway proteins, including CRK/CRKL. Thermal proteome profiling provides an unbiased measure of drug-target engagement and facilitates identification of markers for drug efficacy and toxicity.

ModifiComb, a New Proteomic Tool for Mapping Substoichiometric Post-translational Modifications, Finding Novel Types of Modifications, and Fingerprinting Complex Protein Mixtures

A major challenge in proteomics is to fully identify and characterize the post-translational modification (PTM) patterns present at any given time in cells, tissues, and organisms. Here we present a fast and reliable method (“ModifiComb”) for mapping hundreds types of PTMs at a time, including novel and unexpected PTMs. The high mass accuracy of Fourier transform mass spectrometry provides in many cases unique elemental composition of the PTM through the difference ΔM between the molecular masses of the modified and unmodified peptides, whereas the retention time difference ΔRT between their elution in reversed-phase liquid chromatography provides an additional dimension for PTM identification. Abundant sequence information obtained with complementary fragmentation techniques using ion-neutral collisions and electron capture often locates the modification to a single residue. The (ΔM, ΔRT) maps are representative of the proteome and its overall modification state and may be used for database-independent organism identification, comparative proteomic studies, and biomarker discovery. Examples of newly found modifications include +12.000 Da (+C atom) incorporation into proline residues of peptides from proline-rich proteins found in human saliva. This modification is hypothesized to increase the known activity of the peptide.

Improving Protein Identification Using Complementary Fragmentation Techniques in Fourier Transform Mass Spectrometry

Identification of proteins by MS/MS is performed by matching experimental mass spectra against calculated spectra of all possible peptides in a protein data base. The search engine assigns each spectrum a score indicating how well the experimental data complies with the expected one; a higher score means increased confidence in the identification. One problem is the false-positive identifications, which arise from incomplete data as well as from the presence of misleading ions in experimental mass spectra due to gas-phase reactions, stray ions, contaminants, and electronic noise. We employed a novel technique of reduction of false positives that is based on a combined use of orthogonal fragmentation techniques electron capture dissociation (ECD) and collisionally activated dissociation (CAD). Since ECD and CAD exhibit many complementary properties, their combined use greatly increased the analysis specificity, which was further strengthened by the high mass accuracy (≈1 ppm) afforded by Fourier transform mass spectrometry. The utility of this approach is demonstrated on a whole cell lysate from Escherichia coli. Analysis was made using the data-dependent acquisition mode. Extraction of complementary sequence information was performed prior to data base search using in-house written software. Only masses involved in complementary pairs in the MS/MS spectrum from the same or orthogonal fragmentation techniques were submitted to the data base search. ECD/CAD identified twice as many proteins at a fixed statistically significant confidence level with on average a 64% higher Mascot score. The confidence in protein identification was hereby increased by more than 1 order of magnitude. The combined ECD/CAD searches were on average 20% faster than CAD-only searches. A specially developed test with scrambled MS/MS data revealed that the amount of false-positive identifications was dramatically reduced by the combined use of CAD and ECD.

Electron capture/transfer versus collisionally activated/induced dissociations: Solo or duet?

New ion fragmentation technologies—electron capture dissociation (ECD) and electron-transfer dissociation (ETD)—are based on interaction of multiply charged polypeptides with either free electrons (ECD) or anionic species (ETD). After initial difficulties, these ECD/ETD (ExD) technologies are now being increasingly implemented in high-throughput proteomics work. This critical analysis presents arguments for the combined use of ExD with the conventional low-energy collisional excitation CID/CAD (CxD). It is argued that the database search, a key technology in MS/MS-based proteomics, is vulnerable with respect to the incomplete sequence information obtainable with either of the techniques, peptide MS/MS homology being a major complicating factor. De novo sequencing is viewed as the only adequate answer to this challenge and it can be achieved only with combined use of ExD and CxD. The payoff in the form of additional sequence information is projected to exceed the costs of such implementation. The greatest impact of combining ExD and CxD is expected in high-resolution instruments.

High-Resolution Enabled TMT 8-plexing

Isobaric mass tag-based quantitative proteomics strategies such as iTRAQ and TMT utilize reporter ions in the low-mass range of tandem MS spectra for relative quantification. The number of samples that can be compared in a single experiment (multiplexing) is limited by the number of different reporter ions that can be generated by differential stable isotope incorporation ((15)N, (13)C) across the reporter and the mass balancing parts of the reagents. Here, we demonstrate that a higher multiplexing rate can be achieved by utilizing the 6 mDa mass difference between (15)N- and (13)C-containing reporter fragments, in combination with high-resolution mass spectrometry. Two variants of the TMT127 and TMT129 reagents are available; these are distinguished by the position and the nature of the incorporated stable isotope in the reporter portions of the labels (TMT127L, (12)C(8)H(16)(15)N(1)(+); TMT127H, (12)C(7)(13)C(1)H(16)(14)N(1)(+); TMT129L, (12)C(6)(13)C(2)H(16)(15)N(1)(+); and TMT129H, (12)C(5)(13)C(3)H(16)(14)N(1)(+)). We demonstrate that these variants can be baseline-resolved in Orbitrap Elite higher-energy collision-induced dissociation spectra recorded with a 96 ms transient enabling comparable dynamic range, precision, and accuracy of quantification as 1 Da spaced reporter ions. The increased multiplexing rate enabled determination of inhibitor potencies in chemoproteomic kinase assays covering a wider range of compound concentrations in a single experiment, compared to conventional 6-plex TMT-based assays.

A selective inhibitor reveals PI3Kγ dependence of T H 17 cell differentiation

We devised a high-throughput chemoproteomics method that enabled multiplexed screening of 16,000 compounds against native protein and lipid kinases in cell extracts. Optimization of one chemical series resulted in CZC24832, which is to our knowledge the first selective inhibitor of phosphoinositide 3-kinase γ (PI3Kγ) with efficacy in in vitro and in vivo models of inflammation. Extensive target- and cell-based profiling of CZC24832 revealed regulation of interleukin-17-producing T helper cell (T(H)17) differentiation by PI3Kγ, thus reinforcing selective inhibition of PI3Kγ as a potential treatment for inflammatory and autoimmune diseases.