Generating and navigating proteome maps using mass spectrometry.

Proteomes, the ensembles of all proteins expressed by cells or tissues, are typically analysed by mass spectrometry. Recent technical and computational advances have greatly increased the fraction of a proteome that can be identified and quantified in a single study. Current mass spectrometry-based proteomic strategies have the potential to reproducibly, accurately, quantitatively and comprehensively measure any protein or whole proteomes from cells and tissues at different states. Achieving these goals will require complete proteome maps and analytical strategies that use these maps as prior information and will greatly enhance the impact of proteomics on biological and clinical research.

Protein and gene model inference based on statistical modeling in k-partite graphs.

One of the major goals of proteomics is the comprehensive and accurate description of a proteome. Shotgun proteomics, the method of choice for the analysis of complex protein mixtures, requires that experimentally observed peptides are mapped back to the proteins they were derived from. This process is also known as protein inference. We present Markovian Inference of Proteins and Gene Models (MIPGEM), a statistical model based on clearly stated assumptions to address the problem of protein and gene model inference for shotgun proteomics data. In particular, we are dealing with dependencies among peptides and proteins using a Markovian assumption on k-partite graphs. We are also addressing the problems of shared peptides and ambiguous proteins by scoring the encoding gene models. Empirical results on two control datasets with synthetic mixtures of proteins and on complex protein samples of Saccharomyces cerevisiae, Drosophila melanogaster, and Arabidopsis thaliana suggest that the results with MIPGEM are competitive with existing tools for protein inference.

Deterministic protein inference for shotgun proteomics data provides new insights into Arabidopsis pollen development and function.

Pollen, the male gametophyte of flowering plants, represents an ideal biological system to study developmental processes, such as cell polarity, tip growth, and morphogenesis. Upon hydration, the metabolically quiescent pollen rapidly switches to an active state, exhibiting extremely fast growth. This rapid switch requires relevant proteins to be stored in the mature pollen, where they have to retain functionality in a desiccated environment. Using a shotgun proteomics approach, we unambiguously identified approximately 3500 proteins in Arabidopsis pollen, including 537 proteins that were not identified in genetic or transcriptomic studies. To generate this comprehensive reference data set, which extends the previously reported pollen proteome by a factor of 13, we developed a novel deterministic peptide classification scheme for protein inference. This generally applicable approach considers the gene model-protein sequence-protein accession relationships. It allowed us to classify and eliminate ambiguities inherently associated with any shotgun proteomics data set, to report a conservative list of protein identifications, and to seamlessly integrate data from previous transcriptomics studies. Manual validation of proteins unambiguously identified by a single, information-rich peptide enabled us to significantly reduce the false discovery rate, while keeping valuable identifications of shorter and lower abundant proteins. Bioinformatic analyses revealed a higher stability of pollen proteins compared to those of other tissues and implied a protein family of previously unknown function in vesicle trafficking. Interestingly, the pollen proteome is most similar to that of seeds, indicating physiological similarities between these developmentally distinct tissues.

Identification and functional characterization of N-terminally acetylated proteins in Drosophila melanogaster.

Protein modifications play a major role for most biological processes in living organisms. Amino-terminal acetylation of proteins is a common modification found throughout the tree of life: the N-terminus of a nascent polypeptide chain becomes co-translationally acetylated, often after the removal of the initiating methionine residue. While the enzymes and protein complexes involved in these processes have been extensively studied, only little is known about the biological function of such N-terminal modification events. To identify common principles of N-terminal acetylation, we analyzed the amino-terminal peptides from proteins extracted from Drosophila Kc167 cells. We detected more than 1,200 mature protein N-termini and could show that N-terminal acetylation occurs in insects with a similar frequency as in humans. As the sole true determinant for N-terminal acetylation we could extract the (X)PX rule that indicates the prevention of acetylation under all circumstances. We could show that this rule can be used to genetically engineer a protein to study the biological relevance of the presence or absence of an acetyl group, thereby generating a generic assay to probe the functional importance of N-terminal acetylation. We applied the assay by expressing mutated proteins as transgenes in cell lines and in flies. Here, we present a straightforward strategy to systematically study the functional relevance of N-terminal acetylations in cells and whole organisms. Since the (X)PX rule seems to be of general validity in lower as well as higher eukaryotes, we propose that it can be used to study the function of N-terminal acetylation in all species.

An integrated proteomics and transcriptomics reference data set provides new insights into the Bradyrhizobium japonicum bacteroid metabolism in soybean root nodules.

Bradyrhizobium japonicum, a gram-negative soil bacterium that establishes an N(2)-fixing symbiosis with its legume host soybean (Glycine max), has been used as a symbiosis model system. Using a sensitive geLC-MS/MS proteomics approach, we report the identification of 2315 B. japonicum strain USDA110 proteins (27.8% of the theoretical proteome) that are expressed 21 days post infection in symbiosis with soybean cultivated in growth chambers, substantially expanding the previously known symbiosis proteome. Integration of transcriptomics data generated under the same conditions (2780 expressed genes) allowed us to compile a comprehensive expression profile of B. japonicum during soybean symbiosis, which comprises 3587 genes/proteins (43% of the predicted B. japonicum genes/proteins). Analysis of this data set revealed both the biases and the complementarity of these global profiling technologies. A functional classification and pathway analysis showed that most of the proteins involved in carbon and nitrogen metabolism are expressed, including a complete set of tricarboxylic acid cycle enzymes, several gluconeogenesis and pentose phosphate pathway enzymes, as well as several proteins that were previously not considered to be present during symbiosis. Congruent results were obtained for B. japonicum bacteroids harvested from soybeans grown under field conditions.

Visualizing regulatory interactions in metabolic networks.

BACKGROUND: Direct visualization of data sets in the context of biochemical network drawings is one of the most appealing approaches in the field of data evaluation within systems biology. One important type of information that is very helpful in interpreting and understanding metabolic networks has been overlooked so far. Here we focus on the representation of this type of information given by the strength of regulatory interactions between metabolite pools and reaction steps. RESULTS: The visualization of such interactions in a given metabolic network is based on a novel concept defining the regulatory strength (RS) of effectors regulating certain reaction steps. It is applicable to any mechanistic reaction kinetic formula. The RS values are measures for the strength of an up- or down-regulation of a reaction step compared with the completely non-inhibited or non-activated state, respectively. One numerical RS value is associated to any effector edge contained in the network. The RS is approximately interpretable on a percentage scale where 100% means the maximal possible inhibition or activation, respectively, and 0% means the absence of a regulatory interaction. If many effectors influence a certain reaction step, the respective percentages indicate the proportion in which the different effectors contribute to the total regulation of the reaction step. The benefits of the proposed method are demonstrated with a complex example system of a dynamic E. coli network. CONCLUSION: The presented visualization approach is suitable for an intuitive interpretation of simulation data of metabolic networks under dynamic as well as steady-state conditions. Huge amounts of simulation data can be analyzed in a quick and comprehensive way. An extended time-resolved graphical network presentation provides a series of information about regulatory interaction within the biological system under investigation.

Improved prediction of peptide detectability for targeted proteomics using a rank-based algorithm and organism-specific data.

The in silico prediction of the best-observable "proteotypic" peptides in mass spectrometry-based workflows is a challenging problem. Being able to accurately predict such peptides would enable the informed selection of proteotypic peptides for targeted quantification of previously observed and non-observed proteins for any organism, with a significant impact for clinical proteomics and systems biology studies. Current prediction algorithms rely on physicochemical parameters in combination with positive and negative training sets to identify those peptide properties that most profoundly affect their general detectability. Here we present PeptideRank, an approach that uses learning to rank algorithm for peptide detectability prediction from shotgun proteomics data, and that eliminates the need to select a negative dataset for the training step. A large number of different peptide properties are used to train ranking models in order to predict a ranking of the best-observable peptides within a protein. Empirical evaluation with rank accuracy metrics showed that PeptideRank complements existing prediction algorithms. Our results indicate that the best performance is achieved when it is trained on organism-specific shotgun proteomics data, and that PeptideRank is most accurate for short to medium-sized and abundant proteins, without any loss in prediction accuracy for the important class of membrane proteins. BIOLOGICAL SIGNIFICANCE: Targeted proteomics approaches have been gaining a lot of momentum and hold immense potential for systems biology studies and clinical proteomics. However, since only very few complete proteomes have been reported to date, for a considerable fraction of a proteome there is no experimental proteomics evidence that would allow to guide the selection of the best-suited proteotypic peptides (PTPs), i.e. peptides that are specific to a given proteoform and that are repeatedly observed in a mass spectrometer. We describe a novel, rank-based approach for the prediction of the best-suited PTPs for targeted proteomics applications. By building on methods developed in the field of information retrieval (e.g. web search engines like Google's PageRank), we circumvent the delicate step of selecting positive and negative training sets and at the same time also more closely reflect the experimentalist´s need for selecting e.g. the 5 most promising peptides for targeting a protein of interest. This approach allows to predict PTPs for not yet observed proteins or for organisms without prior experimental proteomics data such as many non-model organisms.

A mighty small heart: the cardiac proteome of adult Drosophila melanogaster.

Drosophila melanogaster is emerging as a powerful model system for the study of cardiac disease. Establishing peptide and protein maps of the Drosophila heart is central to implementation of protein network studies that will allow us to assess the hallmarks of Drosophila heart pathogenesis and gauge the degree of conservation with human disease mechanisms on a systems level. Using a gel-LC-MS/MS approach, we identified 1228 protein clusters from 145 dissected adult fly hearts. Contractile, cytostructural and mitochondrial proteins were most abundant consistent with electron micrographs of the Drosophila cardiac tube. Functional/Ontological enrichment analysis further showed that proteins involved in glycolysis, Ca(2+)-binding, redox, and G-protein signaling, among other processes, are also over-represented. Comparison with a mouse heart proteome revealed conservation at the level of molecular function, biological processes and cellular components. The subsisting peptidome encompassed 5169 distinct heart-associated peptides, of which 1293 (25%) had not been identified in a recent Drosophila peptide compendium. PeptideClassifier analysis was further used to map peptides to specific gene-models. 1872 peptides provide valuable information about protein isoform groups whereas a further 3112 uniquely identify specific protein isoforms and may be used as a heart-associated peptide resource for quantitative proteomic approaches based on multiple-reaction monitoring. In summary, identification of excitation-contraction protein landmarks, orthologues of proteins associated with cardiovascular defects, and conservation of protein ontologies, provides testimony to the heart-like character of the Drosophila cardiac tube and to the utility of proteomics as a complement to the power of genetics in this growing model of human heart disease.