IMP GPI Biosynthesis Report IMP-Bioinformatics
GPI Biosynthesis Main Page  
GPI Site Motif   GPI Site Prediction   Home Page B.E.  

GPI Anchor Biosynthesis Report: TRYP10.0.001820_9 (PIG-P family, Trypanosoma brucei)




BLASTP 2.1.1 [Aug-8-2000]

Reference:
Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schäffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs",  Nucleic Acids Res. 25:3389-3402.

Query= 
         (148 letters)

Database: nr
           887,402 sequences; 277,845,442 total letters

Searching..................................................

Converged !!!


Results of PSI-Blast iteration 2

Distribution of 15 Blast Hits on the Query Sequence


Sequences with E-value BETTER than threshold (0.002)
Sequences producing significant alignments:
sp|P57054|DSR5_HUMAN Down syndrome critical region protein 5 (Do... 181 8e-45
pir||JC7302 Down syndrome critical region protein 5 beta - human 180 1e-44
dbj|BAA95636.1| (AB035745) DSCR5b [Homo sapiens] >gi|9963906|gb|... 179 2e-44
ref|NP_062416.1| (NM_019543) Down syndrome critical region prote... 177 8e-44
pir||JC7301 Down syndrome critical region protein 5 alpha - huma... 177 8e-44
ref|NP_176323.1| (NM_104809) hypothetical protein [Arabidopsis t... 169 2e-41
ref|NP_057514.1| (NM_016430) Down syndrome critical region prote... 153 1e-36
dbj|BAA95512.1| (AP001727) Down syndrome critical region protein... 153 1e-36
ref|NP_496845.1| (NM_064444) BR-1 protein-like [Caenorhabditis e... 125 4e-28
pir||T27023 hypothetical protein Y48E1B.2a - Caenorhabditis elegans 123 1e-27
ref|NP_181476.1| (NM_129502) hypothetical protein [Arabidopsis t... 118 5e-26
emb|CAB16583.2| (Z99295) human downs syndrome critical region-li... 117 8e-26
pir||T38153 hypothetical protein SPAC22A12.13 - fission yeast (... 117 8e-26
Sequences with E-value WORSE than threshold

ref|NP_010725.1| (NC_001136) Protein required for cell viability... 39 0.041
gb|AAF56563.1| (AE003754) CG14550 gene product [Drosophila melan... 36 0.35
Alignments
>sp|P57054|DSR5_HUMAN Down syndrome critical region protein 5 (Down syndrome critical
           region protein C)
 dbj|BAA96872.1| (AB037163) DSCR5b [Homo sapiens]
          Length = 158

 Score =  181 bits (455), Expect = 8e-45
 Identities = 37/141 (26%), Positives = 68/141 (47%), Gaps = 9/141 (6%)

Query: 4   RCTKPLGATPATPTAPPTGRKGHQVAINGFITWILVTMSLVVYFLWAFIPTSFLDMVLAS 63
           R +K  G       +P       + AI GF+ ++      ++Y +WAFIP S+L+ +  +
Sbjct: 17  RLSKAPGKMVENSPSPLP-----ERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLT 71

Query: 64  YYPDKYWAVAIPAILVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVD 123
           Y+P KYWAVA+P  L++ +V    + F + +  T PL     +       + +       
Sbjct: 72  YWPQKYWAVALPVYLLIAIVIGYVLLFGINMMSTSPLDSIHTITDNYAKNQQQKKYQ--- 128

Query: 124 VEDGVPPITDIPVSVASRLLF 144
            E+ +P + DI +S  +++ F
Sbjct: 129 -EEAIPALRDISISEVNQMFF 148
>pir||JC7302 Down syndrome critical region protein 5 beta - human
          Length = 158

 Score =  180 bits (453), Expect = 1e-44
 Identities = 38/141 (26%), Positives = 68/141 (47%), Gaps = 9/141 (6%)

Query: 4   RCTKPLGATPATPTAPPTGRKGHQVAINGFITWILVTMSLVVYFLWAFIPTSFLDMVLAS 63
           R +K  G       +P       + AI GF+ ++      ++Y +WAFIP S+L+ +  +
Sbjct: 17  RLSKAPGKMVENSPSPLP-----ERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLT 71

Query: 64  YYPDKYWAVAIPAILVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVD 123
           Y+P KYWAVA+P  L++ +V    + F L +  T PL     +       + +       
Sbjct: 72  YWPQKYWAVALPVYLLIAIVIGYVLLFGLNMMSTSPLDSIHTITDNYAKNQQQKKYQ--- 128

Query: 124 VEDGVPPITDIPVSVASRLLF 144
            E+ +P + DI +S  +++ F
Sbjct: 129 -EEAIPALRDISISEVNQMFF 148
>dbj|BAA95636.1| (AB035745) DSCR5b [Homo sapiens]
 gb|AAG09757.1|AF237812_1 (AF237812) NPD010 [Homo sapiens]
          Length = 158

 Score =  179 bits (452), Expect = 2e-44
 Identities = 37/141 (26%), Positives = 68/141 (47%), Gaps = 9/141 (6%)

Query: 4   RCTKPLGATPATPTAPPTGRKGHQVAINGFITWILVTMSLVVYFLWAFIPTSFLDMVLAS 63
           R +K  G       +P       + AI GF+ ++      ++Y +WAFIP S+L+ +  +
Sbjct: 17  RLSKAPGKMVENSPSPLP-----ERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLT 71

Query: 64  YYPDKYWAVAIPAILVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVD 123
           Y+P KYWAVA+P  L++ +V    + F + +  T PL     +       + +       
Sbjct: 72  YWPQKYWAVALPVYLLIAIVIGYVLLFGINMMSTSPLDSIHTITDNYAKNQQQKKYQ--- 128

Query: 124 VEDGVPPITDIPVSVASRLLF 144
            E+ +P + DI +S  +++ F
Sbjct: 129 -EEAIPALRDISISEVNQMFF 148
>ref|NP_062416.1| (NM_019543) Down syndrome critical region protein c [Mus musculus]
 gb|AAF32294.1|AF216306_1 (AF216306) DCRC-1 [Mus musculus]
 gb|AAF32295.1|AF216307_1 (AF216307) DCRC-2 [Mus musculus]
 gb|AAF32296.1|AF216308_1 (AF216308) DCRC-3 [Mus musculus]
 gb|AAF32297.1|AF216309_1 (AF216309) DCRC-4 [Mus musculus]
 dbj|BAB23204.1| (AK004169) Down syndrome critical region gene c~data source:MGD,
           source key:MGI:1860433, evidence:ISS~putative [Mus
           musculus]
          Length = 132

 Score =  177 bits (446), Expect = 8e-44
 Identities = 35/127 (27%), Positives = 65/127 (50%), Gaps = 9/127 (7%)

Query: 18  APPTGRKGHQVAINGFITWILVTMSLVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAI 77
           +P       + AI GF+ ++      ++Y +WAF+P S+L+ +  +Y+P KYWAVA+P  
Sbjct: 7   SPLP-----ERAIYGFVLFLSSQFGFILYLVWAFVPESWLNSLGLTYWPQKYWAVALPVY 61

Query: 78  LVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVDVEDGVPPITDIPVS 137
           L++T+V    + F + +  T PL     +       +          ED +P + D+P+S
Sbjct: 62  LLITVVIGYVLLFGINMMSTSPLDSIHTITDNYAKNQQRKNYQ----EDAIPALRDVPIS 117

Query: 138 VASRLLF 144
             +++ F
Sbjct: 118 EVNKMFF 124
>pir||JC7301 Down syndrome critical region protein 5 alpha - human
 dbj|BAA95633.1| (AB035742) DSCR5a1 [Homo sapiens]
 dbj|BAA95634.1| (AB035743) DSCR5a2 [Homo sapiens]
 dbj|BAA95635.1| (AB035744) DSCR5a3 [Homo sapiens]
 dbj|BAA96871.1| (AB037162) DSCR5a [Homo sapiens]
 dbj|BAB12395.1| (AB039659) PIG-P [Homo sapiens]
 gb|AAH05180.1|AAH05180 (BC005180) Unknown (protein for MGC:12168) [Homo sapiens]
 gb|AAH11007.1|AAH11007 (BC011007) Unknown (protein for MGC:13327) [Homo sapiens]
          Length = 134

 Score =  177 bits (446), Expect = 8e-44
 Identities = 34/127 (26%), Positives = 64/127 (49%), Gaps = 9/127 (7%)

Query: 18  APPTGRKGHQVAINGFITWILVTMSLVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAI 77
           +P       + AI GF+ ++      ++Y +WAFIP S+L+ +  +Y+P KYWAVA+P  
Sbjct: 7   SPLP-----ERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLTYWPQKYWAVALPVY 61

Query: 78  LVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVDVEDGVPPITDIPVS 137
           L++ +V    + F + +  T PL     +       + +        E+ +P + DI +S
Sbjct: 62  LLIAIVIGYVLLFGINMMSTSPLDSIHTITDNYAKNQQQKKYQ----EEAIPALRDISIS 117

Query: 138 VASRLLF 144
             +++ F
Sbjct: 118 EVNQMFF 124
>ref|NP_176323.1| (NM_104809) hypothetical protein [Arabidopsis thaliana]
 sp|O64792|DSR5_ARATH Down syndrome critical region protein 5 homolog
 gb|AAC13913.1|AAC13913 (AC004255) T1F9.23 [Arabidopsis thaliana]
          Length = 137

 Score =  169 bits (426), Expect = 2e-41
 Identities = 30/116 (25%), Positives = 59/116 (50%), Gaps = 1/116 (0%)

Query: 30  INGFITWILVTMSLVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAILVMTMVYYLTVH 89
           + GF+  I + ++ V++ +W ++P  FL+ +   YYP KYWA+A+P   ++T++  L  +
Sbjct: 16  VYGFVGSISIVVATVIFLIWGYVPDKFLESIGIYYYPSKYWAMAMPMYSMVTLLVALVFY 75

Query: 90  FLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVDVED-GVPPITDIPVSVASRLLF 144
             L    T   T    +           +  + + ED  + PI+DI ++  + L+F
Sbjct: 76  IGLNFMSTSKPTSLNTLFDDYSREDVNFLPLMKNGEDRPIDPISDIDITRINDLMF 131
>ref|NP_057514.1| (NM_016430) Down syndrome critical region protein 5; Down syndrome
           critical region protein C short form; Down syndrome
           critical region protein C [Homo sapiens]
 gb|AAF32289.1|AF216305_1 (AF216305) Down syndrome critical region protein C short form [Homo
           sapiens]
 dbj|BAA96873.1| (AB037164) DSCR5c [Homo sapiens]
          Length = 108

 Score =  153 bits (385), Expect = 1e-36
 Identities = 30/102 (29%), Positives = 55/102 (53%), Gaps = 4/102 (3%)

Query: 43  LVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAILVMTMVYYLTVHFLLMLYRTDPLTD 102
           +V+Y +WAFIP S+L+ +  +Y+P KYWAVA+P  L++ +V    + F + +  T PL  
Sbjct: 1   MVLYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAVVIGYVLLFGINMMSTSPLDS 60

Query: 103 GFCVAQTNGAVRHESIENLVDVEDGVPPITDIPVSVASRLLF 144
              +       + +        E+ +P + DI +S  +++ F
Sbjct: 61  IHTITDNYAKNQQQKKYQ----EEAIPALRDISISEVNQMFF 98
>dbj|BAA95512.1| (AP001727) Down syndrome critical region protein C [Homo sapiens]
          Length = 108

 Score =  153 bits (384), Expect = 1e-36
 Identities = 30/102 (29%), Positives = 55/102 (53%), Gaps = 4/102 (3%)

Query: 43  LVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAILVMTMVYYLTVHFLLMLYRTDPLTD 102
           +V+Y +WAFIP S+L+ +  +Y+P KYWAVA+P  L++ +V    + F + +  T PL  
Sbjct: 1   MVLYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVIGYVLLFGINMMSTSPLDS 60

Query: 103 GFCVAQTNGAVRHESIENLVDVEDGVPPITDIPVSVASRLLF 144
              +       + +        E+ +P + DI +S  +++ F
Sbjct: 61  IHTITDNYAKNQQQKKYQ----EEAIPALRDISISEVNQMFF 98
>ref|NP_496845.1| (NM_064444) BR-1 protein-like [Caenorhabditis elegans]
 emb|CAB07689.2| (Z93393) contains similarity to Pfam domain: PF01926 (GTPase of
          unknown function), E-value=0.00012,~cDNA EST
          EMBL:AU111138 comes from this gene~cDNA EST
          EMBL:AU114976 comes from this gene [Caenorhabditis
          elegans]
          Length = 890

 Score =  125 bits (313), Expect = 4e-28
 Identities = 25/88 (28%), Positives = 35/88 (39%), Gaps = 6/88 (6%)

Query: 6  TKPLGATPATPTAPPTGRKGH------QVAINGFITWILVTMSLVVYFLWAFIPTSFLDM 59
          T+    T  TP  P    + H         I GF  +I+     V+Y +WA  P   L  
Sbjct: 4  TESSSTTKKTPKTPLPSEEIHLPGPHPARGIYGFALYIVSWTLFVIYLIWAITPVPILYR 63

Query: 60 VLASYYPDKYWAVAIPAILVMTMVYYLT 87
          +  +Y P K WA+AI          Y+T
Sbjct: 64 LGITYIPSKLWALAIGIFFPTAACLYVT 91
>pir||T27023 hypothetical protein Y48E1B.2a - Caenorhabditis elegans
          Length = 1768

 Score =  123 bits (308), Expect = 1e-27
 Identities = 25/88 (28%), Positives = 35/88 (39%), Gaps = 6/88 (6%)

Query: 6  TKPLGATPATPTAPPTGRKGH------QVAINGFITWILVTMSLVVYFLWAFIPTSFLDM 59
          T+    T  TP  P    + H         I GF  +I+     V+Y +WA  P   L  
Sbjct: 4  TESSSTTKKTPKTPLPSEEIHLPGPHPARGIYGFALYIVSWTLFVIYLIWAITPVPILYR 63

Query: 60 VLASYYPDKYWAVAIPAILVMTMVYYLT 87
          +  +Y P K WA+AI          Y+T
Sbjct: 64 LGITYIPSKLWALAIGIFFPTAACLYVT 91
>ref|NP_181476.1| (NM_129502) hypothetical protein [Arabidopsis thaliana]
 pir||T00554 hypothetical protein F12L6.10 - Arabidopsis thaliana
 gb|AAC27835.1| (AC004218) hypothetical protein [Arabidopsis thaliana]
          Length = 773

 Score =  118 bits (294), Expect = 5e-26
 Identities = 17/72 (23%), Positives = 32/72 (43%), Gaps = 1/72 (1%)

Query: 74  IPAILVMTMVYYLTVHFLLMLYRTDPLTDGFCVAQTNGAVRHESIENLVDVED-GVPPIT 132
           +P  L++T++  L  +  L    T P T    +         E    + + ED  + PI+
Sbjct: 1   MPTYLMVTLMLALAFYIGLNFIATPPPTSLNTLFDEYSREPGEFAAGMEEGEDRPIDPIS 60

Query: 133 DIPVSVASRLLF 144
           DI ++  + L+F
Sbjct: 61  DIDITRINDLMF 72
>emb|CAB16583.2| (Z99295) human downs syndrome critical region-like
           [Schizosaccharomyces pombe]
          Length = 120

 Score =  117 bits (293), Expect = 8e-26
 Identities = 27/90 (30%), Positives = 43/90 (47%), Gaps = 7/90 (7%)

Query: 31  NGFITWILVTMSLVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAILVMTMVYYLTVHF 90
            GF+ +++  +   VY +WA  P   L      YY  ++WA+AIP  L + ++Y   +H 
Sbjct: 10  YGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY---IHV 66

Query: 91  LLMLYRTDPLTDGF----CVAQTNGAVRHE 116
           +L  Y T+ LT  F    C+      V  E
Sbjct: 67  VLNAYNTEVLTKPFSSLECIVDQYALVGEE 96
>pir||T38153 hypothetical protein SPAC22A12.13 - fission yeast
           (Schizosaccharomyces pombe)
          Length = 157

 Score =  117 bits (292), Expect = 8e-26
 Identities = 27/90 (30%), Positives = 43/90 (47%), Gaps = 7/90 (7%)

Query: 31  NGFITWILVTMSLVVYFLWAFIPTSFLDMVLASYYPDKYWAVAIPAILVMTMVYYLTVHF 90
            GF+ +++  +   VY +WA  P   L      YY  ++WA+AIP  L + ++Y   +H 
Sbjct: 47  YGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY---IHV 103

Query: 91  LLMLYRTDPLTDGF----CVAQTNGAVRHE 116
           +L  Y T+ LT  F    C+      V  E
Sbjct: 104 VLNAYNTEVLTKPFSSLECIVDQYALVGEE 133
>ref|NP_010725.1| (NC_001136) Protein required for cell viability; Ydr437wp
           [Saccharomyces cerevisiae]
 pir||S69717 hypothetical protein YDR437w - yeast (Saccharomyces cerevisiae)
 gb|AAB64854.1| (U33007) Ydr437wp; CAI: 0.10 [Saccharomyces cerevisiae]
          Length = 140

 Score = 39.3 bits (91), Expect = 0.041
 Identities = 20/90 (22%), Positives = 43/90 (47%), Gaps = 16/90 (17%)

Query: 31  NGFITWILVTMSLVVYFLWAFIPTS--------FLDMVLASYYPDKYWAVAIPAILVMTM 82
             F  ++++T +LV+  +W+ +P+S        F++ +    +P + W + + +I++M M
Sbjct: 7   YWFSQYMIITSTLVLTIIWSILPSSLGEAAPKQFINTL-LDIFPQRRWIITLESIMLMGM 65

Query: 83  VYYLTVHFLLMLYR----TDPLTDGFCVAQ 108
              L  +  L++Y     T PL     V  
Sbjct: 66  ---LCTYIGLLMYNEDTLTPPLDSLSTVTD 92
>gb|AAF56563.1| (AE003754) CG14550 gene product [Drosophila melanogaster]
          Length = 138

 Score = 36.1 bits (83), Expect = 0.35
 Identities = 18/71 (25%), Positives = 32/71 (44%), Gaps = 6/71 (8%)

Query: 40  TMSLVVYFLWAFIPTSFLDMVLAS-YYPDKYWAVAIPAILVM-TMVYYLTVHFLLMLYRT 97
           T+  +VY  WA +P  F        Y PDKY+AV +P ++++    +   ++  + L  T
Sbjct: 4   TVLFLVYVTWALLPVEF----GLHSYLPDKYFAVFVPFLVLVFAWFFAFLIYPAINLSMT 59

Query: 98  DPLTDGFCVAQ 108
             +     V  
Sbjct: 60  VDVDSIASVVD 70
CPU time:    42.65 user secs.	    1.74 sys. secs	   44.39 total secs.

  Database: nr
    Posted date:  Apr 21, 2002  2:19 PM
  Number of letters in database: 277,845,442
  Number of sequences in database:  887,402
  
Lambda     K      H
   0.314    0.169    0.530 

Gapped
Lambda     K      H
   0.270   0.0592    0.230 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 99959590
Number of Sequences: 887402
Number of extensions: 4470056
Number of successful extensions: 24975
Number of sequences better than 10.0: 50
Number of HSP's better than 10.0 without gapping: 25
Number of HSP's successfully gapped in prelim test: 25
Number of HSP's that attempted gapping in prelim test: 24929
Number of HSP's gapped (non-prelim): 52
length of query: 148
length of database: 277,845,442
effective HSP length: 48
effective length of query: 100
effective length of database: 235,250,146
effective search space: 23525014600
effective search space used: 23525014600
T: 11
A: 40
X1: 15 ( 6.8 bits)
X2: 38 (14.8 bits)
X3: 64 (24.9 bits)
S1: 40 (20.7 bits)
S2: 70 (31.3 bits)