Improved Proteomics-Based Drug Mechanism-of-Action Studies Using 16-Plex Isobaric Mass Tags.

Multiplexed quantitative proteomics enabled complex workflows to study the mechanisms by which small molecule drugs interact with the proteome such as thermal proteome profiling (TPP) or multiplexed proteome dynamics profiling (mPDP). TPP measures changes in protein thermal stability in response to drug treatment and thus informs on direct targets and downstream regulation events, while the mPDP approach enables the discovery of regulated protein synthesis and degradation events caused by small molecules and other perturbations. The isobaric mass tags available for multiplexed proteomics have thus far limited the efficiency and sensitivity by which such experiments could be performed. Here we evaluate a recent generation of 16-plex isobaric mass tags and demonstrate the sensitive and time efficient identification of Staurosporine targets in HepG2 cell extracts by recording full thermal denaturation/aggregation profiles of vehicle and compound treated samples in a single mass spectrometry experiment. In 2D-TPP experiments, isothermal titration over seven concentrations per temperature enabled comprehensive selectivity profiling of Staurosporine with EC50 values for kinase targets tightly matching to the kinobeads gold standard assay. Finally, we demonstrate time and condition-based multiplexing of dynamic SILAC labeling experiments to delineate proteome-wide effects of the molecular glue Indisulam on synthesis and degradation rates.

Enhancing Open Modification Searches via a Combined Approach Facilitated by Ursgal.

The identification of peptide sequences and their post-translational modifications (PTMs) is a crucial step in the analysis of bottom-up proteomics data. The recent development of open modification search (OMS) engines allows virtually all PTMs to be searched for. This not only increases the number of spectra that can be matched to peptides but also greatly advances the understanding of the biological roles of PTMs through the identification, and the thereby facilitated quantification, of peptidoforms (peptide sequences and their potential PTMs). Whereas the benefits of combining results from multiple protein database search engines have been previously established, similar approaches for OMS results have been missing so far. Here we compare and combine results from three different OMS engines, demonstrating an increase in peptide spectrum matches of 8-18%. The unification of search results furthermore allows for the combined downstream processing of search results, including the mapping to potential PTMs. Finally, we test for the ability of OMS engines to identify glycosylated peptides. The implementation of these engines in the Python framework Ursgal facilitates the straightforward application of the OMS with unified parameters and results files, thereby enabling yet unmatched high-throughput, large-scale data analysis.

SugarPy facilitates the universal, discovery-driven analysis of intact glycopeptides.

MOTIVATION: Protein glycosylation is a complex post-translational modification with crucial cellular functions in all domains of life. Currently, large-scale glycoproteomics approaches rely on glycan database dependent algorithms and are thus unsuitable for discovery-driven analyses of glycoproteomes. RESULTS: Therefore, we devised SugarPy, a glycan database independent Python module, and validated it on the glycoproteome of human breast milk. We further demonstrated its applicability by analyzing glycoproteomes with uncommon glycans stemming from the green alga Chlamydomonas reinhardtii and the archaeon Haloferax volcanii. SugarPy also facilitated the novel characterization of glycoproteins from the red alga Cyanidioschyzon merolae. AVAILABILITY AND IMPLEMENTATION: The source code is freely available on GitHub (https://github.com/SugarPy/SugarPy), and its implementation in Python ensures support for all operating systems. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Archaeal Proteome Project advances knowledge about archaeal cell biology through comprehensive proteomics.

While many aspects of archaeal cell biology remain relatively unexplored, systems biology approaches like mass spectrometry (MS) based proteomics offer an opportunity for rapid advances. Unfortunately, the enormous amount of MS data generated often remains incompletely analyzed due to a lack of sophisticated bioinformatic tools and field-specific biological expertise for data interpretation. Here we present the initiation of the Archaeal Proteome Project (ArcPP), a community-based effort to comprehensively analyze archaeal proteomes. Starting with the model archaeon Haloferax volcanii, we reanalyze MS datasets from various strains and culture conditions. Optimized peptide spectrum matching, with strict control of false discovery rates, facilitates identifying > 72% of the reference proteome, with a median protein sequence coverage of 51%. These analyses, together with expert knowledge in diverse aspects of cell biology, provide meaningful insights into processes such as N-terminal protein maturation, N-glycosylation, and metabolism. Altogether, ArcPP serves as an invaluable blueprint for comprehensive prokaryotic proteomics.

Nano LC-MS using capillary columns enables accurate quantification of modified ribonucleosides at low femtomol levels.

Post-transcriptional chemical modifications of (t)RNA molecules are crucial in fundamental biological processes, such as translation. Despite their biological importance and accumulating evidence linking them to various human diseases, technical challenges have limited their detection and accurate quantification. Here, we present a sensitive capillary nanoflow liquid chromatography mass spectrometry (nLC-MS) pipeline for quantitative high-resolution analysis of ribonucleoside modifications from complex biological samples. We evaluated two porous graphitic carbon (PGC) materials and one end-capped C18 reference material as stationary phases for reversed-phase separation. We found that these matrices have complementing retention and separation characteristics, including the capability to separate structural isomers. PGC and C18 matrices yielded excellent signal-to-noise ratios in nLC-MS while differing in the separation capability and sensitivity for various nucleosides. This emphasizes the need for tailored LC-MS setups for optimally detecting as many nucleoside modifications as possible. Detection ranges spanning up to six orders of magnitude enable the analysis of individual ribonucleosides down to femtomol concentrations. Furthermore, normalizing the obtained signal intensities to a stable isotope labeled spike-in enabled direct comparison of ribonucleoside levels between different samples. In conclusion, capillary columns coupled to nLC-MS constitute a powerful and sensitive tool for quantitative analysis of modified ribonucleosides in complex biological samples. This setup will be invaluable for further unraveling the intriguing and multifaceted biological roles of RNA modifications.

N-Glycoproteomic Characterization of Mannosidase and Xylosyltransferase Mutant Strains of Chlamydomonasreinhardtii.

At present, only little is known about the enzymatic machinery required for N-glycosylation in Chlamydomonas reinhardtii, leading to the formation of N-glycans harboring Xyl and methylated Man. This machinery possesses new enzymatic features, as C. reinhardtii N-glycans are independent of ß1,2-N-acetylglucosaminyltransferase I. Here we have performed comparative N-glycoproteomic analyses of insertional mutants of mannosidase 1A (IM Man1A ) and xylosyltransferase 1A (IM XylT1A ). The disruption of man1A affected methylation of Man and the addition of terminal Xyl. The absence of XylT1A led to shorter N-glycans compared to the wild type. The use of a IM Man1A xIM XylT1A double mutant revealed that the absence of Man1A suppressed the IM XylT1A phenotype, indicating that the increased N-glycan trimming is regulated by core ß1,2-Xyl and is dependent on Man1A activity. These data point toward an enzymatic cascade in the N-glycosylation pathway of C. reinhardtii with interlinked roles of Man1A and XylT1A. The results described herein represent the first step toward a functional characterization of the enzymatic N-glycosylation machinery in C. reinhardtii.

Identification of methylated GnTI-dependent N-glycans in Botryococcus brauni.

In contrast to mammals and vascular plants, microalgae show a high diversity in the N-glycan structures of complex N-glycoproteins. Although homologues for ß1,2-N-acetylglucosaminyltransferase I (GnTI), a key enzyme in the formation of complex N-glycans, have been identified in several algal species, GnTI-dependent N-glycans have not been detected so far. We have performed an N-glycoproteomic analysis of the hydrocarbon oils accumulating green microalgae Botryococcus braunii. Thereby, the analysis of intact N-glycopeptides allowed the determination of N-glycan compositions. Furthermore, insights into the role of N-glycosylation in B. braunii were gained from functional annotation of the identified N-glycoproteins. In total, 517 unique N-glycosylated peptides have been identified, including intact N-glycopeptides that harbored N-acetylhexosamine (HexNAc) at the nonreducing end. Surprisingly, these GnTI-dependent N-glycans were also found to be modified with (di)methylated hexose. The identification of GnTI-dependent N-glycans in combination with N-glycan methylation in B. braunii revealed an uncommon type of N-glycan processing in this microalgae.

pyQms enables universal and accurate quantification of mass spectrometry data.

Quantitative mass spectrometry (MS) is a key technique in many research areas (1), including proteomics, metabolomics, glycomics, and lipidomics. Because all of the corresponding molecules can be described by chemical formulas, universal quantification tools are highly desirable. Here, we present pyQms, an open-source software for accurate quantification of all types of molecules measurable by MS. pyQms uses isotope pattern matching that offers an accurate quality assessment of all quantifications and the ability to directly incorporate mass spectrometer accuracy. pyQms is, due to its universal design, applicable to every research field, labeling strategy, and acquisition technique. This opens ultimate flexibility for researchers to design experiments employing innovative and hitherto unexplored labeling strategies. Importantly, pyQms performs very well to accurately quantify partially labeled proteomes in large scale and high throughput, the most challenging task for a quantification algorithm.

Comparative transcriptome and proteome analysis reveals a global impact of the nitrogen regulators AreA and AreB on secondary metabolism in Fusarium fujikuroi.

The biosynthesis of multiple secondary metabolites in the phytopathogenic ascomycete Fusarium fujikuroi is strongly affected by nitrogen availability. Here, we present the first genome-wide transcriptome and proteome analysis that compared the wild type and deletion mutants of the two major nitrogen regulators AreA and AreB. We show that AreB acts not simply as an antagonist of AreA counteracting the expression of AreA target genes as suggested based on the yeast model. Both GATA transcription factors affect a large and diverse set of common as well as specific target genes and proteins, acting as activators and repressors. We demonstrate that AreA and AreB are not only involved in fungal nitrogen metabolism, but also in the control of several complex cellular processes like carbon metabolism, transport and secondary metabolism. We show that both GATA transcription factors can be considered as master regulators of secondary metabolism as they affect the expression of more than half of the 47 putative secondary metabolite clusters identified in the genome of F. fujikuroi. While AreA acts as a positive regulator of many clusters under nitrogen-limiting conditions, AreB is able to activate and repress gene clusters (e.g. bikaverin) under nitrogen limitation and sufficiency. In addition, ChIP analyses revealed that loss of AreA or AreB causes histone modifications at some of the regulated gene clusters.

Temporal quantitative phosphoproteomics of ADP stimulation reveals novel central nodes in platelet activation and inhibition.

Adenosine diphosphate (ADP) enhances platelet activation by virtually any other stimulant to complete aggregation. It binds specifically to the G-protein-coupled membrane receptors P2Y1 and P2Y12, stimulating intracellular signaling cascades, leading to integrin αIIbß3 activation, a process antagonized by endothelial prostacyclin. P2Y12 inhibitors are among the most successful antiplatelet drugs, however, show remarkable variability in efficacy. We reasoned whether a more detailed molecular understanding of ADP-induced protein phosphorylation could identify (1) critical hubs in platelet signaling toward aggregation and (2) novel molecular targets for antiplatelet treatment strategies. We applied quantitative temporal phosphoproteomics to study ADP-mediated signaling at unprecedented molecular resolution. Furthermore, to mimic the antagonistic efficacy of endothelial-derived prostacyclin, we determined how Iloprost reverses ADP-mediated signaling events. We provide temporal profiles of 4797 phosphopeptides, 608 of which showed significant regulation. Regulated proteins are implicated in well-known activating functions such as degranulation and cytoskeletal reorganization, but also in less well-understood pathways, involving ubiquitin ligases and GTPase exchange factors/GTPase-activating proteins (GEF/GAP). Our data demonstrate that ADP-triggered phosphorylation occurs predominantly within the first 10 seconds, with many short rather than sustained changes. For a set of phosphorylation sites (eg, PDE3ASer312, CALDAG-GEFISer587, ENSASer109), we demonstrate an inverse regulation by ADP and Iloprost, suggesting that these are central modulators of platelet homeostasis. This study demonstrates an extensive spectrum of human platelet protein phosphorylation in response to ADP and Iloprost, which inversely overlap and represent major activating and inhibitory pathways.

Phosphorylation of GENOMES UNCOUPLED 4 Alters Stimulation of Mg Chelatase Activity in Angiosperms.

GENOMES UNCOUPLED 4 (GUN4) is a positive regulator of light-dependent chlorophyll biosynthesis. GUN4 activates Mg chelatase (MgCh) that catalyzes the insertion of an Mg2+ ion into protoporphyrin IX. We show that Arabidopsis (Arabidopsis thaliana) GUN4 is phosphorylated at Ser 264 (S264), the penultimate amino acid residue at the C terminus. While GUN4 is preferentially phosphorylated in darkness, phosphorylation is reduced upon accumulation of Mg porphyrins. Expression of a phosphomimicking GUN4(S264D) results in an incomplete complementation of the white gun4-2 null mutant and a chlorotic phenotype comparable to gun4 knockdown mutants. Phosphorylated GUN4 has a reduced stimulatory effect on MgCh in vitro and in vivo but retains its protein stability and tetrapyrrole binding capacity. Analysis of GUN4 found in oxygenic photosynthetic organisms reveals the evolution of a C-terminal extension, which harbors the phosphorylation site of GUN4 expressed in angiosperms. Homologs of GUN4 from Synechocystis and Chlamydomonas lack the conserved phosphorylation site found in a C-terminal extension of angiosperm GUN4. Biochemical studies proved the importance of the C-terminal extension for MgCh stimulation and inactivation of GUN4 by phosphorylation in angiosperms. An additional mechanism regulating MgCh activity is proposed. In conjunction with the dark repression of 5-aminolevulinic acid synthesis, GUN4 phosphorylation minimizes the flow of intermediates into the Mg branch of the tetrapyrrole metabolic pathway for chlorophyll biosynthesis.

Calredoxin represents a novel type of calcium-dependent sensor-responder connected to redox regulation in the chloroplast.

Calcium (Ca(2+)) and redox signalling play important roles in acclimation processes from archaea to eukaryotic organisms. Herein we characterized a unique protein from Chlamydomonas reinhardtii that has the competence to integrate Ca(2+)- and redox-related signalling. This protein, designated as calredoxin (CRX), combines four Ca(2+)-binding EF-hands and a thioredoxin (TRX) domain. A crystal structure of CRX, at 1.6 Å resolution, revealed an unusual calmodulin-fold of the Ca(2+)-binding EF-hands, which is functionally linked via an inter-domain communication path with the enzymatically active TRX domain. CRX is chloroplast-localized and interacted with a chloroplast 2-Cys peroxiredoxin (PRX1). Ca(2+)-binding to CRX is critical for its TRX activity and for efficient binding and reduction of PRX1. Thereby, CRX represents a new class of Ca(2+)-dependent 'sensor-responder' proteins. Genetically engineered Chlamydomonas strains with strongly diminished amounts of CRX revealed altered photosynthetic electron transfer and were affected in oxidative stress response underpinning a function of CRX in stress acclimation.

Ursgal, Universal Python Module Combining Common Bottom-Up Proteomics Tools for Large-Scale Analysis.

Proteomics data integration has become a broad field with a variety of programs offering innovative algorithms to analyze increasing amounts of data. Unfortunately, this software diversity leads to many problems as soon as the data is analyzed using more than one algorithm for the same task. Although it was shown that the combination of multiple peptide identification algorithms yields more robust results, it is only recently that unified approaches are emerging; however, workflows that, for example, aim to optimize search parameters or that employ cascaded style searches can only be made accessible if data analysis becomes not only unified but also and most importantly scriptable. Here we introduce Ursgal, a Python interface to many commonly used bottom-up proteomics tools and to additional auxiliary programs. Complex workflows can thus be composed using the Python scripting language using a few lines of code. Ursgal is easily extensible, and we have made several database search engines (X!Tandem, OMSSA, MS-GF+, Myrimatch, MS Amanda), statistical postprocessing algorithms (qvality, Percolator), and one algorithm that combines statistically postprocessed outputs from multiple search engines ("combined FDR") accessible as an interface in Python. Furthermore, we have implemented a new algorithm ("combined PEP") that combines multiple search engines employing elements of "combined FDR", PeptideShaker, and Bayes' theorem.

STATE TRANSITION7-Dependent Phosphorylation Is Modulated by Changing Environmental Conditions, and Its Absence Triggers Remodeling of Photosynthetic Protein Complexes.

In plants and algae, the serine/threonine kinase STN7/STT7, orthologous protein kinases in Chlamydomonas reinhardtii and Arabidopsis (Arabidopsis thaliana), respectively, is an important regulator in acclimation to changing light environments. In this work, we assessed STT7-dependent protein phosphorylation under high light in C. reinhardtii, known to fully induce the expression of light-harvesting complex stress-related protein3 (LHCSR3) and a nonphotochemical quenching mechanism, in relationship to anoxia where the activity of cyclic electron flow is stimulated. Our quantitative proteomics data revealed numerous unique STT7 protein substrates and STT7-dependent protein phosphorylation variations that were reliant on the environmental condition. These results indicate that STT7-dependent phosphorylation is modulated by the environment and point to an intricate chloroplast phosphorylation network responding in a highly sensitive and dynamic manner to environmental cues and alterations in kinase function. Functionally, the absence of the STT7 kinase triggered changes in protein expression and photoinhibition of photosystem I (PSI) and resulted in the remodeling of photosynthetic complexes. This remodeling initiated a pronounced association of LHCSR3 with PSI-light harvesting complex I (LHCI)-ferredoxin-NADPH oxidoreductase supercomplexes. Lack of STT7 kinase strongly diminished PSII-LHCII supercomplexes, while PSII core complex phosphorylation and accumulation were significantly enhanced. In conclusion, our study provides strong evidence that the regulation of protein phosphorylation is critical for driving successful acclimation to high light and anoxic growth environments and gives new insights into acclimation strategies to these environmental conditions.

An evolutionary approach uncovers a diverse response of tRNA 2-thiolation to elevated temperatures in yeast.

Chemical modifications of transfer RNA (tRNA) molecules are evolutionarily well conserved and critical for translation and tRNA structure. Little is known how these nucleoside modifications respond to physiological stress. Using mass spectrometry and complementary methods, we defined tRNA modification levels in six yeast species in response to elevated temperatures. We show that 2-thiolation of uridine at position 34 (s(2)U34) is impaired at temperatures exceeding 30°C in the commonly used Saccharomyces cerevisiae laboratory strains S288C and W303, and in Saccharomyces bayanus. Upon stress relief, thiolation levels recover and we find no evidence that modified tRNA or s(2)U34 nucleosides are actively removed. Our results suggest that loss of 2-thiolation follows accumulation of newly synthesized tRNA that lack s(2)U34 modification due to temperature sensitivity of the URM1 pathway in S. cerevisiae and S. bayanus. Furthermore, our analysis of the tRNA modification pattern in selected yeast species revealed two alternative phenotypes. Most strains moderately increase their tRNA modification levels in response to heat, possibly constituting a common adaptation to high temperatures. However, an overall reduction of nucleoside modifications was observed exclusively in S288C. This surprising finding emphasizes the importance of studies that utilize the power of evolutionary biology, and highlights the need for future systematic studies on tRNA modifications in additional model organisms.

Proton Gradient Regulation5-Like1-Mediated Cyclic Electron Flow Is Crucial for Acclimation to Anoxia and Complementary to Nonphotochemical Quenching in Stress Adaptation.

To investigate the functional importance of Proton Gradient Regulation5-Like1 (PGRL1) for photosynthetic performances in the moss Physcomitrella patens, we generated a pgrl1 knockout mutant. Functional analysis revealed diminished nonphotochemical quenching (NPQ) as well as decreased capacity for cyclic electron flow (CEF) in pgrl1. Under anoxia, where CEF is induced, quantitative proteomics evidenced severe down-regulation of photosystems but up-regulation of the chloroplast NADH dehydrogenase complex, plastocyanin, and Ca2+ sensors in the mutant, indicating that the absence of PGRL1 triggered a mechanism compensatory for diminished CEF. On the other hand, proteins required for NPQ, such as light-harvesting complex stress-related protein1 (LHCSR1), violaxanthin de-epoxidase, and PSII subunit S, remained stable. To further investigate the interrelation between CEF and NPQ, we generated a pgrl1 npq4 double mutant in the green alga Chlamydomonas reinhardtii lacking both PGRL1 and LHCSR3 expression. Phenotypic comparative analyses of this double mutant, together with the single knockout strains and with the P. patens pgrl1, demonstrated that PGRL1 is crucial for acclimation to high light and anoxia in both organisms. Moreover, the data generated for the C. reinhardtii double mutant clearly showed a complementary role of PGRL1 and LHCSR3 in managing high light stress response. We conclude that both proteins are needed for photoprotection and for survival under low oxygen, underpinning a tight link between CEF and NPQ in oxygenic photosynthesis. Given the complementarity of the energy-dependent component of NPQ (qE) and PGRL1-mediated CEF, we suggest that PGRL1 is a capacitor linked to the evolution of the PSII subunit S-dependent qE in terrestrial plants.

The Vac14-interaction network is linked to regulators of the endolysosomal and autophagic pathway.

The scaffold protein Vac14 acts in a complex with the lipid kinase PIKfyve and its counteracting phosphatase FIG4, regulating the interconversion of phosphatidylinositol-3-phosphate to phosphatidylinositol-3,5-bisphosphate. Dysfunctional Vac14 mutants, a deficiency of one of the Vac14 complex components, or inhibition of PIKfyve enzymatic activity results in the formation of large vacuoles in cells. How these vacuoles are generated and which processes are involved are only poorly understood. Here we show that ectopic overexpression of wild-type Vac14 as well as of the PIKfyve-binding deficient Vac14 L156R mutant causes vacuoles. Vac14-dependent vacuoles and PIKfyve inhibitor-dependent vacuoles resulted in elevated levels of late endosomal, lysosomal, and autophagy-associated proteins. However, only late endosomal marker proteins were bound to the membranes of these enlarged vacuoles. In order to decipher the linkage between the Vac14 complex and regulators of the endolysosomal pathway, a protein affinity approach combined with multidimensional protein identification technology was conducted, and novel molecular links were unraveled. We found and verified the interaction of Rab9 and the Rab7 GAP TBC1D15 with Vac14. The identified Rab-related interaction partners support the theory that the regulation of vesicular transport processes and phosphatidylinositol-modifying enzymes are tightly interconnected.

The interplay of light and oxygen in the reactive oxygen stress response of Chlamydomonas reinhardtii dissected by quantitative mass spectrometry.

Light and oxygen are factors that are very much entangled in the reactive oxygen species (ROS) stress response network in plants, algae and cyanobacteria. The first obligatory step in understanding the ROS network is to separate these responses. In this study, a LC-MS/MS based quantitative proteomic approach was used to dissect the responses of Chlamydomonas reinhardtii to ROS, light and oxygen employing an interlinked experimental setup. Application of novel bioinformatics tools allow high quality retention time alignment to be performed on all LC-MS/MS runs increasing confidence in protein quantification, overall sequence coverage and coverage of all treatments measured. Finally advanced hierarchical clustering yielded 30 communities of co-regulated proteins permitting separation of ROS related effects from pure light effects (induction and repression). A community termed redox(II) was identified that shows additive effects of light and oxygen with light as the first obligatory step. Another community termed 4-down was identified that shows repression as an effect of light but only in the absence of oxygen indicating ROS regulation, for example, possibly via product feedback inhibition because no ROS damage is occurring. In summary the data demonstrate the importance of separating light, O2 and ROS responses to define marker genes for ROS responses. As revealed in this study, an excellent candidate is DHAR with strong ROS dependent induction profiles.

pyGCluster, a novel hierarchical clustering approach.

SUMMARY: pyGCluster is a clustering algorithm focusing on noise injection for subsequent cluster validation. The reproducibility of a large amount of clusters obtained with agglomerative hierarchical clustering is assessed. Furthermore, a multitude of different distance-linkage combinations are evaluated. Finally, highly reproducible clusters are meta-clustered into communities. Graphical illustration of the results as node and expression maps is implemented. AVAILABILITY AND IMPLEMENTATION: pyGCluster requires Python 2.7, it is freely available at http://pyGCluster.github.io and published under MIT license. Dependencies are NumPy, SciPy and optionally fastcluster and rpy2. CONTACT: christan@fufezan.net SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online and at http://pyGCluster.github.io.