Algorithm for identification of fusion proteins via mass spectrometry

Identification of fusion proteins has contributed significantly to our understanding of cancer progression, yielding important predictive markers and therapeutic targets. While fusion proteins can be potentially identified by mass spectrometry, all previously found fusion proteins were identified using genomic (rather than mass spectrometry) technologies. This lack of MS/MS applications in studies of fusion proteins is caused by the lack of computational tools that are able to interpret mass spectra from peptides covering unknown fusion breakpoints (fusion peptides). Indeed, the number of potential fusion peptides is so large that the existing MS/MS database search tools become impractical even in the case of small genomes. We explore computational approaches to identifying fusion peptides, propose an algorithm for solving the fusion peptide identification problem, and analyze the performance of this algorithm on simulated data. We further illustrate how this approach can be modified for human exons prediction.     

Prediction of Antimicrobial Activity of Synthetic Peptides by a Decision Tree Model

Antimicrobial resistance is a persistent problem in the public health sphere. However, recent attempts to find effective substitutes to combat infections have been directed at identifying natural antimicrobial peptides in order to circumvent resistance to commercial antibiotics. This study describes the development of synthetic peptides with antimicrobial activity, created in silico by site-directed mutation modeling using wild-type peptides as scaffolds for these mutations. Fragments of antimicrobial peptides were used for modeling with molecular modeling computational tools. To analyze these peptides, a decision tree model, which indicated the action range of peptides on the types of microorganisms on which they can exercise biological activity, was created. The decision tree model was processed using physicochemistry properties from known antimicrobial peptides available at the Antimicrobial Peptide Database (APD). The two most promising peptides were synthesized, and antimicrobial assays showed inhibitory activity against Gram-positive and Gram-negative bacteria. Colossomin C and colossomin D were the most inhibitory peptides at 5 μg/ml against Staphylococcus aureus and Escherichia coli. The methods described in this work and the results obtained are useful for the identification and development of new compounds with antimicrobial activity through the use of computational tools.     

Methods and protocols for prediction of immunogenic epitopes

T-cell recognition of peptide/major histocompatibility complex (MHC) is a prerequisite for cellular immunity. Recently, there has been an influx of bioinformatics tools to facilitate the identification of T-cell epitopes to specific MHC alleles. This article examines existing computational strategies for the study of peptide/MHC interactions. The most important bioinformatics tools and methods with relevance to the study of peptide/MHC interactions have been reviewed. We have also provided guidelines for predicting antigenic peptides based on the availability of existing experimental data.     

Peptide identification by tandem mass spectrometry with alternate fragmentation modes

The high-throughput nature of proteomics mass spectrometry is enabled by a productive combination of data acquisition protocols and the computational tools used to interpret the resulting spectra. One of the key components in mainstream protocols is the generation of tandem mass (MS/MS) spectra by peptide fragmentation using collision induced dissociation, the approach currently used in the large majority of proteomics experiments to routinely identify hundreds to thousands of proteins from single mass spectrometry runs. Complementary to these, alternative peptide fragmentation methods such as electron capture/transfer dissociation and higher-energy collision dissociation have consistently achieved significant improvements in the identification of certain classes of peptides, proteins, and post-translational modifications. Recognizing these advantages, mass spectrometry instruments now conveniently support fine-tuned methods that automatically alternate between peptide fragmentation modes for either different types of peptides or for acquisition of multiple MS/MS spectra from each peptide. But although these developments have the potential to substantially improve peptide identification, their routine application requires corresponding adjustments to the software tools and procedures used for automated downstream processing. This review discusses the computational implications of alternative and alternate modes of MS/MS peptide fragmentation and addresses some practical aspects of using such protocols for identification of peptides and post-translational modifications.     

Peptide identification by database search of mixture tandem mass spectra

In high-throughput proteomics the development of computational methods and novel experimental strategies often rely on each other. In certain areas, mass spectrometry methods for data acquisition are ahead of computational methods to interpret the resulting tandem mass spectra. Particularly, although there are numerous situations in which a mixture tandem mass spectrum can contain fragment ions from two or more peptides, nearly all database search tools still make the assumption that each tandem mass spectrum comes from one peptide. Common examples include mixture spectra from co-eluting peptides in complex samples, spectra generated from data-independent acquisition methods, and spectra from peptides with complex post-translational modifications. We propose a new database search tool (MixDB) that is able to identify mixture tandem mass spectra from more than one peptide. We show that peptides can be reliably identified with up to 95% accuracy from mixture spectra while considering only a 0.01% of all possible peptide pairs (four orders of magnitude speedup). Comparison with current database search methods indicates that our approach has better or comparable sensitivity and precision at identifying single-peptide spectra while simultaneously being able to identify 38% more peptides from mixture spectra at significantly higher precision.     

Cycloquest: identification of cyclopeptides via database search of their mass spectra against genome databases

Hundreds of ribosomally synthesized cyclopeptides have been isolated from all domains of life, the vast majority having been reported in the last 15 years. Studies of cyclic peptides have highlighted their exceptional potential both as stable drug scaffolds and as biomedicines in their own right. Despite this, computational techniques for cyclopeptide identification are still in their infancy, with many such peptides remaining uncharacterized. Tandem mass spectrometry has occupied a niche role in cyclopeptide identification, taking over from traditional techniques such as nuclear magnetic resonance spectroscopy (NMR). MS/MS studies require only picogram quantities of peptide (compared to milligrams for NMR studies) and are applicable to complex samples, abolishing the requirement for time-consuming chromatographic purification. While database search tools such as Sequest and Mascot have become standard tools for the MS/MS identification of linear peptides, they are not applicable to cyclopeptides, due to the parent mass shift resulting from cyclization and different fragmentation patterns of cyclic peptides. In this paper, we describe the development of a novel database search methodology to aid in the identification of cyclopeptides by mass spectrometry and evaluate its utility in identifying two peptide rings from Helianthus annuus, a bacterial cannibalism factor from Bacillus subtilis, and a θ-defensin from Rhesus macaque.     

Computational prediction of proteotypic peptides

Evaluation of: Mallick P, Schirle M, Chen SS et al. Computational prediction of proteotypic peptides for quantitative proteomics. Nat. Biotechnol. 25(1), 125–131 (2007).

Mass spectrometry, the driving analytical force behind proteomics, is primarily used to identify and quantify as many proteins in a complex biological mixture as possible. While there are many ways to prepare samples, one aspect that is common to a vast majority of bottom-up proteomic studies is the digestion of proteins into tryptic peptides prior to their analysis by mass spectrometry. As correctly highlighted by Mallick and colleagues, only a few peptides are repeatedly and consistently identified for any given protein within a complex mixture. While the existence of these proteotypic peptides (to borrow the authors’ terminology) is well known in the proteomics community, there has never been an empirical method to recognize which peptides may be proteotypic for a given protein. In this study, the investigators discovered over 16,000 proteotypic peptides from a collection of over 600,000 peptide identifications obtained from four different analytical platforms. The study examined a number of physicochemical parameters of these peptides to determine which properties were most relevant in defining a proteotypic peptide. These characteristic properties were then used to develop computational tools to predict proteotypic peptides for any given protein within an organism.     

In silico methods for predicting T-cell epitopes: Dr Jekyll or Mr Hyde?

In silico tools offer an attractive alternative strategy to the cumbersome experimental approaches to identify T-cell epitopes. These computational tools have metamorphosed over the years into complex algorithms that attempt to efficiently predict the binding of a plethora of peptides to HLA alleles. In recent years, the scientific community has embraced these techniques to reduce the burden of wet-laboratory experimentation. Although there are some splendid examples of the utility of these methods, there are also evidences where they fall short and remain inconsistent. Hence, are these computational tools ‘Dr Jekyll’ or ‘Mr Hyde’ to the researcher, who wishes to utilize them intrepidly? This article reviews the progress and pitfalls of the in silico tools that identify T-cell epitopes.     

Proteomic analysis of purified Newcastle disease virus particles

Background

Recent advances in liquid chromatography-mass spectrometry (LC-MS) technology have led to more effective approaches for measuring changes in peptide/protein abundances in biological samples. Label-free LC-MS methods have been used for extraction of quantitative information and for detection of differentially abundant peptides/proteins. However, difference detection by analysis of data derived from label-free LC-MS methods requires various preprocessing steps including filtering, baseline correction, peak detection, alignment, and normalization. Although several specialized tools have been developed to analyze LC-MS data, determining the most appropriate computational pipeline remains challenging partly due to lack of established gold standards.

Results

The work in this paper is an initial study to develop a simple model with "presence" or "absence" condition using spike-in experiments and to be able to identify these "true differences" using available software tools. In addition to the preprocessing pipelines, choosing appropriate statistical tests and determining critical values are important. We observe that individual statistical tests could lead to different results due to different assumptions and employed metrics. It is therefore preferable to incorporate several statistical tests for either exploration or confirmation purpose.

Conclusions

The LC-MS data from our spike-in experiment can be used for developing and optimizing LC-MS data preprocessing algorithms and to evaluate workflows implemented in existing software tools. Our current work is a stepping stone towards optimizing LC-MS data acquisition and testing the accuracy and validity of computational tools for difference detection in future studies that will be focused on spiking peptides of diverse physicochemical properties in different concentrations to better represent biomarker discovery of differentially abundant peptides/proteins.     

TOPP--the OpenMS proteomics pipeline

MOTIVATION: Experimental techniques in proteomics have seen rapid development over the last few years. Volume and complexity of the data have both been growing at a similar rate. Accordingly, data management and analysis are one of the major challenges in proteomics. Flexible algorithms are required to handle changing experimental setups and to assist in developing and validating new methods. In order to facilitate these studies, it would be desirable to have a flexible 'toolbox' of versatile and user-friendly applications allowing for rapid construction of computational workflows in proteomics. RESULTS: We describe a set of tools for proteomics data analysis-TOPP, The OpenMS Proteomics Pipeline. TOPP provides a set of computational tools which can be easily combined into analysis pipelines even by non-experts and can be used in proteomics workflows. These applications range from useful utilities (file format conversion, peak picking) over wrapper applications for known applications (e.g. Mascot) to completely new algorithmic techniques for data reduction and data analysis. We anticipate that TOPP will greatly facilitate rapid prototyping of proteomics data evaluation pipelines. As such, we describe the basic concepts and the current abilities of TOPP and illustrate these concepts in the context of two example applications: the identification of peptides from a raw dataset through database search and the complex analysis of a standard addition experiment for the absolute quantitation of biomarkers. The latter example demonstrates TOPP's ability to construct flexible analysis pipelines in support of complex experimental setups. AVAILABILITY: The TOPP components are available as open-source software under the lesser GNU public license (LGPL). Source code is available from the project website at www.OpenMS.de     

Predicting HLA Class I Non-Permissive Amino Acid Residues Substitutions

Prediction of peptide binding to human leukocyte antigen (HLA) molecules is essential to a wide range of clinical entities from vaccine design to stem cell transplant compatibility. Here we present a new structure-based methodology that applies robust computational tools to model peptide-HLA (p-HLA) binding interactions. The method leverages the structural conservation observed in p-HLA complexes to significantly reduce the search space and calculate the system's binding free energy. This approach is benchmarked against existing p-HLA complexes and the prediction performance is measured against a library of experimentally validated peptides. The effect on binding activity across a large set of high-affinity peptides is used to investigate amino acid mismatches reported as high-risk factors in hematopoietic stem cell transplantation.     

In Silico Prediction of Mutant HIV-1 Proteases Cleaving a Target Sequence

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.     

In silico prediction of immunogenic T cell epitopes for HLA-DQ8

Background HLA-DQ alleles are involved in the pathogenesis of hypersensitivity reactions, with HLA-DQ8 associated with several human autoimmune disorders. Limited success has been achieved using sequence-based computational techniques for predicting HLA-DQ8-restricted T cell epitopes while accuracy and efficiency of recently developed structure-based models need to be improved. Results We describe a combined structure-based prediction approach for DQ8-restricted T cell epitope prediction using a recently developed fast and accurate docking protocol, pDOCK, and molecular surface electrostatic potential (MSEP)-based clustering of pMHC binding interfaces. The prediction model was rigorously trained, tested and validated using experimentally verified DQ8 binding and non-binding peptides. High prediction accuracy (average area under the ROC curve, average AROC>0.94) is validated against experimental data. Our model also predicts all binding registers correctly and known T cell activators with 77% accuracy. We also studied the patterns of DQ8-binding peptides and reassure the existence of epitopes not conforming to binding motifs. Conclusions We have developed a model that can be successfully applied as a generic protocol for easy in silico identification of potential immunogenic T cell epitopes. The current model is therefore applicable for screening vaccine candidates irrespective of sequence motifs. We have also illustrated efficient discrimination of different categories of binders from non-binders as well as different categories of pMHC agonists from non-agonists, while accurately predicting the binding registers of DQ8-restricted peptides. This combined approach provides a set of sensitive and specific computational tools to facilitate high-throughput screening of peptides for immunotherapeutic applications such as controlling allergic and autoimmune responses.     

Current progress in computational metabolomics

Being a relatively new addition to the 'omics' field, metabolomics is still evolving its own computational infrastructure and assessing its own computational needs. Due to its strong emphasis on chemical information and because of the importance of linking that chemical data to biological consequences, metabolomics must combine elements of traditional bioinformatics with traditional cheminformatics. This is a significant challenge as these two fields have evolved quite separately and require very different computational tools and skill sets. This review is intended to familiarize readers with the field of metabolomics and to outline the needs, the challenges and the recent progress being made in four areas of computational metabolomics: (i) metabolomics databases; (ii) metabolomics LIMS; (iii) spectral analysis tools for metabolomics and (iv) metabolic modeling.     

The AutoDock suite at 30

The AutoDock suite provides a comprehensive toolset for computational ligand docking and drug design and development. The suite builds on 30 years of methods development, including empirical free energy force fields, docking engines, methods for site prediction, and interactive tools for visualization and analysis. Specialized tools are available for challenging systems, including covalent inhibitors, peptides, compounds with macrocycles, systems where ordered hydration plays a key role, and systems with substantial receptor flexibility. All methods in the AutoDock suite are freely available for use and reuse, which has engendered the continued growth of a diverse community of primary users and third‐party developers.     

tmRDB (tmRNA database)

The tmRNA database (tmRDB) is maintained at the University of Texas Health Science Center at Tyler, Texas, and accessible on the World Wide Web at the URL http://psyche.uthct.edu/dbs/tmRDB/tmRDB.++ +html. Mirror sites are located at Auburn University, Auburn, Alabama (http://www.ag.auburn.edu/mirror/tmRDB/) and the Institute of Biological Sciences, Aarhus, Denmark (http://www.bioinf.au. dk/tmRDB/). The tmRDB provides information and citation links about tmRNA, a molecule that combines functions of tRNA and mRNA in trans-translation. tmRNA is likely to be present in all bacteria and has been found in algae chloroplasts, the cyanelle of Cyanophora paradoxa and the mitochondrion of the flagellate Reclinomonas americana. This release adds 26 new sequences and corresponding predicted tmRNA-encoded tag peptides for a total of 86 tmRNAs, ordered alphabetically and phylogenetically. Secondary structures and three-dimensional models in PDB format for representative molecules are being made available. tmRNA alignments prove individual base pairs and are generated manually assisted by computational tools. The alignments with their corresponding structural annotation can be obtained in various formats, including a new column format designed to improve and simplify computational usability of the data.     

Computationally designed peptide inhibitors of protein-protein interactions in membranes

We recently reported a computational method (CHAMP) for designing sequence-specific peptides that bind to the membrane-embedded portions of transmembrane proteins. We successfully applied this method to design membrane-spanning peptides targeting the transmembrane domains of the alpha IIb subunit of integrin alpha IIbbeta 3. Previously, we demonstrated that these CHAMP peptides bind specifically with reasonable affinity to isolated transmembrane helices of the targeted transmembrane region. These peptides also induced integrin alpha IIbbeta 3 activation due to disruption of the helix-helix interactions between the transmembrane domains of the alpha IIb and beta 3 subunits. In this paper, we show the direct interaction of the designed anti-alpha IIb CHAMP peptide with isolated full-length integrin alpha IIbbeta 3 in detergent micelles. Further, the behavior of the designed peptides in phospholipid bilayers is essentially identical to their behavior in detergent micelles. In particular, the peptides assume a membrane-spanning alpha-helical conformation that does not disrupt bilayer integrity. The activity and selectivity of the CHAMP peptides were further explored in platelets, comfirming that anti-alpha IIb activates wild-type alpha IIbbeta 3 in whole cells as a result of its disruption of the protein-protein interactions between the alpha and beta subunits in the transmembrane regions. These results demonstrate that CHAMP is a successful chemical biology approach that can provide specific tools for probing the transmembrane domains of proteins.     

Signal peptide prediction based on analysis of experimentally verified cleavage sites

A number of computational tools are available for detecting signal peptides, but their abilities to locate the signal peptide cleavage sites vary significantly and are often less than satisfactory. We characterized a set of 270 secreted recombinant human proteins by automated Edman analysis and used the verified cleavage sites to evaluate the success rate of a number of computational prediction programs. An examination of the frequency of amino acid in the N-terminal region of the data set showed a preference of proline and glutamine but a bias against tyrosine. The data set was compared to the SWISS-PROT database and revealed a high percentage of discrepancies with cleavage site annotations that were computationally generated. The best program for predicting signal sequences was found to be SignalP 2.0-NN with an accuracy of 78.1% for cleavage site recognition. The new data set can be utilized for refining prediction algorithms, and we have built an improved version of profile hidden Markov model for signal peptides based on the new data.