Flexible and accessible workflows for improved proteogenomic analysis using the Galaxy framework.

Proteogenomics combines large-scale genomic and transcriptomic data with mass-spectrometry-based proteomic data to discover novel protein sequence variants and improve genome annotation. In contrast with conventional proteomic applications, proteogenomic analysis requires a number of additional data processing steps. Ideally, these required steps would be integrated and automated via a single software platform offering accessibility for wet-bench researchers as well as flexibility for user-specific customization and integration of new software tools as they emerge. Toward this end, we have extended the Galaxy bioinformatics framework to facilitate proteogenomic analysis. Using analysis of whole human saliva as an example, we demonstrate Galaxy's flexibility through the creation of a modular workflow incorporating both established and customized software tools that improve depth and quality of proteogenomic results. Our customized Galaxy-based software includes automated, batch-mode BLASTP searching and a Peptide Sequence Match Evaluator tool, both useful for evaluating the veracity of putative novel peptide identifications. Our complex workflow (approximately 140 steps) can be easily shared using built-in Galaxy functions, enabling their use and customization by others. Our results provide a blueprint for the establishment of the Galaxy framework as an ideal solution for the emerging field of proteogenomics.

Insulin signalling to mTOR mediated by the Akt/PKB substrate PRAS40.

Insulin stimulates protein synthesis and cell growth by activation of the protein kinases Akt (also known as protein kinase B, PKB) and mammalian target of rapamycin (mTOR). It was reported that Akt activates mTOR by phosphorylation and inhibition of tuberous sclerosis complex 2 (TSC2). However, in recent studies the physiological requirement of Akt phosphorylation of TSC2 for mTOR activation has been questioned. Here, we identify PRAS40 (proline-rich Akt/PKB substrate 40 kDa) as a novel mTOR binding partner that mediates Akt signals to mTOR. PRAS40 binds the mTOR kinase domain and its interaction with mTOR is induced under conditions that inhibit mTOR signalling, such as nutrient or serum deprivation or mitochondrial metabolic inhibition. Binding of PRAS40 inhibits mTOR activity and suppresses constitutive activation of mTOR in cells lacking TSC2. PRAS40 silencing inactivates insulin-receptor substrate-1 (IRS-1) and Akt, and uncouples the response of mTOR to Akt signals. Furthermore, PRAS40 phosphorylation by Akt and association with 14-3-3, a cytosolic anchor protein, are crucial for insulin to stimulate mTOR. These findings identify PRAS40 as an important regulator of insulin sensitivity of the Akt-mTOR pathway and a potential target for the treatment of cancers, insulin resistance and hamartoma syndromes.

Deep metaproteomic analysis of human salivary supernatant.

The human salivary proteome is extremely complex, including proteins from salivary glands, serum, and oral microbes. Much has been learned about the host component, but little is known about the microbial component. Here we report a metaproteomic analysis of salivary supernatant pooled from six healthy subjects. For deep interrogation of the salivary proteome, we combined protein dynamic range compression (DRC), multidimensional peptide fractionation, and high-mass accuracy MS/MS with a novel two-step peptide identification method using a database of human proteins plus those translated from oral microbe genomes. Peptides were identified from 124 microbial species as well as uncultured phylotypes such as TM7. Streptococcus, Rothia, Actinomyces, Prevotella, Neisseria, Veilonella, Lactobacillus, Selenomonas, Pseudomonas, Staphylococcus, and Campylobacter were abundant among the 65 genera from 12 phyla represented. Taxonomic diversity in our study was broadly consistent with metagenomic studies of saliva. Proteins mapped to 20 KEGG pathways, with carbohydrate metabolism, amino acid metabolism, energy metabolism, translation, membrane transport, and signal transduction most represented. The communities sampled appear to be actively engaged in glycolysis and protein synthesis. This first deep metaproteomic catalog from human salivary supernatant provides a baseline for future studies of shifts in microbial diversity and protein activities potentially associated with oral disease.

Workflow for analysis of high mass accuracy salivary data set using MaxQuant and ProteinPilot search algorithm.

LTQ Orbitrap data analyzed with ProteinPilot can be further improved by MaxQuant raw data processing, which utilizes precursor-level high mass accuracy data for peak processing and MGF creation. In particular, ProteinPilot results from MaxQuant-processed peaklists for Orbitrap data sets resulted in improved spectral utilization due to an improved peaklist quality with higher precision and high precursor mass accuracy (HPMA). The output and postsearch analysis tools of both workflows were utilized for previously unexplored features of a three-dimensional fractionated and hexapeptide library (ProteoMiner) treated whole saliva data set comprising 200 fractions. ProteinPilot's ability to simultaneously predict multiple modifications showed an advantage from ProteoMiner treatment for modified peptide identification. We demonstrate that complementary approaches in the analysis pipeline provide comprehensive results for the whole saliva data set acquired on an LTQ Orbitrap. Overall our results establish a workflow for improved protein identification from high mass accuracy data.

Quantitative nuclear proteomics identifies mTOR regulation of DNA damage response.

Cellular nutritional and energy status regulates a wide range of nuclear processes important for cell growth, survival, and metabolic homeostasis. Mammalian target of rapamycin (mTOR) plays a key role in the cellular responses to nutrients. However, the nuclear processes governed by mTOR have not been clearly defined. Using isobaric peptide tagging coupled with linear ion trap mass spectrometry, we performed quantitative proteomics analysis to identify nuclear processes in human cells under control of mTOR. Within 3 h of inhibiting mTOR with rapamycin in HeLa cells, we observed down-regulation of nuclear abundance of many proteins involved in translation and RNA modification. Unexpectedly, mTOR inhibition also down-regulated several proteins functioning in chromosomal integrity and up-regulated those involved in DNA damage responses (DDRs) such as 53BP1. Consistent with these proteomic changes and DDR activation, mTOR inhibition enhanced interaction between 53BP1 and p53 and increased phosphorylation of ataxia telangiectasia mutated (ATM) kinase substrates. ATM substrate phosphorylation was also induced by inhibiting protein synthesis and suppressed by inhibiting proteasomal activity, suggesting that mTOR inhibition reduces steady-state (abundance) levels of proteins that function in cellular pathways of DDR activation. Finally, rapamycin-induced changes led to increased survival after radiation exposure in HeLa cells. These findings reveal a novel functional link between mTOR and DDR pathways in the nucleus potentially operating as a survival mechanism against unfavorable growth conditions.

Large-scale phosphoproteomics analysis of whole saliva reveals a distinct phosphorylation pattern.

In-depth knowledge of bodily fluid phosphoproteomes, such as whole saliva, is limited. To better understand the whole saliva phosphoproteome, we generated a large-scale catalog of phosphorylated proteins. To circumvent the wide dynamic range of phosphoprotein abundance in whole saliva, we combined dynamic range compression using hexapeptide beads, strong cation exchange HPLC peptide fractionation, and immobilized metal affinity chromatography prior to mass spectrometry. In total, 217 unique phosphopeptides sites were identified representing 85 distinct phosphoproteins at 2.3% global FDR. From these peptides, 129 distinct phosphorylation sites were identified of which 57 were previously known, but only 11 of which had been previously identified in whole saliva. Cellular localization analysis revealed salivary phosphoproteins had a distribution similar to all known salivary proteins, but with less relative representation in "extracellular" and "plasma membrane" categories compared to salivary glycoproteins. Sequence alignment showed that phosphorylation occurred at acidic-directed kinase, proline-directed, and basophilic motifs. This differs from plasma phosphoproteins, which predominantly occur at Golgi casein kinase recognized sequences. Collectively, these results suggest diverse functions for salivary phosphoproteins and multiple kinases involved in their processing and secretion. In all, this study should lay groundwork for future elucidation of the functions of salivary protein phosphorylation.

Hexapeptide libraries for enhanced protein PTM identification and relative abundance profiling in whole human saliva.

Dynamic range compression (DRC) by hexapeptide libraries increases MS/MS-based identification of lower-abundance proteins in complex mixtures. However, two unanswered questions impede fully realizing DRC's potential in shotgun proteomics. First, does DRC enhance identification of post-translationally modified proteins? Second, can DRC be incorporated into a workflow enabling relative protein abundance profiling? We sought to answer both questions analyzing human whole saliva. Addressing question one, we coupled DRC with covalent glycopeptide enrichment and MS/MS. With DRC we identified ∼2 times more N-linked glycoproteins and their glycosylation sites than without DRC, dramatically increasing the known salivary glycoprotein catalog. Addressing question two, we compared differentially stable isotope-labeled saliva samples pooled from healthy and metastatic breast cancer women using a multidimensional peptide fractionation-based workflow, analyzing in parallel one sample portion with DRC and one portion without. Our workflow categorizes proteins with higher absolute abundance, whose relative abundance ratios are altered by DRC, from proteins of lower absolute abundance detected only after DRC. Within each of these salivary protein categories, we identified novel abundance changes putatively associated with breast cancer, demonstrating feasibility and benefits of DRC for relative abundance profiling. Collectively, our results bring us closer to realizing the full potential of DRC for proteomic studies.

Performance of preclinical models in predicting drug-induced liver injury in humans: a systematic review.

Drug-induced liver injury (DILI) causes one in three market withdrawals due to adverse drug reactions, causing preventable human suffering and massive financial loss. We applied evidence-based methods to investigate the role of preclinical studies in predicting human DILI using two anti-diabetic drugs from the same class, but with different toxicological profiles: troglitazone (withdrawn from US market due to DILI) and rosiglitazone (remains on US market). Evidence Stream 1: A systematic literature review of in vivo studies on rosiglitazone or troglitazone was conducted (PROSPERO registration CRD42018112353). Evidence Stream 2: in vitro data on troglitazone and rosiglitazone were retrieved from the US EPA ToxCast database. Evidence Stream 3: troglitazone- and rosiglitazone-related DILI cases were retrieved from WHO Vigibase. All three evidence stream analyses were conducted according to evidence-based methodologies and performed according to pre-registered protocols. Evidence Stream 1: 9288 references were identified, with 42 studies included in analysis. No reported biomarker for either drug indicated a strong hazard signal in either preclinical animal or human studies. All included studies had substantial limitations, resulting in "low" or "very low" certainty in findings. Evidence Stream 2: Troglitazone was active in twice as many in vitro assays (129) as rosiglitazone (60), indicating a strong signal for more off-target effects. Evidence Stream 3: We observed a fivefold difference in both all adverse events and liver-related adverse events reported, and an eightfold difference in fatalities for troglitazone, compared to rosiglitazone. In summary, published animal and human trials failed to predict troglitazone's potential to cause severe liver injury in a wider patient population, while in vitro data showed marked differences in the two drugs' off-target activities, offering a new paradigm for reducing drug attrition in late development and in the market. This investigation concludes that death and disability due to adverse drug reactions may be prevented if mechanistic information is deployed at early stages of drug development by pharmaceutical companies and is considered by regulators as a part of regulatory submissions.

Novel In Situ Collection of Tumor Interstitial Fluid from a Head and Neck Squamous Carcinoma Reveals a Unique Proteome with Diagnostic Potential.

INTRODUCTION: Tumors lack normal drainage of secreted fluids and consequently build up tumor interstitial fluid (TIF). Unlike other bodily fluids, TIF likely contains a high proportion of tumor-specific proteins with potential as biomarkers. METHODS: Here, we evaluated a novel technique using a unique ultrafiltration catheter for in situ collection of TIF and used it to generate the first catalog of TIF proteins from a head and neck squamous cell carcinoma (HNSCC). To maximize proteomic coverage, TIF was immunodepleted for high abundance proteins and digested with trypsin, and peptides were fractionated in three dimensions prior to mass spectrometry. RESULTS: We identified 525 proteins with high confidence. The HNSCC TIF proteome was distinct compared to proteomes of other bodily fluids. It contained a relatively high proportion of proteins annotated by Gene Ontology as "extracellular" compared to other secreted fluid and cellular proteomes, indicating minimal cell lysis from our in situ collection technique. Several proteins identified are putative biomarkers of HNSCC, supporting our catalog's value as a source of potential biomarkers. CONCLUSIONS: In all, we demonstrate a reliable new technique for in situ TIF collection and provide the first HNSCC TIF protein catalog with value as a guide for others seeking to develop tumor biomarkers. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12014-010-9050-3) contains supplementary material, which is available to authorized users.

A dynamic range compression and three-dimensional peptide fractionation analysis platform expands proteome coverage and the diagnostic potential of whole saliva.

Comprehensive identification of proteins in whole human saliva is critical for appreciating its full diagnostic potential. However, this is challenged by the large dynamic range of protein abundance within the fluid. To address this problem, we used an analysis platform that coupled hexapeptide libraries for dynamic range compression (DRC) with three-dimensional (3D) peptide fractionation. Our approach identified 2340 proteins in whole saliva and represents the largest saliva proteomic dataset generated using a single analysis platform. Three-dimensional peptide fractionation involving sequential steps of preparative isoelectric focusing (IEF), strong cation exchange, and capillary reversed-phase liquid chromatography was essential for maximizing gains from DRC. Compared to saliva not treated with hexapeptide libraries, DRC substantially increased identified proteins across physicochemical and functional categories. Approximately 20% of total salivary proteins are also seen in plasma, and proteins in both fluids show comparable functional diversity and disease-linkage. However, for a subset of diseases, saliva has higher apparent diagnostic potential. These results expand the potential for whole saliva in health monitoring/diagnostics and provide a general platform for improving proteomic coverage of complex biological samples.

Advancing mammalian cell culture engineering using genome-scale technologies.

Mammalian cell-derived protein therapeutic production has changed the landscape of human healthcare in the past two decades. The importance of protein therapeutics has motivated the search for more cost-effective and efficient cell lines capable of producing high quality protein products. The factors contributing to optimal producer cell lines are often complex, and not simply conferred by one gene or gene product, which makes an understanding of system-wide properties for better engineering of optimized cell lines essential. Genome-scale technologies (genomics, transcriptomics and proteomics) enable such engineering studies. However, the use of these technologies in cell culture engineering is still in its infancy. Here, we summarize current knowledge of cell properties important for the design of efficient protein-producing mammalian cell lines, and highlight relevant studies to-date that use genome-scale technologies in these cell systems. We also provide a focused review of relevant alternative and emerging technologies, which have seen limited use in cell culture engineering, but hold great potential for significant advancements in protein therapeutic production.

Caenorhabditis elegans OSR-1 regulates behavioral and physiological responses to hyperosmotic environments.

The molecular mechanisms that enable multicellular organisms to sense and modulate their responses to hyperosmotic environments are poorly understood. Here, we employ Caenorhabditis elegans to characterize the response of a multicellular organism to osmotic stress and establish a genetic screen to isolate mutants that are osmotic stress resistant (OSR). In this study, we describe the cloning of a novel gene, osr-1, and demonstrate that it regulates osmosensation, adaptation, and survival in hyperosmotic environments. Whereas wild-type animals exposed to hyperosmotic conditions rapidly lose body volume, motility, and viability, osr-1(rm1) mutant animals maintain normal body volume, motility, and viability even upon chronic exposures to high osmolarity environments. In addition, osr-1(rm1) animals are specifically resistant to osmotic stress and are distinct from previously characterized osmotic avoidance defective (OSM) and general stress resistance age-1(hx546) mutants. OSR-1 is expressed in the hypodermis and intestine, and expression of OSR-1 in hypodermal cells rescues the osr-1(rm1) phenotypes. Genetic epistasis analysis indicates that OSR-1 regulates survival under osmotic stress via CaMKII and a conserved p38 MAP kinase signaling cascade and regulates osmotic avoidance and resistance to acute dehydration likely by distinct mechanisms. We suggest that OSR-1 plays a central role in integrating stress detection and adaptation responses by invoking multiple signaling pathways to promote survival under hyperosmotic environments.

In-gel peptide IEF sample preparation for LC/MS analysis.

The technique of proteolytically digesting a sample and identifying its protein components by liquid chromatography followed by mass spectrometry (LC-MS) is a widely used analytical tool. Prior fractionation by isoelectric focusing (IEF) may be performed to increase the depth of proteome coverage. Here, we describe a method for in-gel IEF separation of a proteolytic digest that utilizes commercially available immobilized pH gradient (IPG) strips and a widely used IEF instrument.

2-D Western blotting for evaluation of antibodies developed for detection of host cell protein.

Recombinant proteins generated for therapeutic use must be substantially free of residual host cell protein (HCP). The presence of host cell protein (HCP) is usually assayed by ELISA using a polyclonal antibody mixture raised against a population of proteins derived from the host cell background. This antibody should recognize as high a proportion as possible of the potential HCPs in a given sample. A recommended method for evaluating the assay involves two-dimensional electrophoretic separation followed by Western blotting.We present here a method using commercial anti-HCP antibody and samples derived from Chinese Hamster Ovary (CHO) cells. The 2-D electrophoresis procedure gives highly reproducible spot patterns and entire procedure can be completed in less than 2 days. Software analysis enables the straightforward generation of percent coverage values for the antibody when used to probe HCP-containing samples.

Genetic interactions among ZDS1,2, CDC37, and protein kinase CK2 in Saccharomyces cerevisiae.

We report here the identification of the homologous gene pair ZDS1,2 as multicopy suppressors of a temperature-sensitive allele (cka2-13(ts)) of the CKA2 gene encoding the alpha' catalytic subunit of protein kinase CK2. Overexpression of ZDS1,2 suppressed the temperature sensitivity, geldanamycin (GA) sensitivity, slow growth, and flocculation of multiple cka2 alleles and enhanced CK2 activity in vivo toward a known physiological substrate, Fpr3. Consistent with the existence of a recently described positive feedback loop between CK2 and Cdc37, overexpression of ZDS1,2 also suppressed the temperature sensitivity, abnormal morphology, and GA sensitivity of a CK2 phosphorylation-deficient mutant of CDC37, cdc37-S14A, as well as the GA sensitivity of a cdc37-1 allele. A likely basis for all of these effects is our observation that ZDS1,2 overexpression enhances Cdc37 protein levels. Activation of the positive feedback loop between CK2 and Cdc37 likely contributes to the pleiotropic nature of ZDS1,2, as both CK2 and Cdc37 regulate diverse cellular functions.

Three-dimensional peptide fractionation for highly sensitive nanoscale LC-based shotgun proteomic analysis of complex protein mixtures.

To sensitively analyze complex protein mixtures by mass spectrometry-based shotgun proteomics, researchers have employed platforms that couple orthogonal peptide fractionation methods using nanoscale HPLC. Commonly used platforms have coupled either strong cation exchange (SCX) HPLC or preparative isoelectric focusing (IEF) with nanoscale reversed-phase (nanoRP) HPLC fractionation of peptides. Coupling two dimensions of peptide fractionation, prior to mass spectrometric analysis, increases the sensitivity for identifying low abundance proteins. However, the large dynamic range of protein abundance and high level of complexity of protein mixtures derived from many biological sources, such as bodily fluids, require additional steps of peptide fractionation. To address this shortcoming, we have developed a platform combining three dimensions of peptide fractionation as follows: (1) preparative IEF; (2) SCX HPLC; and (3) nanoRP HPLC. This platform significantly increases the sensitivity of shotgun proteomic analysis in complex protein mixtures. Here, we describe the implementation of this three-dimensional peptide fractionation platform for proteomic studies of complex mixtures.

Hsf1 activation inhibits rapamycin resistance and TOR signaling in yeast revealed by combined proteomic and genetic analysis.

TOR kinases integrate environmental and nutritional signals to regulate cell growth in eukaryotic organisms. Here, we describe results from a study combining quantitative proteomics and comparative expression analysis in the budding yeast, S. cerevisiae, to gain insights into TOR function and regulation. We profiled protein abundance changes under conditions of TOR inhibition by rapamycin treatment, and compared this data to existing expression information for corresponding gene products measured under a variety of conditions in yeast. Among proteins showing abundance changes upon rapamycin treatment, almost 90% of them demonstrated homodirectional (i.e., in similar direction) transcriptomic changes under conditions of heat/oxidative stress. Because the known downstream responses regulated by Tor1/2 did not fully explain the extent of overlap between these two conditions, we tested for novel connections between the major regulators of heat/oxidative stress response and the TOR pathway. Specifically, we hypothesized that activation of regulator(s) of heat/oxidative stress responses phenocopied TOR inhibition and sought to identify these putative TOR inhibitor(s). Among the stress regulators tested, we found that cells (hsf1-R206S, F256S and ssa1-3 ssa2-2) constitutively activated for heat shock transcription factor 1, Hsf1, inhibited rapamycin resistance. Further analysis of the hsf1-R206S, F256S allele revealed that these cells also displayed multiple phenotypes consistent with reduced TOR signaling. Among the multiple Hsf1 targets elevated in hsf1-R206S, F256S cells, deletion of PIR3 and YRO2 suppressed the TOR-regulated phenotypes. In contrast to our observations in cells activated for Hsf1, constitutive activation of other regulators of heat/oxidative stress responses, such as Msn2/4 and Hyr1, did not inhibit TOR signaling. Thus, we propose that activated Hsf1 inhibits rapamycin resistance and TOR signaling via elevated expression of specific target genes in S. cerevisiae. Additionally, these results highlight the value of comparative expression analyses between large-scale proteomic and transcriptomic datasets to reveal new regulatory connections.

Proteomic mapping of 4-hydroxynonenal protein modification sites by solid-phase hydrazide chemistry and mass spectrometry.

The modification of proteins by the cytotoxic, reactive aldehyde 4-hydroxynonenal (HNE) is known to alter protein function and impair cellular mechanisms. In order to identify susceptible amino acid sites of HNE modification within complex biological mixtures by microcapillary liquid chromatography and linear ion trap tandem mass spectrometry, we have developed a solid-phase capture and release strategy that utilizes reversible hydrazide chemistry to enrich HNE-modified peptides. To maximize the detection of fragment ions diagnostic of HNE modification, both neutral loss-dependent acquisition of MS/MS/MS spectra and the pulsed Q dissociation operation mode were employed. When the solid-phase hydrazide enrichment strategy was applied to a yeast lysate treated with HNE, 125 distinct amino acid sites of HNE modification were mapped on 67 different proteins. The endogenous susceptibility of many of these proteins to HNE modification was demonstrated by analyzing HNE-treated yeast cell cultures with a complementary biotin hydrazide enrichment strategy. Further analysis revealed that the majority of amino acid sites susceptible to HNE modification were histidine residues, with most of these sites being flanked by basic amino acid residues, and predicted to be solvent exposed. These results demonstrate the effectiveness of this novel strategy as a general platform for proteome-scale identification of amino acid sites susceptible to HNE modification from within complex mixtures.

Preparative peptide isoelectric focusing as a tool for improving the identification of lysine-acetylated peptides from complex mixtures.

Protein sequence database searching of tandem mass spectrometry data is commonly employed to identify post-translational modifications (PTMs) to peptides in global proteomic studies. In these studies, the accurate identification of these modified peptides relies on strategies to ensure high-confidence results from sequence database searching in which differential mass shift parameters are employed to identify PTMs to specific amino acids. Using lysine acetylation as an example PTM, we have observed that the inclusion of differential modification information in sequence database searching dramatically increases the potential for false-positive sequence matches to modified peptides, making the confident identification of true sequence matches difficult. In a proof-of-principle study of whole cell yeast lysates, we demonstrate the combination of preparative isoelectric focusing using free-flow electrophoresis, and an adjusted peptide isoelectric point prediction algorithm, as an effective means to increase the confidence of lysine-acetylated peptide identification. These results demonstrate the potential utility of this general strategy for improving the identification of PTMs which cause a shift to the intrinsic isoelectric point of peptides.

iTRAQ reagent-based quantitative proteomic analysis on a linear ion trap mass spectrometer.

For proteomic analysis using tandem mass spectrometry, linear ion trap instruments provide unsurpassed sensitivity but unreliably detect low mass peptide fragments, precluding their use with iTRAQ reagent-labeled samples. Although the popular LTQ linear ion trap supports analyzing iTRAQ reagent-labeled peptides via pulsed Q dissociation, PQD, its effectiveness remains questionable. Using a standard mixture, we found careful tuning of relative collision energy necessary for fragmenting iTRAQ reagent-labeled peptides, and increasing microscan acquisition and repeat count improves quantification but identifies somewhat fewer peptides. We developed software to calculate abundance ratios via summing reporter ion intensities across spectra matching to each protein, thereby providing maximized accuracy. Testing found that results closely corresponded between analysis using optimized LTQ-PQD settings plus our software and using a Qstar instrument. Thus, we demonstrate the effectiveness of LTQ-PQD analyzing iTRAQ reagent-labeled peptides, and provide guidelines for successful quantitative proteomic studies.