Proteomic Analyses Reveal Divergent Ubiquitylation Site Patterns in Murine Tissues

Posttranslational modifications of proteins increase the complexity of the cellular proteome and enable rapid regulation of protein functions in response to environmental changes. Protein ubiquitylation is a central regulatory posttranslational modification that controls numerous biological processes including proteasomal degradation of proteins, DNA damage repair and innate immune responses. Here we combine high-resolution mass spectrometry with single-step immunoenrichment of di-glycine modified peptides for mapping of endogenous putative ubiquitylation sites in murine tissues. We identify more than 20,000 unique ubiquitylation sites on proteins involved in diverse biological processes. Our data reveals that ubiquitylation regulates core signaling pathways common for each of the studied tissues. In addition, we discover that ubiquitylation regulates tissue-specific signaling networks. Many tissue-specific ubiquitylation sites were obtained from brain highlighting the complexity and unique physiology of this organ. We further demonstrate that different di-glycine-lysine-specific monoclonal antibodies exhibit sequence preferences, and that their complementary use increases the depth of ubiquitylation site analysis, thereby providing a more unbiased view of protein ubiquitylation.Ubiquitin is a small 76-amino-acid protein that is conjugated to the ε-amino group of lysines in a highly orchestrated enzymatic cascade involving ubiquitin activating (E1), ubiquitin conjugating (E2), and ubiquitin ligase (E3) enzymes (1). Ubiquitylation is involved in the regulation of diverse cellular processes including protein degradation (2, 3, 4), DNA damage repair (5, 6), DNA replication (7), cell surface receptor endocytosis, and innate immune signaling (8, 9). Deregulation of protein ubiquitylation is implicated in the development of cancer and neurodegenerative diseases (10, 11). Inhibitors targeting the ubiquitin proteasome system are used in the treatment of hematologic malignancies such as multiple myeloma (12, 13).Recent developments in the mass spectrometry (MS)-based proteomics have greatly expedited proteome-wide analysis of posttranslational modifications (PTMs) (14–17). Large-scale mapping of ubiquitylation sites by mass spectrometry is based on the identification of the di-glycine remnant that results from trypsin digestion of ubiquitylated proteins and remains attached to ubiquitylated lysines (18). Recently, two monoclonal antibodies were developed that specifically recognize di-glycine remnant modified peptides enabling their efficient enrichment from complex peptide mixtures (19, 20). These antibodies have been used to identify thousands of endogenous ubiquitylation sites in human cells, and to quantify site-specific changes in ubiquitylation in response to different cellular perturbations (20–22). It should be noted that the di-glycine remnant is not specific for proteins modified by ubiquitin but also proteins modified by NEDD8 or ISG15 generate an identical di-glycine remnant on modified lysines making it impossible to distinguish between these modifications by mass spectrometry. However, expression of NEDD8 in mouse tissues was shown to be developmentally down-regulated (23), and ISG15 expression in bovine tissues is low in the absence of interferon stimulation (24). In cell culture experiments it was shown that a great majority of sites identified using di-glycine-lysine-specific antibodies stems from ubiquitylated peptides (20).The rates of cell proliferation and protein turnover in mammals vary dramatically between different tissues. Immortalized cell lines, often derived from cancer, are selected for high proliferation rates and fail to represent the complex conditions in tissues. Tissue proteomics can help to gain a more comprehensive understanding of physiological processes in multicellular organisms. Analysis of tissue proteome and PTMs can provide important insights into tissue-specific processes and signaling networks that regulate these processes (25–32). In addition, development of mass spectrometric methods for analysis of PTMs in diseased tissues might lead to the identification of biomarkers.In this study, we combined high-resolution mass spectrometry with immunoenrichment of di-glycine modified peptides to investigate endogenous ubiquitylation sites in murine tissues. We identified more than 20,000 ubiquitylation sites from five different murine tissues and report the largest ubiquitylation dataset obtained from mammalian tissues to date. Furthermore, we compared the performance of the two monoclonal di-glycine-lysine-specific antibodies available for enrichment of ubiquitylated peptides, and reveal their relative preferences for different amino acids flanking ubiquitylation sites.     

High-definition De Novo Sequencing of Crustacean Hyperglycemic Hormone (CHH)-family Neuropeptides

A complete understanding of the biological functions of large signaling peptides (>4 kDa) requires comprehensive characterization of their amino acid sequences and post-translational modifications, which presents significant analytical challenges. In the past decade, there has been great success with mass spectrometry-based de novo sequencing of small neuropeptides. However, these approaches are less applicable to larger neuropeptides because of the inefficient fragmentation of peptides larger than 4 kDa and their lower endogenous abundance. The conventional proteomics approach focuses on large-scale determination of protein identities via database searching, lacking the ability for in-depth elucidation of individual amino acid residues. Here, we present a multifaceted MS approach for identification and characterization of large crustacean hyperglycemic hormone (CHH)-family neuropeptides, a class of peptide hormones that play central roles in the regulation of many important physiological processes of crustaceans. Six crustacean CHH-family neuropeptides (8–9.5 kDa), including two novel peptides with extensive disulfide linkages and PTMs, were fully sequenced without reference to genomic databases. High-definition de novo sequencing was achieved by a combination of bottom-up, off-line top-down, and on-line top-down tandem MS methods. Statistical evaluation indicated that these methods provided complementary information for sequence interpretation and increased the local identification confidence of each amino acid. Further investigations by MALDI imaging MS mapped the spatial distribution and colocalization patterns of various CHH-family neuropeptides in the neuroendocrine organs, revealing that two CHH-subfamilies are involved in distinct signaling pathways.Neuropeptides and hormones comprise a diverse class of signaling molecules involved in numerous essential physiological processes, including analgesia, reward, food intake, learning and memory (1). Disorders of the neurosecretory and neuroendocrine systems influence many pathological processes. For example, obesity results from failure of energy homeostasis in association with endocrine alterations (2, 3). Previous work from our lab used crustaceans as model organisms found that multiple neuropeptides were implicated in control of food intake, including RFamides, tachykinin related peptides, RYamides, and pyrokinins (4–6).Crustacean hyperglycemic hormone (CHH)¹ family neuropeptides play a central role in energy homeostasis of crustaceans (7–17). Hyperglycemic response of the CHHs was first reported after injection of crude eyestalk extract in crustaceans. Based on their preprohormone organization, the CHH family can be grouped into two sub-families: subfamily-I containing CHH, and subfamily-II containing molt-inhibiting hormone (MIH) and mandibular organ-inhibiting hormone (MOIH). The preprohormones of the subfamily-I have a CHH precursor related peptide (CPRP) that is cleaved off during processing; and preprohormones of the subfamily-II lack the CPRP (9). Uncovering their physiological functions will provide new insights into neuroendocrine regulation of energy homeostasis.Characterization of CHH-family neuropeptides is challenging. They are comprised of more than 70 amino acids and often contain multiple post-translational modifications (PTMs) and complex disulfide bridge connections (7). In addition, physiological concentrations of these peptide hormones are typically below picomolar level, and most crustacean species do not have available genome and proteome databases to assist MS-based sequencing.MS-based neuropeptidomics provides a powerful tool for rapid discovery and analysis of a large number of endogenous peptides from the brain and the central nervous system. Our group and others have greatly expanded the peptidomes of many model organisms (3, 18–33). For example, we have discovered more than 200 neuropeptides with several neuropeptide families consisting of as many as 20–40 members in a simple crustacean model system (5, 6, 25–31, 34). However, a majority of these neuropeptides are small peptides with 5–15 amino acid residues long, leaving a gap of identifying larger signaling peptides from organisms without sequenced genome. The observed lack of larger size peptide hormones can be attributed to the lack of effective de novo sequencing strategies for neuropeptides larger than 4 kDa, which are inherently more difficult to fragment using conventional techniques (34–37). Although classical proteomics studies examine larger proteins, these tools are limited to identification based on database searching with one or more peptides matching without complete amino acid sequence coverage (36, 38).Large populations of neuropeptides from 4–10 kDa exist in the nervous systems of both vertebrates and invertebrates (9, 39, 40). Understanding their functional roles requires sufficient molecular knowledge and a unique analytical approach. Therefore, developing effective and reliable methods for de novo sequencing of large neuropeptides at the individual amino acid residue level is an urgent gap to fill in neurobiology. In this study, we present a multifaceted MS strategy aimed at high-definition de novo sequencing and comprehensive characterization of the CHH-family neuropeptides in crustacean central nervous system. The high-definition de novo sequencing was achieved by a combination of three methods: (1) enzymatic digestion and LC-tandem mass spectrometry (MS/MS) bottom-up analysis to generate detailed sequences of proteolytic peptides; (2) off-line LC fractionation and subsequent top-down MS/MS to obtain high-quality fragmentation maps of intact peptides; and (3) on-line LC coupled to top-down MS/MS to allow rapid sequence analysis of low abundance peptides. Combining the three methods overcomes the limitations of each, and thus offers complementary and high-confidence determination of amino acid residues. We report the complete sequence analysis of six CHH-family neuropeptides including the discovery of two novel peptides. With the accurate molecular information, MALDI imaging and ion mobility MS were conducted for the first time to explore their anatomical distribution and biochemical properties.     

The Effect of Proteasome Inhibition on the Generation of the Human Leukocyte Antigen (HLA) Peptidome

The Major histocompatibility complex (MHC) class I peptidome is thought to be generated mostly through proteasomal degradation of cellular proteins, a notion that is based on the alterations in presentation of selected peptides following proteasome inhibition. We evaluated the effects of proteasome inhibitors, epoxomicin and bortezomib, on human cultured cancer cells. Because the inhibitors did not reduce the level of presentation of the cell surface human leukocyte antigen (HLA) molecules, we followed their effects on the rates of synthesis of both HLA peptidome and proteome of the cells, using dynamic stable isotope labeling in tissue culture (dynamic-SILAC). The inhibitors reduced the rates of synthesis of most cellular proteins and HLA peptides, yet the synthesis rates of some of the proteins and HLA peptides was not decreased by the inhibitors and of some even increased. Therefore, we concluded that the inhibitors affected the production of the HLA peptidome in a complex manner, including modulation of the synthesis rates of the source proteins of the HLA peptides, in addition to their effect on their degradation. The collected data may suggest that the current reliance on proteasome inhibition may overestimate the centrality of the proteasome in the generation of the MHC peptidome. It is therefore suggested that the relative contribution of the proteasomal and nonproteasomal pathways to the production of the MHC peptidome should be revaluated in accordance with the inhibitors effects on the synthesis rates of the source proteins of the MHC peptides.The repertoires and levels of peptides, presented by the major histocompatibility complex (MHC)¹ class I molecules at the cells'' surface, are modulated by multiple factors. These include the rates of synthesis and degradation of their source proteins, the transport efficacy of the peptides through the transporter associated with antigen processing (TAP) into the endoplasmic reticulum (ER), their subsequent processing and loading onto the MHC molecules within the ER, and the rates of transport of the MHC molecules with their peptide cargo to the cell surface. The off-rates of the presented peptides, the residence time of the MHC complexes at the cell surface, and their retrograde transport back into the cytoplasm, influence, as well, the presented peptidomes (reviewed in (1)). Even though significant portions of the MHC class I peptidomes are thought to be derived from newly synthesized proteins, including misfolded proteins, defective ribosome products (DRiPs), and short lived proteins (SLiPs), most of the MHC peptidome is assumed to originate from long-lived proteins, which completed their functional cellular roles or became defective (retirees), (reviewed in (2)).The main protease, supplying the MHC peptidome production pipeline, is thought to be the proteasome (3). It is also responsible for generation of the final carboxyl termini of the MHC peptides (4), (reviewed in (5)). The final trimming of the n-termini of the peptides, within the endoplasmic reticulum (ER), is thought to be performed by amino peptidases, such as ERAP1/ERAAP, which fit the peptides into their binding groove on the MHC molecules (6) (reviewed in (7)). Nonproteasomal proteolytic pathways were also suggested as possible contributors to the MHC peptidome, including proteolysis by the ER resident Signal peptide peptidase (8, 9), the cytoplasmic proteases Insulin degrading enzyme (10), Tripeptidyl peptidase (11–13), and a number of proteases within the endolysosome pathway (reviewed recently in (14–17)). In contrast to the mostly cytoplasmic and ER production of the MHC class I peptidome, the class II peptidome is produced in a special compartment, associated with the endolysosome pathway (18–20). This pathway is also thought to participate in the cross presentation of class I peptides, derived from proteins up-taken by professional antigen presenting cells (21), (reviewed in (15–17, 22)).The centrality of the proteasomes in the generation of the MHC peptidome was deduced mostly from the observed change in presentation levels of small numbers of selected peptides, following proteasome inhibition (3, 23). Even the location of some of the genes encoding the catalytic subunits of the immunoproteasome (LMP2 and LMP7) (24) within the MHC class II genomic locus, was suggested to support the involvement of the proteasome in the generation of the MHC class I peptidome (25). Similar conclusions were deduced from alterations in peptide presentation, following expression of the catalytic subunits of the immunoproteasome (26), (reviewed in (5)). Yet, although most of the reports indicated reductions in presentation of selected peptides by proteasome inhibition (3, 27–29), others have observed only limited, and sometimes even opposite effects (23, 30–32).The matter is further complicated by the indirect effects of proteasome inhibition used for such studies on the arrest of protein synthesis by the cells (33–35), on the transport rates of the MHC molecules to the cell surface, and on their retrograde transport back to the vesicular system (36) (reviewed in (37)). Proteasome inhibition likely causes shortage of free ubiquitin, reduced supply of free amino acids, and induces an ER unfolded protein response (UPR), which signals the cells to block most (but not all) cellular protein synthesis (reviewed in (38)). Because a significant portion of the MHC peptidome originates from degradation of DRiPs and SLiPs (reviewed in (2)), arrest of new protein synthesis should influence the presentation of their derived MHC peptides. Taken together, these arguments may suggest that merely following the changes in the presentation levels of the MHC molecules, or even of specific MHC peptides, after proteasome inhibition, does not provide the full picture for deducing the relative contribution of the proteasomal pathway to the production of the MHC peptidome (reviewed in (7)).MHC peptidome analysis can be performed relatively easily by modern capillary chromatography combined with mass spectrometry (reviewed in (39)). The peptides are recovered from immunoaffinity purified MHC molecules after detergent solubilization of the cells (40, 41), from soluble MHC molecules secreted to the cells'' growth medium (42, 43) or from patients'' plasma (44). The purified peptides pools are resolved by capillary chromatography and the individual peptides are identified and quantified by tandem mass spectrometry (40), (reviewed in (45–47)). In cultured cells, quantitative analysis can also be followed by metabolic incorporation of stable isotope labeled amino acids (SILAC) (48). Furthermore, the rates of de novo synthesis of both MHC peptides and their proteins of origin can be followed using the dynamic-SILAC proteomics approach (49) with its further adaptation to HLA peptidomics (50–52).This study attempts to define the relative contribution of the proteasomes to the production of HLA class I peptidome by simultaneously following the effects of proteasome inhibitors, epoxomicin and bortezomib (Velcade), on the rates of de novo synthesis of both the HLA class I peptidome and the cellular proteome of the same MCF-7 human breast cancer cultured cells. The proteasome inhibitors did not reduce the levels of HLA presentations, yet affected the rates of production of both the HLA peptidome and cellular proteome, mostly decreasing, but also increasing some of the synthesis rates of the HLA peptides and cellular proteins. Thus, we suggest that the degree of contribution of the proteasomal pathway to the production of the HLA-I peptidome should be re-evaluated in accordance with their effects on the entire HLA class-I peptidome, while taking into consideration the inhibitors'' effects on the synthesis (and degradation) rates of the source proteins of each of the studied HLA peptides.     

Endogenous Peptide Discovery of the Rat Circadian Clock: A FOCUSED STUDY OF THE SUPRACHIASMATIC NUCLEUS BY ULTRAHIGH PERFORMANCE TANDEM MASS SPECTROMETRY*

Understanding how a small brain region, the suprachiasmatic nucleus (SCN), can synchronize the body''s circadian rhythms is an ongoing research area. This important time-keeping system requires a complex suite of peptide hormones and transmitters that remain incompletely characterized. Here, capillary liquid chromatography and FTMS have been coupled with tailored software for the analysis of endogenous peptides present in the SCN of the rat brain. After ex vivo processing of brain slices, peptide extraction, identification, and characterization from tandem FTMS data with <5-ppm mass accuracy produced a hyperconfident list of 102 endogenous peptides, including 33 previously unidentified peptides, and 12 peptides that were post-translationally modified with amidation, phosphorylation, pyroglutamylation, or acetylation. This characterization of endogenous peptides from the SCN will aid in understanding the molecular mechanisms that mediate rhythmic behaviors in mammals.Central nervous system neuropeptides function in cell-to-cell signaling and are involved in many physiological processes such as circadian rhythms, pain, hunger, feeding, and body weight regulation (1–4). Neuropeptides are produced from larger protein precursors by the selective action of endopeptidases, which cleave at mono- or dibasic sites and then remove the C-terminal basic residues (1, 2). Some neuropeptides undergo functionally important post-translational modifications (PTMs),1 including amidation, phosphorylation, pyroglutamylation, or acetylation. These aspects of peptide synthesis impact the properties of neuropeptides, further expanding their diverse physiological implications. Therefore, unveiling new peptides and unreported peptide properties is critical to advancing our understanding of nervous system function.Historically, the analysis of neuropeptides was performed by Edman degradation in which the N-terminal amino acid is sequentially removed. However, analysis by this method is slow and does not allow for sequencing of the peptides containing N-terminal PTMs (5). Immunological techniques, such as radioimmunoassay and immunohistochemistry, are used for measuring relative peptide levels and spatial localization, but these methods only detect peptide sequences with known structure (6). More direct, high throughput methods of analyzing brain regions can be used.Mass spectrometry, a rapid and sensitive method that has been used for the analysis of complex biological samples, can detect and identify the precise forms of neuropeptides without prior knowledge of peptide identity, with these approaches making up the field of peptidomics (7–12). The direct tissue and single neuron analysis by MALDI MS has enabled the discovery of hundreds of neuropeptides in the last decade, and the neuronal homogenate analysis by fractionation and subsequent ESI or MALDI MS has yielded an equivalent number of new brain peptides (5). Several recent peptidome studies, including the work by Dowell et al. (10), have used the specificity of FTMS for peptide discovery (10, 13–15). Here, we combine the ability to fragment ions at ultrahigh mass accuracy (16) with a software pipeline designed for neuropeptide discovery. We use nanocapillary reversed-phase LC coupled to 12 Tesla FTMS for the analysis of peptides present in the suprachiasmatic nucleus (SCN) of rat brain.A relatively small, paired brain nucleus located at the base of the hypothalamus directly above the optic chiasm, the SCN contains a biological clock that generates circadian rhythms in behaviors and homeostatic functions (17, 18). The SCN comprises ∼10,000 cellular clocks that are integrated as a tissue level clock which, in turn, orchestrates circadian rhythms throughout the brain and body. It is sensitive to incoming signals from the light-sensing retina and other brain regions, which cause temporal adjustments that align the SCN appropriately with changes in environmental or behavioral state. Previous physiological studies have implicated peptides as critical synchronizers of normal SCN function as well as mediators of SCN inputs, internal signal processing, and outputs; however, only a small number of peptides have been identified and explored in the SCN, leaving unresolved many circadian mechanisms that may involve peptide function.Most peptide expression in the SCN has only been studied through indirect antibody-based techniques (19–29), although we recently used MS approaches to characterize several peptides detected in SCN releasates (30). Previous studies indicate that the SCN expresses a rich diversity of peptides relative to other brain regions studied with the same techniques. Previously used immunohistochemical approaches are not only inadequate for comprehensively evaluating PTMs and alternate isoforms of known peptides but are also incapable of exhaustively examining the full peptide complement of this complex biological network of peptidergic inputs and intrinsic components. A comprehensive study of SCN peptidomics is required that utilizes high resolution strategies for directly analyzing the peptide content of the neuronal networks comprising the SCN.In our study, the SCN was obtained from ex vivo coronal brain slices via tissue punch and subjected to multistage peptide extraction. The SCN tissue extract was analyzed by FTMS/MS, and the high resolution MS and MS/MS data were processed using ProSightPC 2.0 (16), which allows the identification and characterization of peptides or proteins from high mass accuracy MS/MS data. In addition, the Sequence Gazer included in ProSightPC was used for manually determining PTMs (31, 32). As a result, a total of 102 endogenous peptides were identified, including 33 that were previously unidentified, and 12 PTMs (including amidation, phosphorylation, pyroglutamylation, and acetylation) were found. The present study is the first comprehensive peptidomics study for identifying peptides present within the mammalian SCN. In fact, this is one of the first peptidome studies to work with discrete brain nuclei as opposed to larger brain structures and follows up on our recent report using LC-ion trap for analysis of the peptides in the supraoptic nucleus (33); here, the use of FTMS allows a greater range of PTMs to be confirmed and allows higher confidence in the peptide assignments. This information on the peptides in the SCN will serve as a basis to more exhaustively explore the extent that previously unreported SCN neuropeptides may function in SCN regulation of mammalian circadian physiology.     

Development of a Novel Cross-linking Strategy for Fast and Accurate Identification of Cross-linked Peptides of Protein Complexes

Knowledge of elaborate structures of protein complexes is fundamental for understanding their functions and regulations. Although cross-linking coupled with mass spectrometry (MS) has been presented as a feasible strategy for structural elucidation of large multisubunit protein complexes, this method has proven challenging because of technical difficulties in unambiguous identification of cross-linked peptides and determination of cross-linked sites by MS analysis. In this work, we developed a novel cross-linking strategy using a newly designed MS-cleavable cross-linker, disuccinimidyl sulfoxide (DSSO). DSSO contains two symmetric collision-induced dissociation (CID)-cleavable sites that allow effective identification of DSSO-cross-linked peptides based on their distinct fragmentation patterns unique to cross-linking types (i.e. interlink, intralink, and dead end). The CID-induced separation of interlinked peptides in MS/MS permits MS³ analysis of single peptide chain fragment ions with defined modifications (due to DSSO remnants) for easy interpretation and unambiguous identification using existing database searching tools. Integration of data analyses from three generated data sets (MS, MS/MS, and MS³) allows high confidence identification of DSSO cross-linked peptides. The efficacy of the newly developed DSSO-based cross-linking strategy was demonstrated using model peptides and proteins. In addition, this method was successfully used for structural characterization of the yeast 20 S proteasome complex. In total, 13 non-redundant interlinked peptides of the 20 S proteasome were identified, representing the first application of an MS-cleavable cross-linker for the characterization of a multisubunit protein complex. Given its effectiveness and simplicity, this cross-linking strategy can find a broad range of applications in elucidating the structural topology of proteins and protein complexes.Proteins form stable and dynamic multisubunit complexes under different physiological conditions to maintain cell viability and normal cell homeostasis. Detailed knowledge of protein interactions and protein complex structures is fundamental to understanding how individual proteins function within a complex and how the complex functions as a whole. However, structural elucidation of large multisubunit protein complexes has been difficult because of a lack of technologies that can effectively handle their dynamic and heterogeneous nature. Traditional methods such as nuclear magnetic resonance (NMR) analysis and x-ray crystallography can yield detailed information on protein structures; however, NMR spectroscopy requires large quantities of pure protein in a specific solvent, whereas x-ray crystallography is often limited by the crystallization process.In recent years, chemical cross-linking coupled with mass spectrometry (MS) has become a powerful method for studying protein interactions (1–3). Chemical cross-linking stabilizes protein interactions through the formation of covalent bonds and allows the detection of stable, weak, and/or transient protein-protein interactions in native cells or tissues (4–9). In addition to capturing protein interacting partners, many studies have shown that chemical cross-linking can yield low resolution structural information about the constraints within a molecule (2, 3, 10) or protein complex (11–13). The application of chemical cross-linking, enzymatic digestion, and subsequent mass spectrometric and computational analyses for the elucidation of three-dimensional protein structures offers distinct advantages over traditional methods because of its speed, sensitivity, and versatility. Identification of cross-linked peptides provides distance constraints that aid in constructing the structural topology of proteins and/or protein complexes. Although this approach has been successful, effective detection and accurate identification of cross-linked peptides as well as unambiguous assignment of cross-linked sites remain extremely challenging due to their low abundance and complicated fragmentation behavior in MS analysis (2, 3, 10, 14). Therefore, new reagents and methods are urgently needed to allow unambiguous identification of cross-linked products and to improve the speed and accuracy of data analysis to facilitate its application in structural elucidation of large protein complexes.A number of approaches have been developed to facilitate MS detection of low abundance cross-linked peptides from complex mixtures. These include selective enrichment using affinity purification with biotinylated cross-linkers (15–17) and click chemistry with alkyne-tagged (18) or azide-tagged (19, 20) cross-linkers. In addition, Staudinger ligation has recently been shown to be effective for selective enrichment of azide-tagged cross-linked peptides (21). Apart from enrichment, detection of cross-linked peptides can be achieved by isotope-labeled (22–24), fluorescently labeled (25), and mass tag-labeled cross-linking reagents (16, 26). These methods can identify cross-linked peptides with MS analysis, but interpretation of the data generated from interlinked peptides (two peptides connected with the cross-link) by automated database searching remains difficult. Several bioinformatics tools have thus been developed to interpret MS/MS data and determine interlinked peptide sequences from complex mixtures (12, 14, 27–32). Although promising, further developments are still needed to make such data analyses as robust and reliable as analyzing MS/MS data of single peptide sequences using existing database searching tools (e.g. Protein Prospector, Mascot, or SEQUEST).Various types of cleavable cross-linkers with distinct chemical properties have been developed to facilitate MS identification and characterization of cross-linked peptides. These include UV photocleavable (33), chemical cleavable (19), isotopically coded cleavable (24), and MS-cleavable reagents (16, 26, 34–38). MS-cleavable cross-linkers have received considerable attention because the resulting cross-linked products can be identified based on their characteristic fragmentation behavior observed during MS analysis. Gas-phase cleavage sites result in the detection of a “reporter” ion (26), single peptide chain fragment ions (35–38), or both reporter and fragment ions (16, 34). In each case, further structural characterization of the peptide product ions generated during the cleavage reaction can be accomplished by subsequent MSⁿ1 analysis. Among these linkers, the “fixed charge” sulfonium ion-containing cross-linker developed by Lu et al. (37) appears to be the most attractive as it allows specific and selective fragmentation of cross-linked peptides regardless of their charge and amino acid composition based on their studies with model peptides.Despite the availability of multiple types of cleavable cross-linkers, most of the applications have been limited to the study of model peptides and single proteins. Additionally, complicated synthesis and fragmentation patterns have impeded most of the known MS-cleavable cross-linkers from wide adaptation by the community. Here we describe the design and characterization of a novel and simple MS-cleavable cross-linker, DSSO, and its application to model peptides and proteins and the yeast 20 S proteasome complex. In combination with new software developed for data integration, we were able to identify DSSO-cross-linked peptides from complex peptide mixtures with speed and accuracy. Given its effectiveness and simplicity, we anticipate a broader application of this MS-cleavable cross-linker in the study of structural topology of other protein complexes using cross-linking and mass spectrometry.     

Impact of Temperature Dependent Sampling Procedures in Proteomics and Peptidomics – A Characterization of the Liver and Pancreas Post Mortem Degradome

Little is known about the nature of post mortem degradation of proteins and peptides on a global level, the so-called degradome. This is especially true for nonneural tissues. Degradome properties in relation to sampling procedures on different tissues are of great importance for the studies of, for instance, post translational modifications and/or the establishment of clinical biobanks. Here, snap freezing of fresh (<2 min post mortem time) mouse liver and pancreas tissue is compared with rapid heat stabilization with regard to effects on the proteome (using two-dimensional differential in-gel electrophoresis) and peptidome (using label free liquid chromatography). We report several proteins and peptides that exhibit heightened degradation sensitivity, for instance superoxide dismutase in liver, and peptidyl-prolyl cis-trans isomerase and insulin C-peptides in pancreas. Tissue sampling based on snap freezing produces a greater amount of degradation products and lower levels of endogenous peptides than rapid heat stabilization. We also demonstrate that solely snap freezing related degradation can be attenuated by subsequent heat stabilization. We conclude that tissue sampling involving a rapid heat stabilization step is preferable to freezing with regard to proteomic and peptidomic sample quality.The evolving maturation of the field of proteomics has, in the same way as in genomics, highlighted the need of better sampling procedures and sample preparation methodologies to minimize the effect of post mortem alterations. The aspect of sample quality is not new in any way and is relevant in most biomedical fields but has only lately started to receive adequate attention. The main factors influencing sample quality is storage temperature of the body until tissue removal (foremost a problem in clinical settings and extraction of less accessible tissue samples from model organisms) and post mortem interval (PMI)¹ (1–3). Post mortem degradation in during PMI is a well known compromising problem when studying endogenous peptides (2, 3) and has also been proven to affect the results of polypeptide (here defined as proteins larger than 10 kDa) studies (3–8). PMI degradation has mainly been studied on human or mouse brain tissue, using two-dimensional electrophoresis (2-DE), SDS-PAGE, and immunoblotting (1, 3–12). There are also a few proteomic studies on muscle tissue degradation in livestock (13–16).We and others have previously explored the effect of focused microwave irradiation with regard to sample quality, demonstrating that this method is more reliable than snap freezing in liquid nitrogen, especially with regard to post-translational modification (PTM) stability (2, 3, 17–20). An alternative method based on cryostat dissection with subsequent heat treatment through boiling has also been reported to improve endogenous peptide sample quality (21). Besides focused microwave irradiation, which is specifically used for rodent brain tissue sampling, we have also demonstrated the efficiency of rapid heat stabilization through conductivity with regard to sample degradation (3, 22). Although somewhat constrained by its dependence on how quickly the tissue is harvested from the body, the latter procedure has the added advantage that it can be used on any type of tissue and species, fresh as well as frozen. This study will compare effects of sampling procedures on the liver and pancreas degradome following rapid heat stabilization, the more traditional snap freezing, or the combination of snap freezing with subsequent heat stabilization.To summarize, this study investigated the effects of post mortem degradation in pancreas and liver. Both tissues are well studied because of their multiple functions in the body and their involvement in different diseases such as diabetes or hepatocarcinoma. Pancreas is especially interesting in this context as it displays endocrine secretion of peptides, and exocrine secretion of digestive enzymes, the later making it a protease rich tissue. We used both two-dimensional difference in gel electrophoresis (2D-DIGE) and label free liquid chromatography mass spectrometry (LC-MS) based differential peptide display (2, 18), the later to better investigate changes in small molecular fragment that are not easily detectable by gel-based methods. 2D-DIGE is an unrivaled methodology to characterize alterations in isoform patterns, which is an important aspect considering that post-translational modifications (PTMs) such as phosphorylations are especially sensitive to post mortem influence within a few minutes PMI (3). The peptidomics approach has been used in several studies to point out early post mortem changes and protein degradation that tissue undergo following sampling and is therefore a well-suited method (3, 18, 22).     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Sampling From the Proteome to the Human Leukocyte Antigen-DR (HLA-DR) Ligandome Proceeds Via High Specificity

Comprehensive analysis of the complex nature of the Human Leukocyte Antigen (HLA) class II ligandome is of utmost importance to understand the basis for CD4⁺ T cell mediated immunity and tolerance. Here, we implemented important improvements in the analysis of the repertoire of HLA-DR-presented peptides, using hybrid mass spectrometry-based peptide fragmentation techniques on a ligandome sample isolated from matured human monocyte-derived dendritic cells (DC). The reported data set constitutes nearly 14 thousand unique high-confident peptides, i.e. the largest single inventory of human DC derived HLA-DR ligands to date. From a technical viewpoint the most prominent finding is that no single peptide fragmentation technique could elucidate the majority of HLA-DR ligands, because of the wide range of physical chemical properties displayed by the HLA-DR ligandome. Our in-depth profiling allowed us to reveal a strikingly poor correlation between the source proteins identified in the HLA class II ligandome and the DC cellular proteome. Important selective sieving from the sampled proteome to the ligandome was evidenced by specificity in the sequences of the core regions both at their N- and C- termini, hence not only reflecting binding motifs but also dominant protease activity associated to the endolysosomal compartments. Moreover, we demonstrate that the HLA-DR ligandome reflects a surface representation of cell-compartments specific for biological events linked to the maturation of monocytes into antigen presenting cells. Our results present new perspectives into the complex nature of the HLA class II system and will aid future immunological studies in characterizing the full breadth of potential CD4⁺ T cell epitopes relevant in health and disease.Human Leukocyte Antigen (HLA)¹ class II molecules on professional antigen presenting cells such as dendritic cells (DC) expose peptide fragments derived from exogenous and endogenous proteins to be screened by CD4⁺ T cells (1, 2). The activation and recruitment of CD4⁺ T cells recognizing disease-related peptide antigens is critical for the development of efficient antipathogen or antitumor immunity. Furthermore, the presentation of self-peptides and their interaction with CD4⁺ T cells is essential to maintain immunological tolerance and homeostasis (3). Knowledge of the nature of HLA class II-presented peptides on DC is of great importance to understand the rules of antigen processing and peptide binding motifs (4), whereas the identity of disease-related antigens may provide new knowledge on immunogenicity and leads for the development of vaccines and immunotherapy (5, 6).Mass spectrometry (MS) has proven effective for the analysis HLA class II-presented peptides (4, 7, 8). MS-based ligandome studies have demonstrated that HLA class II molecules predominantly present peptides derived from exogenous proteins that entered the cells by endocytosis and endogenous proteins that are associated with the endo-lysosomal compartments (4). Yet proteins residing in the cytosol, nucleus or mitochondria can also be presented by HLA class II molecules, primarily through autophagy (9–11). Multiple studies have mapped the HLA class II ligandome of antigen presenting cells in the context of infectious pathogens (12), autoimmune diseases (13–17) or cancer (14, 18, 19), or those that are essential for self-tolerance in the human thymus (3, 20). Notwithstanding these efforts, and certainly not in line with the extensive knowledge on the HLA class I ligandome (21), the nature of the HLA class II-presented peptide repertoire and particular its relationship to the cellular source proteome remains poorly understood.To advance our knowledge on the HLA-DR ligandome on activated DC without having to deal with limitations in cell yield from peripheral human blood (12, 21, 22) or tissue isolates (3), we explored the use of MUTZ-3 cells. This cell line has been used as a model of human monocyte-derived DCs. MUTZ-3 cells can be matured to act as antigen presenting cells and express then high levels of HLA class II molecules, and can be propagated in vitro to large cell densities (23–25). We also evaluated the performance of complementary and hybrid MS fragmentation techniques electron-transfer dissociation (ETD), electron-transfer/higher-energy collision dissociation (EThcD) (26), and higher-energy collision dissociation (HCD) to sequence and identify the HLA class II ligandome. Together this workflow allowed for the identification of an unprecedented large set of about 14 thousand unique peptide sequences presented by DC derived HLA-DR molecules, providing an in-depth view of the complexity of the HLA class II ligandome, revealing underlying features of antigen processing and surface-presentation to CD4⁺ T cells.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Proteomics Identification of Drosophila Small Interfering RNA-associated Factors

The Drosophila melanogaster RNA-induced silencing complex (RISC) forms a large ribonucleoprotein particle on small interfering RNAs (siRNAs) and catalyzes target mRNA cleavage during RNA interference (RNAi). Dicer-2, R2D2, Loquacious, and Argonaute-2 are examples of RISC-associated factors that are involved in RNAi. Holo-RISC is an ∼80 S small interfering ribonucleoprotein, which suggests that there are many additional proteins that participate in the RNAi pathway. In this study, we used siRNA affinity capture combined with mass spectrometry to identify novel components of the Drosophila RNAi machinery. Our study identified both established RISC components and novel siRNA-associated factors, many of which contain domains that are consistent with potential roles in RNAi. Functional analysis of these novel siRNA-associated proteins suggests that these factors may play an important role in RNAi.Small RNAs can regulate gene expression through a collection of mechanisms broadly termed RNA silencing. Small RNA-mediated silencing mechanisms occur in most species (1–5). The ability to silence the expression of specific genes using small RNAs via RNA interference (RNAi)¹ has greatly facilitated our understanding of gene function in eukaryotes. In addition, small RNA-mediated gene silencing has therapeutic potential and holds promise for the treatment of specific diseases (6). Understanding the mechanism of RNAi and identifying the components of the RNAi machinery are essential for harnessing its full potential in both genome-wide screens and therapeutic applications.Recently, high throughput sequencing technology has revealed the presence of endogenous siRNAs in plant, fly, worm, and mammalian cells (7–16). These endogenous siRNAs target transposable element RNAs, pseudogene RNAs, and protein-coding mRNAs (17). Therefore, the endogenous siRNA pathway seems to have evolved as a mechanism of cellular defense against selfish genetic elements. The roles of these siRNAs in development and cell physiology are poorly understood.Drosophila melanogaster is a well characterized model system for studying RNAi. In Drosophila, long double-stranded RNAs (dsRNAs) are processed by the endonuclease Dicer-2 into 21-nucleotide siRNAs (18). After processing, these siRNAs form an initiator complex with Dicer-2 and the dsRNA-binding domain (dsRBD)-containing protein R2D2 (19–23). This R2D2-Dicer-2 Initiator (RDI) complex transitions to a larger siRNP called the RISC loading complex (21, 22, 24, 25) and then to pre-RISC (26). Subsequently, pre-RISC matures into holo-RISC, which includes the catalytic activity necessary for target mRNA cleavage (21, 25, 27). The endonuclease subunit responsible for target cleavage in holo-RISC is Argonaute-2 (Ago2) (28, 29), which uses the guide strand of the siRNA duplex to target complementary mRNA sequences for cleavage and degradation.Studies of the RDI complex strongly suggest that it includes no other proteins besides Dicer-2 and R2D2 (22). Additional proteins such as Ago2 are present in pre-RISC and holo-RISC, but nonetheless the complete compositions of the RISC loading complex, pre-RISC, and holo-RISC are unknown. Furthermore, holo-RISC sediments at ∼80 S during sucrose gradient centrifugation (30). These observations indicate that additional protein factors associate with siRNAs. In this study, we identified siRNA-binding proteins from Drosophila embryo extracts. Target cleavage assays and immunoblotting of our siRNA affinity-selected proteins suggest that we purified active holo-RISC components. Proteomics analysis of the affinity matrix revealed both established and novel siRNA-associated proteins. Functional analyses of a subset of these factors suggest that they play important roles in RNAi.     

N-Terminal-oriented Proteogenomics of the Marine Bacterium Roseobacter Denitrificans Och114 using N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP) Labeling and Diagonal Chromatography

Given the ease of whole genome sequencing with next-generation sequencers, structural and functional gene annotation is now purely based on automated prediction. However, errors in gene structure are frequent, the correct determination of start codons being one of the main concerns. Here, we combine protein N termini derivatization using (N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP Ac-OSu) as a labeling reagent with the COmbined FRActional DIagonal Chromatography (COFRADIC) sorting method to enrich labeled N-terminal peptides for mass spectrometry detection. Protein digestion was performed in parallel with three proteases to obtain a reliable automatic validation of protein N termini. The analysis of these N-terminal enriched fractions by high-resolution tandem mass spectrometry allowed the annotation refinement of 534 proteins of the model marine bacterium Roseobacter denitrificans OCh114. This study is especially efficient regarding mass spectrometry analytical time. From the 534 validated N termini, 480 confirmed existing gene annotations, 41 highlighted erroneous start codon annotations, five revealed totally new mis-annotated genes; the mass spectrometry data also suggested the existence of multiple start sites for eight different genes, a result that challenges the current view of protein translation initiation. Finally, we identified several proteins for which classical genome homology-driven annotation was inconsistent, questioning the validity of automatic annotation pipelines and emphasizing the need for complementary proteomic data. All data have been deposited to the ProteomeXchange with identifier PXD000337.Recent developments in mass spectrometry and bioinformatics have established proteomics as a common and powerful technique for identifying and quantifying proteins at a very broad scale, but also for characterizing their post-translational modifications and interaction networks (1, 2). In addition to the avalanche of proteomic data currently being reported, many genome sequences are established using next-generation sequencing, fostering proteomic investigations of new cellular models. Proteogenomics is a relatively recent field in which high-throughput proteomic data is used to verify coding regions within model genomes to refine the annotation of their sequences (2–8). Because genome annotation is now fully automated, the need for accurate annotation for model organisms with experimental data is crucial. Many projects related to genome re-annotation of microorganisms with the help of proteomics have been recently reported, such as for Mycoplasma pneumoniae (9), Rhodopseudomonas palustris (10), Shewanella oneidensis (11), Thermococcus gammatolerans (12), Deinococcus deserti (13), Salmonella thyphimurium (14), Mycobacterium tuberculosis (15, 16), Shigella flexneri (17), Ruegeria pomeroyi (18), and Candida glabrata (19), as well as for higher organisms such as Anopheles gambiae (20) and Arabidopsis thaliana (4, 5).The most frequently reported problem in automatic annotation systems is the correct identification of the translational start codon (21–23). The error rate depends on the primary annotation system, but also on the organism, as reported for Halobacterium salinarum and Natromonas pharaonis (24), Deinococcus deserti (21), and Ruegeria pomeroyi (18), where the error rate is estimated above 10%. Identification of a correct translational start site is essential for the genetic and biochemical analysis of a protein because errors can seriously impact subsequent biological studies. If the N terminus is not correctly identified, the protein will be considered in either a truncated or extended form, leading to errors in bioinformatic analyses (e.g. during the prediction of its molecular weight, isoelectric point, cellular localization) and major difficulties during its experimental characterization. For example, a truncated protein may be heterologously produced as an unfolded polypeptide recalcitrant to structure determination (25). Moreover, N-terminal modifications, which are poorly documented in annotation databases, may occur (26, 27).Unfortunately, the poor polypeptide sequence coverage obtained for the numerous low abundance proteins in current shotgun MS/MS proteomic studies implies that the overall detection of N-terminal peptides obtained in proteogenomic studies is relatively low. Different methods for establishing the most extensive list of protein N termini, grouped under the so-called “N-terminomics” theme, have been proposed to selectively enrich or improve the detection of these peptides (2, 28, 29). Large N-terminome studies have recently been reported based on resin-assisted enrichment of N-terminal peptides (30) or terminal amine isotopic labeling of substrates (TAILS) coupled to depletion of internal peptides with a water-soluble aldehyde-functionalized polymer (31–35). Among the numerous N-terminal-oriented methods (2), specific labeling of the N terminus of intact proteins with N-tris(2,4,6-trimethoxyphenyl)phosphonium acetyl succinamide (TMPP-Ac-OSu)¹ has proven reliable (21, 36–39). TMPP-derivatized N-terminal peptides have interesting properties for further LC-MS/MS mass spectrometry: (1) an increase in hydrophobicity because of the trimethoxyphenyl moiety added to the peptides, increasing their retention times in reverse phase chromatography, (2) improvement of their ionization because of the introduction of a positively charged group, and (3) a much simpler fragmentation pattern in tandem mass spectrometry. Other reported approaches rely on acetylation, followed by trypsin digestion, and then biotinylation of free amino groups (40); guanidination of lysine lateral chains followed by N-biotinylation of the N termini and trypsin digestion (41); or reductive amination of all free amino groups with formaldehyde preceeding trypsin digestion (42). Recently, we applied the TMPP method to the proteome of the Deinococcus deserti bacterium isolated from upper sand layers of the Sahara desert (13). This method enabled the detection of N-terminal peptides allowing the confirmation of 278 translation initiation codons, the correction of 73 translation starts, and the identification of non-canonical translation initiation codons (21). However, most TMPP-labeled N-terminal peptides are hidden among the more abundant internal peptides generated after proteolysis of a complex proteome, precluding their detection. This results in disproportionately fewer N-terminal validations, that is, 5 and 8% of total polypeptides coded in the theoretical proteomes of Mycobacterium smegmatis (37) and Deinococcus deserti (21) with a total of 342 and 278 validations, respectively.An interesting chromatographic method to fractionate peptide mixtures for gel-free high-throughput proteome analysis has been developed over the last years and applied to various topics (43, 44). This technique, known as COmbined FRActional DIagonal Chromatography (COFRADIC), uses a double chromatographic separation with a chemical reaction in between to change the physico-chemical properties of the extraneous peptides to be resolved from the peptides of interest. Its previous applications include the separation of methionine-containing peptides (43), N-terminal peptide enrichment (45, 46), sulfur amino acid-containing peptides (47), and phosphorylated peptides (48). COFRADIC was identified as the best method for identification of N-terminal peptides of two archaea, resulting in the identification of 240 polypeptides (9% of the theoretical proteome) for Halobacterium salinarum and 220 (8%) for Natronomonas pharaonis (24).Taking advantage of both the specificity of TMPP labeling, the resolving power of COFRADIC for enrichment, and the increase in information through the use of multiple proteases, we performed the proteogenomic analysis of a marine bacterium from the Roseobacter clade, namely Roseobacter denitrificans OCh114. This novel approach allowed us to validate and correct 534 unique proteins (13% of the theoretical proteome) with TMPP-labeled N-terminal signatures obtained using high-resolution tandem mass spectrometry. We corrected 41 annotations and detected five new open reading frames in the R. denitrificans genome. We further identified eight distinct proteins showing direct evidence for multiple start sites.     

Development and Evaluation of Normalization Methods for Label-free Relative Quantification of Endogenous Peptides

The performances of 10 different normalization methods on data of endogenous brain peptides produced with label-free nano-LC-MS were evaluated. Data sets originating from three different species (mouse, rat, and Japanese quail), each consisting of 35–45 individual LC-MS analyses, were used in the study. Each sample set contained both technical and biological replicates, and the LC-MS analyses were performed in a randomized block fashion. Peptides in all three data sets were found to display LC-MS analysis order-dependent bias. Global normalization methods will only to some extent correct this type of bias. Only the novel normalization procedure RegrRun (linear regression followed by analysis order normalization) corrected for this type of bias. The RegrRun procedure performed the best of the normalization methods tested and decreased the median S.D. by 43% on average compared with raw data. This method also produced the smallest fraction of peptides with interblock differences while producing the largest fraction of differentially expressed peaks between treatment groups in all three data sets. Linear regression normalization (Regr) performed second best and decreased median S.D. by 38% on average compared with raw data. All other examined methods reduced median S.D. by 20–30% on average compared with raw data.Peptidomics is defined as the analysis of the peptide content within an organism, tissue, or cell (1–3). The proteome and peptidome have common features, but there are also prominent differences. Proteomics generally identifies proteins by using the information of biologically inactive peptides derived from tryptic digestion, whereas peptidomics tries to identify endogenous peptides using single peptide sequence information only (4). Endogenous neuropeptides are peptides used for intracellular signaling that can act as neurotransmitters or neuromodulators in the nervous system. These polypeptides of 3–100 amino acids can be abundantly produced in large neural populations or in trace levels from single neurons (5) and are often generated through the cleavage of precursor proteins. However, unwanted peptides can also be created through post-mortem induced proteolysis (6). The later aspect complicates the technical analysis of neuropeptides as post-mortem conditions increase the number of degradation peptides. The possibility to detect, identify, and quantify lowly expressed neuropeptides using label-free LC-MS techniques has improved with the development of new sample preparation techniques including rapid heating of the tissue, which prevents protein degradation and inhibition of post-mortem proteolytic activity (7, 8).It has been suggested by us (4, 5) and others (9) that comparing the peptidome between samples of e.g. diseased and normal tissue may lead to the discovery of biologically relevant peptides of certain pathological or pharmacological events. However, differences in relative peptide abundance measurements may not only originate from biological differences but also from systematic bias and noise. To reduce the effects of experimentally induced variability it is common to normalize the raw data. This is a concept well known in the area of genomics studies using gene expression microarrays (10–12). As a consequence, many methods developed for microarray data have also been adapted for normalizing peptide data produced with LC-MS techniques (10–16). Normally the underlying assumption for applying these techniques is that the total or mean/median peak abundances should be equal across different experiments, in this case between LC-MS analyses. Global normalization methods refer to cases where all peak abundances are used to determine a single normalization factor between experiments (13, 15, 16), a subset of peaks assumed to be similarly abundant between experiments (16) is used, or spiked-in peptides are used as internal standards. In a study by Callister et al. (14), normalization methods for tryptic LC-FTICR-MS peptide data were compared. The authors concluded that global or iterative linear regression works best in most cases but also recommended that the best procedure should be selected for each data set individually. Methods used for normalizing LC-MS data have been reviewed previously (14, 17, 18), but to our knowledge only Callister et al. (14) have used small data sets to systematically evaluate such methods. None of these studies have targeted data of endogenous peptides.In this study, the effects of 10 different normalization methods were evaluated on data produced by a nano-LC system coupled to an electrospray Q-TOF or linear trap quadrupole (LTQ)¹ mass spectrometer. Normalization methods that originally were developed for gene expression data were used, and one novel method, linear regression followed by analysis order normalization (RegrRun), is presented. The normalization methods were evaluated using three data sets of endogenous brain peptides originating from three different species (mouse, rat, and Japanese quail), each consisting of 35–45 individual LC-MS analyses. Each data set contained both technical and biological replicates.