A data analysis framework for combining multiple batches increases the power of isobaric proteomics experiments.

We present a framework for the analysis of multiplexed mass spectrometry proteomics data that reduces estimation error when combining multiple isobaric batches. Variations in the number and quality of observations have long complicated the analysis of isobaric proteomics data. Here we show that the power to detect statistical associations is substantially improved by utilizing models that directly account for known sources of variation in the number and quality of observations that occur across batches.In a multibatch benchmarking experiment, our open-source software (msTrawler) increases the power to detect changes, especially in the range of less than twofold changes, while simultaneously increasing quantitative proteome coverage by utilizing more low-signal observations. Further analyses of previously published multiplexed datasets of 4 and 23 batches highlight both increased power and the ability to navigate complex missing data patterns without relying on unverifiable imputations or discarding reliable measurements.

Discovery of reactive peptide inhibitors of human papillomavirus oncoprotein E6.

Human papillomavirus (HPV) infections account for nearly all cervical cancer cases, which is the fourth most common cancer in women worldwide. High-risk variants, including HPV16, drive tumorigenesis in part by promoting the degradation of the tumor suppressor p53. This degradation is mediated by the HPV early protein 6 (E6), which recruits the E3 ubiquitin ligase E6AP and redirects its activity towards ubiquitinating p53. Targeting the protein interaction interface between HPV E6 and E6AP is a promising modality to mitigate HPV-mediated degradation of p53. In this study, we designed a covalent peptide inhibitor, termed reactide, that mimics the E6AP LXXLL binding motif by selectively targeting cysteine 58 in HPV16 E6 with quantitative conversion. This reactide provides a starting point in the development of covalent peptidomimetic inhibitors for intervention against HPV-driven cancers.

High-Throughput Proteome Profiling of Plasma and Native Plasma Complexes Using Native Chromatography.

We describe a high-throughput method for co-fractionation mass spectrometry (CF-MS) profiling for native plasma protein profiling. CF-MS allows the profiling of endogenous protein complexes between samples. Proteins often interact with other proteins and form macromolecular complexes that are different in disease states as well as cell states and cell types. This protocol describes an example for the sample preparation of 954 individual size exclusion chromatography (SEC) fractions, derived from 18 plasma samples that were separated into 53 fractions. Eighteen plasma samples were chosen based on the TMTpro multiplexing, but this methodology can be adapted for fewer or larger numbers of samples as appropriate. Our automated sample preparation method allows for high-throughput native plasma profiling, and we provide detailed methods for both a label-free and an isobaric labeling approach, discuss the merits of each approach, and detail the advantages of combining these strategies for comprehensive native plasma proteome profiling.

Conditional Fragment Ion Probabilities Improve Database Searching for Nonmonoisotopic Precursors.

Stochastic, intensity-based precursor isolation can result in isotopically enriched fragment ions. This problem is exacerbated for large peptides and stable isotope labeling experiments using deuterium or 15N. For stable isotope labeling experiments, incomplete and ubiquitous labeling strategies result in the isolation of peptide ions composed of many distinct structural isomers. Unfortunately, existing proteomics search algorithms do not account for this variability in isotopic incorporation, and thus often yield poor peptide and protein identification rates. We sought to resolve this shortcoming by deriving the expected isotopic distributions of each fragment ion and incorporating them into the theoretical mass spectra used for peptide-spectrum-matching. We adapted the Comet search platform to integrate a modified spectral prediction algorithm we term Conditional fragment Ion Distribution Search (CIDS). Comet-CIDS uses a traditional database searching strategy, but for each candidate peptide we compute the isotopic distribution of each fragment to better match the observed m/z distributions. Evaluating previously generated D2O and 15N labeled data sets, we found that Comet-CIDS identified more confident peptide spectral matches and higher protein sequence coverage compared to traditional theoretical spectra generation, with the magnitude of improvement largely determined by the amount of labeling in the sample.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology.

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

Triple-threat quantitative multiplexed plasma proteomics analysis on immune complex disease MRL-lpr mice.

Systemic lupus erythematosus is a common autoimmune inflammatory disease which is associated with increases in autoantibodies and immune complexes that deposit in the kidney. The MRL-lpr mouse is a common mouse model used for the study of lupus and immune complex glomerulonephritis but very little is known about the plasma proteome changes in this model. We performed in-depth quantitative proteome profiling on MRL-lpr and control (strain MpJ) mice to investigate the changes in the proteome, immunoglobulins and their glycoproteome as well as protein and immune complexes. Methodologies used included immunohistochemistry, immunoglobulin isotyping, multiplexed proteome profiling, immunoglobulin immunoprecipitation with glycoproteome profiling, and size exclusion chromatography (SEC) profiling to enable a comprehensive proteome profiling of proteins and protein complexes. We also used a novel native multiplexed plasma proteome profiling (NativeMP3) method that relies on native enrichment of plasma proteins enabling ultra-deep single shot profiling where we identified 922 plasma proteins at 1% false discovery rate (FDR) in a single shot mass spectrometry run. We observed many large plasma protein differences between the MRL-lpr and control strain including differences in the immunoglobulins, immunoglobulins against specific antigens, chemokines, and proteases as well as changes in protein complexes such as the immunoproteasome.

Novel insights from a multiomics dissection of the Hayflick limit.

The process wherein dividing cells exhaust proliferative capacity and enter into replicative senescence has become a prominent model for cellular aging in vitro. Despite decades of study, this cellular state is not fully understood in culture and even much less so during aging. Here, we revisit Leonard Hayflick's original observation of replicative senescence in WI-38 human lung fibroblasts equipped with a battery of modern techniques including RNA-seq, single-cell RNA-seq, proteomics, metabolomics, and ATAC-seq. We find evidence that the transition to a senescent state manifests early, increases gradually, and corresponds to a concomitant global increase in DNA accessibility in nucleolar and lamin associated domains. Furthermore, we demonstrate that senescent WI-38 cells acquire a striking resemblance to myofibroblasts in a process similar to the epithelial to mesenchymal transition (EMT) that is regulated by t YAP1/TEAD1 and TGF-ß2. Lastly, we show that verteporfin inhibition of YAP1/TEAD1 activity in aged WI-38 cells robustly attenuates this gene expression program.

Automated 16-Plex Plasma Proteomics with Real-Time Search and Ion Mobility Mass Spectrometry Enables Large-Scale Profiling in Naked Mole-Rats and Mice.

Performing large-scale plasma proteome profiling is challenging due to limitations imposed by lengthy preparation and instrument time. We present a fully automated multiplexed proteome profiling platform (AutoMP3) using the Hamilton Vantage liquid handling robot capable of preparing hundreds to thousands of samples. To maximize protein depth in single-shot runs, we combined 16-plex Tandem Mass Tags (TMTpro) with high-field asymmetric waveform ion mobility spectrometry (FAIMS Pro) and real-time search (RTS). We quantified over 40 proteins/min/sample, doubling the previously published rates. We applied AutoMP3 to investigate the naked mole-rat plasma proteome both as a function of the circadian cycle and in response to ultraviolet (UV) treatment. In keeping with the lack of synchronized circadian rhythms in naked mole-rats, we find few circadian patterns in plasma proteins over the course of 48 h. Furthermore, we quantify many disparate changes between mice and naked mole-rats at both 48 h and one week after UV exposure. These species differences in plasma protein temporal responses could contribute to the pronounced cancer resistance observed in naked mole-rats. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [1] partner repository with the dataset identifier PXD022891.

Systematic identification of cancer cell vulnerabilities to natural killer cell-mediated immune surveillance.

Only a subset of cancer patients respond to T-cell checkpoint inhibitors, highlighting the need for alternative immunotherapeutics. We performed CRISPR-Cas9 screens in a leukemia cell line to identify perturbations that enhance natural killer effector functions. Our screens defined critical components of the tumor-immune synapse and highlighted the importance of cancer cell interferon-γ signaling in modulating NK activity. Surprisingly, disrupting the ubiquitin ligase substrate adaptor DCAF15 strongly sensitized cancer cells to NK-mediated clearance. DCAF15 disruption induced an inflamed state in leukemic cells, including increased expression of lymphocyte costimulatory molecules. Proteomic and biochemical analysis revealed that cohesin complex members were endogenous client substrates of DCAF15. Genetic disruption of DCAF15 was phenocopied by treatment with indisulam, an anticancer drug that functions through DCAF15 engagement. In AML patients, reduced DCAF15 expression was associated with improved survival. These findings suggest that DCAF15 inhibition may have useful immunomodulatory properties in the treatment of myeloid neoplasms.

eIF2B activator prevents neurological defects caused by a chronic integrated stress response.

The integrated stress response (ISR) attenuates the rate of protein synthesis while inducing expression of stress proteins in cells. Various insults activate kinases that phosphorylate the GTPase eIF2 leading to inhibition of its exchange factor eIF2B. Vanishing White Matter (VWM) is a neurological disease caused by eIF2B mutations that, like phosphorylated eIF2, reduce its activity. We show that introduction of a human VWM mutation into mice leads to persistent ISR induction in the central nervous system. ISR activation precedes myelin loss and development of motor deficits. Remarkably, long-term treatment with a small molecule eIF2B activator, 2BAct, prevents all measures of pathology and normalizes the transcriptome and proteome of VWM mice. 2BAct stimulates the remaining activity of mutant eIF2B complex in vivo, abrogating the maladaptive stress response. Thus, 2BAct-like molecules may provide a promising therapeutic approach for VWM and provide relief from chronic ISR induction in a variety of disease contexts.

Proteasomal control of cytokinin synthesis protects Mycobacterium tuberculosis against nitric oxide.

One of several roles of the Mycobacterium tuberculosis proteasome is to defend against host-produced nitric oxide (NO), a free radical that can damage numerous biological macromolecules. Mutations that inactivate proteasomal degradation in Mycobacterium tuberculosis result in bacteria that are hypersensitive to NO and attenuated for growth in vivo, but it was not known why. To elucidate the link between proteasome function, NO resistance, and pathogenesis, we screened for suppressors of NO hypersensitivity in a mycobacterial proteasome ATPase mutant and identified mutations in Rv1205. We determined that Rv1205 encodes a pupylated proteasome substrate. Rv1205 is a homolog of the plant enzyme LONELY GUY, which catalyzes the production of hormones called cytokinins. Remarkably, we report that an obligate human pathogen secretes several cytokinins. Finally, we determined that the Rv1205-dependent accumulation of cytokinin breakdown products is likely responsible for the sensitization of Mycobacterium tuberculosis proteasome-associated mutants to NO.

Effects of MEK inhibitors GSK1120212 and PD0325901 in vivo using 10-plex quantitative proteomics and phosphoproteomics.

Multiplexed isobaric tag based quantitative proteomics and phosphoproteomics strategies can comprehensively analyze drug treatments effects on biological systems. Given the role of mitogen-activated protein/extracellular signal-regulated kinase (MEK) signaling in cancer and mitogen-activated protein kinase (MAPK)-dependent diseases, we sought to determine if this pathway could be inhibited safely by examining the downstream molecular consequences. We used a series of tandem mass tag 10-plex experiments to analyze the effect of two MEK inhibitors (GSK1120212 and PD0325901) on three tissues (kidney, liver, and pancreas) from nine mice. We quantified ∼ 6000 proteins in each tissue, but significant protein-level alterations were minimal with inhibitor treatment. Of particular interest was kidney tissue, as edema is an adverse effect of these inhibitors. From kidney tissue, we enriched phosphopeptides using titanium dioxide (TiO2 ) and quantified 10 562 phosphorylation events. Further analysis by phosphotyrosine peptide immunoprecipitation quantified an additional 592 phosphorylation events. Phosphorylation motif analysis revealed that the inhibitors decreased phosphorylation levels of proline-x-serine-proline (PxSP) and serine-proline (SP) sites, consistent with extracellular-signal-regulated kinase (ERK) inhibition. The MEK inhibitors had the greatest decrease on the phosphorylation of two proteins, Barttin and Slc12a3, which have roles in ion transport and fluid balance. Further studies will provide insight into the effect of these MEK inhibitors with respect to edema and other adverse events in mouse models and human patients.

An ERK/Cdk5 axis controls the diabetogenic actions of PPARγ.

Obesity-linked insulin resistance is a major precursor to the development of type 2 diabetes. Previous work has shown that phosphorylation of PPARγ (peroxisome proliferator-activated receptor γ) at serine 273 by cyclin-dependent kinase 5 (Cdk5) stimulates diabetogenic gene expression in adipose tissues. Inhibition of this modification is a key therapeutic mechanism for anti-diabetic drugs that bind PPARγ, such as the thiazolidinediones and PPARγ partial agonists or non-agonists. For a better understanding of the importance of this obesity-linked PPARγ phosphorylation, we created mice that ablated Cdk5 specifically in adipose tissues. These mice have both a paradoxical increase in PPARγ phosphorylation at serine 273 and worsened insulin resistance. Unbiased proteomic studies show that extracellular signal-regulated kinase (ERK) kinases are activated in these knockout animals. Here we show that ERK directly phosphorylates serine 273 of PPARγ in a robust manner and that Cdk5 suppresses ERKs through direct action on a novel site in MAP kinase/ERK kinase (MEK). Importantly, pharmacological inhibition of MEK and ERK markedly improves insulin resistance in both obese wild-type and ob/ob mice, and also completely reverses the deleterious effects of the Cdk5 ablation. These data show that an ERK/Cdk5 axis controls PPARγ function and suggest that MEK/ERK inhibitors may hold promise for the treatment of type 2 diabetes.

Assessing subunit dependency of the Plasmodium proteasome using small molecule inhibitors and active site probes.

The ubiquitin-proteasome system (UPS) is a potential pathway for therapeutic intervention for pathogens such as Plasmodium, the causative agent of malaria. However, due to the essential nature of this proteolytic pathway, proteasome inhibitors must avoid inhibition of the host enzyme complex to prevent toxic side effects. The Plasmodium proteasome is poorly characterized, making rational design of inhibitors that induce selective parasite killing difficult. In this study, we developed a chemical probe that labels all catalytic sites of the Plasmodium proteasome. Using this probe, we identified several subunit selective small molecule inhibitors of the parasite enzyme complex. Treatment with an inhibitor that is specific for the ß5 subunit during blood stage schizogony led to a dramatic decrease in parasite replication while short-term inhibition of the ß2 subunit did not affect viability. Interestingly, coinhibition of both the ß2 and ß5 catalytic subunits resulted in enhanced parasite killing at all stages of the blood stage life cycle and reduced parasite levels in vivo to barely detectable levels. Parasite killing was achieved with overall low host toxicity, something that has not been possible with existing proteasome inhibitors. Our results highlight differences in the subunit dependency of the parasite and human proteasome, thus providing a strategy for development of potent antimalarial drugs with overall low host toxicity.

iPSC-derived neurons from GBA1-associated Parkinson's disease patients show autophagic defects and impaired calcium homeostasis.

Mutations in the acid ß-glucocerebrosidase (GBA1) gene, responsible for the lysosomal storage disorder Gaucher's disease (GD), are the strongest genetic risk factor for Parkinson's disease (PD) known to date. Here we generate induced pluripotent stem cells from subjects with GD and PD harbouring GBA1 mutations, and differentiate them into midbrain dopaminergic neurons followed by enrichment using fluorescence-activated cell sorting. Neurons show a reduction in glucocerebrosidase activity and protein levels, increase in glucosylceramide and α-synuclein levels as well as autophagic and lysosomal defects. Quantitative proteomic profiling reveals an increase of the neuronal calcium-binding protein 2 (NECAB2) in diseased neurons. Mutant neurons show a dysregulation of calcium homeostasis and increased vulnerability to stress responses involving elevation of cytosolic calcium. Importantly, correction of the mutations rescues such pathological phenotypes. These findings provide evidence for a link between GBA1 mutations and complex changes in the autophagic/lysosomal system and intracellular calcium homeostasis, which underlie vulnerability to neurodegeneration.

Carboxyl-ester lipase maturity-onset diabetes of the young is associated with development of pancreatic cysts and upregulated MAPK signaling in secretin-stimulated duodenal fluid.

Carboxyl-ester lipase (CEL) maturity-onset diabetes of the young (MODY) is a monogenic form of diabetes and pancreatic exocrine dysfunction due to mutations in the CEL gene encoding CEL. The pathogenic mechanism for diabetes development is unknown. Since CEL is expressed mainly in pancreatic acinar cells, we asked whether we could find structural pancreatic changes in CEL-MODY subjects during the course of diabetes development. Furthermore, we hypothesized that the diseased pancreas releases proteins that are detectable in pancreatic fluid and potentially reflect activation or inactivation of disease-specific pathways. We therefore investigated nondiabetic and diabetic CEL-mutation carriers by pancreatic imaging studies and secretin-stimulated duodenal juice sampling. The secretin-stimulated duodenal juice was studied using cytokine assays, mass spectrometry (MS) proteomics, and multiplexed MS-based measurement of kinase activities. We identified multiple pancreatic cysts in all eight diabetic mutation carriers but not in any of the four nondiabetic mutation carriers or the six healthy controls. Furthermore, we identified upregulated mitogen-activated protein kinase (MAPK) target proteins and MAPK-driven cytokines and increased MAPK activity in the secretin-stimulated duodenal juice. These findings show that subjects with CEL-MODY develop multiple pancreatic cysts by the time they develop diabetes and that upregulated MAPK signaling in the pancreatic secretome may reflect the pathophysiological development of pancreatic cysts and diabetes.

Increasing throughput in targeted proteomics assays: 54-plex quantitation in a single mass spectrometry run.

Targeted proteomics assays such as those measuring end points in activity assays are sensitive and specific but often lack in throughput. In an effort to significantly increase throughput, a comparison was made between the traditional approach which utilizes an internal standard and the multiplexing approach which relies on isobaric tagging. A kinase activity assay was used for proof of concept, and experiments included three biological replicates for every condition. Results from the two approaches were highly similar with the multiplexing showing greater throughput. Two novel 6-plex isobaric tags were added for a total of three 6-plex experiments (18-plex) in a single run. Next, three mass variants of the target peptide were labeled with the three isobaric tags giving nine 6-plex reactions for 54-plex quantitation in a single run. Since the multiplexing approach allows all samples to be combined prior to purification and acquisition, the 54-plex approach resulted in a significant reduction in purification resources (time, reagents, etc.) and a ~50-fold improvement in acquisition throughput. We demonstrate the 54-plex assay in several ways including measuring inhibition of PKA activity in MCF7 cell lysates for a panel of nine compounds.

Absence of diabetes and pancreatic exocrine dysfunction in a transgenic model of carboxyl-ester lipase-MODY (maturity-onset diabetes of the young).

BACKGROUND: CEL-MODY is a monogenic form of diabetes with exocrine pancreatic insufficiency caused by mutations in CARBOXYL-ESTER LIPASE (CEL). The pathogenic processes underlying CEL-MODY are poorly understood, and the global knockout mouse model of the CEL gene (CELKO) did not recapitulate the disease. We therefore aimed to create and phenotype a mouse model specifically over-expressing mutated CEL in the pancreas. METHODS: We established a monotransgenic floxed (flanking LOX sequences) mouse line carrying the human CEL mutation c.1686delT and crossed it with an elastase-Cre mouse to derive a bitransgenic mouse line with pancreas-specific over-expression of CEL carrying this disease-associated mutation (TgCEL). Following confirmation of murine pancreatic expression of the human transgene by real-time quantitative PCR, we phenotyped the mouse model fed a normal chow and compared it with mice fed a 60% high fat diet (HFD) as well as the effects of short-term and long-term cerulein exposure. RESULTS: Pancreatic exocrine function was normal in TgCEL mice on normal chow as assessed by serum lipid and lipid-soluble vitamin levels, fecal elastase and fecal fat absorption, and the normoglycemic mice exhibited normal pancreatic morphology. On 60% HFD, the mice gained weight to the same extent as controls, had normal pancreatic exocrine function and comparable glucose tolerance even after resuming normal diet and follow up up to 22 months of age. The cerulein-exposed TgCEL mice gained weight and remained glucose tolerant, and there were no detectable mutation-specific differences in serum amylase, islet hormones or the extent of pancreatic tissue inflammation. CONCLUSIONS: In this murine model of human CEL-MODY diabetes, we did not detect mutation-specific endocrine or exocrine pancreatic phenotypes, in response to altered diets or exposure to cerulein.

Correlation profiling for determining kinase-substrate relationships.

Arguably the most important issue in phosphorylation studies is the determination of kinase-substrate relationships. Kinase Activity-Abundance Correlation (KAAC) profiling is a technique that can be used to suggest the kinase responsible for a particular phosphorylation event. The method involves the separation of a lysate at the protein level, after which fractions are analyzed to produce an activity profile for each substrate peptide. This activity profile is correlated with the kinase abundance profiles (obtained using shotgun proteomics) to obtain candidate kinases responsible for phosphorylation of the substrate peptide. We demonstrate this approach by determining potential kinase-substrate pairs for six peptides, selected from the literature, that were shown to be upregulated during mitosis. Finally, for a subset of fractions we evaluated the use of stable isotope labeling to improve quantification.

Mycobacterium tuberculosis prokaryotic ubiquitin-like protein-deconjugating enzyme is an unusual aspartate amidase.

Deamidase of Pup (Dop), the prokaryotic ubiquitin-like protein (Pup)-deconjugating enzyme, is critical for the full virulence of Mycobacterium tuberculosis and is unique to bacteria, providing an ideal target for the development of selective chemotherapies. We used a combination of genetics and chemical biology to characterize the mechanism of depupylation. We identified an aspartate as a potential nucleophile in the active site of Dop, suggesting a novel protease activity to target for inhibitor development.