Proteomics Informatics Databases, data repositories and - PowerPoint PPT Presentation

Proteomics Informatics – Databases, data repositories and standardization (Week 7)

Protein Sequence Databases

RefSeq Distinguishing Features of the RefSeq collection include: • non-redundancy • explicitly linked nucleotide and protein sequences • updates to reflect current knowledge of sequence data and biology • data validation and format consistency • ongoing curation by NCBI staff and collaborators, with reviewed records indicated http://www.ncbi.nlm.nih.gov/books/NBK21091/

Ensembl • genome information for sequenced chordate genomes. • evidenced-based gene sets for all supported species • large-scale whole genome multiple species alignments across vertebrates • variation data resources for 17 species and regulation annotations based on ENCODE and other data sets. http://www.ensembl.org/

UniProt The mission of UniProt is to provide the scientific community with a comprehensive, high-quality and freely accessible resource of protein sequence and functional information. http://www.uniprot.org/

Species-Centric Consortia For some organisms, there are consortia that provide high-quality databases: Yeast (http://yeastgenome.org/) Fly (http://flybase.org/) Arabidopsis (http://arabidopsis.org/)

FASTA RefSeq: >gi|168693669|ref|NP_001108231.1| zinc finger protein 683 [Homo sapiens] MKEESAAQLGCCHRPMALGGTGGSLSPSLDFQLFRGDQVFSACRPLPDMVDAHGPSCASWLCPLPLAPGRSALLACLQDL DLNLCTPQPAPLGTDLQGLQEDALSMKHEPPGLQASSTDDKKFTVKYPQNKDKLGKQPERAGEGAPCPAFSSHNSSSPPP LQNRKSPSPLAFCPCPPVNSISKELPFLLHAFYPGYPLLLPPPHLFTYGALPSDQCPHLLMLPQDPSYPTMAMPSLLMMV NELGHPSARWETLLPYPGAFQASGQALPSQARNPGAGAAPTDSPGLERGGMASPAKRVPLSSQTGTAALPYPLKKKNGKI LYECNICGKSFGQLSNLKVHLRVHSGERPFQCALCQKSFTQLAHLQKHHLVHTGERPHKCSVCHKRFSSSSNLKTHLRLH SGARPFQCSVCRSRFTQHIHLKLHHRLHAPQPCGLVHTQLPLASLACLAQWHQGALDLMAVASEKHMGYDIDEVKVSSTS QGKARAVSLSSAGTPLVMGQDQNN Ensembl: >ENSMUSP00000131420 pep:known supercontig:NCBIM37:NT_166407:104574:105272 gene:ENSMUSG00000092057 transcript:ENSMUST00000167991 MFSLMKKRRRKSSSNTLRNIVGCRISHCWKEGNEPVTQWKAIVLGQLPTNPSLYLVKYDGIDSIYGQELYSDDRILNLKVL PPIVVFPQVRDAHLARALVGRAVQQKFERKDGSEVNWRGVVLAQVPIMKDLFYITYKKDPALYAYQLLDDYKEGNLHMIPD TPPAEERSGGDSDVLIGNWVQYTRKDGSKKFGKVVYQVLDNPSVFFIKFHGDIHIYVYTMVPKILEVEKS UniProt: >sp|Q16695|H31T_HUMAN Histone H3.1t OS=Homo sapiens GN=HIST3H3 PE=1 SV=3 MARTKQTARKSTGGKAPRKQLATKVARKSAPATGGVKKPHRYRPGTVALREIRRYQKSTELLIRKLPFQRLMREIAQDFK TDLRFQSSAVMALQEACESYLVGLFEDTNLCVIHAKRVTIMPKDIQLARRIRGERA http://en.wikipedia.org/wiki/FASTA_format

PEFF - PSI Extended Fasta Format >sp:P06748 \ID=NPM_HUMAN \Pname=(Nucleophosmin) (NPM) (Nucleolar phosphoprotein B23) (Numatrin) (Nucleolar protein NO38) \NcbiTaxId=9606 \ModRes=(125|MOD:00046)(199|MOD:00047) \Length=294 >sp:P00761 \ID=TRYP_PIG \Pname=(Trypsin precursor) (EC 3.4.21.4) \NcbiTaxId=9823 \Variant=(20|20|V) \Processed=(1|8|PROPEP)(9|231|CHAIN) \Length=231 http://www.psidev.info/node/363

Sample-specific protein sequence databases Samples Peptides MS Protein DB Identified and quantified peptides and proteins

Sample-specific protein sequence databases Next-generation sequencing Samples of the genome and transcriptome Peptides MS Sample-specific Protein DB Identified and quantified peptides and proteins

Sample-specific protein sequence databases Next-generation sequencing Samples of the genome and transcriptome Peptides MS Sample-specific Protein DB Novel Expression Somatic and germ-line mutations TCGA G AGCTG Identified and quantified TCGA G AGCTG Exon 1 Exon 2 TCGA G AGCTG TCGA G AGCTG peptides and proteins TCGA G AGCTG Exon 1 TCGATAGCTG Alternative Splicing Gene Fusions Gene Y Gene X Gene X Gene Y Exon 2 Exon 1 Exon 2 Exon 1 Exon 1 Exon 2 Exon 3

Proteomics and Transcriptomics of Breast Tumors Xenograft Primary tumor Breast tumor RNA-Seq Illumina HiSeq Sample Loading Sample Collection Chamber Chamber 1 2 3 4 5 6 7 8 Stacking Gel Resolving Gel MS/MS ABI 5600 Triple TOF ---250,000 ---150,000 ---100,000 ---75,000 ---50,000 ---37,000 ---25,000 ---20,000 ---15,000 ---10,000

Germline and Somatic Variants 1 Basal germline Basal somatic Luminal germline 0.1 Luminal somatic Frequency 0.01 0.001 0.0001 0.00001 0 5 10 15 Number of amino acid changes The frequency of proteins as a function of the number of amino acid changes due to germline and somatic variants for the basal and luminal breast tumor xenografts

Alternative Splicing 100000 Model Exons Alternative splicing Novel expression 10000 Number of junctions 1000 100 10 1 1 10 100 1000 10000 Number of reads The number of exon/exon junctions as a function of the number of RNA-Seq reads for the basal breast tumor xenograft.

Protein identification using sample-specific sequence databases Protein DB Germline Tumor 362 variants genome sequence + germline Somatic 9 / somatic variants variants Potentially Tumor novel RNA-Seq 1114 peptides Spans splice Potentially 70 site novel peptides

Data Repositories

ProteomeExchange http://www.proteomeexchange.org/

PRIDE http://www.ebi.ac.uk/pride/

PeptideAtlas http://www.peptideatlas.org/

The Global Proteome Machine Databases (GPMDB) http://gpmdb.thegpm.org

Comparison with GPMDB Most proteins show very reproducible peptide patterns

Comparison with GPMDB Query Spectrum Best match In GPMDB Second best match In GPMDB

GPMDB usage last month

GPMDB usage last month

GPMDB Data Crowdsourcing Any lab performs experiments Raw data sent to public repository (TRANCHE, PRIDE) Data imported by GPMDB Data analyzed & accepted/rejected Accepted information loaded into public collection General community uses information and inspects data

Information for including a data set in GPMDB 1. MS/MS data (required) 1. MS raw data files 2. ASCII files: mzXML, mzML, MGF, DTA, etc. 3. Analysis files: DAT, MSF, BIOML 2. Sample Information (supply if possible) 1. Species : human, yeast 2. Cell/tissue type & subcellular localization 3. Reagents: urea, formic acid, etc. 4. Quantitation: SILAC, iTRAQ 5. Proteolysis agent: trypsin, Lys-C 3. Project information (suggested) 1. Project name 2. Contact information

How to characterize the evidence in GPMDB for a protein? High confidence Medium confidence Low confidence No observation

Statistical model for 212 observations of TP53 Star End N -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 Skew Kurt t 214 248 539 0.15 0.18 0.22 0.17 0.15 0.07 0.03 0.01 0.01 0.00 -0.01 -2.01 249 267 1010 0.04 0.09 0.13 0.16 0.16 0.14 0.13 0.06 0.04 0.05 -0.08 -1.89 182 196 832 0.09 0.15 0.20 0.19 0.18 0.13 0.05 0.01 0.00 0.00 -0.12 -1.84 250 267 4 0.25 0.00 0.25 0.00 0.25 0.00 0.00 0.00 0.00 0.25 0.48 -2.28 1 24 269 0.10 0.12 0.12 0.17 0.12 0.12 0.14 0.04 0.04 0.03 -0.33 -0.88 24 65 51 0.22 0.22 0.20 0.14 0.06 0.00 0.04 0.08 0.02 0.04 0.47 -1.62 66 101 334 0.09 0.08 0.11 0.11 0.09 0.11 0.09 0.13 0.08 0.12 0.10 -1.21 249 273 60 0.02 0.00 0.20 0.10 0.13 0.25 0.20 0.07 0.03 0.00 0.45 -1.36 214 242 10 0.00 0.10 0.00 0.00 0.00 0.00 0.30 0.20 0.20 0.20 0.54 -1.39 214 239 32 0.03 0.06 0.16 0.16 0.09 0.22 0.09 0.16 0.00 0.03 0.20 -0.99 111 120 117 0.09 0.20 0.15 0.26 0.29 0.01 0.00 0.00 0.00 0.00 0.62 -1.36 251 267 16 0.00 0.00 0.13 0.25 0.19 0.13 0.13 0.13 0.06 0.00 0.24 -0.60 214 241 14 0.00 0.00 0.00 0.07 0.29 0.21 0.07 0.29 0.00 0.07 0.87 -0.97 159 174 100 0.30 0.25 0.31 0.03 0.07 0.03 0.01 0.00 0.00 0.00 0.99 -1.07 68 101 10 0.00 0.00 0.00 0.00 0.00 0.20 0.10 0.10 0.30 0.30 0.86 -0.91 235 248 30 0.00 0.03 0.00 0.00 0.30 0.20 0.23 0.13 0.03 0.07 0.81 -0.82

Statistical model for observations of DNAH2

Statistical model for observations of GRAP2

DNA Repair

DNA Repair

TP53BP1:p, tumor protein p53 binding protein 1

TP53BP1:p, tumor protein p53 binding protein 1

Sequence Annotations

TP53BP1:p, tumor protein p53 binding protein 1

TP53BP1:p, tumor protein p53 binding protein 1

Peptide observations, catalase Peptide Sequence Observations FSTVAGESGSADTVR 2633 FNTANDDNVTQVR 2432 AFYVNVLNEEQR 1722 LVNANGEAVYCK 1701 GPLLVQDVVFTDEMAHFDR 1637 LSQEDPDYGIR 1560 LFAYPDTHR 1499 NLSVEDAAR 1400 FYTEDGNWDLVGNNTPIFFIR 1386 ADVLTTGAGNPVGDK 1338

Peptide frequency ( ω ), catalase ω Peptide Sequence FSTVAGESGSADTVR 0.08 FNTANDDNVTQVR 0.07 AFYVNVLNEEQR 0.05 LVNANGEAVYCK 0.05 GPLLVQDVVFTDEMAHFDR 0.05 LSQEDPDYGIR 0.04 LFAYPDTHR 0.04 NLSVEDAAR 0.04 FYTEDGNWDLVGNNTPIFFIR 0.04 ADVLTTGAGNPVGDK 0.04

Global frequency of observation ( ω ), catalase 0.08 0.06 ω 0.04 0.02 0.00 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Peptide sequences

Omega ( Ω ) value for a protein identification For any set peptides observed in an experiment assigned to a particular protein ( 1 to j ): ∑ Ω = ω ( ) protein j j Ω protein ≤ ( ) 1

Protein Ω’s for a set of identifications Ω ( z= 2) Ω ( z= 3) Protein ID SERPINB1 0.88 0.82 SNRPD1 0.88 0.59 CFL1 0.81 0.87 SNRPE 0.8 0.81 PPIA 0.79 0.64 CSTA 0.79 0.36 PFN1 0.76 0.61 CAT 0.71 0.78 GLRX 0.66 0.8 CALM1 0.62 0.76 FABP5 0.57 0.17

Retention Time Distribution