diff options
Diffstat (limited to 'Master/texmf-dist/doc/latex/texshade/AQP1.phd')
| -rw-r--r-- | Master/texmf-dist/doc/latex/texshade/AQP1.phd | 1907 |
1 files changed, 1907 insertions, 0 deletions
diff --git a/Master/texmf-dist/doc/latex/texshade/AQP1.phd b/Master/texmf-dist/doc/latex/texshade/AQP1.phd new file mode 100644 index 00000000000..4f6fa61c328 --- /dev/null +++ b/Master/texmf-dist/doc/latex/texshade/AQP1.phd @@ -0,0 +1,1907 @@ + + +From phd@EMBL-Heidelberg.de Wed Nov 25 10:24:25 1998 +Date: Tue, 24 Nov 1998 17:45:25 +0100 +From: Protein Prediction <phd@EMBL-Heidelberg.de> +To: eric.beitz@uni-tuebingen.de +Subject: PredictProtein + +The following information has been received by the server: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +________________________________________________________________________________ + +reference predict_h25873 (Tue Nov 24 17:43:21 MET 1998) +from eric.beitz@uni-tuebingen.de +password(###) +resp MAIL +orig HTML +prediction of: -secondary structure (PHDsec)-solvent accessibility (PHDacc)- +return msf format +# no description +MASEIKKKLFWRAVVAEFLAMTLFVFISIGSALGFNYPLERNQTLVQDNVKVSLAFGLSIATLAQSVGHISGAHSNPAVT +LGLLLSCQISILRAVMYIIAQCVGAIVASAILSGITSSLLENSLGRNDLARGVNSGQGLGIEIIGTLQLVLCVLATTDRR +RRDLGGSAPLAIGLSVALGHLLAIDYTGCGINPARSFGSAVLTRNFSNHWIFWVGPFIGSALAVLIYDFILAPRSSDFTD +RMKVWTSGQVEEYDLDADDINSRVEMKPK + +________________________________________________________________________________ + +Result of PROSITE search (Amos Bairoch): +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +please quote: A Bairoch, P Bucher & K Hofmann: The PROSITE database, +its status in 1997. Nucl. Acids Res., 1997, 25, 217-221. + +________________________________________________________________________________ + +-------------------------------------------------------- + +-------------------------------------------------------- + +Pattern-ID: ASN_GLYCOSYLATION PS00001 PDOC00001 +Pattern-DE: N-glycosylation site +Pattern: N[^P][ST][^P] + 42 NQTL + 250 NFSN + +Pattern-ID: GLYCOSAMINOGLYCAN PS00002 PDOC00002 +Pattern-DE: Glycosaminoglycan attachment site +Pattern: SG.G + 135 SGQG + +Pattern-ID: PKC_PHOSPHO_SITE PS00005 PDOC00005 +Pattern-DE: Protein kinase C phosphorylation site +Pattern: [ST].[RK] + 157 TDR + 398 TDR + +Pattern-ID: CK2_PHOSPHO_SITE PS00006 PDOC00006 +Pattern-DE: Casein kinase II phosphorylation site +Pattern: [ST].{2}[DE] + 118 SLLE + 383 SRVE + +Pattern-ID: MYRISTYL PS00008 PDOC00008 +Pattern-DE: N-myristoylation site +Pattern: G[^EDRKHPFYW].{2}[STAGCN][^P] + 30 GSALGF + 92 GLSIAT + 179 GLLLSC + 288 GAIVAS + 407 GITSSL + 544 GVNSGQ + 722 GLSVAL + 917 GINPAR + 1141 GSALAV + +Pattern-ID: PROKAR_LIPOPROTEIN PS00013 PDOC00013 +Pattern-DE: Prokaryotic membrane lipoprotein lipid attachment site +Pattern: [^DERK]{6}[LIVMFWSTAG]{2}[LIVMFYSTAGCQ][AGS]C + 77 PAVTLGLLLSC + +Pattern-ID: MIP PS00221 PDOC00193 +Pattern-DE: MIP family signature +Pattern: [HNQA].NP[STA][LIVMF][ST][LIVMF][GSTAFY] + 74 HSNPAVTLG + +________________________________________________________________________________ + +Result of ProDom domain search (Corpet, Gouzy, Kahn): +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +- please quote: ELL Sonnhammer & D Kahn, Prot. Sci., 1994, 3, 482-492 + +________________________________________________________________________________ + +--- ------------------------------------------------------------ +--- Results from running BLAST against PRODOM domains +--- +--- PLEASE quote: +--- F Corpet, J Gouzy, D Kahn (1998). The ProDom database +--- of protein domain families. Nucleic Ac Res 26:323-326. +--- +--- BEGIN of BLASTP output +BLASTP 1.4.7 [16-Oct-94] [Build 17:06:52 Oct 31 1994] + +Reference: Altschul, Stephen F., Warren Gish, Webb Miller, Eugene W. Myers, +and David J. Lipman (1990). Basic local alignment search tool. J. Mol. Biol. +215:403-10. + +Query= prot (#) ppOld, no description /home/phd/server/work/predict_h25873 + (269 letters) + +Database: /home/phd/ut/prodom/prodom_34_2 + 53,597 sequences; 6,740,067 total letters. +Searching..................................................done + + Smallest + Sum + High Probability +Sequences producing High-scoring Segment Pairs: Score P(N) N + + 390 p34.2 (45) MIP(6) AQP1(4) GLPF(4) // PROTEIN INTRIN... 270 2.0e-32 1 + 45663 p34.2 (1) AQPZ_ECOLI // AQUAPORIN Z. 90 3.2e-13 2 + 45611 p34.2 (1) AQP2_HUMAN // AQUAPORIN-CD (AQP-CD) (WAT... 136 6.0e-13 1 + 304 p34.2 (61) AQP2(10) GLPF(6) MIP(5) // PROTEIN CHANN... 121 9.2e-11 1 + 45607 p34.2 (1) PMIP_NICAL // POLLEN-SPECIFIC MEMBRANE I... 80 1.2e-07 2 + 45606 p34.2 (1) BIB_DROME // NEUROGENIC PROTEIN BIG BRAIN. 80 1.2e-05 2 + 2027 p34.2 (15) GLPF(9) AQP3(2) // PROTEIN FACILITATOR ... 60 3.4e-05 2 + 45615 p34.2 (1) GLPF_STRPN // GLYCEROL UPTAKE FACILITATO... 63 0.024 1 + 45638 p34.2 (1) AQP5_HUMAN // AQUAPORIN 5. 61 0.044 1 + +>390 p34.2 (45) MIP(6) AQP1(4) GLPF(4) // PROTEIN INTRINSIC CHANNEL WATER + AQUAPORIN TONOPLAST MEMBRANE FOR PLASMA LENS + Length = 88 + + Score = 270 (125.3 bits), Expect = 2.0e-32, P = 2.0e-32 + Identities = 47/67 (70%), Positives = 56/67 (83%) + +Query: 156 TTDRRRRDLGGSAPLAIGLSVALGHLLAIDYTGCGINPARSFGSAVLTRNFSNHWIFWVG 215 + T D+RR +GGSAPL IG SVALGHL+ I YTGCG+NPARSFG AV+T NF+NHW++WVG +Sbjct: 22 TDDKRRGSVGGSAPLPIGFSVALGHLIGIPYTGCGMNPARSFGPAVVTGNFTNHWVYWVG 81 + +Query: 216 PFIGSAL 222 + P IG+ L +Sbjct: 82 PIIGAVL 88 + + Score = 95 (44.1 bits), Expect = 2.3e-06, P = 2.3e-06 + Identities = 20/33 (60%), Positives = 23/33 (69%) + +Query: 136 GQGLGIEIIGTLQLVLCVLATTDRRRRDLGGSA 168 + GQ L +EIIGT QLV CV ATTD +RR G + +Sbjct: 1 GQNLVVEIIGTFQLVYCVFATTDDKRRGSVGGS 33 + +>45663 p34.2 (1) AQPZ_ECOLI // AQUAPORIN Z. + Length = 96 + + Score = 90 (41.8 bits), Expect = 3.2e-13, Sum P(2) = 3.2e-13 + Identities = 18/36 (50%), Positives = 25/36 (69%) + +Query: 166 GSAPLAIGLSVALGHLLAIDYTGCGINPARSFGSAV 201 + G AP+AIGL++ L HL++I T +NPARS A+ +Sbjct: 25 GFAPIAIGLALTLIHLISIPVTNTSVNPARSTAVAI 60 + + Score = 63 (29.2 bits), Expect = 3.2e-13, Sum P(2) = 3.2e-13 + Identities = 11/25 (44%), Positives = 14/25 (56%) + +Query: 210 WIFWVGPFIGSALAVLIYDFILAPR 234 + W FWV P +G + LIY +L R +Sbjct: 71 WFFWVVPIVGGIIGGLIYRTLLEKR 95 + +>45611 p34.2 (1) AQP2_HUMAN // AQUAPORIN-CD (AQP-CD) (WATER CHANNEL PROTEIN FOR + RENAL COLLECTING DUCT) (ADH WATER CHANNEL) (AQUAPORIN 2) (COLLECTING DUCT + WATER CHANNEL PROTEIN) (WCH-CD). + Length = 49 + + Score = 136 (63.1 bits), Expect = 6.0e-13, P = 6.0e-13 + Identities = 23/42 (54%), Positives = 34/42 (80%) + +Query: 50 VKVSLAFGLSIATLAQSVGHISGAHSNPAVTLGLLLSCQISI 91 + +++++AFGL I TL Q++GHISGAH NPAVT+ L+ C +S+ +Sbjct: 8 LQIAMAFGLGIGTLVQALGHISGAHINPAVTVACLVGCHVSV 49 + +>304 p34.2 (61) AQP2(10) GLPF(6) MIP(5) // PROTEIN CHANNEL WATER AQUAPORIN + INTRINSIC DUCT COLLECTING FOR TONOPLAST WCH-CD + Length = 43 + + Score = 121 (56.1 bits), Expect = 9.2e-11, P = 9.2e-11 + Identities = 24/43 (55%), Positives = 31/43 (72%) + +Query: 70 ISGAHSNPAVTLGLLLSCQISILRAVMYIIAQCVGAIVASAIL 112 + ISG H NPAVT+GLL+ + LRAV YI AQ +GA+ +A+L +Sbjct: 1 ISGGHINPAVTIGLLIGGRFPFLRAVFYIAAQLLGAVAGAALL 43 + +>45607 p34.2 (1) PMIP_NICAL // POLLEN-SPECIFIC MEMBRANE INTEGRAL PROTEIN. + Length = 69 + + Score = 80 (37.1 bits), Expect = 1.2e-07, Sum P(2) = 1.2e-07 + Identities = 17/54 (31%), Positives = 32/54 (59%) + +Query: 149 LVLCVLATTDRRRRDLGGSAPLAIGLSVALGHLLAIDYTGCGINPARSFGSAVL 202 + L++ V++ R +G A +A+G+++ L +A +G +NPARS G A++ +Sbjct: 13 LLMFVISGVATDDRAIGQVAGIAVGMTITLNVFVAGPISGASMNPARSIGPAIV 66 + + Score = 34 (15.8 bits), Expect = 1.2e-07, Sum P(2) = 1.2e-07 + Identities = 8/18 (44%), Positives = 11/18 (61%) + +Query: 136 GQGLGIEIIGTLQLVLCV 153 + GQ L IEII + L+ + +Sbjct: 1 GQSLAIEIIISFLLMFVI 18 + +>45606 p34.2 (1) BIB_DROME // NEUROGENIC PROTEIN BIG BRAIN. + Length = 119 + + Score = 80 (37.1 bits), Expect = 1.2e-05, Sum P(2) = 1.2e-05 + Identities = 15/34 (44%), Positives = 24/34 (70%) + +Query: 1 MASEIKKKLFWRAVVAEFLAMTLFVFISIGSALG 34 + M +EI+ FWR++++E LA ++VFI G+A G +Sbjct: 55 MQAEIRTLEFWRSIISECLASFMYVFIVCGAAAG 88 + + Score = 39 (18.1 bits), Expect = 1.2e-05, Sum P(2) = 1.2e-05 + Identities = 9/17 (52%), Positives = 12/17 (70%) + +Query: 53 SLAFGLSIATLAQSVGH 69 + +LA GL++ATL Q H +Sbjct: 103 ALASGLAMATLTQCFLH 119 + +>2027 p34.2 (15) GLPF(9) AQP3(2) // PROTEIN FACILITATOR GLYCEROL UPTAKE + AQUAPORIN DIFFUSION UPTAKE/EFFLUX PEPX 5'REGION ORF1 + Length = 55 + + Score = 60 (27.8 bits), Expect = 3.4e-05, Sum P(2) = 3.4e-05 + Identities = 17/46 (36%), Positives = 20/46 (43%) + +Query: 156 TTDRRRRDLGGSAPLAIGLSVALGHLLAIDYTGCGINPARSFGSAV 201 + T D GG PL +G V + TG INPAR FG + +Sbjct: 10 TDDGNNVPSGGLHPLMVGFLVMGIGMSLGGTTGYAINPARDFGPRI 55 + + Score = 37 (17.2 bits), Expect = 3.4e-05, Sum P(2) = 3.4e-05 + Identities = 7/10 (70%), Positives = 8/10 (80%) + +Query: 149 LVLCVLATTD 158 + L+ CVLA TD +Sbjct: 2 LIACVLALTD 11 + +>45615 p34.2 (1) GLPF_STRPN // GLYCEROL UPTAKE FACILITATOR PROTEIN. + Length = 26 + + Score = 63 (29.2 bits), Expect = 0.025, P = 0.024 + Identities = 13/23 (56%), Positives = 18/23 (78%) + +Query: 205 NFSNHWIFWVGPFIGSALAVLIY 227 + ++S WI VGP IG+ALAVL++ +Sbjct: 1 DWSYAWIPVVGPVIGAALAVLVF 23 + +>45638 p34.2 (1) AQP5_HUMAN // AQUAPORIN 5. + Length = 27 + + Score = 61 (28.3 bits), Expect = 0.045, P = 0.044 + Identities = 11/19 (57%), Positives = 18/19 (94%) + +Query: 50 VKVSLAFGLSIATLAQSVG 68 + ++++LAFGL+I TLAQ++G +Sbjct: 8 LQIALAFGLAIGTLAQALG 26 + +Parameters: + E=0.1 + B=500 + + V=500 + -ctxfactor=1.00 + + Query ----- As Used ----- ----- Computed ---- + Frame MatID Matrix name Lambda K H Lambda K H + +0 0 BLOSUM62 0.322 0.138 0.394 same same same + + Query + Frame MatID Length Eff.Length E S W T X E2 S2 + +0 0 269 269 0.10 69 3 11 22 0.22 33 + +Statistics: + Query Expected Observed HSPs HSPs + Frame MatID High Score High Score Reportable Reported + +0 0 59 (27.4 bits) 270 (125.3 bits) 14 14 + + Query Neighborhd Word Excluded Failed Successful Overlaps + Frame MatID Words Hits Hits Extensions Extensions Excluded + +0 0 5349 3124825 609708 2510548 4569 2 + + Database: /home/phd/ut/prodom/prodom_34_2 + Release date: unknown + Posted date: 12:24 PM MET DST May 06, 1998 + # of letters in database: 6,740,067 + # of sequences in database: 53,597 + # of database sequences satisfying E: 9 + No. of states in DFA: 564 (111 KB) + Total size of DFA: 226 KB (256 KB) + Time to generate neighborhood: 0.03u 0.00s 0.03t Real: 00:00:00 + Time to search database: 9.80u 0.03s 9.83t Real: 00:00:10 + Total cpu time: 9.90u 0.06s 9.96t Real: 00:00:10 +--- END of BLASTP output +--- ------------------------------------------------------------ +--- +--- Again: these results were obtained based on the domain data- +--- base collected by Daniel Kahn and his coworkers in Toulouse. +--- +--- PLEASE quote: +--- F Corpet, J Gouzy, D Kahn (1998). The ProDom database +--- of protein domain families. Nucleic Ac Res 26:323-326. +--- +--- The general WWW page is on: +---- --------------------------------------- +--- http://www.toulouse.inra.fr/prodom.html +---- --------------------------------------- +--- +--- For WWW graphic interfaces to PRODOM, in particular for your +--- protein family, follow the following links (each line is ONE +--- single link for your protein!!): +--- +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=390 ==> multiple alignment, consensus, PDB and PROSITE links of domain 390 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=390 ==> graphical output of all proteins having domain 390 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45663 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45663 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45663 ==> graphical output of all proteins having domain 45663 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45611 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45611 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45611 ==> graphical output of all proteins having domain 45611 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=304 ==> multiple alignment, consensus, PDB and PROSITE links of domain 304 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=304 ==> graphical output of all proteins having domain 304 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45607 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45607 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45607 ==> graphical output of all proteins having domain 45607 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45606 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45606 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45606 ==> graphical output of all proteins having domain 45606 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=2027 ==> multiple alignment, consensus, PDB and PROSITE links of domain 2027 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=2027 ==> graphical output of all proteins having domain 2027 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45615 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45615 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45615 ==> graphical output of all proteins having domain 45615 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom1=45638 ==> multiple alignment, consensus, PDB and PROSITE links of domain 45638 +http://www.toulouse.inra.fr/prodom/cgi-bin/ReqProdomII.pl?id_dom2=45638 ==> graphical output of all proteins having domain 45638 +--- +--- NOTE: if you want to use the link, make sure the entire line +--- is pasted as URL into your browser! +--- +--- END of PRODOM +--- ------------------------------------------------------------ + +________________________________________________________________________________ + +--- Database used for sequence comparison: +--- SEQBASE RELEASE 34.0 OF EMBL/SWISS-PROT WITH 59021 SEQUENCES + +The alignment that has been used as input to the network is: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +________________________________________________________________________________ + +--- ------------------------------------------------------------ +--- MAXHOM multiple sequence alignment +--- ------------------------------------------------------------ +--- +--- MAXHOM ALIGNMENT HEADER: ABBREVIATIONS FOR SUMMARY +--- ID : identifier of aligned (homologous) protein +--- STRID : PDB identifier (only for known structures) +--- PIDE : percentage of pairwise sequence identity +--- WSIM : percentage of weighted similarity +--- LALI : number of residues aligned +--- NGAP : number of insertions and deletions (indels) +--- LGAP : number of residues in all indels +--- LSEQ2 : length of aligned sequence +--- ACCNUM : SwissProt accession number +--- NAME : one-line description of aligned protein +--- +--- MAXHOM ALIGNMENT HEADER: SUMMARY +ID STRID IDE WSIM LALI NGAP LGAP LEN2 ACCNUM NAME +aqp1_rat 100 100 269 0 0 269 P29975 PROXIMAL TUBULE) (AQUAPOR +aqp1_mouse 98 99 269 0 0 269 Q02013 PROXIMAL TUBULE) (AQUAPOR +aqp1_human 93 97 269 0 0 269 P29972 PROXIMAL TUBULE) (AQUAPOR +aqp1_bovin 90 95 269 1 2 271 P47865 PROXIMAL TUBULE) (AQUAPOR +aqp1_sheep 90 94 269 2 3 272 P56401 PROXIMAL TUBULE) (AQUAPOR +aqpa_ranes 78 89 268 2 5 272 P50501 AQUAPORIN FA-CHIP. +aqp2_dasno 49 73 109 1 7 109 P79164 PROTEIN) (WCH-CD) (FRAGME +aqp2_bovin 49 73 109 1 7 109 P79099 PROTEIN) (WCH-CD) (FRAGME +aqp2_canfa 48 72 109 1 7 109 P79144 PROTEIN) (WCH-CD) (FRAGME +aqp2_rabit 48 73 109 1 7 109 P79213 PROTEIN) (WCH-CD) (FRAGME +aqp2_elema 47 72 109 1 7 109 P79168 PROTEIN) (WCH-CD) (FRAGME +aqp2_horse 47 72 109 1 7 109 P79165 PROTEIN) (WCH-CD) (FRAGME +aqp2_proha 47 73 109 1 7 109 P79229 PROTEIN) (WCH-CD) (FRAGME +mip_rat 46 73 259 1 7 261 P09011 LENS FIBER MAJOR INTRINSI +aqp2_oryaf 46 72 109 1 7 109 P79200 PROTEIN) (WCH-CD) (FRAGME +mip_mouse 46 73 261 1 7 263 P51180 LENS FIBER MAJOR INTRINSI +mip_ranpi 45 73 261 1 7 263 Q06019 LENS FIBER MAJOR INTRINSI +mip_bovin 45 73 261 1 7 263 P06624 LENS FIBER MAJOR INTRINSI +mip_human 45 73 261 1 7 263 P30301 LENS FIBER MAJOR INTRINSI +mip_chick 45 72 110 1 1 112 P28238 LENS FIBER MAJOR INTRINSI +aqp5_rat 44 71 262 2 8 265 P47864 AQUAPORIN 5. +aqp5_human 44 71 262 2 8 265 P55064 AQUAPORIN 5. +aqp2_human 44 72 261 2 8 271 P41181 PROTEIN) (WCH-CD). +aqp4_human 43 70 266 2 5 323 P55087 AQUAPORIN 4 (WCH4) (MERCU +aqp4_rat 43 70 266 2 5 323 P47863 AQUAPORIN 4 (WCH4) (MERCU +aqp4_mouse 43 69 265 3 6 322 P55088 AQUAPORIN 4 (WCH4) (MERCU +aqp2_rat 42 71 261 2 8 271 P34080 PROTEIN) (WCH-CD). +aqp2_mouse 42 71 261 2 8 271 P56402 PROTEIN) (WCH-CD). +wc2a_arath 42 67 248 4 12 287 P43286 PLASMA MEMBRANE INTRINSIC +aqp6_human 42 68 260 2 9 282 Q13520 AQUAPORIN 6 (AQUAPORIN-2 +wc2c_arath 41 66 248 4 12 285 P30302 INTRINSIC PROTEIN) (WSI-T +wc2b_arath 41 66 248 4 12 285 P43287 PLASMA MEMBRANE INTRINSIC +wc1c_arath 41 65 238 4 10 286 Q08733 (TMP-B). +wc1b_arath 41 65 238 4 10 286 Q06611 (TMP-A). +tipw_lyces 40 65 237 4 10 286 Q08451 (RIPENING-ASSOCIATED MEMB +wc1a_arath 40 64 238 4 10 286 P43285 PLASMA MEMBRANE INTRINSIC +tipw_pea 40 64 237 4 11 289 P25794 RESPONSIVE PROTEIN 7A). +tipa_arath 38 64 250 3 9 268 P26587 TONOPLAST INTRINSIC PROTE +aqua_atrca 38 64 246 4 10 282 P42767 AQUAPORIN. +dip_antma 38 65 242 2 4 250 P33560 PROBABLE TONOPLAST INTRIN +aqpz_ecoli 37 59 220 4 17 231 P48838 AQUAPORIN Z (BACTERIAL NO +tip2_tobac 37 64 242 2 4 250 P24422 TONOPLAST INTRINSIC PROTE +tip1_tobac 37 64 242 2 4 250 P21653 TONOPLAST INTRINSIC PROTE +tipg_arath 33 62 241 2 4 251 P25818 TONOPLAST INTRINSIC PROTE +bib_drome 33 60 260 4 10 700 P23645 NEUROGENIC PROTEIN BIG BR +tipr_arath 33 62 243 2 4 253 P21652 TONOPLAST INTRINSIC PROTE +tipa_phavu 33 62 246 2 4 256 P23958 TONOPLAST INTRINSIC PROTE +tipg_orysa 32 62 240 2 5 250 P50156 TONOPLAST INTRINSIC PROTE +--- +--- MAXHOM ALIGNMENT: IN MSF FORMAT +MSF of: /home/phd/server/work/predict_h25873-22040.hssp from: 1 to: 269 + /home/phd/server/work/predict_h25873-22040.msfRet MSF: 269 Type: P 24-Nov-98 17:44:5 Check: 3448 .. + + Name: predict_h258 Len: 269 Check: 8331 Weight: 1.00 + Name: aqp1_rat Len: 269 Check: 8331 Weight: 1.00 + Name: aqp1_mouse Len: 269 Check: 7552 Weight: 1.00 + Name: aqp1_human Len: 269 Check: 6501 Weight: 1.00 + Name: aqp1_bovin Len: 269 Check: 7067 Weight: 1.00 + Name: aqp1_sheep Len: 269 Check: 7582 Weight: 1.00 + Name: aqpa_ranes Len: 269 Check: 4844 Weight: 1.00 + Name: aqp2_dasno Len: 269 Check: 8933 Weight: 1.00 + Name: aqp2_bovin Len: 269 Check: 9649 Weight: 1.00 + Name: aqp2_canfa Len: 269 Check: 8990 Weight: 1.00 + Name: aqp2_rabit Len: 269 Check: 8787 Weight: 1.00 + Name: aqp2_elema Len: 269 Check: 9381 Weight: 1.00 + Name: aqp2_horse Len: 269 Check: 8993 Weight: 1.00 + Name: aqp2_proha Len: 269 Check: 8855 Weight: 1.00 + Name: mip_rat Len: 269 Check: 9773 Weight: 1.00 + Name: aqp2_oryaf Len: 269 Check: 8554 Weight: 1.00 + Name: mip_mouse Len: 269 Check: 9723 Weight: 1.00 + Name: mip_ranpi Len: 269 Check: 5937 Weight: 1.00 + Name: mip_bovin Len: 269 Check: 1430 Weight: 1.00 + Name: mip_human Len: 269 Check: 372 Weight: 1.00 + Name: mip_chick Len: 269 Check: 4658 Weight: 1.00 + Name: aqp5_rat Len: 269 Check: 9033 Weight: 1.00 + Name: aqp5_human Len: 269 Check: 6547 Weight: 1.00 + Name: aqp2_human Len: 269 Check: 6209 Weight: 1.00 + Name: aqp4_human Len: 269 Check: 2589 Weight: 1.00 + Name: aqp4_rat Len: 269 Check: 4412 Weight: 1.00 + Name: aqp4_mouse Len: 269 Check: 2845 Weight: 1.00 + Name: aqp2_rat Len: 269 Check: 5748 Weight: 1.00 + Name: aqp2_mouse Len: 269 Check: 6526 Weight: 1.00 + Name: wc2a_arath Len: 269 Check: 4866 Weight: 1.00 + Name: aqp6_human Len: 269 Check: 9404 Weight: 1.00 + Name: wc2c_arath Len: 269 Check: 6187 Weight: 1.00 + Name: wc2b_arath Len: 269 Check: 7328 Weight: 1.00 + Name: wc1c_arath Len: 269 Check: 8575 Weight: 1.00 + Name: wc1b_arath Len: 269 Check: 9544 Weight: 1.00 + Name: tipw_lyces Len: 269 Check: 9283 Weight: 1.00 + Name: wc1a_arath Len: 269 Check: 598 Weight: 1.00 + Name: tipw_pea Len: 269 Check: 9253 Weight: 1.00 + Name: tipa_arath Len: 269 Check: 6544 Weight: 1.00 + Name: aqua_atrca Len: 269 Check: 2848 Weight: 1.00 + Name: dip_antma Len: 269 Check: 9619 Weight: 1.00 + Name: aqpz_ecoli Len: 269 Check: 5641 Weight: 1.00 + Name: tip2_tobac Len: 269 Check: 490 Weight: 1.00 + Name: tip1_tobac Len: 269 Check: 622 Weight: 1.00 + Name: tipg_arath Len: 269 Check: 3231 Weight: 1.00 + Name: bib_drome Len: 269 Check: 7687 Weight: 1.00 + Name: tipr_arath Len: 269 Check: 4476 Weight: 1.00 + Name: tipa_phavu Len: 269 Check: 5563 Weight: 1.00 + Name: tipg_orysa Len: 269 Check: 3537 Weight: 1.00 + +// + + 1 50 +predict_h258 MASEIKKKLF WRAVVAEFLA MTLFVFISIG SALGFNYPLE RNQTLVQDNV +aqp1_rat MASEIKKKLF WRAVVAEFLA MTLFVFISIG SALGFNYPLE RNQTLVQDNV +aqp1_mouse MASEIKKKLF WRAVVAEFLA MTLFVFISIG SALGFNYPLE RNQTLVQDNV +aqp1_human MASEFKKKLF WRAVVAEFLA TTLFVFISIG SALGFKYPVG NNQTAVQDNV +aqp1_bovin MASEFKKKLF WRAVVAEFLA MILFIFISIG SALGFHYPIK SNQTtvQDNV +aqp1_sheep MASEFKKKLF WRAVVAEFLA MILFIFISIG SALGFHYPIK SNQTtvQDNV +aqpa_ranes MASEFKKKAF WRAVIAEFLA MILFVFISIG AALGFNFPIE EKANQtqDIV +aqp2_dasno ......SVAF SRAVLAEFLA TLIFVFFGLG SALSWPQALP S.......VL +aqp2_bovin ......SIAF SRAVLAEFLA TLLFVFFGLG SALNWPQALP S.......VL +aqp2_canfa ......SVAF SRAVFAEFLA TLLFVFFGLG SALNWPQALP S.......VL +aqp2_rabit ......SIAF SRAVFAEFLA TLLFVFFGLG SALNWPSALP S.......TL +aqp2_elema ......SIAF SRAVFSEFLA TLLFVFFGLG SALNWPQALP S.......VL +aqp2_horse ......SIAF SRAVLAEFLA TLLFVFFGLG SALNWPQAMP S.......VL +aqp2_proha ......SIAF SRAVLSEFLA TLLFVFFGLG SALNWPQALP S.......VL +mip_rat ...ELRSASF WRAIFAEFFA TLFYVFFGLG SSLRWA.... ...PGPLHVL +aqp2_oryaf ......SIAF SKAVFSEFLA TLLFVFFGLG SALNWPQALP S.......GL +mip_mouse .MWELRSASF WRAIFAEFFA TLFYVFFGLG ASLRWA.... ...PGPLHVL +mip_ranpi .MWEFRSFSF WRAVFAEFFG TMFYVFFGLG ASLKWAAGPA .......NVL +mip_bovin .MWELRSASF WRAICAEFFA SLFYVFFGLG ASLRWA.... ...PGPLHVL +mip_human .MWELRSASF WRAIFAEFFA TLFYVFFGLG SSLRWA.... ...PGPLHVL +mip_chick .......... .......... .......... .......... .......... +aqp5_rat MKKEVCSLAF FKAVFAEFLA TLIFVFFGLG SALKWPSALP T.......IL +aqp5_human MKKEVCSVAF LKAVFAEFLA TLIFVFFGLG SALKWPSALP T.......IL +aqp2_human .MWELRSIAF SRAVFAEFLA TLLFVFFGLG SALNWPQALP S.......VL +aqp4_human AFKGVWTQAF WKAVTAEFLA MLIFVLLSLG STINWG...G TEKPLPVDMV +aqp4_rat AFKGVWTQAF WKAVTAEFLA MLIFVLLSVG STINWG...G SENPLPVDMV +aqp4_mouse AFKGVWTQAF WKAVSAEFLA TLIFVL.GVG STINWG...G SENPLPVDMV +aqp2_rat .MWELRSIAF SRAVLAEFLA TLLFVFFGLG SALQWASSPP S.......VL +aqp2_mouse .MWELRSIAY CRAVLAEFLA TLLFVFFGLG SALQWASSPP S.......VL +wc2a_arath DGAELKKWSF YRAVIAEFVA TLLFLYITVL TVIGYKIQSD TDAGGVdgIL +aqp6_human MLACRLWKAI SRALFAEFLA TGLYVFFGVG SVMRWPTALP S.......VL +wc2c_arath DAEELTKWSL YRAVIAEFVA TLLFLYVTVL TVIGYKIQSD TKAGGVdgIL +wc2b_arath DADELTKWSL YRAVIAEFVA TLLFLYITVL TVIGYKIQSD TKAGGVdgIL +wc1c_arath EPGELSSWSF YRAGIAEFIA TFLFLYITVL TVMGVKRA.. PNMCASVGIQ +wc1b_arath EPGELASWSF WRAGIAEFIA TFLFLYITVL TVMGVKR..S PNMCASVGIQ +tipw_lyces EPGELSSWSF YRAGIAEFMA TFLFLYITIL TVMGLKRSDS LCSSV..GIQ +wc1a_arath EPGELSSWSF WRAGIAEFIA TFLFLYITVL TVMGVKR..S PNMCASVGIQ +tipw_pea EPSELTSWSF YRAGIAEFIA TFLFLYITVL TVMGVVRESS KCKTV..GIQ +tipa_arath RADEATHPDS IRATLAEFLS TFVFVFAAEG SILSLDKLYW EHAAHAGTni +aqua_atrca DMGELKLWSF WRAAIAEFIA TLLFLYITVA TVIGYKKETD PCASVGL..L +dip_antma SIGDSFSVAS IKAYVAEFIA TLLFVFAGVG SAIAYNKLTS DAALDPAGLV +aqpz_ecoli .........M FRKLAAECFG TFWLVFGGCG SAVLAAGFPE ....LGIGFA +tip2_tobac SIGDSFSVGS LKAYVAEFIA TLLFVFAGVG SAIAYNKLTA DAALDPAGLV +tip1_tobac SIGDSFSVGS LKAYVAEFIA TLLFVFAGVG SAIAYNKLTA DAALDPAGLV +tipg_arath RPDEATRPDA LKAALAEFIS TLIFVVAGSG SGMAFNKLTE NGATTPSGLV +bib_drome MQAEIRTLEF WRSIISECLA SFMYVFIVCG AAAGVGVGAS VSSVL....L +tipr_arath RPDEATRPDA LKAALAEFIS TLIFVVAGSG SGMAFNKLTE NGATTPSGLV +tipa_phavu RTDEATHPDS MRASLAEFAS TFIFVFAGEG SGLALVKIYQ DSAFSAGELL +tipg_orysa SHQEVYHPGA LKAALAEFIS TLIFVFAGQG SGMAFSKLTG GGATTPAGLI + + 51 100 +predict_h258 KVSLAFGLSI ATLAQSVGHI SGAHSNPAVT LGLLLSCQIS ILRAVMYIIA +aqp1_rat KVSLAFGLSI ATLAQSVGHI SGAHSNPAVT LGLLLSCQIS ILRAVMYIIA +aqp1_mouse KVSLAFGLSI ATLAQSVGHI SGAHLNPAVT LGLLLSCQIS ILRAVMYIIA +aqp1_human KVSLAFGLSI ATLAQSVGHI SGAHLNPAVT LGLLLSCQIS IFRALMYIIA +aqp1_bovin KVSLAFGLSI ATLAQSVGHI SGAHLNPAVT LGLLLSCQIS VLRAIMYIIA +aqp1_sheep KVSLAFGLSI ATLAQSVGHI SGAHLNPAVT LGLLLSCQIS ILRAIMYIIA +aqpa_ranes KVSLAFGISI ATMAQSVGHV SGAHLNPAVT LGCLLSCQIS ILKAVMYIIA +aqp2_dasno QIALAFGLAI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +aqp2_bovin QIAMAFGLAI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAVFYVAA +aqp2_canfa QIAMAFGLGI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +aqp2_rabit QIAMAFGLGI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +aqp2_elema QIAMAFGLAI GTLVQTLGHI SGAHINPAVT VACLVGCHVS FLRATFYLAA +aqp2_horse QIAMAFGLAI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +aqp2_proha QIAMAFGLAI GTLVQTLGHI SGAHINPAVT IACLVGCHVS FLRALFYLAA +mip_rat QVALAFGLAL ATLVQTVGHI SGAHVNPAVT FAFLVGSQMS LLRAFCYIAA +aqp2_oryaf QIAMAFGLAI GTLVQTLGHI SGAHINPAVT VACLVGCHVS FLRAIFYVAA +mip_mouse QVALAFGLAL ATLVQTVGHI SGAHVNPAVT FAFLVGSQMS LLRAFCYIAA +mip_ranpi VIALAFGLVL ATMVQSIGHV SGAHINPAVT FAFLIGSQMS LFRAIFYIAA +mip_bovin QVALAFGLAL ATLVQAVGHI SGAHVNPAVT FAFLVGSQMS LLRAICYMVA +mip_human QVAMAFGLAL ATLVQSVGHI SGAHVNPAVT FAFLVGSQMS LLRAFCYMAA +mip_chick .......... .......... .......... .......... .......... +aqp5_rat QISIAFGLAI GTLAQALGPV SGGHINPAIT LALLIGNQIS LLRAVFYVAA +aqp5_human QIALAFGLAI GTLAQALGPV SGGHINPAIT LALLVGNQIS LLRAFFYVAA +aqp2_human QIAMAFGLGI GTLVQALGHI SGAHINPAVT VACLVGCHVS VLRAAFYVAA +aqp4_human LISLCFGLSI ATMVQCFGHI SGGHINPAVT VAMVCTRKIS IAKSVFYIAA +aqp4_rat LISLCFGLSI ATMVQCFGHI SGGHINPAVT VAMVCTRKIS IAKSVFYITA +aqp4_mouse LISLCFGLSI ATMVQCLGHI SGGHINPAVT VAMVCTRKIS IAKSVFYIIA +aqp2_rat QIAVAFGLGI GILVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +aqp2_mouse QIAVAFGLGI GTLVQALGHV SGAHINPAVT VACLVGCHVS FLRAAFYVAA +wc2a_arath GIAWAFGGMI FILVYCTAGI SGGHINPAVT FGLFLARKVS LPRALLYIIA +aqp6_human QIAITFNLVT AMAVQVTWKT SGAHANPAVT LAFLVGSHIS LPRAVAYVAA +wc2c_arath GIAWAFGGMI FILVYCTAGI SGGHINPAVT FGLFLARKVS LIRAVLYMVA +wc2b_arath GIAWAFGGMI FILVYCTAGI SGGHINPAVT FGLFLARKVS LIRAVLYMVA +wc1c_arath GIAWAFGGMI FALVYCTAGI SGGHINPAVT FGLFLARKLS LTRAVFYIVM +wc1b_arath GIAWAFGGMI FALVYCTAGI SGGHINPAVT FGLFLARKLS LTRAVYYIVM +tipw_lyces GVAWAFGGMI FALVYCTAGI SGGHINPAVT FGLFLARKLS LTRAVFYMVM +wc1a_arath GIAWAFGGMI FALVYCTAGI SGGHINPAVT FGLFLARKLS LTRALYYIVM +tipw_pea GIAWAFGGMI FALVYCTAGI SGGHINPAVT FGLFLARKLS LTRAIFYMVM +tipa_arath LVALAHAFAL FAAVSAAINV SGGHVNPAVT FGALVGGRVT AIRAIYYWIA +aqua_atrca GIAWSFGGMI FVLVYCTAGI SGGHINPAVT FGLFLARKVS LLRALVYMIA +dip_antma AVAVAHAFAL FVGVSMAANV SGGHLNPAVT LGLAVGGNIT ILTGLFYWIA +aqpz_ecoli GVALAFGLTV LTMAFAVGHI SGGHFNPAVT IGLWAGGRFP AKEVVGYVIA +tip2_tobac AVAVAHAFAL FVGVSIAANI SGGHLNPAVT LGLAVGGNIT ILTGFFYWIA +tip1_tobac AVAVAHAFAL FVGVSIAANI SGGHLNPAVT LGLAVGGNIT ILTGFFYWIA +tipg_arath AAAVAHAFGL FVAVSVGANI SGGHVNPAVT FGAFIGGNIT LLRGILYWIA +bib_drome ATALASGLAM ATLTQCFLHI SGAHINPAVT LALCVVRSIS PIRAAMYITA +tipr_arath AAAVAHAFGL FVAVSVGANI SGGHVNPAVT FGAFIGGNIT LLRGILYWIA +tipa_phavu ALALAHAFAL FAAVSASMHV SGGHVNPAVS FGALIGGRIS VIRAVYYWIA +tipg_orysa AAAVAHAFAL FVAVSVGANI SGGHVNPAVT FGAFVGGNIT LFRGLLYWIA + + 101 150 +predict_h258 QCVGAIVASA ILSGITSSLL ENSLGRNDLA RGVNSGQGLG IEIIGTLQLV +aqp1_rat QCVGAIVASA ILSGITSSLL ENSLGRNDLA RGVNSGQGLG IEIIGTLQLV +aqp1_mouse QCVGAIVATA ILSGITSSLV DNSLGRNDLA HGVNSGQGLG IEIIGTLQLV +aqp1_human QCVGAIVATA ILSGITSSLT GNSLGRNDLA DGVNSGQGLG IEIIGTLQLV +aqp1_bovin QCVGAIVATA ILSGITSSLP DNSLGLNALA PGVNSGQGLG IEIIGTLQLV +aqp1_sheep QCVGAIVATV ILSGITSSLP DNSLGLNALA PGVNSGQGLG IEIIGTLQLV +aqpa_ranes QCLGAVVATA ILSGITSGLE NNSLGLNGLS PGVSAGQGLG VEILVTFQLV +aqp2_dasno QLLGAVAGAA ILHEITPPDV RG........ .......... .......... +aqp2_bovin QLLGAVAGAA LLHEITPPAI RG........ .......... .......... +aqp2_canfa QLLGAVAGAA LLHEITPPHV RG........ .......... .......... +aqp2_rabit QLLGAVAGAA LLHEITPAEV RG........ .......... .......... +aqp2_elema QLLGAVAGAA LLHELTPPDI RG........ .......... .......... +aqp2_horse QLLGAVAGAA LLHEITPPDI RR........ .......... .......... +aqp2_proha QLLGAVAGAA LLHELTPPDI RG........ .......... .......... +mip_rat QLLGAVAGAA VLYSVTPPAV RGNLALNTLH AGVSVGQATT VEIFLTLQFV +aqp2_oryaf QLLGAVAGAA LLHELTPPDI RG........ .......... .......... +mip_mouse QLLGAVAGAA VLYSVTPPAV RGNLALNTLH TGVSVGQATT VEIFLTLQFV +mip_ranpi QLLGAVAGAA VLYGVTPAAI RGNLALNTLH PGVSLGQATT VEIFLTLQFV +mip_bovin QLLGAVAGAA VLYSVTPPAV RGNLALNTLH PGVSVGQATI VEIFLTLQFV +mip_human QLLGAVAGAA VLYSVTPPAV RGNLALNTLH PAVSVGQATT VEIFLTLQFV +mip_chick .......... .......... .......... .......... .......... +aqp5_rat QLVGAIAGAG ILYWLAPLNA RGNLAVNALN NNTTPGKAMV VELILTFQLA +aqp5_human QLVGAIAGAG ILYGVAPLNA RGNLAVNALN NNTTQGQAMV VELILTFQLA +aqp2_human QLLGAVAGAA LLHEITPADI RGDLAVNALS NSTTAGQAVT VELFLTLQLV +aqp4_human QCLGAIIGAG ILYLVTPPSV VGGLGVTMVH GNLTAGHGLL VELIITFQLV +aqp4_rat QCLGAIIGAG ILYLVTPPSV VGGLGVTTVH GNLTAGHGLL VELIITFQLV +aqp4_mouse QCLGAIIGAG ILYLVTPPSV VGGLGVTTVH GNLTAGHGLL VELIITFQLV +aqp2_rat QLLGAVAGAA ILHEITPVEI RGDLAVNALH NNATAGQAVT VELFLTMQLV +aqp2_mouse QLLGAVAGAA ILHEITPVEI RGDLAVNALH NNATAGQAVT VELFLTMQLV +wc2a_arath QCLGAICGVG FVKAFQSSYY TRYGGgnSLA DGYSTGTGLA AEIIGTFVLV +aqp6_human QLVGATVGAA LLYGVMPGDI RETLGINVVR NSVSTGQAVA VELLLTLQLV +wc2c_arath QCLGAICGVG FVKAFQSSHY VNYGGgnFLA DGYNTGTGLA AEIIGTFVLV +wc2b_arath QCLGAICGVG FRQSFQSSYY DRYGGgnSLA DGYNTGTGLA AEIIGTFVLV +wc1c_arath QCLGAICGAG VVKGFQPNPY QtgGGANTVA HGYTKGSGLG AEIIGTFVLV +wc1b_arath QCLGAICGAG VVKGFQPKQY QagGGANTIA HGYTKGSGLG AEIIGTFVLV +tipw_lyces QCLGAICGAG VVKGFMVGPY QrgGGANVVN PGYTKGDGLG AEIIGTFVLV +wc1a_arath QCLGAICGAG VVKGFQPKQY QagGGANTVA HGYTKGSGLG AEIIGTFVLV +tipw_pea QVLGAICGAG VVKGFEGKQR FGDLNgnFVA PGYTKGDGLG AEIVGTFILV +tipa_arath QLLGAILACL LLRLTTNGMR PVGFR...LA SGVGAVNGLV LEIILTFGLV +aqua_atrca QCAGAICGVG LVKAFMKGPY NqgGGANSVA LGYNKGTAFG AELIGTFVLV +dip_antma QCLGSTVACL LLKFVTNGL. ..SVPTHGVA AGMDAIQGVV MEIIITFALV +aqpz_ecoli QVVGGIVAAA LLYLIASGKT GFDAAASGFA sgYSMLSALV VELVLSAGFL +tip2_tobac QLLGSTVACL LLKYVTNGL. ..AVPTHGVA AGLNGFQGVV MEIIITFALV +tip1_tobac QLLGSTVACL LLKYVTNGL. ..AVPTHGVA AGLNGLQGVV MEIIITFALV +tipg_arath QLLGSVVACL ILKFATGGLA VPAFG...LS AGVGVLNAFV FEIVMTFGLV +bib_drome QCGGGIAGAA LLYGVTVPGY QGNLQAasHS AALAAWERFG VEFILTSLVV +tipr_arath QLLGSVVACL ILKFATGGLA VPPFG...LS AGVGVLNAFV FEIVMTFGLV +tipa_phavu QLLGSIVAAL VLRLVTNNMR PSGF...HVS PGVGVGHMFI LEVVMTFGLM +tipg_orysa QLLGSTVACF LLRFSTGGLA TGTFGL.... TGVSVWEALV LEIVMTFGLV + + 151 200 +predict_h258 LCVLATTDRR RRDLGGSAPL AIGLSVALGH LLAIDYTGCG INPARSFGSA +aqp1_rat LCVLATTDRR RRDLGGSAPL AIGLSVALGH LLAIDYTGCG INPARSFGSA +aqp1_mouse LCVLATTDRR RRDLGGSAPL AIGLSVALGH LLAIDYTGCG INPARSFGSA +aqp1_human LCVLATTDRR RRDLGGSAPL AIGLSVALGH LLAIDYTGCG INPARSFGSA +aqp1_bovin LCVLATTDRR RRDLGGSGPL AIGFSVALGH LLAIDYTGCG INPARSFGSS +aqp1_sheep LCVLATTDRR RrdLGDSGPL AIGFSVALGH LLAIDYTGCG INPARSFGSS +aqpa_ranes LCVVAVTDRR RHDVSGSVPL AIGLSVALGH LIAIDYTGCG MNPARSFGSA +aqp2_dasno .......... .......... .......... .......... .......... +aqp2_bovin .......... .......... .......... .......... .......... +aqp2_canfa .......... .......... .......... .......... .......... +aqp2_rabit .......... .......... .......... .......... .......... +aqp2_elema .......... .......... .......... .......... .......... +aqp2_horse .......... .......... .......... .......... .......... +aqp2_proha .......... .......... .......... .......... .......... +mip_rat LCIFATYDER RNGRMGSVAL AVGFSLTLGH LFGMYYTGAG MNPARSFAPA +aqp2_oryaf .......... .......... .......... .......... .......... +mip_mouse LCIFATYDER RNGRMGSVAL AVGFSLTLGH LFGMYYTGAG MNPARSFAPA +mip_ranpi LCIFATYDER RNGRLGSVSL AIGFSLTLGH LFGLYYTGAS MNPARSFAPA +mip_bovin LCIFATYDER RNGRLGSVAL AVGFSLTLGH LFGMYYTGAG MNPARSFAPA +mip_human LCIFATYDER RNGQLGSVAL AVGFSLALGH LFGMYYTGAG MNPARSFAPA +mip_chick ........DR HDGRPGSAAL PVGFSLALGH LFGIPFTGAG MNPARSFAPA +aqp5_rat LCIFSSTDSR RTSPVGSPAL SIGLSVTLGH LVGIYFTGCS MNPARSFGPA +aqp5_human LCIFASTDSR RTSPVGSPAL SIGLSVTLGH LVGIYFTGCS MNPARSFGPA +aqp2_human LCIFASTDER RGENPGTPAL SIGFSVALGH LLGIHYTGCS MNPARSLAPA +aqp4_human FTIFASCDSK RTDVTGSIAL AIGFSVAIGH LFAINYTGAS MNPARSFGPA +aqp4_rat FTIFASCDSK RTDVTGSVAL AIGFSVAIGH LFAINYTGAS MNPARSFGPA +aqp4_mouse FTVFASCDSK RTDVTGSIAL AIGFSVAIGH LFAINYTGAS MNPARSFGPA +aqp2_rat LCIFASTDER RGDNLGSPAL SIGFSVTLGH LLGIYFTGCS MNPARSLAPA +aqp2_mouse LCIFASTDER RSDNLGSPAL SIGFSVTLGH LLGIYFTGCS MNPARSLAPA +wc2a_arath YTVFSATDPK RSavPVLAPL PIGFAVFMVH LATIPITGTG INPARSFGAA +aqp6_human LCVFASTDSR QTS..GSPAT MIGISWALGH LIGILFTGCS MNPARSFGPA +wc2c_arath YTVFSATDPK RNavPVLAPL PIGFAVFMVH LATIPITGTG INPARSFGAA +wc2b_arath YTVFSATDPK RNavPVLAPL PIGFAVFMVH LATIPITGTG INPARSFGAS +wc1c_arath YTVFSATDAK RSavPILAPL PIGFAVFLVH LATIPITGTG INPARSLGAA +wc1b_arath YTVFSATDAK RNavPILAPL PIGFAVFLVH LATIPITGTG INPARSLGAA +tipw_lyces YTVFSATDAK RNavPILAPL PIGFAVFLVH LATIPITGTG INPARSLGAA +wc1a_arath YTVFSATDAK RNavPILAPL PIGFAVFLVH LATIPITATG INPARSLGAA +tipw_pea YTVFSATDAK RSavPILAPL PIGFAVFLVH LATIPITGTG INPARSLGAA +tipa_arath YVVYStiDPK RGSLGIIAPL AIGLIVGANI LVGGPFSGAS MNPARAFGPA +aqua_atrca YTVFSATDPK RSavPILAPL PIGFAVFMVH LATIPITGTG INPARSFGAA +dip_antma YTVYAtaDPK KGSLGVIAPI AIGFIVGANI LAAGPFSGGS MNPARSFGPA +aqpz_ecoli LVIHGATDKF APA..GFAPI AIGLALTLIH LISIPVTNTS VNPARSTAVA +tip2_tobac YTVYAtaDPK KGSLGTIAPI AIGFIVGANI LAAGPFSGGS MNPARSFGPA +tip1_tobac YTVYAtaDPK KGSLGTIAPI AIGFIVGANI LAAGPFSGGS MNPARSFGPA +tipg_arath YTVYAtiDPK NGSLGTIAPI AIGFIVGANI LAGGAFSGAS MNPAVAFGPA +bib_drome LCYFVSTDPM KKFMGNS.AA SIGCAYSACC FVSMPYLN.. ..PARSLGPS +tipr_arath YTVYAtiDPK NGSLGTIAPI AIGFIVGANI LAGGAFSGAS MNPAVAFGPA +tipa_phavu YTVYGtiDPK RGAVSYIAPL AIGLIVGANI LVGGPFDGAC MNPALAFGPS +tipg_orysa YTVYAtvDPK KGSLGTIAPI AIGFIVGANI LVGGAFDGAS MNPAVSFGPA + + 201 250 +predict_h258 VLTRNFSNHW IFWVGPFIGS ALAVLIYDFI LAPRSSDFTD RMKVWTSGQV +aqp1_rat VLTRNFSNHW IFWVGPFIGS ALAVLIYDFI LAPRSSDFTD RMKVWTSGQV +aqp1_mouse VLTRNFSNHW IFWVGPFIGG ALAVLIYDFI LAPRSSDFTD RMKVWTSGQV +aqp1_human VITHNFSNHW IFWVGPFIGG ALAVLIYDFI LAPRSSDLTD RVKVWTSGQV +aqp1_bovin VITHNFQDHW IFWVGPFIGA ALAVLIYDFI LAPRSSDLTD RVKVWTSGQV +aqp1_sheep VITHNFQDHW IFWVGPFIGA ALAVLIYDFI LAPRSSDLTD RVKVWTSGQV +aqpa_ranes VLTKNFTYHW IFWVGPMIGG AAAAIIYDFI LAPRTSDLTD RMKVWTNGQV +aqp2_dasno .......... .......... .......... .......... .......... +aqp2_bovin .......... .......... .......... .......... .......... +aqp2_canfa .......... .......... .......... .......... .......... +aqp2_rabit .......... .......... .......... .......... .......... +aqp2_elema .......... .......... .......... .......... .......... +aqp2_horse .......... .......... .......... .......... .......... +aqp2_proha .......... .......... .......... .......... .......... +mip_rat ILTRNFSNHW VYWVGPIIGG GLGSLLYDFL LFPRLKSVSE RLSILKGARP +aqp2_oryaf .......... .......... .......... .......... .......... +mip_mouse ILTRNFSNHW VYWVGPIIGG GLGSLLYDFL LFPRLKSVSE RLSILKGARP +mip_ranpi VLTRNFTNHW VYWVGPIIGG ALGGLVYDFI LFPRMRGLSE RLSILKGARP +mip_bovin ILTRNFTNHW VYWVGPVIGA GLGSLLYDFL LFPRLKSVSE RLSILKGSRP +mip_human ILTGNFTNHW VYWVGPIIGG GLGSLLYDFL LFPRLKSISE RLSVLKGAKP +mip_chick VITRNFTNHW VFWAGPLLGA ALAALLYELA LCPRARSMAE RLAV.LRGEP +aqp5_rat VVMNRFssHW VFWVGPIVGA MLAAILYFYL LFPSSLSLHD RVAVVKGTYE +aqp5_human VVMNRFsaHW VFWVGPIVGA VLAAILYFYL LFPNSLSLSE RVAIIKGTYE +aqp2_human VVTGKFDDHW VFWIGPLVGA ILGSLLYNYV LFPPAKSLSE RLAVLKGLEp +aqp4_human VIMGNWENHW IYWVGPIIGA VLAGGLYEYV FCPDVEFKRR FKEAFSKaqT +aqp4_rat VIMGNWENHW IYWVGPIIGA VLAGALYEYV FCPDVELKRR LKEAFSKaqT +aqp4_mouse VIMGNWANHW IYWVGPIMGA VLAGALYEYV FCPDVELKRR LKEAFSKaqT +aqp2_rat VVTGKFDDHW VFWIGPLVGA IIGSLLYNYL LFPSAKSLQE RLAVLKGLEp +aqp2_mouse VVTGKFDDHW VFWIGPLVGA IIGSLLYNYL LFPSTKSLQE RLAVLKGLEp +wc2a_arath VIYnpWDDHW IFWVGPFIGA AIAAFYHQFV LRASGSKSLG SFRSAANV.. +aqp6_human IIIGKFTVHW VFWVGPLMGA LLASLIYNFV LFPDTKTLAQ RLAILTGTVE +wc2c_arath VIFnpWDDHW IFWVGPFIGA TIAAFYHQFV LRASGSKSLG SFRSAANV.. +wc2b_arath VIYnpWDDHW IFWVGPFIGA AIAAFYHQFV LRASGSKSLG SFRSAANV.. +wc1c_arath IIYnaWDDHW IFWVGPFIGA ALAALYHQLV IRAIPFKSRS .......... +wc1b_arath IIFnaWDDHW VFWVGPFIGA ALAALYHVIV IRAIPFKSRS .......... +tipw_lyces IIYnaWNDHW IFWVGPMIGA ALAAIYHQII IRAMPFHRS. .......... +wc1a_arath IIYnsWDDHW VFWVGPFIGA ALAALYHVVV IRAIPFKSRS .......... +tipw_pea IVFngWNDHW IFWVGPFIGA ALAALYHQVV IRAIPFKSK. .......... +tipa_arath LVGWRWHDHW IYWVGPFIGS ALAALIYEYM VIPTEPPTHH AHGVHQPLAP +aqua_atrca VIyrVWDDHW IFWVGPFVGA LAAAAYHQYV LRAAAIKALG SFRSNPTN.. +dip_antma VASGDFSQNW IYWAGPLIGG ALAGFIYGDV FITAHAPLPT SEDYA..... +aqpz_ecoli IFQgaLEQLW FFWVVPIVGG IIGGLIYRTL LEKRD..... .......... +tip2_tobac VVAGDFSQNW IYWAGPLIGG GLAGFIYGDV FIGCHTPLPT SEDYA..... +tip1_tobac VVAGDFSQNW IYWAGPLIGG GLAGFIYGDV FIGCHTPLPT SEDYA..... +tipg_arath VVSWTWTNHW VYWAGPLVGG GIAGLIYEVF FINTTHEQLP TTDY...... +bib_drome FVLNKWDSHW VYWFGPLVGG MASGLVYEYI FNSRNRNLRH NKGSIDNDSS +tipr_arath VVSWTWTNHW VYWAGPLVGG GIAGLIYEVF FINTTHTSSS NHRLLN.... +tipa_phavu LVGWQWHQHW IFWVGPLLGA ALAALVYEYA VIPIEPPPHH HQPLATEDY. +tipg_orysa LVSWSWESQW VYWVGPLIGG GLAGVIYEVL FISHTHEQLP TTDY...... + + 251 269 +predict_h258 EEYDLDADDI NSRVEMKPK +aqp1_rat EEYDLDADDI NSRVEMKPK +aqp1_mouse EEYDLDADDI NSRVEMKPK +aqp1_human EEYDLDADDI NSRVEMKPK +aqp1_bovin EEYDLDADDI NSRVEMKPK +aqp1_sheep EEYDLDADDI NSRVEMKPK +aqpa_ranes EEYELDGDD. NTRVEMKPK +aqp2_dasno .......... ......... +aqp2_bovin .......... ......... +aqp2_canfa .......... ......... +aqp2_rabit .......... ......... +aqp2_elema .......... ......... +aqp2_horse .......... ......... +aqp2_proha .......... ......... +mip_rat SDSNGQPEGT GEPVELKTQ +aqp2_oryaf .......... ......... +mip_mouse SDSNGQPEGT GEPVELKTQ +mip_ranpi AEPEGQQEAT GEPIELKTQ +mip_bovin SESNGQPEVT GEPVELKTQ +mip_human DVSNGQPEVT GEPVELNTQ +mip_chick PAAAPPPEPP AEPLELKTQ +aqp5_rat PEEDWEDHRE ERKKTIELT +aqp5_human PDEDWEEQRE ERKKTMELT +aqp2_human tDWEEREVRR RQSVELHSP +aqp4_human KGSYMEVEDN RSQVETDDL +aqp4_rat KGSYMEVEDN RSQVETEDL +aqp4_mouse KGSYMEVEDN RSQVETEDL +aqp2_rat tDWEEREVRR RQSVELHSP +aqp2_mouse tDWEEREVRR RQSVELHSP +wc2a_arath .......... ......... +aqp6_human VGTGARAGAE PLKKESQPG +wc2c_arath .......... ......... +wc2b_arath .......... ......... +wc1c_arath .......... ......... +wc1b_arath .......... ......... +tipw_lyces .......... ......... +wc1a_arath .......... ......... +tipw_pea .......... ......... +tipa_arath EDY....... ......... +aqua_atrca .......... ......... +dip_antma .......... ......... +aqpz_ecoli .......... ......... +tip2_tobac .......... ......... +tip1_tobac .......... ......... +tipg_arath .......... ......... +bib_drome SIHSEDELNY DMDMEKPNK +tipr_arath .......... ......... +tipa_phavu .......... ......... +tipg_orysa .......... ......... + +________________________________________________________________________________ + + Prediction of: + + - secondary structure, by PHDsec + - solvent accessibility, by PHDacc + - and helical transmembrane regions, by PHDhtm + + PHD: Profile fed neural network systems from HeiDelberg + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Author: Burkhard Rost + EMBL, Heidelberg, FRG + Meyerhofstrasse 1, 69 117 Heidelberg + Internet: Predict-Help@EMBL-Heidelberg.DE + + All rights reserved. + + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + Secondary structure prediction by PHDsec: + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Author: Burkhard Rost + EMBL, Heidelberg, FRG + Meyerhofstrasse 1, 69 117 Heidelberg + Internet: Rost@EMBL-Heidelberg.DE + + All rights reserved. + +About the network method +~~~~~~~~~~~~~~~~~~~~~~~ + +The network procedure is described in detail in: +1) Rost, Burkhard; Sander, Chris: + Prediction of protein structure at better than 70% accuracy. + J. Mol. Biol., 1993, 232, 584-599. + +A brief description is given in: + Rost, Burkhard; Sander, Chris: + Improved prediction of protein secondary structure by use of se- + quence profiles and neural networks. + Proc. Natl. Acad. Sci. U.S.A., 1993, 90, 7558-7562. + +The PHD mail server is described in: +2) Rost, Burkhard; Sander, Chris; Schneider, Reinhard: + PHD - an automatic mail server for protein secondary structure + prediction. + CABIOS, 1994, 10, 53-60. + +The latest improvement steps (up to 72%) are explained in: +3) Rost, Burkhard; Sander, Chris: + Combining evolutionary information and neural networks to predict + protein secondary structure. + Proteins, 1994, 19, 55-72. + +To be quoted for publications of PHD output: + Papers 1-3 for the prediction of secondary structure and the pre- + diction server. + +About the input to the network +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The prediction is performed by a system of neural networks. +The input is a multiple sequence alignment. It is taken from an HSSP +file (produced by the program MaxHom: + Sander, Chris & Schneider, Reinhard: Database of Homology-Derived + Structures and the Structural Meaning of Sequence Alignment. + Proteins, 1991, 9, 56-68. + +For optimal results the alignment should contain sequences with varying +degrees of sequence similarity relative to the input protein. +The following is an ideal situation: + ++-----------------+----------------------+ +| sequence: | sequence identity | ++-----------------+----------------------+ +| target sequence | 100 % | +| aligned seq. 1 | 90 % | +| aligned seq. 2 | 80 % | +| ... | ... | +| aligned seq. 7 | 30 % | ++-----------------+----------------------+ + +Estimated Accuracy of Prediction +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +A careful cross validation test on some 250 protein chains (in total +about 55,000 residues) with less than 25% pairwise sequence identity +gave the following results: + +++================++-----------------------------------------+ +|| Qtotal = 72.1% || ("overall three state accuracy") | +++================++-----------------------------------------+ + ++----------------------------+-----------------------------+ +| Qhelix (% of observed)=70% | Qhelix (% of predicted)=77% | +| Qstrand(% of observed)=62% | Qstrand(% of predicted)=64% | +| Qloop (% of observed)=79% | Qloop (% of predicted)=72% | ++----------------------------+-----------------------------+ +.......................................................................... + +These percentages are defined by: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +| number of correctly predicted residues +|Qtotal = --------------------------------------- (*100) +| number of all residues +| +| no of res correctly predicted to be in helix +|Qhelix (% of obs) = -------------------------------------------- (*100) +| no of all res observed to be in helix +| +| +| no of res correctly predicted to be in helix +|Qhelix (% of pred)= -------------------------------------------- (*100) +| no of all residues predicted to be in helix + +.......................................................................... + +Averaging over single chains +~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The most reasonable way to compute the overall accuracies is the above +quoted percentage of correctly predicted residues. However, since the +user is mainly interested in the expected performance of the prediction +for a particular protein, the mean value when averaging over protein +chains might be of help as well. Computing first the three state +accuracy for each protein chain, and then averaging over 250 chains +yields the following average: + ++-------------------------------====--+ +| Qtotal/averaged over chains = 72.2% | ++-------------------------------====--+ +| standard deviation = 9.3% | ++-------------------------------------+ + +.......................................................................... + +Further measures of performance +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +Matthews correlation coefficient: + ++---------------------------------------------+ +| Chelix = 0.63, Cstrand = 0.53, Cloop = 0.52 | ++---------------------------------------------+ +.......................................................................... + +Average length of predicted secondary structure segments: + +. +------------+----------+ +. | predicted | observed | ++-----------+------------+----------+ +| Lhelix = | 10.3 | 9.3 | +| Lstrand = | 5.0 | 5.3 | +| Lloop = | 7.2 | 5.9 | ++-----------+------------+----------+ +.......................................................................... + +The accuracy matrix in detail: + ++---------------------------------------+ +| number of residues with H, E, L | ++---------+------+------+------+--------+ +| |net H |net E |net L |sum obs | ++---------+------+------+------+--------+ +| obs H |12447 | 1255 | 3990 | 17692 | +| obs E | 949 | 7493 | 3750 | 12192 | +| obs L | 2604 | 2875 |19962 | 25441 | ++---------+------+------+------+--------+ +| sum Net |16000 |11623 |27702 | 55325 | ++---------+------+------+------+--------+ + +Note: This table is to be read in the following manner: + 12447 of all residues predicted to be in helix, were observed to + be in helix, 949 however belong to observed strands, 2604 to + observed loop regions. The term "observed" refers to the DSSP + assignment of secondary structure calculated from 3D coordinates + of experimentally determined structures (Dictionary of Secondary + Structure of Proteins: Kabsch & Sander (1983) Biopolymers, 22, + 2577-2637). + +Position-specific reliability index +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The network predicts the three secondary structure types using real +numbers from the output units. The prediction is assigned by choosing +the maximal unit ("winner takes all"). However, the real numbers +contain additional information. +E.g. the difference between the maximal and the second largest output +unit can be used to derive a "reliability index". This index is given +for each residue along with the prediction. The index is scaled to +have values between 0 (lowest reliability), and 9 (highest). +The accuracies (Qtot) to be expected for residues with values above a +particular value of the index are given below as well as the fraction +of such residues (%res).: + ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| index| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | +| %res |100.0| 99.2| 90.4| 80.9| 71.6| 62.5| 52.8| 42.3| 29.8| 14.1| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| | | | | | | | | | | | +| Qtot | 72.1| 72.3| 74.8| 77.7| 80.3| 82.9| 85.7| 88.5| 91.1| 94.2| +| | | | | | | | | | | | ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| H%obs| 70.4| 70.6| 73.7| 77.1| 80.1| 83.1| 86.0| 89.3| 92.5| 96.4| +| E%obs| 61.5| 61.7| 63.7| 66.6| 69.1| 71.7| 74.6| 77.0| 77.8| 68.1| +| | | | | | | | | | | | +| H%prd| 77.8| 78.0| 80.0| 82.6| 84.7| 86.9| 89.2| 91.3| 93.1| 95.4| +| E%prd| 64.5| 64.7| 67.8| 71.0| 74.2| 77.6| 81.4| 85.1| 89.8| 93.5| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ + +The above table gives the cumulative results, e.g. 62.5% of all +residues have a reliability of at least 5. The overall three-state +accuracy for this subset of almost two thirds of all residues is 82.9%. +For this subset, e.g., 83.1% of the observed helices are correctly +predicted, and 86.9% of all residues predicted to be in helix are +correct. + +.......................................................................... + +The following table gives the non-cumulative quantities, i.e. the +values per reliability index range. These numbers answer the question: +how reliable is the prediction for all residues labeled with the +particular index i. + ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| index| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | +| %res | 8.8| 9.5| 9.3| 9.1| 9.7| 10.5| 12.5| 15.7| 14.1| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| | | | | | | | | | | +| Qtot | 46.6| 50.6| 57.7| 62.6| 67.9| 74.2| 82.2| 88.3| 94.2| +| | | | | | | | | | | ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| H%obs| 36.8| 42.3| 49.5| 55.2| 61.7| 69.9| 78.8| 87.4| 96.4| +| E%obs| 44.7| 44.5| 52.1| 55.4| 60.9| 68.0| 75.9| 81.0| 68.1| +| | | | | | | | | | | +| H%prd| 49.9| 52.5| 60.3| 64.2| 69.2| 77.5| 85.4| 89.9| 95.4| +| E%prd| 41.7| 47.1| 53.6| 57.0| 64.0| 71.6| 78.8| 88.8| 93.5| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+ + +For example, for residues with Relindex = 5 64% of all predicted betha- +strand residues are correctly identified. + + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + Solvent accessibility prediction by PHDacc: + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Author: Burkhard Rost + EMBL, Heidelberg, FRG + Meyerhofstrasse 1, 69 117 Heidelberg + Internet: Rost@EMBL-Heidelberg.DE + + All rights reserved. + +About the network method +~~~~~~~~~~~~~~~~~~~~~~~ + +The network for prediction of secondary structure is described in +detail in: + Rost, Burkhard; Sander, Chris: + Prediction of protein structure at better than 70% accuracy. + J. Mol. Biol., 1993, 232, 584-599. + +The analysis of the prediction of solvent exposure is given in: + Rost, Burkhard; Sander, Chris: + Conservation and prediction of solvent accessibility in protein + families. Proteins, 1994, 20, 216-226. + +To be quoted for publications of PHD exposure prediction: + Both papers quoted above. + +Definition of accessibility +~~~~~~~~~~~~~~~~~~~~~~~~~~ + +For training the residue solvent accessibility the DSSP (Dictionary of +Secondary Structure of Proteins; Kabsch & Sander (1983) Biopolymers, 22, +2577-2637) values of accessible surface area have been used. The +prediction provides values for the relative solvent accessibility. The +normalisation is the following: + +| ACCESSIBILITY (from DSSP in Angstrom) +|RELATIVE_ACCESSIBILITY = ------------------------------------- * 100 +| MAXIMAL_ACC (amino acid type i) + +where MAXIMAL_ACC (i) is the maximal accessibility of amino acid type i. +The maximal values are: + ++----+----+----+----+----+----+----+----+----+----+----+----+ +| A | B | C | D | E | F | G | H | I | K | L | M | +| 106| 160| 135| 163| 194| 197| 84| 184| 169| 205| 164| 188| ++----+----+----+----+----+----+----+----+----+----+----+----+ +| N | P | Q | R | S | T | V | W | X | Y | Z | +| 157| 136| 198| 248| 130| 142| 142| 227| 180| 222| 196| ++----+----+----+----+----+----+----+----+----+----+----+ + +Notation: one letter code for amino acid, B stands for D or N; Z stands + for E or Q; and X stands for undetermined. + +The relative solvent accessibility can be used to estimate the number +of water molecules (W) in contact with the residue: + +W = ACCESSIBILITY /10 + +The prediction is given in 10 states for relative accessibility, with + +RELATIVE_ACCESSIBILITY = (PREDICTED_ACC * PREDICTED_ACC) + +where PREDICTED_ACC = 0 - 9. + +Estimated Accuracy of Prediction +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +A careful cross validation test on some 238 protein chains (in total +about 62,000 residues) with less than 25% pairwise sequence identity +gave the following results: + +Correlation +........... + +The correlation between observed and predicted solvent accessibility +is: + +----------- +corr = 0.53 +----------- + +This value ought to be compared to the worst and best case prediction +scenario: random prediction (corr = 0.0) and homology modelling +(corr = 0.66). (Note: homology modelling yields a relative accurate +prediction in 3D if, and only if, a significantly identical sequence +has a known 3D structure.) + +3-state accuracy +................ + +Often the relative accessibility is projected onto, e.g., 3 states: + b = buried (here defined as < 9% relative accessibility), + i = intermediate ( 9% <= rel. acc. < 36% ), + e = exposed ( rel. acc. >= 36% ). + +A projection onto 3 states or 2 states (buried/exposed) enables the +compilation of a 3- and 2-state prediction accuracy. PHD reaches an +overall 3-state accuracy of: + Q3 = 57.5% +(compared to 35% for random prediction and 70% for homology modelling). + +In detail: + ++-----------------------------------+-------------------------+ +| Qburied (% of observed)=77% | Qb (% of predicted)=60% | +| Qintermediate (% of observed)= 9% | Qi (% of predicted)=44% | +| Qexposed (% of observed)=78% | Qe (% of predicted)=56% | ++-----------------------------------+-------------------------+ + +10-state accuracy +................. + +The network predicts relative solvent accessibility in 10 states, with +state i (i = 0-9) corresponding to a relative solvent accessibility of +i*i %. The 10-state accuracy of the network is: + + Q10 = 24.5% + +.......................................................................... + +These percentages are defined by: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +| number of correctly predicted residues +|Q3 = --------------------------------------- (*100) +| number of all residues +| +| no of res. correctly predicted to be buried +|Qburied (% of obs) = ------------------------------------------- (*100) +| no of all res. observed to be buried +| +| +| no of res. correctly predicted to be buried +|Qburied (% of pred)= ------------------------------------------- (*100) +| no of all residues predicted to be buried + +.......................................................................... + +Averaging over single chains +~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The most reasonable way to compute the overall accuracies is the above +quoted percentage of correctly predicted residues. However, since the +user is mainly interested in the expected performance of the prediction +for a particular protein, the mean value when averaging over protein +chains might be of help as well. Computing first the correlation +between observed and predicted accessibility for each protein chan, and +then averaging over all 238 chains yields the following average: + ++-------------------------------====--+ +| corr/averaged over chains = 0.53 | ++-------------------------------====--+ +| standard deviation = 0.11 | ++-------------------------------------+ + +.......................................................................... + +Further details of performance accuracy +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The accuracy matrix in detail: +.............................. + +-------+----------------------------------------------------+----------- +\ PHD | 0 1 2 3 4 5 6 7 8 9 | SUM %obs +-------+----------------------------------------------------+----------- +OBS 0 | 8611 140 8 44 82 169 772 334 27 0 | 10187 16.6 +OBS 1 | 4367 164 0 50 106 231 738 346 44 3 | 6049 9.8 +OBS 2 | 3194 168 1 68 125 303 951 513 42 7 | 5372 8.7 +OBS 3 | 2760 159 8 80 136 327 1246 746 58 19 | 5539 9.0 +OBS 4 | 2312 144 2 72 166 396 1615 1245 124 19 | 6095 9.9 +OBS 5 | 1873 96 3 84 138 425 1979 1834 187 27 | 6646 10.8 +OBS 6 | 1387 67 1 60 80 278 2237 2627 231 51 | 7019 11.4 +OBS 7 | 1082 35 0 32 56 225 1871 3107 302 60 | 6770 11.0 +OBS 8 | 660 25 0 27 43 136 1206 2374 325 87 | 4883 7.9 +OBS 9 | 325 20 2 27 29 74 648 1159 366 214 | 2864 4.7 +-------+----------------------------------------------------+----------- +SUM |26571 1018 25 544 961 2564 13263 14285 1706 487 | +-------+----------------------------------------------------+----------- + +Note: This table is to be read in the following manner: + 8611 of all residues predicted to be in exposed by 0%, were + observed with 0% relative accessibility. However, 325 of all + residues predicted to have 0% are observed as completely exposed + (obs = 9 -> rel. acc. >= 81%). The term "observed" refers to the + DSSP compilation of area of solvent accessibility calculated from + 3D coordinates of experimentally determined structures (Diction- + ary of Secondary Structure of Proteins: Kabsch & Sander (1983) + Biopolymers, 22, 2577-2637). + +Accuracy for each amino acid: +............................. + ++---+------------------------------+-----+-------+------+ +|AA | Q3 b%o b%p i%o i%p e%o e%p | Q10 | corr | N | ++---+------------------------------+-----+-------+------+ +| A | 59.0 87 60 2 38 66 57 | 31 | 0.530 | 5054 | +| C | 62.0 91 67 5 39 25 21 | 34 | 0.244 | 893 | +| D | 56.5 21 45 6 49 94 57 | 20 | 0.321 | 3536 | +| E | 60.8 9 40 3 41 98 61 | 21 | 0.347 | 3743 | +| F | 63.3 94 67 9 46 29 37 | 27 | 0.366 | 2436 | +| G | 52.1 75 51 1 31 67 53 | 22 | 0.405 | 4787 | +| H | 50.9 63 53 23 45 71 50 | 18 | 0.442 | 1366 | +| I | 64.9 95 68 6 41 30 38 | 34 | 0.360 | 3437 | +| K | 66.6 2 11 2 37 98 67 | 23 | 0.267 | 3652 | +| L | 61.6 93 65 8 44 31 40 | 31 | 0.368 | 5016 | +| M | 60.1 92 64 5 39 45 44 | 29 | 0.452 | 1371 | +| N | 55.5 45 45 8 38 87 59 | 17 | 0.410 | 2923 | +| P | 53.0 48 48 9 39 83 56 | 18 | 0.364 | 2920 | +| Q | 54.3 27 44 7 44 92 56 | 20 | 0.344 | 2225 | +| R | 49.9 15 47 36 47 76 51 | 18 | 0.372 | 2765 | +| S | 55.6 69 53 3 51 81 56 | 22 | 0.464 | 3981 | +| T | 51.8 61 51 8 38 78 53 | 21 | 0.432 | 3740 | +| V | 61.1 93 65 5 40 39 42 | 34 | 0.418 | 4156 | +| W | 56.2 85 62 20 49 29 27 | 21 | 0.318 | 891 | +| Y | 49.7 73 52 33 49 36 38 | 19 | 0.359 | 2301 | ++---+------------------------------+-----+-------+------+ + +Abbreviations: + +AA: amino acid in one-letter code +b%o, i%o, e%o: = Qburied, Qintermediate, Qexposed (% of observed), + i.e. percentage of correct prediction in each state, see above +b%p, i%p, e%p: = Qburied, Qintermediate, Qexposed (% of predicted), + i.e. probability of correct prediction in each state, see above +b%o: = Qburied (% of observed), see above +Q10: percentage of correctly predicted residues in each of the 10 + states of predicted relative accessibility. +corr: correlation between predicted and observed rel. acc. +N: number of residues in data set + +Accuracy for different secondary structure: +........................................... + ++--------+------------------------------+----+-------+-------+ +| type | Q3 b%o b%p i%o i%p e%o e%p |Q10 | corr | N | ++--------+------------------------------+----+-------+-------+ +| helix | 59.5 79 64 8 44 80 56 | 27 | 0.574 | 20100 | +| strand | 61.3 84 73 9 46 69 37 | 35 | 0.524 | 13356 | +| loop | 54.4 64 43 11 44 78 61 | 18 | 0.442 | 27968 | ++--------+------------------------------+----+-------+-------+ + +Abbreviations as before. + +Position-specific reliability index +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The network predicts the 10 states for relative accessibility using real +numbers from the output units. The prediction is assigned by choosing +the maximal unit ("winner takes all"). However, the real numbers +contain additional information. +E.g. the difference between the maximal and the second largest output +unit (with the constraint that the second largest output is compiled +among all units at least 2 positions off the maximal unit) can be used +to derive a "reliability index". This index is given for each residue +along with the prediction. The index is scaled to have values between +0 (lowest reliability), and 9 (highest). +The accuracies (Q3, corr, asf.) to be expected for residues with values +above a particular value of the index are given below as well as the +fraction of such residues (%res).: + ++---+------------------------------+----+-------+-------+ +|RI | Q3 b%o b%p i%o i%p e%o e%p |Q10 | corr | %res | ++---+------------------------------+----+-------+-------+ +| 0 | 57.5 77 60 9 44 78 56 | 24 | 0.535 | 100.0 | +| 1 | 59.1 76 63 9 45 82 57 | 25 | 0.560 | 91.2 | +| 2 | 61.7 79 66 4 47 87 58 | 27 | 0.594 | 77.1 | +| 3 | 66.6 87 70 1 51 89 63 | 30 | 0.650 | 57.1 | +| 4 | 70.0 89 72 0 83 91 67 | 32 | 0.686 | 45.8 | +| 5 | 72.9 92 75 0 0 93 70 | 34 | 0.722 | 35.6 | +| 6 | 76.3 95 77 0 0 93 75 | 36 | 0.769 | 24.7 | +| 7 | 79.0 97 79 0 0 93 78 | 39 | 0.803 | 16.0 | +| 8 | 80.9 98 80 0 0 91 81 | 43 | 0.824 | 9.6 | +| 9 | 81.2 99 80 0 0 88 83 | 45 | 0.828 | 5.9 | ++---+------------------------------+----+-------+-------+ + +Abbreviations as before. + +The above table gives the cumulative results, e.g. 45.8% of all +residues have a reliability of at least 4. The correlation for this +most reliably predicted half of the residues is 0.686, i.e. a value +comparable to what could be expected if homology modelling were +possible. For this subset of 45.8% of all residues, 89% of the buried +residues are correctly predicted, and 72% of all residues predicted to +be buried are correct. + +.......................................................................... + +The following table gives the non-cumulative quantities, i.e. the +values per reliability index range. These numbers answer the question: +how reliable is the prediction for all residues labeled with the +particular index i. + ++---+------------------------------+----+-------+-------+ +|RI | Q3 b%o b%p i%o i%p e%o e%p |Q10 | corr | %res | ++---+------------------------------+----+-------+-------+ +| 0 | 40.9 79 40 16 41 21 40 | 14 | 0.175 | 8.8 | +| 1 | 45.4 61 46 28 44 48 44 | 17 | 0.278 | 14.1 | +| 2 | 47.4 53 52 10 46 80 44 | 19 | 0.343 | 19.9 | +| 3 | 52.9 75 59 4 50 77 47 | 23 | 0.439 | 11.4 | +| 4 | 60.0 81 63 0 83 84 56 | 25 | 0.547 | 10.1 | +| 5 | 65.2 82 70 0 0 93 62 | 28 | 0.607 | 10.9 | +| 6 | 71.3 90 72 0 0 94 70 | 31 | 0.692 | 8.8 | +| 7 | 76.0 94 76 0 0 95 75 | 34 | 0.762 | 6.3 | +| 8 | 80.5 97 81 0 0 94 79 | 39 | 0.808 | 3.8 | +| 9 | 81.2 99 80 0 0 88 83 | 45 | 0.828 | 5.9 | ++---+------------------------------+----+-------+-------+ + +For example, for residues with RI = 4 83% of all predicted intermediate +residues are correctly predicted as such. + + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + Prediction of helical transmembrane segments by PHDhtm: + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Author: Burkhard Rost + EMBL, Heidelberg, FRG + Meyerhofstrasse 1, 69 117 Heidelberg + Internet: Rost@EMBL-Heidelberg.DE + + All rights reserved. + +About the network method +~~~~~~~~~~~~~~~~~~~~~~~ + +The PHD mail server is described in: + Rost, Burkhard; Sander, Chris; Schneider, Reinhard: + PHD - an automatic mail server for protein secondary structure + prediction. + CABIOS, 1994, 10, 53-60. + +To be quoted for publications of PHDhtm output: + Rost, Burkhard; Casadio, Rita; Fariselli, Piero; Sander, Chris: + Prediction of helical transmembrane segments at 95% accuracy. + Protein Science, 1995, 4, 521-533. + +Estimated Accuracy of Prediction +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +A cross validation test on 69 helical trans-membrane proteins (in total +about 30,000 residues) with less than 25% pairwise sequence identity +gave the following results: + +++================++-----------------------------------------+ +|| Qtotal = 94.7% || ("overall two state accuracy") | +++================++-----------------------------------------+ + ++----------------------------+-----------------------------+ +| Qhelix (% of observed)=92% | Qhelix (% of predicted)=83% | +| Qloop (% of observed)=96% | Qloop (% of predicted)=97% | ++----------------------------+-----------------------------+ + +.......................................................................... + +These percentages are defined by: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +| number of correctly predicted residues +|Qtotal = --------------------------------------- (*100) +| number of all residues +| +| no of res correctly predicted to be in helix +|Qhelix (% of obs) = -------------------------------------------- (*100) +| no of all res observed to be in helix +| +| +| no of res correctly predicted to be in helix +|Qhelix (% of pred)= -------------------------------------------- (*100) +| no of all residues predicted to be in helix + +.......................................................................... + +Further measures of performance +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +Matthews correlation coefficient: + ++---------------------------------------------+ +| Chelix = 0.84, Cloop = 0.84 | ++---------------------------------------------+ +.......................................................................... + +Average length of predicted secondary structure segments: + +| +------------+----------+ +| | predicted | observed | ++-----------+------------+----------+ +| Lhelix = | 24.6 | 22.2 | ++-----------+------------+----------+ +.......................................................................... + +The accuracy matrix in detail: + ++---------------------------------+ +| number of residues with H, L | ++---------+------+-------+--------+ +| |net H | net L |sum obs | ++---------+------+-------+--------+ +| obs H | 5214 | 492 | 5706 | +| obs L | 1050 | 22423 | 23473 | ++---------+------+-------+--------+ +| sum Net | 6264 | 22915 | 29179 | ++---------+------+-------+--------+ + +Note: This table is to be read in the following manner: + 5214 of all residues predicted to be in a helical trans-membrane + region, were observed to be in the lipid bilayer, 1050 however + were observed either inside or outside of the protein, i.e. in + loop (or non-membrane) regions. The term "observed" refers to DSSP + assignment of secondary structure calculated from 3D coordinates + of experimentally determined structures (Dictionary of Secondary + Structure of Proteins: Kabsch & Sander (1983) Biopolymers, 22, + 2577-2637) where these were available. For all other proteins, + the assignment of trans-membrane segments has been taken from the + Swissprot data bank (Bairoch, A.; Boeckmann, B.: The SWISS-PROT + protein sequence data bank. Nucl. Acids Res. 20: 2019-2022, 1992). + +.......................................................................... + +Overlap between predicted and observed segments: + ++-----------------+---------------+----------------+ +| segment overlap | % of observed | % of predicted | +| Sov helix | 95.6% | 95.5% | +| Sov loop | 83.6% | 97.2% | ++-----------------+---------------+----------------+ +| Sov total | 86.0% | 96.8% | ++-----------------+---------------+----------------+ + + Definition of Sov in: Rost et al., JMB, 1994, 235, 13-26. + + As helical trans-membrane segments are longer than globular heli- + ces, correctly predicted segments can easily be made out. PHDhtm + misses 5 out of 258 observed segments, predicts 6 where non is + observed and 3 times the predicted helical segment overlaps two + observed regions. Thus, in total more than 95% of all segments + are correctly predicted. + +.......................................................................... + +Entropy of prediction (information measure): + ++-----------------+ +| I = 0.64 | ++-----------------+ + + (For comparison: homology modelling of globular proteins in three + states: I=0.62.) + Definition of Sov in: Rost et al., JMB, 1994, 235, 13-26. + +Position-specific reliability index +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +The network predicts two states: helical trans-membrane region and rest +using two output units. The prediction is assigned by choosing the ma- +ximal unit ("winner takes all"). However, the real numbers of the out- +put units contain additional information. +E.g. the difference between the two output units can be used to derive +a "reliability index". This index is given for each residue along with +the prediction. The index is scaled to have values between 0 (lowest +reliability), and 9 (highest). +The accuracies (Qtot) to be expected for residues with values above a +particular value of the index are given below as well as the fraction +of such residues (%res).: + ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| index| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | +| %res |100.0| 98.8| 97.3| 95.9| 94.1| 92.3| 89.9| 86.2| 75.0| 66.8| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| | | | | | | | | | | | +| Qtot | 94.7| 95.2| 95.6| 96.2| 96.7| 97.2| 97.7| 98.4| 99.4| 99.8| +| | | | | | | | | | | | ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ +| H%obs| 91.8| 92.9| 93.8| 94.4| 95.0| 95.7| 96.2| 96.8| 95.5| 78.7| +| L%obs| 95.3| 95.7| 96.1| 96.6| 97.0| 97.5| 98.1| 98.8| 99.7|100.0| +| | | | | | | | | | | | +| H%prd| 82.7| 83.8| 85.0| 86.7| 88.1| 89.7| 91.4| 93.8| 96.3| 97.1| +| L%prd| 97.9| 98.3| 98.5| 98.7| 98.8| 99.0| 99.2| 99.4| 99.7| 99.9| ++------+-----+-----+-----+-----+-----+-----+-----+-----+-----+-----+ + +The above table gives the cumulative results, e.g. 92.3% of all +residues have a reliability of at least 5. The overall two-state +accuracy for this subset is 97.2%. For this subset, e.g., 95.7% of +the observed helical trans-membrane residues are correctly predicted, +and 89.7% of all residues predicted to be in helical trans-membrane +segment are correct. + +The resulting network (PHD) prediction is: +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +________________________________________________________________________________ + + PHD: Profile fed neural network systems from HeiDelberg + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Prediction of: + secondary structure, by PHDsec + solvent accessibility, by PHDacc + and helical transmembrane regions, by PHDhtm + + Author: + Burkhard Rost + EMBL, 69012 Heidelberg, Germany + Internet: Rost@EMBL-Heidelberg.DE + + All rights reserved. + + The network systems are described in: + + PHDsec: B Rost & C Sander: JMB, 1993, 232, 584-599. + B Rost & C Sander: Proteins, 1994, 19, 55-72. + PHDacc: B Rost & C Sander: Proteins, 1994, 20, 216-226. + PHDhtm: B Rost et al.: Prot. Science, 1995, 4, 521-533. + + Some statistics + ~~~~~~~~~~~~~~~ + + Percentage of amino acids: + +--------------+--------+--------+--------+--------+--------+ + | AA: | L | A | S | G | I | + | % of AA: | 13.0 | 10.0 | 9.7 | 8.9 | 8.6 | + +--------------+--------+--------+--------+--------+--------+ + | AA: | V | R | T | F | D | + | % of AA: | 7.8 | 5.2 | 4.5 | 4.5 | 4.5 | + +--------------+--------+--------+--------+--------+--------+ + | AA: | N | Q | E | P | K | + | % of AA: | 4.1 | 3.0 | 3.0 | 2.6 | 2.6 | + +--------------+--------+--------+--------+--------+--------+ + | AA: | Y | M | W | H | C | + | % of AA: | 1.9 | 1.9 | 1.5 | 1.5 | 1.5 | + +--------------+--------+--------+--------+--------+--------+ + + Percentage of secondary structure predicted: + +--------------+--------+--------+--------+ + | SecStr: | H | E | L | + | % Predicted: | 43.9 | 16.7 | 39.4 | + +--------------+--------+--------+--------+ + + According to the following classes: + all-alpha: %H>45 and %E< 5; all-beta : %H<5 and %E>45 + alpha-beta : %H>30 and %E>20; mixed: rest, + this means that the predicted class is: mixed class + + PHD output for your protein + ~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + Tue Nov 24 17:44:57 1998 + Jury on: 10 different architectures (version 5.94_317 ). + Note: differently trained architectures, i.e., different versions can + result in different predictions. + + About the protein + ~~~~~~~~~~~~~~~~~ + + HEADER /home/phd/server/work/predict_h25873-220 + COMPND + SOURCE + AUTHOR + SEQLENGTH 269 + NCHAIN 1 chain(s) in predict_h25873-22040 data set + NALIGN 48 + (=number of aligned sequences in HSSP file) + + Abbreviations: PHDsec + ~~~~~~~~~~~~~~~~~~~~~ + + sequence: + AA : amino acid sequence + secondary structure: + HEL: H=helix, E=extended (sheet), blank=other (loop) + PHD: Profile network prediction HeiDelberg + Rel: Reliability index of prediction (0-9) + detail: + prH: 'probability' for assigning helix + prE: 'probability' for assigning strand + prL: 'probability' for assigning loop + note: the 'probabilites' are scaled to the interval 0-9, e.g., + prH=5 means, that the first output node is 0.5-0.6 + subset: + SUB: a subset of the prediction, for all residues with an expected + average accuracy > 82% (tables in header) + note: for this subset the following symbols are used: + L: is loop (for which above " " is used) + ".": means that no prediction is made for this residue, as the + reliability is: Rel < 5 + + Abbreviations: PHDacc + ~~~~~~~~~~~~~~~~~~~~~ + + SS : secondary structure + HEL: H=helix, E=extended (sheet), blank=other (loop) + solvent accessibility: + 3st: relative solvent accessibility (acc) in 3 states: + b = 0-9%, i = 9-36%, e = 36-100%. + PHD: Profile network prediction HeiDelberg + Rel: Reliability index of prediction (0-9) + O_3: observed relative acc. in 3 states: B, I, E + note: for convenience a blank is used intermediate (i). + P_3: predicted relative accessibility in 3 states + 10st:relative accessibility in 10 states: + = n corresponds to a relative acc. of n*n % + subset: + SUB: a subset of the prediction, for all residues with an expected + average correlation > 0.69 (tables in header) + note: for this subset the following symbols are used: + "I": is intermediate (for which above " " is used) + ".": means that no prediction is made for this residue, as the + reliability is: Rel < 4 + + Abbreviations: PHDhtm + ~~~~~~~~~~~~~~~~~~~~~ + + secondary structure: + HL: T=helical transmembrane region, blank=other (loop) + PHD: Profile network prediction HeiDelberg + PHDF:filtered prediction, i.e., too long transmembrane segments + are split, too short ones are deleted + Rel: Reliability index of prediction (0-9) + detail: + prH: 'probability' for assigning helical transmembrane region + prL: 'probability' for assigning loop + note: the 'probabilites' are scaled to the interval 0-9, e.g., + prH=5 means, that the first output node is 0.5-0.6 + subset: + SUB: a subset of the prediction, for all residues with an expected + average accuracy > 82% (tables in header) + note: for this subset the following symbols are used: + L: is loop (for which above " " is used) + ".": means that no prediction is made for this residue, as the + reliability is: Rel < 5 + + protein: predict length 269 + + ....,....1....,....2....,....3....,....4....,....5....,....6 + AA |MASEIKKKLFWRAVVAEFLAMTLFVFISIGSALGFNYPLERNQTLVQDNVKVSLAFGLSI| + PHD sec | HHHHHHHHHHHHHHHHHHHHHHHHHHEE HHHHHHHHHHHHH| + Rel sec |998443148899999999999998997676530312469989998623353579999999| + detail: + prH sec |001223468899999999999998888777653112210000000145566788999999| + prE sec |000011000000000000000001001111233542100000000000323211000000| + prL sec |998665420100000000000000000011112244578988998753100000000000| + subset: SUB sec |LLL.....HHHHHHHHHHHHHHHHHHHHHHH......LLLLLLLLL...H.HHHHHHHHH| + + ACCESSIBILITY + 3st: P_3 acc |eeeebee bbb bbbbbbbbbbbbbbbbbbbbbebeee eeeeeeeeebbbbbbbbbbbb| + 10st: PHD acc |997706650005000000000000000000000607775779776677000000000000| + Rel acc |735421110541467608662789996343122133420454330023453975664547| + subset: SUB acc |e.ee.....bb.bbbb.bbb.bbbbbb.b.......e..eee......bb.bbbbbbbbb| + ....,....7....,....8....,....9....,....10...,....11...,....12 + AA |ATLAQSVGHISGAHSNPAVTLGLLLSCQISILRAVMYIIAQCVGAIVASAILSGITSSLL| + PHD sec |HHHHHHHHHE HHHHEHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH | + Rel sec |999996412122653167703135552356779999999999999999999998467213| + detail: + prH sec |998986544334223477843456665567779999999999999999999998611343| + prE sec |001001123420010000145432101221110000000000000000000000000000| + prL sec |000001232245765521000000123210000000000000000000000000278555| + subset: SUB sec |HHHHHH......LL..HHH....HHH..HHHHHHHHHHHHHHHHHHHHHHHHHH.LL...| + + ACCESSIBILITY + 3st: P_3 acc |bbbbebbbebbbbbb bbbbbbbbbbbebbbbbbbbbbbbbbbbbbbbbbbbeebbeeeb| + 10st: PHD acc |000060006000000500000000000600000000000000000000000067006760| + Rel acc |456515321655013144869663400154551757478936465465467713401400| + subset: SUB acc |bbbb.b...bbb....bbbbbbb.b...bbbb.bbbbbbb.bbbbbbbbbbb..b..e..| + ....,....13...,....14...,....15...,....16...,....17...,....18 + AA |ENSLGRNDLARGVNSGQGLGIEIIGTLQLVLCVLATTDRRRRDLGGSAPLAIGLSVALGH| + PHD sec | HHH EEEEEEEEEEEEEEEEEEE E E HHHHHH| + Rel sec |359985212134223651899898866789799875436658889963211351457756| + detail: + prH sec |320002345432332111000000000000100000221120000000001113567767| + prE sec |100000000000011014899888877789789886100000000013544222221111| + prL sec |568986543466545763100000011100000112567768889975454564210111| + subset: SUB sec |.LLLLL.........LL.EEEEEEEEEEEEEEEEEE..LLLLLLLLL.....L..HHHHH| + + ACCESSIBILITY + 3st: P_3 acc |eeebbbebbbeebeebeebbbbbbbbbbbbbbbbbbbeeeeeeeebbbbbbbbbbbbbbb| + 10st: PHD acc |677000600077076077000000000000000000077767767000000000000000| + Rel acc |133100124043040233247198656399879530035414413123255869586654| + subset: SUB acc |........b.e..e.....bb.bbbbb.bbbbbb....ee.ee......bbbbbbbbbbb| + ....,....19...,....20...,....21...,....22...,....23...,....24 + AA |LLAIDYTGCGINPARSFGSAVLTRNFSNHWIFWVGPFIGSALAVLIYDFILAPRSSDFTD| + PHD sec |HEEEE E HHHEEEE EEEEEE HHHHHHHHHHHHHEEEEE | + Rel sec |321341126989622145152653534229996251699999999973147525556642| + detail: + prH sec |521100000000145432463121122000000114789999999875421111121124| + prE sec |244564431000000000015765121358997510000000000013467642110000| + prL sec |233234457889754567411012655530002364200000000010010136667765| + subset: SUB sec |........LLLLL....H.H.EE.L....EEEE.L.HHHHHHHHHHH...EE.LLLLL..| + + ACCESSIBILITY + 3st: P_3 acc |bbbbebbbbbbebb bbbbbbbbeebeebbbbbbbbbbbbbbbbbbbbbbbbeeeee ee| + 10st: PHD acc |000060000006005000000007606600000000000000000000000076777577| + Rel acc |754424240102242141047612131118967874356346635751777031345044| + subset: SUB acc |bbbb.b.b.....b..b..bbb.......bbbbbbb.bb.bbb.bbb.bbb....ee.ee| + ....,....25...,....26...,....27...,....28...,....29...,....30 + AA |RMKVWTSGQVEEYDLDADDINSRVEMKPK| + PHD sec |HHHHHH | + Rel sec |66775259975467555457776422699| + detail: + prH sec |77887520012221222221111100000| + prE sec |00000000000000000000001233200| + prL sec |11112379987678777678887655799| + subset: SUB sec |HHHHH.LLLLL.LLLLL.LLLLL...LLL| + + ACCESSIBILITY + 3st: P_3 acc |ebebbeeeeeeeeeeeeeeeeeebeeeee| + 10st: PHD acc |60700787677777677777767067789| + Rel acc |10411563134335144444514212559| + subset: SUB acc |..e..ee...e..e.eeeeee.e...eee| + + PHDhtm Helical transmembrane prediction + note: PHDacc and PHDsec are reliable for water- + soluble globular proteins, only. Thus, + please take the predictions above with + particular caution wherever transmembrane + helices are predicted by PHDhtm! + + PHDhtm +--- +--- PhdTopology REFINEMENT AND TOPOLOGY PREDICTION: SYMBOLS +--- AA : amino acid in one-letter code +--- PHD htm : HTM's predicted by the PHD neural network +--- system (T=HTM, ' '=not HTM) +--- Rel htm : Reliability index of prediction (0-9, 0 is low) +--- detail : Neural network output in detail +--- prH htm : 'Probability' for assigning a helical trans- +--- membrane region (HTM) +--- prL htm : 'Probability' for assigning a non-HTM region +--- note: 'Probabilites' are scaled to the interval +--- 0-9, e.g., prH=5 means, that the first +--- output node is 0.5-0.6 +--- subset : Subset of more reliable predictions +--- SUB htm : All residues for which the expected average +--- accuracy is > 82% (tables in header). +--- note: for this subset the following symbols are used: +--- L: is loop (for which above ' ' is used) +--- '.': means that no prediction is made for this, +--- residue as the reliability is: Rel < 5 +--- other : predictions derived based on PHDhtm +--- PHDFhtm : filtered prediction, i.e., too long HTM's are +--- split, too short ones are deleted +--- PHDRhtm : refinement of neural network output +--- PHDThtm : topology prediction based on refined model +--- symbols used: +--- i: intra-cytoplasmic +--- T: transmembrane region +--- o: extra-cytoplasmic +--- +--- PhdTopology REFINEMENT AND TOPOLOGY PREDICTION + ....,....1....,....2....,....3....,....4....,....5....,....6 + AA |MASEIKKKLFWRAVVAEFLAMTLFVFISIGSALGFNYPLERNQTLVQDNVKVSLAFGLSI| + PHD htm | TTTTTTTTTTTTTTTTTTT TTTTTTTTTTTT| + detail: | | + prH htm |000000000001136788999999999988875321110000000123678889999988| + prL htm |999999999998863211000000000011124678889999999876321110000011| + other: | | + PHDFhtm | TTTTTTTTTTTTTTTTTTT TTTTTTTTTTT| + PHDRhtm | TTTTTTTTTTTTTTTTTT TTTTTTTTTTT| + PHDThtm |iiiiiiiiiiiiiiTTTTTTTTTTTTTTTTTToooooooooooooooooTTTTTTTTTTT| + subset: | | + SUB htm |............................................................| + ....,....7....,....8....,....9....,....10...,....11...,....12 + AA |ATLAQSVGHISGAHSNPAVTLGLLLSCQISILRAVMYIIAQCVGAIVASAILSGITSSLL| + PHD htm |TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT | + detail: | | + prH htm |888888877777666677788888888888888888888888888888888876543211| + prL htm |111111122222333322211111111111111111111111111111111123456788| + other: | | + PHDFhtm |TTTTTTTTTTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT | + PHDRhtm |TTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTT | + PHDThtm |TTTTTTTTiiiiiiiiiiiiiTTTTTTTTTTTTTTTTTTTTTTTTToooooooooooooo| + subset: | | + SUB htm |............................................................| + ....,....13...,....14...,....15...,....16...,....17...,....18 + AA |ENSLGRNDLARGVNSGQGLGIEIIGTLQLVLCVLATTDRRRRDLGGSAPLAIGLSVALGH| + PHD htm | TTTTTTTTTTTTTTTTTTT TTTTTTTTTTTTT| + detail: | | + prH htm |000000000001234567788888999988887643211111111235788899998888| + prL htm |999999999998765432211111000011112356788888888764211100001111| + other: | | + PHDFhtm | TTTTTTTTTTTTTTTTTTT TTTTTTTTTTTTT| + PHDRhtm | TTTTTTTTTTTTTTTTTT TTTTTTTTTTTT| + PHDThtm |ooooooooooooooooTTTTTTTTTTTTTTTTTTiiiiiiiiiiiiiiTTTTTTTTTTTT| + subset: | | + SUB htm |............................................................| + ....,....19...,....20...,....21...,....22...,....23...,....24 + AA |LLAIDYTGCGINPARSFGSAVLTRNFSNHWIFWVGPFIGSALAVLIYDFILAPRSSDFTD| + PHD htm |TTTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT | + detail: | | + prH htm |888887765443432233334566777777788888888888888888887542100000| + prL htm |111112234556567766665433222222211111111111111111112457899999| + other: | | + PHDFhtm |TTTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT | + PHDRhtm |TTTTTT TTTTTTTTTTTTTTTTTTT | + PHDThtm |TTTTTToooooooooooooooooooooooooTTTTTTTTTTTTTTTTTTTiiiiiiiiii| + subset: | | + SUB htm |............................................................| + ....,....25...,....26...,....27...,....28...,....29...,....30 + AA |RMKVWTSGQVEEYDLDADDINSRVEMKPK| + PHD htm | | + detail: | | + prH htm |00000000000000000000000000000| + prL htm |99999999999999999999999999999| + other: | | + PHDFhtm | | + PHDRhtm | | + PHDThtm |iiiiiiiiiiiiiiiiiiiiiiiiiiiii| + subset: | | + SUB htm |.............................| +--- +--- PhdTopology REFINEMENT AND TOPOLOGY PREDICTION END +--- + +________________________________________________________________________________ + +________________________________________________________________________________ + +----------------------------------------------------------------------------- +--- PredictProtein: NEWS from January, 1997 --- +--- --- +--- Dear user, --- +--- --- +--- as of January 1, 1997, EMBL has effectively decided to not --- +--- support the PredictProtein service by personal resources. I do --- +--- maintain the program, so to speak, in my private time. However, --- +--- my contract obliges me to do science, instead. Unfortunately, --- +--- the computer environment at EMBL is at the same time starting --- +--- to become increasingly unstable. Consequence of these two re- --- +--- cent developments is that the PredictProtein service is not as --- +--- stable as it was. --- +--- --- +--- I apologise for the problems this may cause. In particular, --- +--- I apologise for my inability to reply to the 20-30 daily, per- --- +--- sonal mails, and suggest to re-submit requests after 24 hours! --- +--- --- +--- Hoping that I shall find a more convenient solution for the --- +--- future of the PredictProtein I remain with my best regards, --- +--- --- +--- Burkhard Rost --- +----------------------------------------------------------------------------- +--- PredictProtein: NEWS from April, 1998 --- +--- --- +-------------------------------- --- +--- MOVING PredictProtein --- +--- There appears to be light on the horizon! PP will may be having --- +--- many hickups over the next months (as I shall leave EMBL). How- --- +--- ever, the server seems to have a fair chance of survival thanks --- +--- to a major support that is being raised by Columbia University, --- +--- New York, U.S.A.). I hope that this will settle the issue for --- +--- the years to come ... --- +-------------------------------- --- +--- WARNING --- +--- After a major rewriting of most of the PP code over the last, --- +--- I am afraid that not all errors have been traced by me, yet. --- +--- Thus, please have mercy and report any bug you'll encounter! --- +--- THANKS, Burkhard Rost --- +-------------------------------- --- +--- NEW PREDICTION DEFAULTS --- +--- * Coiled-coil regions: now by default the program COILS written by --- +--- Andrei Lupas is run on your sequence. An output is returned if a --- +--- coiled-coil region has been detected. --- +--- * Functional sequence motifs: now by default the PROSITE database --- +--- written by Amos Bairoch, Philip Bucher and Kay Hofmann is scanned --- +--- for sequence motifs. An output is returned if any motif has been --- +--- detected. --- +-------------------------------- --- +--- see http://www.embl-heidelberg.de/predictprotein/ppNews.html --- +--- for a description of the following new options. --- +--- NEW INPUT OPTION --- +--- * Your input sequence(s) in FASTA-list format ("# FASTA list ") --- +--- NEW OUTPUT OPTIONS --- +--- * Return also BLASTP output ("return blast") --- +--- * Return prediction additionally in RDB format ("return phd rdb") --- +--- * Return topits hssp ("return topits hssp") --- +--- * Return topits strip ("return topits strip") --- +--- * Return topits own ("return topits own") --- +--- * Return no coils ("return no coils") --- +--- * Return no prosite ("return no prosite") --- +----------------------------------------------------------------------------- |