Using a simple, environment friendly proteome extraction (TOP), we were able to optimize the analysis of clinical samples. Using our TOP method we analyzed a clinical cohort of microsatellite stable (MSS) and unstable (MSI-H) colorectal carcinoma (CRC). We identified a tumor cell specific, STAT1-centered, immune signature expressed by the MSI-H tumor cells. We then showed that long, but not short, exposure to Interferon- induces a similar signature in vitro. We identified 10 different temporal protein expression patterns, classifying the Interferon- protein temporal regulation in CRC. Our data sheds light on the changes that tumor cells undergo under long-term immunological pressure in vivo, the importance of STAT proteins in specific biological scenarios. The data generated could help find novel clinical biomarkers and therapeutic approaches.

A key point in achieving accurate intact glycopeptide identification is the definition of the glycan composition file that is used to match experimental with theoretical masses by a glycoproteomics search engine. At present, these files are mainly built from searching the literature and/or querying data sources focused on posttranslational modifications. Most glycoproteomics search engines include a default composition file that is readily used when processing MS data. We introduce here a glycan composition visualizing and comparative tool associated with the GlyConnect database and called GlyConnect Compozitor. It offers a web interface through which the database can be queried to bring out contextual information relative to a set of glycan compositions. The tool takes advantage of compositions being related to one another through shared monosaccharide counts and outputs interactive graphs summarizing information searched in the database. These results provide a guide for selecting or deselecting compositions in a file in order to reflect the context of a study as closely as possible. They also confirm the consistency of a set of compositions based on the content of the GlyConnect database. As part of the tool collection of the Glycomics@ExPASy initiative, Compozitor is hosted at https://glyconnect.expasy.org/compozitor/ where it can be run as a web application. It is also directly accessible from the GlyConnect database.

The Discoidin, CUB, and LCCL domain-containing protein (DCBLD) family consists of two type-I transmembrane scaffolding receptors, DCBLD1 and DCBLD2, which play important roles in development and cancer. The nonreceptor tyrosine kinases FYN and ABL are known to drive phosphorylation of tyrosine residues in YXXP motifs within the intracellular domains of DCBLD family members, which leads to the recruitment of the Src homology 2 (SH2) domain of the adaptors CT10 regulator of kinase (CRK) and CRK-like (CRKL). We previously characterized the FYN- and ABL-driven phosphorylation of DCBLD family YXXP motifs. However, we have identified additional FYN- and ABL-dependent phosphorylation sites on DCBLD1 and DCBLD2. This suggests that beyond CRK and CRKL, additional DCBLD interactors may be regulated by FYN and ABL activity. Here, we report a quantitative proteomics approach in which we map the FYN- and ABL-regulated interactomes of DCBLD family members. We found FYN and ABL regulated the binding of several signaling molecules to DCBLD1 and DCBLD2, including members of the 14-3-3 family of adaptors. Biochemical investigation of the DCBLD2/14-3-3 interaction revealed ABL-induced binding of 14-3-3 family members directly to DCBLD2.

Glutathionylation is an important posttranslational modification that protects proteins from further oxidative damage as well as influencing protein structure and activity. In the present study, we demonstrate that the cysteine-42 residue in protein arginine N-methyltransferase 5 (PRMT5) is glutathionylated in aged mice or in cells that have been exposed to oxidative stress. Deglutathionylation of this protein is catalyzed by glutaredoxin-1 (Grx1). Using mutagenesis and subsequent biochemical analyses, we show that glutathionylation decreased the binding affinity of PRMT5 with methylosome protein-50 (MEP50) and reduced the methyltransferase activity of PRMT5. Furthermore, overexpression of PRMT5-C42A mutant caused a significant increase in histone methylation in HEK293T and A549 cells and promoted cell growth, whereas overexpression of the PRMT5-C42D mutant, a mimic of glutathionylated PRMT5, inhibited cell proliferation. Taken together, our results demonstrate a new mechanism of regulation of PRMT5 methyltransferases activity and suggest that PRMT5 glutathionylation is partly responsible for reactive oxygen species-mediated cell growth inhibition.

Studies in the yeast Saccharomyces cerevisiae have helped define mechanisms underlying the activity of the ubiquitin–proteasome system (UPS), uncover the proteasome assembly pathway, and link the UPS to the maintenance of cellular homeostasis. However, the spectrum of UPS substrates is incompletely defined, even though multiple techniques—including MS—have been used. Therefore, we developed a substrate trapping proteomics workflow to identify previously unknown UPS substrates. We first generated a yeast strain with an epitope tagged proteasome subunit to which a proteasome inhibitor could be applied. Parallel experiments utilized inhibitor insensitive strains or strains lacking the tagged subunit. After affinity isolation, enriched proteins were resolved, in-gel digested, and analyzed by high resolution liquid chromatography-tandem MS. A total of 149 proteasome partners were identified, including all 33 proteasome subunits. When we next compared data between inhibitor sensitive and resistant cells, 27 proteasome partners were significantly enriched. Among these proteins were known UPS substrates and proteins that escort ubiquitinated substrates to the proteasome. We also detected Erg25 as a high-confidence partner. Erg25 is a methyl oxidase that converts dimethylzymosterol to zymosterol, a precursor of the plasma membrane sterol, ergosterol. Because Erg25 is a resident of the endoplasmic reticulum (ER) and had not previously been directly characterized as a UPS substrate, we asked whether Erg25 is a target of the ER associated degradation (ERAD) pathway, which most commonly mediates proteasome-dependent destruction of aberrant proteins. As anticipated, Erg25 was ubiquitinated and associated with stalled proteasomes. Further, Erg25 degradation depended on ERAD-associated ubiquitin ligases and was regulated by sterol synthesis. These data expand the cohort of lipid biosynthetic enzymes targeted for ERAD, highlight the role of the UPS in maintaining ER function, and provide a novel tool to uncover other UPS substrates via manipulations of our engineered strain.

Co-fractionation MS (CF-MS) is a technique with potential to characterize endogenous and unmanipulated protein complexes on an unprecedented scale. However this potential has been offset by a lack of guidelines for best-practice CF-MS data collection and analysis. To obtain such guidelines, this study thoroughly evaluates novel and published Saccharomyces cerevisiae CF-MS data sets using very high proteome coverage libraries of yeast gold standard complexes. A new method for identifying gold standard complexes in CF-MS data, Reference Complex Profiling, and the Extending 'Guilt-by-Association' by Degree (EGAD) R package are used for these evaluations, which are verified with concurrent analyses of published human data. By evaluating data collection designs, which involve fractionation of cell lysates, it is found that near-maximum recall of complexes can be achieved with fewer samples than published studies. Distributing sample collection across orthogonal fractionation methods, rather than a single high resolution data set, leads to particularly efficient recall. By evaluating 17 different similarity scoring metrics, which are central to CF-MS data analysis, it is found that two metrics rarely used in past CF-MS studies – Spearman and Kendall correlations – and the recently introduced Co-apex metric frequently maximize recall, whereas a popular metric—Euclidean distance—delivers poor recall. The common practice of integrating external genomic data into CF-MS data analysis is also evaluated, revealing that this practice may improve the precision and recall of known complexes but is generally unsuitable for predicting novel complexes in model organisms. If studying nonmodel organisms using orthologous genomic data, it is found that particular subsets of fractionation profiles (e.g. the lowest abundance quartile) should be excluded to minimize false discovery. These assessments are summarized in a series of universally applicable guidelines for precise, sensitive and efficient CF-MS studies of known complexes, and effective predictions of novel complexes for orthogonal experimental validation.

After ejaculation, mammalian spermatozoa must undergo a process known as capacitation in order to successfully fertilize the oocyte. Several post-translational modifications occur during capacitation, including sialylation, which despite being limited to a few proteins, seems to be essential for proper sperm-oocyte interaction. Regardless of its importance, to date, no single study has ever identified nor quantified which glycoproteins bearing terminal sialic acid (Sia) are altered during capacitation. Here we characterize sialylation during mouse sperm capacitation. Using tandem MS coupled with liquid chromatography (LC–MS/MS), we found 142 nonreductant peptides, with 9 of them showing potential modifications on their sialylated oligosaccharides during capacitation. As such, N-linked sialoglycopeptides from C4b-binding protein, endothelial lipase (EL), serine proteases 39 and 52, testis-expressed protein 101 and zonadhesin were reduced following capacitation. In contrast, mitochondrial aconitate hydratase (aconitase; ACO2), a TCA cycle enzyme, was the only protein to show an increase in Sia content during capacitation. Interestingly, although the loss of Sia within EL (N62) was accompanied by a reduction in its phospholipase A₁ activity, a decrease in the activity of ACO2 (i.e. stereospecific isomerization of citrate to isocitrate) occurred when sialylation increased (N612). The latter was confirmed by N612D recombinant protein tagged with both His and GFP. The replacement of Sia for the negatively charged Aspartic acid in the N612D mutant caused complete loss of aconitase activity compared with the WT. Computer modeling show that N612 sits atop the catalytic site of ACO2. The introduction of Sia causes a large conformational change in the alpha helix, essentially, distorting the active site, leading to complete loss of function. These findings suggest that the switch from oxidative phosphorylation, over to glycolysis that occurs during capacitation may come about through sialylation of ACO2.

Renal Cell Carcinoma (RCC) is one of the most commonly diagnosed cancers worldwide with research efforts dramatically improving understanding of the biology of the disease. To investigate the role of the immune system in treatment-naïve clear cell Renal Cell Carcinoma (ccRCC), we interrogated the immune infiltrate in patient-matched ccRCC tumor samples, benign normal adjacent tissue (NAT) and peripheral blood mononuclear cells (PBMCs isolated from whole blood, focusing our attention on the myeloid cell infiltrate. Using flow cytometric, MS, and ExCYT analysis, we discovered unique myeloid populations in PBMCs across patient samples. Furthermore, normal adjacent tissues and ccRCC tissues contained numerous myeloid populations with a unique signature for both tissues. Enrichment of the immune cell (CD45⁺) fraction and subsequent gene expression analysis revealed a number of myeloid-related genes that were differentially expressed. These data provide evidence, for the first time, of an immunosuppressive and pro-tumorigenic role of myeloid cells in early, clinically localized ccRCC. The identification of a number of immune proteins for therapeutic targeting provides a rationale for investigation into the potential efficacy of earlier intervention with single-agent or combination immunotherapy for ccRCC.

Protein synthesis on the endoplasmic reticulum (ER) requires the dynamic coordination of numerous cellular components. Together, resident ER membrane proteins, cytoplasmic translation factors, and both integral membrane and cytosolic RNA-binding proteins operate in concert with membrane-associated ribosomes to facilitate ER-localized translation. Little is known, however, regarding the spatial organization of ER-localized translation. This question is of growing significance as it is now known that ER-bound ribosomes contribute to secretory, integral membrane, and cytosolic protein synthesis alike. To explore this question, we utilized quantitative proximity proteomics to identify neighboring protein networks for the candidate ribosome interactors SEC61β (subunit of the protein translocase), RPN1 (oligosaccharyltransferase subunit), SEC62 (translocation integral membrane protein), and LRRC59 (ribosome binding integral membrane protein). Biotin labeling time course studies of the four BioID reporters revealed distinct labeling patterns that intensified but only modestly diversified as a function of labeling time, suggesting that the ER membrane is organized into discrete protein interaction domains. Whereas SEC61β and RPN1 reporters identified translocon-associated networks, SEC62 and LRRC59 reporters revealed divergent protein interactomes. Notably, the SEC62 interactome is enriched in redox-linked proteins and ER luminal chaperones, with the latter likely representing proximity to an ER luminal chaperone reflux pathway. In contrast, the LRRC59 interactome is highly enriched in SRP pathway components, translation factors, and ER-localized RNA-binding proteins, uncovering a functional link between LRRC59 and mRNA translation regulation. Importantly, analysis of the LRRC59 interactome by native immunoprecipitation identified similar protein and functional enrichments. Moreover, [³⁵S]-methionine incorporation assays revealed that siRNA silencing of LRRC59 expression reduced steady state translation levels on the ER by ca. 50%, and also impacted steady state translation levels in the cytosol compartment. Collectively, these data reveal a functional domain organization for the ER and identify a key role for LRRC59 in the organization and regulation of local translation.

The EGFR tyrosine kinase inhibitor gefitinib is commonly used for lung cancer patients. However, some patients eventually become resistant to gefitinib and develop progressive disease. Here, we indicate that ecto-ATP synthase, which ectopically translocated from mitochondrial inner membrane to plasma membrane, is considered as a potential therapeutic target for drug-resistant cells. Quantitative multi-omics profiling reveals that ecto-ATP synthase inhibitor mediates CK2-dependent phosphorylation of DNA topoisomerase IIα (topo IIα) at serine 1106 and subsequently increases the expression of long noncoding RNA, GAS5. Additionally, we also determine that downstream of GAS5, p53 pathway, is activated by ecto-ATP synthase inhibitor for regulation of programed cell death. Interestingly, GAS5-proteins interactomic profiling elucidates that GAS5 associates with topo IIα and subsequently enhancing the phosphorylation level of topo IIα. Taken together, our findings suggest that ecto-ATP synthase blockade is an effective therapeutic strategy via regulation of CK2/phospho-topo IIα/GAS5 network in gefitinib-resistant lung cancer cells.

Small admixtures in water, e.g. of metal ions, often act as cell growth regulators. Here we report that enrichment of deuterium content in water, normally found at 8 mm concentration, two-three folds increases cell proliferation and lowers the oxidative stress level as well. Acting as an anti-oxidant, deuterium-enriched water prevents the toxic effect of such oxidative agents as hydrogen peroxide and auranofin. This action is opposite to that of deuterium depletion that is known to suppress cell growth and induce oxidative stress in mitochondria. We thus hypothesize that deuterium may be a natural cell growth regulator that controls mitochondrial oxidation-reduction balance. Because growth acceleration is reduced approximately by half by addition to water a minute amount (0.15%) of ¹⁸O isotope, at least part of the deuterium effect on cell growth can be explained by the isotopic resonance phenomenon. A slight (2-fold) enrichment of deuterium in water accelerates human cell growth. Quantitative MS based proteomics determined changes in protein abundances and redox states and found that deuterium-enriched water acts mainly through decreasing ROS production in mitochondria. This action is opposite to that of deuterium depletion that suppresses cell growth by inducing oxidative stress. Thus deuterium may be a natural cell growth regulator that controls mitochondrial oxidation-reduction balance. The role of isotopic resonance in this effect was validated by further experiments on bacteria.

Amino acid hydroxylation is a common post-translational modification, which generally regulates protein interactions or adds a functional group that can be further modified. Such hydroxylation is currently considered irreversible, necessitating the degradation and re-synthesis of the entire protein to reset the modification. Here we present evidence that the cellular machinery can reverse FIH-mediated asparagine hydroxylation on intact proteins. These data suggest that asparagine hydroxylation is a flexible and dynamic post-translational modification akin to modifications involved in regulating signaling networks, such as phosphorylation, methylation and ubiquitylation.

We performed an in-depth characterization and comparison of the pediatric and adult urinary glycomes using a nanoLC-MS/MS based glycomics method, which included normal healthy pediatric (1–10 years, n = 21) and adult (21–50 years, n = 22) individuals. A total of 116 N-glycan compositions were identified, and 46 of them could be reproducibly quantified. We performed quantitative comparisons of the 46 glycan compositions between different age and sex groups. The results showed significant quantitative changes between the pediatric and adult cohorts. The pediatric urinary N-glycome was found to contain a higher level of high-mannose (HM), asialylated/afucosylated glycans (excluding HM), neutral fucosylated and agalactosylated glycans, and a lower level of trisialylated glycans compared with the adult. We further analyzed gender-associated glycan changes in the pediatric and adult group, respectively. In the pediatric group, there was almost no difference of glycan levels between males and females. In adult, the majority of glycans were more abundant in males than females, except the high-mannose and tetrasialylated glycans. These findings highlight the importance to consider age-matching and adult sex-matching for urinary glycan studies. The identified normal pediatric and adult urinary glycomes can serve as a baseline reference for comparisons to other disease states affected by glycosylation.

Despite the crucial function of the small intestine in nutrient uptake our understanding of the molecular events underlying the digestive function is still rudimentary. Recent studies demonstrated that enterocytes do not direct the entire dietary triacylglycerol toward immediate chylomicron synthesis. Especially after high-fat challenges, parts of the resynthesized triacylglycerol are packaged into cytosolic lipid droplets for transient storage in the endothelial layer of the small intestine. The reason for this temporary storage of triacylglycerol is not completely understood. To utilize lipids from cytosolic lipid droplets for chylomicron synthesis in the endoplasmic reticulum, stored triacylglycerol has to be hydrolyzed either by cytosolic lipolysis or lipophagy. Interestingly, triacylglycerol storage and chylomicron secretion rates are unevenly distributed along the small intestine, with the proximal jejunum exhibiting the highest intermittent storage capacity. We hypothesize that correlating hydrolytic enzyme activities with the reported distribution of triacylglycerol storage and chylomicron secretion in different sections of the small intestine is a promising strategy to determine key enzymes in triacylglycerol remobilization. We employed a serine hydrolase specific activity-based labeling approach in combination with quantitative proteomics to identify and rank hydrolases based on their relative activity in 11 sections of the small intestine. Moreover, we identified several clusters of enzymes showing similar activity distribution along the small intestine. Merging our activity-based results with substrate specificity and subcellular localization known from previous studies, carboxylesterase 2e and arylacetamide deacetylase emerge as promising candidates for triacylglycerol mobilization from cytosolic lipid droplets in enterocytes.

Extracellular vesicles (EVs) secreted by the epididymal epithelium transfer to spermatozoa key proteins that are essential in promoting motility and subsequent fertilization success. Using the domestic cat model, the objectives were to (1) characterize and compare protein content of EVs between segments of the epididymis, and (2) compare EV protein compositions between normo- and teratospermic individuals (producing >60% of abnormal spermatozoa). Epididymal EVs from adult cats were isolated and assessed via liquid chromatography tandem MS. Both male types shared 3008 proteins in total, with 98 and 20 EV proteins unique to normospermic and teratospermic males, respectively. Expression levels of several proteins changed between epididymal segments in both male types. Several proteins in both groups were related to sperm motility (e.g. hexokinase 1, adenylate kinase isoenzyme) and zona pellucida or oolemma binding (e.g. disintegrin and metalloproteinase domain proteins, zona binding proteins 1 and 2). Interestingly, seven cauda-derived EV proteins trended downward in teratospermic compared with normospermic males, which may relate to poor sperm quality. Collective results revealed, for the first time, EV proteins related to sequential sperm maturation with differences observed between normospermic and teratospermic individuals.

Endometrial carcinoma (EC) is the most common gynecologic malignancy in the United States, with limited effective targeted therapies. Endometrial tumors exhibit frequent alterations in protein kinases, yet only a small fraction of the kinome has been therapeutically explored. To identify kinase therapeutic avenues for EC, we profiled the kinome of endometrial tumors and normal endometrial tissues using Multiplexed Inhibitor Beads and Mass Spectrometry (MIB-MS). Our proteomics analysis identified a network of kinases overexpressed in tumors, including Serine/Arginine-Rich Splicing Factor Kinase 1 (SRPK1). Immunohistochemical (IHC) analysis of endometrial tumors confirmed MIB-MS findings and showed SRPK1 protein levels were highly expressed in endometrioid and uterine serous cancer (USC) histological subtypes. Moreover, querying large-scale genomics studies of EC tumors revealed high expression of SRPK1 correlated with poor survival. Loss-of-function studies targeting SRPK1 in an established USC cell line demonstrated SRPK1 was integral for RNA splicing, as well as cell cycle progression and survival under nutrient deficient conditions. Profiling of USC cells identified a compensatory response to SRPK1 inhibition that involved EGFR and the up-regulation of IGF1R and downstream AKT signaling. Co-targeting SRPK1 and EGFR or IGF1R synergistically enhanced growth inhibition in serous and endometrioid cell lines, representing a promising combination therapy for EC.

The absence of the dystrophin protein in Duchenne muscular dystrophy (DMD) results in myofiber fragility and a plethora of downstream secondary pathologies. Although a variety of experimental therapies are in development, achieving effective treatments for DMD remains exceptionally challenging, not least because the pathological consequences of dystrophin loss are incompletely understood. Here we have performed proteome profiling in tibialis anterior muscles from two murine DMD models (mdx and mdx52) at three ages (8, 16, and 80 weeks of age), all n = 3. High-resolution isoelectric focusing liquid chromatography-tandem MS (HiRIEF-LC–MS/MS) was used to quantify the expression of 4974 proteins across all 27 samples. The two dystrophic models were found to be highly similar, whereas multiple proteins were differentially expressed relative to WT (C57BL/6) controls at each age. Furthermore, 1795 proteins were differentially expressed when samples were pooled across ages and dystrophic strains. These included numerous proteins associated with the extracellular matrix and muscle function that have not been reported previously. Pathway analysis revealed multiple perturbed pathways and predicted upstream regulators, which together are indicative of cross-talk between inflammatory, metabolic, and muscle growth pathways (e.g. TNF, INF, NF-B, SIRT1, AMPK, PGC-1α, PPARs, ILK, and AKT/PI3K). Upregulation of CAV3, MVP and PAK1 protein expression was validated in dystrophic muscle by Western blot. Furthermore, MVP was upregulated during, but not required for, the differentiation of C2C12 myoblasts suggesting that this protein may affect muscle regeneration. This study provides novel insights into mutation-independent proteomic signatures characteristic of the dystrophic phenotype and its progression with aging.

Sepsis-induced acute kidney injury (S-AKI) is the most common complication in hospitalized and critically ill patients, highlighted by a rapid decline of kidney function occurring a few hours or days after sepsis onset. Systemic inflammation elicited by microbial infections is believed to lead to kidney damage under immunocompromised conditions. However, although AKI has been recognized as a disease with long-term sequelae, partly because of the associated higher risk of chronic kidney disease (CKD), the understanding of kidney pathophysiology at the molecular level and the global view of dynamic regulations in situ after S-AKI, including the transition to CKD, remains limited. Existing studies of S-AKI mainly focus on deriving sepsis biomarkers from body fluids. In the present study, we constructed a mid-severity septic murine model using cecal ligation and puncture (CLP), and examined the temporal changes to the kidney proteome and phosphoproteome at day 2 and day 7 after CLP surgery, corresponding to S-AKI and the transition to CKD, respectively, by employing an ultrafast and economical filter-based sample processing method combined with the label-free quantitation approach. Collectively, we identified 2,119 proteins and 2950 phosphosites through multi-proteomics analyses. Among them, we identified an array of highly promising candidate marker proteins indicative of disease onset and progression accompanied by immunoblot validations, and further denoted the pathways that are specifically responsive to S-AKI and its transition to CKD, which include regulation of cell metabolism regulation, oxidative stress, and energy consumption in the diseased kidneys. Our data can serve as an enriched resource for the identification of mechanisms and biomarkers for sepsis-induced kidney diseases.

Specific E3 ligases target tumor suppressors for degradation. Inhibition of such E3 ligases may be an important approach to cancer treatment. RNF146 is a RING domain and PARylation-dependent E3 ligase that functions as an activator of the β-catenin/Wnt and YAP/Hippo pathways by targeting the degradation of several tumor suppressors. Tankyrases 1 and 2 (TNKS1/2) are the only known poly-ADP-ribosyltransferases that require RNF146 to degrade their substrates. However, systematic identification of RNF146 substrates have not yet been performed. To uncover substrates of RNF146 that are targeted for degradation, we generated RNF146 knockout cells and TNKS1/2-double knockout cells and performed proteome profiling with label-free quantification as well as transcriptome analysis. We identified 160 potential substrates of RNF146, which included many known substrates of RNF146 and TNKS1/2 and 122 potential TNKS-independent substrates of RNF146. In addition, we validated OTU domain-containing protein 5 and Protein mono-ADP-ribosyltransferase PARP10 as TNKS1/2-independent substrates of RNF146 and SARDH as a novel substrate of TNKS1/2 and RNF146. Our study is the first proteome-wide analysis of potential RNF146 substrates. Together, these findings not only demonstrate that proteome profiling can be a useful general approach for the systemic identification of substrates of E3 ligases but also reveal new substrates of RNF146, which provides a resource for further functional studies.

AAA+ ATPases constitute a large family of proteins that are involved in a plethora of cellular processes including DNA disassembly, protein degradation and protein complex disassembly. They typically form a hexametric ring-shaped structure with six subunits in a (pseudo) 6-fold symmetry. In a subset of AAA+ ATPases that facilitate protein unfolding and degradation, six subunits cooperate to translocate protein substrates through a central pore in the ring. The number and type of nucleotides in an AAA+ ATPase hexamer is inherently linked to the mechanism that underlies cooperation among subunits and couples ATP hydrolysis with substrate translocation. We conducted a native MS study of a monodispersed form of PAN, an archaeal proteasome AAA+ ATPase, to determine the number of nucleotides bound to each hexamer of the WT protein. We utilized ADP and its analogs (TNP-ADP and mant-ADP), and a nonhydrolyzable ATP analog (AMP-PNP) to study nucleotide site occupancy within the PAN hexamer in ADP- and ATP-binding states, respectively. Throughout all experiments we used a Walker A mutant (PAN^K217A) that is impaired in nucleotide binding as an internal standard to mitigate the effects of residual solvation on mass measurement accuracy and to serve as a reference protein to control for nonspecific nucleotide binding. This approach led to the unambiguous finding that a WT PAN hexamer carried – from expression host – six tightly bound ADP molecules that could be exchanged for ADP and ATP analogs. Although the Walker A mutant did not bind ADP analogs, it did bind AMP-PNP, albeit at multiple stoichiometries. We observed variable levels of hexamer dissociation and an appearance of multimeric species with the over-charged molecular ion distributions across repeated experiments. We posit that these phenomena originated during ESI process at the final stages of ESI droplet evolution.

Plasmodium, the malaria parasite, undergoes a complex life cycle alternating between a vertebrate host and a mosquito vector of the genus Anopheles. In red blood cells of the vertebrate host, Plasmodium multiplies asexually or differentiates into gamete precursors, the male and female gametocytes, responsible for parasite transmission. Sexual stage maturation occurs in the midgut of the mosquito vector, where male and female gametes egress from the host erythrocytes to fuse and form a zygote. Gamete egress entails the successive rupture of two membranes surrounding the parasite, the parasitophorous vacuole membrane and the erythrocyte plasma membrane. In this study, we used the rodent model parasite Plasmodium berghei to design a label-free quantitative proteomic approach aimed at identifying gender-related proteins differentially released/secreted by purified mature gametocytes when activated to form gametes. We compared the abundance of molecules secreted by wild type gametocytes of both genders with that of a transgenic line defective in male gamete maturation and egress. This enabled us to provide a comprehensive data set of egress-related molecules and their gender specificity. Using specific antibodies, we validated eleven candidate molecules, predicted as either gender-specific or common to both male and female gametocytes. All of them localize to punctuate, vesicle-like structures that relocate to cell periphery upon activation, but only three of them localize to the gametocyte-specific secretory vesicles named osmiophilic bodies. Our results confirm that the egress process involves a tightly coordinated secretory apparatus that includes different types of vesicles and may put the basis for functional studies aimed at designing novel transmission-blocking molecules.

Mallory-Denk-bodies (MDBs) are hepatic protein aggregates associated with inflammation both clinically and in MDB-inducing models. Similar protein aggregation in neurodegenerative diseases also triggers inflammation and NF-B activation. However, the precise mechanism that links protein aggregation to NF-B-activation and inflammatory response remains unclear. Herein we find that treating primary hepatocytes with MDB-inducing agents (N-methylprotoporphyrin (NMPP), protoporphyrin IX (PPIX), or Zinc-protoporphyrin IX (ZnPP)) elicited an IBα-loss with consequent NF-B activation. Four known mechanisms of IBα-loss i.e. the canonical ubiquitin-dependent proteasomal degradation (UPD), autophagic-lysosomal degradation, calpain degradation and translational inhibition, were all probed and excluded. Immunofluorescence analyses of ZnPP-treated cells coupled with 8 M urea/CHAPS-extraction revealed that this IBα-loss was due to its sequestration along with IBβ into insoluble aggregates, thereby releasing NF-B. Through affinity pulldown, proximity biotinylation by antibody recognition, and other proteomic analyses, we verified that NF-B subunit p65, which stably interacts with IBα under normal conditions, no longer binds to it upon ZnPP-treatment. Additionally, we identified 10 proteins that interact with IBα under baseline conditions, aggregate upon ZnPP-treatment, and maintain the interaction with IBα after ZnPP-treatment, either by cosequestering into insoluble aggregates or through a different mechanism. Of these 10 proteins, the nucleoporins Nup153 and Nup358/RanBP2 were identified through RNA-interference, as mediators of IBα-nuclear import. The concurrent aggregation of IBα, NUP153, and RanBP2 upon ZnPP-treatment, synergistically precluded the nuclear entry of IBα and its consequent binding and termination of NF-B activation. This novel mechanism may account for the protein aggregate-induced inflammation observed in liver diseases, thus identifying novel targets for therapeutic intervention. Because of inherent commonalities this MDB cell model is a bona fide protoporphyric model, making these findings equally relevant to the liver inflammation associated with clinical protoporphyria.

At neuronal synapses, activation of group I metabotropic glutamate receptors (mGluR1/5) triggers a form of long-term depression (mGluR-LTD) that relies on new protein synthesis and the internalization of AMPA-type glutamate receptors. Dysregulation of these processes has been implicated in the development of mental disorders such as autism spectrum disorders and therefore merit a better understanding on a molecular level. Here, to study mGluR-induced signaling pathways, we integrated quantitative phosphoproteomics with the analyses of newly synthesized proteins via bio-orthogonal amino acids (azidohomoalanine) in a pulsed labeling strategy in cultured hippocampal neurons stimulated with DHPG, a specific agonist for group I mGluRs. We identified several kinases with important roles in DHPG-induced mGluR activation, which we confirmed using small molecule kinase inhibitors. Furthermore, changes in the AMPA receptor endocytosis pathway in both protein synthesis and protein phosphorylation were identified, whereby Intersectin-1 was validated as a novel player in this pathway. This study revealed several new insights into the molecular pathways downstream of group I mGluR activation in hippocampal neurons, and provides a rich resource for further analyses.

Huanglongbing (HLB) is the most devastating and widespread citrus disease. All commercial citrus varieties are susceptible to the HLB-associated bacterium, Candidatus Liberibacter asiaticus (CLas), which resides in the phloem. The phloem is part of the plant vascular system and is involved in sugar transport. To investigate the plant response to CLas, we enriched for proteins surrounding the phloem in an HLB susceptible sweet orange variety, Washington navel (Citrus sinensis (L) Osbeck). Quantitative proteomics revealed global changes in the citrus proteome after CLas inoculation. Plant metabolism and translation were suppressed, whereas defense-related proteins such as peroxidases, proteases and protease inhibitors were induced in the vasculature. Transcript accumulation and enzymatic activity of plant peroxidases in CLas infected sweet orange varieties under greenhouse and field conditions were assessed. Although peroxidase transcript accumulation was induced in CLas infected sweet orange varieties, peroxidase enzymatic activity varied. Specific serine proteases were up-regulated in Washington navel in the presence of CLas based on quantitative proteomics. Subsequent activity-based protein profiling revealed increased activity of two serine proteases, and reduced activity of one protease in two C. sinensis sweet orange varieties under greenhouse and field conditions. The observations in the current study highlight global reprogramming of the citrus vascular proteome and differential regulation of enzyme classes in response to CLas infection. These results open an avenue for further investigation of diverse responses to HLB across different environmental conditions and citrus genotypes.

Stroke remains a leading cause of death and disability worldwide. Despite continuous advances, the identification of key molecular signatures in the hyper-acute phase of ischemic stroke is still a primary interest for translational research on stroke diagnosis, prognosis, and treatment. Data integration from high-throughput -omics techniques has become crucial to unraveling key interactions among different molecular elements in complex biological contexts, such as ischemic stroke. Thus, we used advanced data integration methods for a multi-level joint analysis of transcriptomics and proteomics data sets obtained from mouse brains at 2 h after cerebral ischemia. By modeling net-like correlation structures, we identified an integrated network of genes and proteins that are differentially expressed at a very early stage after stroke. We validated 10 of these deregulated elements in acute stroke, and changes in their expression pattern over time after cerebral ischemia were described. Of these, CLDN20, GADD45G, RGS2, BAG5, and CTNND2 were next evaluated as blood biomarkers of cerebral ischemia in mice and human blood samples, which were obtained from stroke patients and patients presenting stroke-mimicking conditions. Our findings indicate that CTNND2 levels in blood might potentially be useful for distinguishing ischemic strokes from stroke-mimicking conditions in the hyper-acute phase of the disease. Furthermore, circulating GADD45G content within the first 6 h after stroke could also play a key role in predicting poor outcomes in stroke patients. For the first time, we have used an integrative biostatistical approach to elucidate key molecules in the initial stages of stroke pathophysiology and highlight new notable molecules that might be further considered as blood biomarkers of ischemic stroke.

This article has been withdrawn by the authors. We discovered an error after this manuscript was published as a Paper in Press. Specifically, we learned that the structures of glycans presented for the PD-L1 peptide were drawn and labeled incorrectly. We wish to withdraw this article and submit a corrected version for review.

The histidine kinase Hik33 plays important roles in mediating cyanobacterial response to divergent types of abiotic stresses including cold, salt, high light (HL), and osmotic stresses. However, how these functions are regulated by Hik33 remains to be addressed. Using a hik33-deficient strain (hik33) of Synechocystis sp. PCC 6803 (Synechocystis) and quantitative proteomics, we found that Hik33 depletion induces differential protein expression highly similar to that induced by divergent types of stresses. This typically includes downregulation of proteins in photosynthesis and carbon assimilation that are necessary for cell propagation, and upregulation of heat shock proteins, chaperons, and proteases that are important for cell survival. This observation indicates that depletion of Hik33 alone mimics divergent types of abiotic stresses, and that Hik33 could be important for preventing abnormal stress response in the normal condition. Moreover, we found the majority of proteins of plasmid origin were significantly upregulated in hik33, though their biological significance remains to be addressed. Together, the systematically characterized Hik33-regulated cyanobacterial proteome, which is largely involved in stress responses, builds the molecular basis for Hik33 as a general regulator of stress responses.

This article has been withdrawn by the authors. This article did not comply with the editorial guidelines of MCP. Specifically, single peptide based protein identifications of 9-19% were included in the analysis and discussed in the results and conclusions. We wish to withdraw this article and resubmit a clarified, corrected manuscript for review.

Glycopeptides in peptide or digested protein samples pose a number of analytical and bioinformatics challenges beyond those posed by unmodified peptides or peptides with smaller posttranslational modifications. Exact structural elucidation of glycans is generally beyond the capability of a single mass spectrometry experiment, so a reasonable level of identification for tandem mass spectrometry, taken by several glycopeptide software tools, is that of peptide sequence and glycan composition, meaning the number of monosaccharides of each distinct mass, for example HexNAc(2)Hex(5) rather than man5. Even at this level, however, glycopeptide analysis poses challenges:  finding glycopeptide spectra when they are a tiny fraction of the total spectra; assigning spectra with unanticipated glycans, not in the initial glycan database; and finding, scoring, and labeling diagnostic peaks in tandem mass spectra.  Here we discuss recent improvements to Byonic, a glycoproteomics search program, that address these three issues. Byonic now supports filtering spectra by m/z peaks, so that the user can limit attention to spectra with diagnostic peaks, for example, at least two out of three of 204.087 for HexNAc, 274.092 for NeuAc (with water loss), and 366.139 for HexNAc-Hex, all within a set mass tolerance, for example, ± 0.01 Daltons. Also new is glycan "wildcard" search, which allows an unspecified mass within a user-set mass range to be applied to N- or O-linked glycans and enables assignment of spectra with unanticipated glycans. Finally the next release of Byonic supports user-specified peak annotations from user-defined posttranslational modifications. We demonstrate the utility of these new software features by finding previously unrecognized glycopeptides in publicly available data, including glycosylated neuropeptides from rat brain.

Sparkling wine is an alcoholic beverage enjoyed around the world. The sensory properties of sparkling wine depend on a complex interplay between the chemical and biochemical components in the final product. Glycoproteins have been linked to positive and negative qualities in sparkling wine, but the glycosylation profiles of sparkling wine have not been previously investigated in detail. We analysed the glyco/proteome of sparkling wines using protein- and glycopeptide-centric approaches. We developed an automated workflow that created ion libraries to analyse Sequential Window Acquisition of all THeoretical mass spectra (SWATH) Data Independent Acquisition (DIA) mass spectrometry data based on glycopeptides identified by Byonic. We applied our workflow to three pairs of experimental sparkling wines to assess the effects of aging on lees and of different yeast strains used in the Liqueur de Tirage for secondary fermentation. We found that aging a cuvée on lees for 24 months compared to 8 months led to a dramatic decrease in overall protein abundance and an enrichment in large glycans at specific sites in some proteins. Secondary fermentation of a Riesling wine with Saccharomyces cerevisiae yeast strain Siha4 produced more yeast proteins and glycoproteins than with S. cerevisiae yeast strain DV10. The abundance and glycosylation profiles of grape glycoproteins were also different between grape varieties. This work represents the first in-depth study into protein- and peptide-specific glycosylation in sparkling wines and describes a quantitative glycoproteomic SWATH/DIA workflow that is broadly applicable to other sample types.

Many cell surface and secreted proteins are modified by the covalent addition of glycans that play an important role in the development of multicellular organisms. These glycan modifications enable communication between cells and the extracellular matrix via interactions with specific glycan-binding lectins and the regulation of receptor-mediated signaling. Aberrant protein glycosylation has been associated with the development of several muscular diseases suggesting essential glycan- and lectin-mediated functions in myogenesis and muscle development but our molecular understanding of the precise glycans, catalytic enzymes and lectins involved remain only partially understood. Here, we quantified dynamic remodeling of the membrane-associated proteome during a time-course of myogenesis in cell culture. We observed wide-spread changes in the abundance of several important lectins and enzymes facilitating glycan biosynthesis. Glycomics-based quantification of released N-linked glycans confirmed remodeling of the glycome consistent with the regulation of glycosyltransferases and glycosidases responsible for their formation including a previously unknown di-galactose-to-sialic acid switch supporting a functional role of these glycoepitopes in myogenesis. Furthermore, dynamic quantitative glycoproteomic analysis with multiplexed stable isotope labelling and analysis of enriched glycopeptides with multiple fragmentation approaches identified glycoproteins modified by these regulated glycans including several integrins and growth factor receptors. Myogenesis was also associated with the regulation of several lectins most notably the up-regulation of galectin-1 (LGALS1). CRISPR/Cas9-mediated deletion of Lgals1 inhibited differentiation and myotube formation suggesting an early functional role of galectin-1 in the myogenic program. Importantly, similar changes in N-glycosylation and the up-regulation of galectin-1 during postnatal skeletal muscle development were observed in mice. Treatment of new-born mice with recombinant adeno-associated viruses to overexpress galectin-1 in the musculature resulted in enhanced muscle mass. Our data form a valuable resource to further understand the glycobiology of myogenesis and will aid the development of intervention strategies to promote healthy muscle development or regeneration.

CD4+ T cell responses are crucial for inducing and maintaining effective anti-cancer immunity, and the identification of human leukocyte antigen class II (HLA-II) cancer-specific epitopes is key to the development of potent cancer immunotherapies. In many tumor types, and especially in glioblastoma (GBM), HLA-II complexes are hardly ever naturally expressed. Hence, little is known about immunogenic HLA-II epitopes in GBM. With stable expression of the class II major histocompatibility complex transactivator (CIITA) coupled to a detailed and sensitive mass spectrometry based immunopeptidomics analysis, we here uncovered a remarkable breadth of the HLA-ligandome in HROG02, HROG17 and RA GBM cell lines. The effect of CIITA expression on the induction of the HLA-II presentation machinery was striking in each of the three cell lines, and it was significantly higher compared to interferon gamma (IFN) treatment. In total, we identified 16,123 unique HLA-I peptides and 32,690 unique HLA-II peptides. In order to genuinely define the identified peptides as true HLA ligands, we carefully characterized their association with the different HLA allotypes. In addition, we identified 138 and 279 HLA-I and HLA-II ligands, respectively, most of which are novel in GBM, derived from known GBM-associated tumor-antigens that have been used as source proteins for a variety of GBM vaccines. Our data further indicate that CIITA-expressing GBM cells acquired an antigen presenting cell-like phenotype as we found that they directly present external proteins as HLA-II ligands. Not only that CIITA-expressing GBM cells are attractive models for antigen discovery endeavors, but also such engineered cells have great therapeutic potential through massive presentation of a diverse antigenic repertoire.

The glycoprotein spike (S) on the surface of SARS-CoV-2 is a determinant for viral invasion and host immune response. Herein, we characterized the site-specific N-glycosylation of S protein at the level of intact glycopeptides. All 22 potential N-glycosites were identified in the S-protein protomer and were found to be preserved among the 753 SARS-CoV-2 genome sequences. The glycosites exhibited glycoform heterogeneity as expected for a human cell-expressed protein subunit. We identified masses that correspond to 157 N-glycans, primarily of the complex type. In contrast, the insect cell-expressed S protein contained 38 N-glycans, completely of the high-mannose type. Our results revealed that the glycan types were highly determined by the differential processing of N-glycans among human and insect cells, regardless of the glycosites’ location. Moreover, the N-glycan compositions were conserved among different sizes of subunits. Our study indicate that the S protein N-glycosylation occurs regularly at each site, albeit the occupied N-glycans were diverse and heterogenous. This N-glycosylation landscape and the differential N-glycan patterns among distinct host cells are expected to shed light on the infection mechanism and present a positive view for the development of vaccines and targeted drugs.

Esophageal squamous cell cancer (ESCC) is an aggressive malignancy with poor therapeutic outcomes. However, the alterations in proteins and post-translational modifications (PTMs) leading to the pathogenesis of ESCC remains unclear. Here, we provide the comprehensive characterization of the proteome, phosphorylome, lysine acetylome and succinylome for ESCC and matched control cells using quantitative proteomic approach. We identify abnormal protein and post-translational modification (PTM) pathways, including significantly downregulated lysine succinylation sites in cancer cells. Focusing on hyposuccinylation, we reveal that this altered PTM was enriched on enzymes of metabolic pathways inextricably linked with cancer metabolism. Importantly, ESCC malignant behaviors such as cell migration are inhibited once the level of succinylation was restored in vitro or in vivo. This effect was further verified by mutations to disrupt succinylation sites in candidate proteins. Meanwhile, we found that succinylation has a negative regulatory effect on histone methylation to promote cancer migration. Finally, hyposuccinylation is confirmed in primary ESCC specimens. Our findings together demonstrate that lysine succinylation may alter ESCC metabolism and migration, providing new insights into the functional significance of PTM in cancer biology.

The choice for adjuvant chemotherapy in stage II colorectal cancer (CRC) is controversial as many patients are cured by surgery alone and it is difficult to identify patients with high-risk of recurrence of the disease. There is a need for better stratification of this group of patients. Mass spectrometry imaging could identify patients at risk. We report here the N-glycosylation signatures of the different cell populations in a group of stage II CRC tissue samples. The cancer cells, compared to normal epithelial cells, have increased levels of sialylation and high-mannose glycans, as well as decreased levels of fucosylation and highly branched N-glycans. When looking at the interface between cancer and its microenvironment, it seems that the cancer N-glycosylation signature spreads into the surrounding stroma at the invasive front of the tumor. This finding was more outspoken in patients with a worse outcome within this sample group.

Paramphistomosis, caused by the rumen fluke, Calicophoron daubneyi, is a parasitic infection of ruminant livestock which has seen a rapid rise in prevalence throughout Western Europe in recent years. Following ingestion of metacercariae (parasite cysts) by the mammalian host, newly-excysted juveniles (NEJs) emerge and invade the duodenal submucosa which causes significant pathology in heavy infections. The immature larvae then migrate upwards, along the gastrointestinal tract, and enter the rumen where they mature and begin to produce eggs. Despite their emergence, and sporadic outbreaks of acute disease, we know little about the molecular mechanisms used by C. daubneyi to establish infection, acquire nutrients and to avoid the host immune response. Here, transcriptome analysis of four intra-mammalian life-cycle stages, integrated with secretome analysis of the NEJ and adult parasites (responsible for acute and chronic disease respectively), revealed how the expression and secretion of selected families of virulence factors and immunomodulators are regulated in accordance with fluke development and migration. Our data show that whilst a family of cathepsins B with varying S2 sub-site residues (indicating distinct substrate specificities) are differentially secreted by NEJs and adult flukes, cathepsins L and F are secreted in low abundance by NEJs only. We found that C. daubneyi has an expanded family of aspartic peptidases, which is up-regulated in adult worms, although they are underrepresented in the secretome. The most abundant proteins in adult fluke secretions were helminth defence molecules (HDMs) that likely establish an immune environment permissive to fluke survival and/or neutralise pathogen-associated molecular patterns (PAMPs) such as bacterial lipopolysaccharide in the microbiome-rich rumen. The distinct collection of molecules secreted by C. daubneyi allowed the development of the first coproantigen-based ELISA for paramphistomosis which, importantly, did not recognise antigens from other helminths commonly found as co-infections with rumen fluke.

Antibodies play essential roles in both diagnostics and therapeutics. Epitope mapping is essential to understand how an antibody works and to protect intellectual property. Given the millions of antibodies for which epitope information is lacking, there is a need for high-throughput epitope mapping. To address this, we developed a strategy, Antibody binding epitope Mapping (AbMap), by combining a phage displayed peptide library with next generation sequencing. Using AbMap, profiles of the peptides bound by 202 antibodies were determined in a single test, and linear epitopes were identified for >50% of the antibodies. Using spike protein (S1 and S2)-enriched antibodies from the convalescent serum of one COVID-19 patient as the input, both linear and conformational epitopes of spike protein specific antibodies were identified. We defined peptide-binding profile of an antibody as the Binding Capacity (BiC). Conceptually, the BiC could serve as a systematic and functional descriptor of any antibody. Requiring at least one order of magnitude less time and money to map linear epitopes than traditional technologies, AbMap allows for high-throughput epitope mapping and creates many possibilities.

Thermal proteome profiling (TPP) allows for the unbiased detection of drug – target protein engagements in vivo. Traditionally, one cell type is used for TPP studies, with the risk of missing important differentially expressed target proteins. The use of whole organisms would circumvent this problem. Zebrafish embryos are amenable to such an approach. Here, we used TPP on whole zebrafish embryo lysate to identify protein targets of napabucasin, a compound that may affect Signal transducer and activator of transcription 3 (Stat3) signaling through an ill-understood mechanism. In zebrafish embryos, napabucasin induced developmental defects consistent with inhibition of Stat3 signaling. TPP profiling showed no distinct shift in Stat3 upon napabucasin treatment, but effects were detected on the oxidoreductase, Pora, which might explain effects on Stat3 signaling. Interestingly, thermal stability of several aldehyde dehydrogenases (Aldhs) was affected. Moreover, napabucasin activated ALDH enzymatic activity in vitro. Aldhs have crucial roles in retinoic acid metabolism and functionally we validated napabucasin-mediated activation of the retinoic acid pathway in zebrafish in vivo. We conclude that TPP profiling in whole zebrafish embryo lysate is feasible and facilitates direct correlation of in vivo effects of small molecule drugs with their protein targets.

The complexity and dynamics of the immensely heterogeneous glycoproteome of the prostate cancer (PCa) tumour micro-environment remain incompletely mapped, a knowledge gap that impedes our molecular-level understanding of the disease. To this end, we have used sensitive glycomics and glycoproteomics to map the protein-, cell- and tumour grade-specific N- and O-glycosylation in surgically-removed PCa tissues spanning five histological grades (n = 10/grade) and tissues from patients with benign prostatic hyperplasia (n = 5). Quantitative glycomics revealed PCa grade-specific alterations of the oligomannosidic-, paucimannosidic- and branched sialylated complex-type N-glycans, and dynamic remodelling of the sialylated core 1- and core 2-type O-glycome. Deep quantitative glycoproteomics identified ~7,400 unique N-glycopeptides from 500 N-glycoproteins and ~500 unique O-glycopeptides from nearly 200 O-glycoproteins. With reference to a recent Tissue and Blood Atlas, our data indicate that paucimannosidic glycans of the PCa tissues arise mainly from immune cell-derived glycoproteins. Further, the grade-specific PCa glycosylation arises primarily from dynamics in the cellular makeup of the PCa tumour microenvironment across grades involving increased oligomannosylation of prostate-derived glycoproteins and decreased bisecting GlcNAcylation of N-glycans carried by the extracellular matrix proteins. Further, elevated expression of several oligosaccharyltransferase subunits and enhanced N-glycoprotein site occupancy were observed associated with PCa progression. Finally, correlations between the protein-specific glycosylation and PCa progression were observed including increased site-specific core 2-type O-glycosylation of collagen VI. In conclusion, integrated glycomics and glycoproteomics have enabled new insight into the complexity and dynamics of the tissue glycoproteome associated with PCa progression generating an important resource to explore the underpinning disease mechanisms.

Malaria elimination is still pending on the development of novel tools that rely on a deep understanding of parasite biology. Proteins of all living cells undergo a myriad number of posttranslational modifications (PTMs) that are critical to multifarious life processes. An extensive proteome-wide dissection revealed a fine PTM map of most proteins in both Plasmodium falciparum, the causative agent of severe malaria, and the infected red blood cells. More than two-thirds of proteins of the parasite and its host cell underwent extensive and dynamic modification throughout the erythrocytic developmental stage. PTMs critically modulate the virulence factors involved in the host-parasite interaction and pathogenesis. Furthermore, P. falciparum stabilized the supporting proteins of erythrocyte origin by selective de-modification. Collectively, our multiple omic analyses, apart from having furthered a deep understanding of the systems biology of P. falciparum and malaria pathogenesis, provide a valuable resource for mining new antimalarial targets.

Biological functions emerge from complex and dynamic networks of protein-protein interactions. Because these protein-protein interaction networks, or interactomes, represent pairwise connections within a hierarchically organized system, it is often useful to identify higher-order associations embedded within them, such as multi-member protein complexes. Graph-based clustering techniques are widely used to accomplish this goal, and dozens of field-specific and general clustering algorithms exist. However, interactomes can be prone to errors, especially when inferred from high-throughput biochemical assays. Therefore, robustness to network-level noise is an important criterion for any clustering algorithm that aims to generate robust, reproducible clusters. Here, we tested the robustness of a range of graph-based clustering algorithms in the presence of noise, including algorithms common across domains and those specific to protein networks. Strikingly, we found that all of the clustering algorithms tested here markedly amplified noise within the underlying protein interaction network. Randomly rewiring only 1% of network edges yielded more than a 50% change in clustering results, indicating that clustering markedly amplified network-level noise. Moreover, we found the impact of network noise on individual clusters was not uniform: some clusters were consistently robust to injected noise while others were not. To assist in assessing this, we developed the clust.perturb R package and Shiny web application to measure the reproducibility of clusters by randomly perturbing the network. We show that clust.perturb results are predictive of real-world cluster stability: poorly reproducible clusters as identified by clust.perturb are significantly less likely to be reclustered across experiments. We conclude that graph-based clustering amplifies noise in protein interaction networks, but quantifying the robustness of a cluster to network noise can separate stable protein complexes from spurious associations.

Protease activity has been associated with pathological processes that can lead to cancer development and progression. However, understanding the pathological unbalance in proteolysis is challenging since changes can occur simultaneously at protease, their inhibitor and substrate levels. Here, we present a pipeline that combines peptidomics, proteomics and peptidase predictions for studying proteolytic events in the saliva of seventy-nine patients and their association with oral squamous cell carcinoma (OSCC) prognosis. Our findings revealed differences in the saliva peptidome of patients with (pN+) or without (pN0) lymph node metastasis and delivered a panel of ten endogenous peptides correlated with poor prognostic factors plus five molecules able to classify pN0 and pN+ patients (ROC-AUC>0.85). In addition, endo- and exopeptidases putatively implicated in the processing of differential peptides were investigated using cancer tissue gene expression data from publicly repositories reinforcing their association with poorer survival rates and prognosis in oral cancer. The dynamics of the OSCC-related proteolysis was further explored via the proteomic profiling of saliva. This revealed that peptidase/endopeptidase inhibitors exhibited reduced levels in the saliva of pN+ patients, as confirmed by SRM-MS, whilst minor changes were detected in the level of saliva proteases. Taken together, our results indicated that proteolytic activity is accentuated in the saliva of OSCC patients with lymph node metastasis and, at least in part, this is modulated by reduced levels of salivary peptidase inhibitors. Therefore, this integrated pipeline provided better comprehension and discovery of molecular features with implications in the oral cancer metastasis prognosis.

Modulation of the host cell is integral to the survival and replication of microbial pathogens. Several intracellular bacterial pathogens deliver bacterial proteins, termed ‘effector proteins’ into the host cell during infection by sophisticated protein translocation systems, which manipulate cellular processes and functions. The functional contribution of individual effectors is poorly characterised, particularly in intracellular bacterial pathogens with large effector protein repertoires. Technical caveats have limited the capacity to study these proteins during a native infection, with many effector proteins having only been demonstrated to be translocated during over-expression of tagged versions. Here we developed a novel strategy to examine effector proteins in the context of infection. We coupled a broad, unbiased proteomics-based screen with organelle purification to study the host-pathogen interactions occurring between the host cell mitochondrion and the Gram-negative, Q fever pathogen Coxiella burnetii. We identify 4 novel mitochondrially-targeted C. burnetii effector proteins, renamed Mitochondrial Coxiella effector protein (Mce) B to E. Examination of the subcellular localisation of ectopically expressed proteins confirmed their mitochondrial localisation, demonstrating the robustness of our approach. Subsequent biochemical analysis and affinity enrichment proteomics of one of these effector proteins, MceC, revealed the protein localises to the inner membrane and can interact with components of the mitochondrial quality control machinery. Our study adapts high-sensitivity proteomics to study intracellular host-pathogen interactions, providing a robust strategy to examine the sub-cellular localisation of effector proteins during native infection. This approach could be applied to a range of pathogens and host cell compartments to provide a rich map of effector dynamics throughout infection.

Hirschsprung disease (HSCR) is a heterogeneous group of neurocristopathy characterized by the absence of the enteric ganglia along a variable length of the intestine. Genetic defects play a major role in the pathogenesis of HSCR while family studies of pathogenic variants in all the known genes (loci) only demonstrate incomplete penetrance and variable expressivity for unknown reasons. Here, we applied large-scale, quantitative proteomics of human colon tissues from 21 patients using iTRAQ method followed by bioinformatics analysis. Selected findings were confirmed by parallel reaction monitoring (PRM) verification. At last the interesting differentially expressed proteins were confirmed by western blot. A total of 5341 proteins in human colon tissues were identified. Among them, 664 proteins with >1.2-fold difference were identified in 6 groups: groups A1 and A2 pooled protein from the ganglionic and aganglionic colon of male, long-segment HSCR patients (L-HSCR, n=7); groups B1 and B2 pooled protein from the ganglionic and aganglionic colon of male, short-segment HSCR patients (S-HSCR, n=7); and groups C1 and C2 pooled protein from the ganglionic and aganglionic colon of female, S-HSCR patients (n=7). Based on these analyses, 49 proteins from 5 pathways were selected for PRM verification, including ribosome, endocytosis, spliceosome, oxidative phosphorylation and cell adhesion. The downregulation of three neuron projection development genes ARF4, KIF5B and RAB8A in the aganglionic part of the colon were verified in 15 paired colon samples using WB. The findings of this study will shed new light on the pathogenesis of HSCR and facilitate the development of therapeutic targets.

Thyroglobulin (Tg) is a secreted iodoglycoprotein serving as the precursor for T3 and T4 hormones. Many characterized Tg gene mutations produce secretion-defective variants resulting in congenital hypothyroidism (CH). Tg processing and secretion is controlled by extensive interactions with chaperone, trafficking, and degradation factors comprising the secretory proteostasis network. While dependencies on individual proteostasis network components are known, the integration of proteostasis pathways mediating Tg protein quality control and the molecular basis of mutant Tg misprocessing remain poorly understood. We employ a multiplexed quantitative affinity purification–mass spectrometry approach to define the Tg proteostasis interactome and changes between WT and several CH-variants. Mutant Tg processing is associated with common imbalances in proteostasis engagement including increased chaperoning, oxidative folding, and engagement by targeting factors for ER-associated degradation (ERAD). Furthermore, we reveal mutation-specific changes in engagement with N-glycosylation components, suggesting distinct requirements for one Tg variant on dual engagement of both oligosaccharyltransferase complex isoforms for degradation. Modulating dysregulated proteostasis components and pathways may serve as a therapeutic strategy to restore Tg secretion and thyroid hormone biosynthesis.

The Rhizobium-legume symbiosis is a beneficial interaction in which the bacterium converts atmospheric nitrogen into ammonia and delivers it to the plant in exchange for carbon compounds. This symbiosis implies the adaptation of bacteria to live inside host plant cells. In this work we apply RP-LC-MS/MS and  iTRAQ techniques to study the proteomic profile of endosymbiotic cells (bacteroids) induced by Rhizobium leguminosarum bv viciae strain UPM791 in legume nodules. Nitrogenase subunits, tricarboxylic acid cycle enzymes, and stress response proteins are amongst the most abundant from over one thousand rhizobial proteins identified in pea (Pisum sativum) bacteroids. Comparative analysis of bacteroids induced in pea and in lentil (Lens culinaris)nodules revealed the existence of a significant host-specific differential response affecting dozens of bacterial proteins, including stress-related proteins, transcriptional regulators, and proteins involved in the carbon and nitrogen metabolisms. A mutant affected in one of these proteins, homologous to a GntR-like transcriptional regulator, showed a symbiotic performance significantly  impaired in symbiosis with pea, but not with lentil plants. Analysis of the proteomes of bacteroids isolated from both hosts also revealed the presence of different sets of plant-derived nodule-specific cysteine rich (NCR) peptides, indicating that the endosymbiotic bacteria find a host-specific cocktail of chemical stressors inside the nodule. By studying variations of the bacterial response to different plant cell environments we will be able to identify specific limitations imposed by the host that might give us clues for the improvement of rhizobial performance.

The early detection of pancreatic ductal adenocarcinoma is a complex clinical obstacle yet is key to improving the overall likelihood of patient survival. Current and prospective carbohydrate biomarkers CA19-9 and sTRA are sufficient for surveilling disease progression yet are not approved for delineating PDAC from other abdominal cancers and non-cancerous pancreatic pathologies. To further understand these glycan epitopes, an imaging mass spectrometry approach was utilized to assess the N-glycome of the human pancreas and pancreatic cancer in a cohort of PDAC patients represented by tissue microarrays and whole tissue sections. Orthogonally, these same tissues were characterized by multi-round immunofluorescence which defined expression of CA19-9 and sTRA as well as other lectins towards carbohydrate epitopes with the potential to improve PDAC diagnosis. These analyses revealed distinct differences not only in N-glycan spatial localization across both healthy and diseased tissues but importantly between different biomarker-categorized tissue samples. Unique sulfated bi-antennary N-glycans were detected specifically in normal pancreatic islets. N-glycans from CA19-9 expressing tissues tended to be bi-, tri- and tetra-antennary structures with both core and terminal fucose residues and bisecting N-acetylglucosamines. These N-glycans were detected in less abundance in sTRA-expressing tumor tissues, which favored tri- and tetra-antennary structures with polylactosamine extensions. Increased sialylation of N-glycans was detected in all tumor tissues. A candidate new biomarker derived from IMS was further explored by fluorescence staining with selected lectins on the same tissues. The lectins confirmed the expression of the epitopes in cancer cells and revealed different tumor-associated staining patterns between glycans with bisecting GlcNAc and those with terminal GlcNAc. Thus, the combination of lectin-IHC and IMS techniques produces more complete information for tumor classification than the individual analyses alone. These findings potentiate the development of early assessment technologies to rapidly and specifically identify PDAC in the clinic that may directly impact patient outcomes.

Aspergillus flavus (A. flavus), a pathogenic fungus, can produce carcinogenic and toxic aflatoxins that are a serious agricultural and medical threat worldwide. Attempts to decipher the aflatoxin biosynthetic pathway have been hampered by the lack of a high-quality genome annotation for A. flavus. To address this gap, we performed a comprehensive proteogenomic analysis using high-accuracy mass spectrometry data for this pathogen. The resulting high-quality dataset confirmed the translation of 8,724 previously-predicted genes, and identified 732 novel proteins, 269 splice variants, 447 single amino acid variants, 188 revised genes. A subset of novel proteins was experimentally validated by RT-PCR and synthetic peptides. Further functional annotation suggested that a number of the identified novel proteins may play roles in aflatoxin biosynthesis and stress responses in A. flavus. This comprehensive strategy also identified a wide range of post-translational modifications (PTMs), including 3,461 modification sites from 1,765 proteins. Functional analysis suggested the involvement of these modified proteins in the regulation of cellular metabolic and aflatoxin biosynthetic pathways. Together, we provided a high quality annotation of A. flavus genome and revealed novel insights into the mechanisms of aflatoxin production and pathogenicity in this pathogen.

The molecular mechanism associated with mammalian meiosis has yet to be fully explored, and one of the main reasons for this lack of exploration is that some meiosis-essential genes are still unknown. The profiling of gene expression during spermatogenesis has been performed in previous studies, yet few studies have aimed to find new functional genes. Since there is a huge gap between the number of genes that are able to be quantified and the number of genes that can be characterized by phenotype screening in one assay, an efficient method to rank quantified genes according to phenotypic relevance is of great importance. We proposed to rank genes by the probability of their function in mammalian meiosis based on global protein abundance using machine learning. Here, nine types of germ cells focusing on continual substages of meiosis prophase I were isolated, and the corresponding proteomes were quantified by high-resolution mass spectrometry. By combining meiotic labels annotated from the MGI mouse knockout database and the spermatogenesis proteomics dataset, a supervised machine learning package, FuncProFinder, was developed to rank meiosis-essential candidates. Of the candidates whose functions were unannotated, four of ten genes with the top prediction scores, Zcwpw1, Tesmin, 1700102P08Rik and Kctd19, were validated as meiosis-essential genes by knockout mouse models. Therefore,  mammalian meiosis-essential genes could be efficiently predicted based on the protein abundance dataset, which provides a paradigm for other functional gene mining from a related abundance dataset.

The histopathological subtype of lung adenocarcinoma (LUAD) is closely associated with prognosis. Micropapillary or solid predominant LUAD tends to relapse after surgery at an early stage, whereas lepidic pattern shows a favorable outcome. However, the molecular mechanism underlying this phenomenon remains unknown. Here, we recruited 31 lepidic predominant LUADs (LR: low-risk subtype group) and 28 micropapillary or solid predominant LUADs (HR: high-risk subtype group). Tissues of these cases were obtained and label-free quantitative proteomic and bioinformatic analyses were performed. Additionally, prognostic impact of targeted proteins was validated using The Cancer Genome Atlas databases (n=492) and tissue microarrays composed of early-stage LUADs (n=228). A total of 192 differentially expressed proteins were identified between tumor tissues of LR and HR and three clusters were identified via hierarchical clustering excluding eight proteins. Cluster 1 (65 proteins) showed a sequential decrease in expression from normal tissues to tumor tissues of LR and then to HR and was predominantly enriched in pathways such as tyrosine metabolism and ECM-receptor interaction, and increased matched mRNA expression of 18 proteins from this cluster predicted favorable prognosis. Cluster 2 (70 proteins) demonstrated a sequential increase in expression from normal tissues to tumor tissues of LR and then to HR and was mainly enriched in pathways such as extracellular organization, DNA replication and cell cycle, and high matched mRNA expression of 25 proteins indicated poor prognosis. Cluster 3 (49 proteins) showed high expression only in LR, with high matched mRNA expression of 20 proteins in this cluster indicating favorable prognosis. Furthermore, high expression of ERO1A and FEN1 at protein level predicted poor prognosis in early-stage LUAD, supporting the mRNA results. In conclusion, we discovered key differentially expressed proteins and pathways between low-risk and high-risk subtypes of early-stage LUAD. Some of these proteins could serve as potential biomarkers in prognostic evaluation.