Difference between revisions of "Sequence-based analyses Gaucher Disease"
(→P04062) |
(→P47863) |
||
(42 intermediate revisions by the same user not shown) | |||
Line 353: | Line 353: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
− | <figtable id="tab:P35462_trans_png"> |
||
+ | File:3pbl_opm.png|Figure 6: Dopamine D3 receptor(P35462,3PBL) showing in OPM. |
||
− | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0" align="left" width="380px" |
||
− | + | File:3pbl_pdbtm.png|Figure 7: Dopamine D3 receptor(P35462,3PBL) showing in PDBTM. |
|
+ | </gallery> |
||
− | |[[File:3pbl_pdbtm.png|thumb|200px|left|Figure 7: Dopamine D3 receptor(P35462,3PBL) showing in PDBTM.]] |
||
− | |} |
||
− | </figtable> |
||
− | <br style="clear:both;"> |
||
− | '''Table 9''': Visualization of transmenbrane helix prediction results for P35462 in OPM and PDBTM. |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
Line 400: | Line 396: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | <gallery widths=200px heights=150px perrow=4 caption=""> |
||
− | <figtable id="tab:Q9YDF8_trans_png"> |
||
+ | File:1orq_opm.png|Figure 8: Potassium channel KvAP(Q9YDF8,1ORQ) showing in OPM. |
||
− | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0" align="left" width="380px" |
||
− | + | File:1ors_opm.png|Figure 9: Potassium channel KvAP(Q9YDF8,1ORS) showing in OPM. |
|
− | + | File:1orq_pdbtm.png|Figure 10: Potassium channel KvAP(Q9YDF8,1ORQ) showing in PDBTM. |
|
− | + | File:1ors_pdbtm.png|Figure 11: Potassium channel KvAP(Q9YDF8,1ORS) showing in PDBTM. |
|
+ | </gallery> |
||
− | |[[File:1ors_pdbtm.png|thumb|200px|left|Figure 11: Potassium channel KvAP(Q9YDF8,1ORS) showing in PDBTM.]] |
||
− | |} |
||
− | </figtable> |
||
− | <br style="clear:both;"> |
||
− | '''Table 11''': Visualization of transmenbrane helix prediction results for Q9YDF8 in OPM and PDBTM. |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
Line 446: | Line 438: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
− | <figtable id="tab:P47863_trans_png"> |
||
+ | File:2d57_opm.png|Figure 12: Aquaporin-4(P47863,2D57) showing in OPM. |
||
− | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0" align="left" width="380px" |
||
− | + | File:2d57_pdbtm.png|Figure 13: Aquaporin-4(P47863,2D57) showing in PDBTM. |
|
+ | </gallery> |
||
− | |[[File:2d57_pdbtm.png|thumb|200px|left|Aquaporin-4(P47863,2D57) showing in PDBTM.]] |
||
− | |} |
||
− | </figtable> |
||
− | |||
− | <br style="clear:both;"> |
||
− | '''Table 13''': Visualization of transmenbrane helix prediction results for P47863 in OPM and PDBTM. |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
Line 461: | Line 448: | ||
The online server of [http://www.cbs.dtu.dk/services/SignalP/ SignalP] with version 4.0 was used to predict signal peptides for our protein [http://www.uniprot.org/uniprot/P04062 P04062] and other three proteins [http://www.uniprot.org/uniprot/P02768 P02768], [http://www.uniprot.org/uniprot/P47863 P47863] and [http://www.uniprot.org/uniprot/P11279 P11279](<xr id="tab:ss_sequences_signa"/>): We also used Polyphobius and the information from uniprot and [http://www.signalpeptide.de/ Signal Peptide Website] to compare and validate the prediction results. |
The online server of [http://www.cbs.dtu.dk/services/SignalP/ SignalP] with version 4.0 was used to predict signal peptides for our protein [http://www.uniprot.org/uniprot/P04062 P04062] and other three proteins [http://www.uniprot.org/uniprot/P02768 P02768], [http://www.uniprot.org/uniprot/P47863 P47863] and [http://www.uniprot.org/uniprot/P11279 P11279](<xr id="tab:ss_sequences_signa"/>): We also used Polyphobius and the information from uniprot and [http://www.signalpeptide.de/ Signal Peptide Website] to compare and validate the prediction results. |
||
+ | |||
+ | SignalP used three different scores, C, S and Y. Two additional scores, S-mean and the D-score, are reported in the output. The algorithm of SignalP employed two different neural networks, one for predicting the actual signal peptide and one for predicting the position of the signal peptidase cleavage site. <br/> |
||
+ | * The S-score is reported for every single amino acid position. If S-score is high, it suggests that the corresponding amino acid is part of a signal peptide. On the other hand, a low score indicate that the amino acid is part of a mature protein. <br/> |
||
+ | * The C-score is so called ''cleavage site'' score which should only be significantly high at the cleavage site. <br/> |
||
+ | * The Y-max is calculated by combining the C-score with the S-score and produces a better cleavage site prediction than the raw C-score alone.<br/> |
||
+ | * The S-mean is the average of the S-score for the length of the predicted signal peptide. |
||
+ | * The D-score is calculated as a weighted average of the S-mean and the Y-max scores.<br/> |
||
+ | |||
+ | Observing the above scores will help us to locate the signal peptides. |
||
+ | |||
<figtable id="tab:ss_sequences_signa"> |
<figtable id="tab:ss_sequences_signa"> |
||
Line 479: | Line 476: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
− | '''Table |
+ | '''Table 11''': The proteins used to predict signal peptides. |
<br style="clear:both;"> |
<br style="clear:both;"> |
||
=== P04062 === |
=== P04062 === |
||
+ | |||
+ | By using the standard option, SignalP did not return a positive result for our protein Glucosylceramidase [http://www.uniprot.org/uniprot/P04062 P04062], even the s score and c score were quite high. After lowering the Cutoff value to 0.35, SignalP returned a finding of signal peptide which was validated by other tools and database(<xr id="tab:ss_sinal_p04062"/>). |
||
+ | # Measure Position Value Cutoff signal peptide? |
||
+ | max. C 40 0.305 |
||
+ | max. Y 40 0.396 |
||
+ | max. S 36 0.684 |
||
+ | mean S 1-39 0.323 |
||
+ | D 1-39 0.367 0.350 YES |
||
+ | Name=sp_P04062_GLCM_HUMAN SP='YES' Cleavage site between pos. 39 and 40: ASG-AR D=0.367 D-cutoff=0.350 Networks=SignalP-TM |
||
<figtable id="tab:ss_sinal_p04062"> |
<figtable id="tab:ss_sinal_p04062"> |
||
Line 497: | Line 503: | ||
|- |
|- |
||
| [http://www.uniprot.org/uniprot/P04062 uniprot] || 1-39 |
| [http://www.uniprot.org/uniprot/P04062 uniprot] || 1-39 |
||
+ | |- |
||
+ | | PolyPhobius || 1-39 |
||
|- |
|- |
||
| style="border-style: solid; border-width: 0 0 1px 0" | |
| style="border-style: solid; border-width: 0 0 1px 0" | |
||
Line 502: | Line 510: | ||
|- |
|- |
||
| style="border-style: solid; border-width: 0 0 0px 0" | signal peptide sequence |
| style="border-style: solid; border-width: 0 0 0px 0" | signal peptide sequence |
||
− | | style="border-style: solid; border-width: 0 0 0px 0" | |
+ | | style="border-style: solid; border-width: 0 0 0px 0" | MEFSSPSREECPKPLSRVSIMAGSLTGLLLLQAVSWASG |
|} |
|} |
||
</figtable> |
</figtable> |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
− | '''Table |
+ | '''Table 12''': Signal peptides prediction of P04062. |
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | |||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
+ | File:sp_P04062_GLCM_HUMAN.png|Figure 14: signal peptides prediction of Glucosylceramidase (P04062) by using SignalP. |
||
+ | File:p04062.png |Figure 15: signal peptides of Glucosylceramidase (P04062) showing in [http://www.signalpeptide.de/index.php?sess=&m=myprotein&s=details&id=26600&listname= Signal Peptide Website]. |
||
+ | </gallery> |
||
=== P02768 === |
=== P02768 === |
||
+ | For the protein Serum albumin ([http://www.uniprot.org/uniprot/P02768 P02768]), SignalP returned positive prediction which was validated by using other tools. |
||
+ | # Measure Position Value Cutoff signal peptide? |
||
+ | max. C 19 0.710 |
||
+ | max. Y 19 0.798 |
||
+ | max. S 2 0.929 |
||
+ | mean S 1-18 0.890 |
||
+ | D 1-18 0.848 0.450 YES |
||
+ | Name=sp_P02768_ALBU_HUMAN SP='YES' Cleavage site between pos. 18 and 19: AYS-RG D=0.848 D-cutoff=0.450 Networks=SignalP-noTM |
||
<figtable id="tab:ss_sinal_p02768"> |
<figtable id="tab:ss_sinal_p02768"> |
||
{| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 1px 0" align="left" width="480px" |
{| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 1px 0" align="left" width="480px" |
||
Line 524: | Line 545: | ||
|- |
|- |
||
| [http://www.uniprot.org/uniprot/p02768 uniprot] || 1-18 |
| [http://www.uniprot.org/uniprot/p02768 uniprot] || 1-18 |
||
+ | |- |
||
+ | | PolyPhobius || 1-18 |
||
|- |
|- |
||
| style="border-style: solid; border-width: 0 0 1px 0" | |
| style="border-style: solid; border-width: 0 0 1px 0" | |
||
Line 534: | Line 557: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
− | '''Table |
+ | '''Table 13''': Signal peptides prediction of P02768. |
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | |||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
+ | File:sp_P02768_ALBU_HUMAN.png|Figure 16: signal peptides prediction of Serum albumin (P02768) by using SignalP. |
||
+ | File:p02768.png |Figure 17: signal peptides of Serum albumin (P02768) showing in [http://www.signalpeptide.de/index.php?sess=&m=myprotein&s=details&id=22229&listname= Signal Peptide Website]. |
||
+ | </gallery> |
||
=== P47863 === |
=== P47863 === |
||
+ | |||
+ | For the protein Aquaporin-4 [http://www.uniprot.org/uniprot/P47863 P47863], as expected from its alternative name "Mercurial-insensitive water channel", SignalP returned very low scores for which it was not meaningful to lower the cutoff value. All the other tools reported no finding either (<xr id="tab:ss_sinal_P47863"/>). |
||
+ | |||
+ | # Measure Position Value Cutoff signal peptide? |
||
+ | max. C 41 0.209 |
||
+ | max. Y 41 0.164 |
||
+ | max. S 56 0.183 |
||
+ | mean S 1-40 0.139 |
||
+ | D 1-40 0.154 0.500 NO |
||
+ | Name=sp_P47863_AQP4_RAT SP='NO' D=0.154 D-cutoff=0.500 Networks=SignalP-TM |
||
+ | |||
+ | <figtable id="tab:ss_sinal_P47863"> |
||
+ | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 1px 0" align="left" width="480px" |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | P47863 |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | position |
||
+ | |- |
||
+ | |- |
||
+ | | [http://www.cbs.dtu.dk/services/SignalP/ SignalP] || - |
||
+ | |- |
||
+ | | [http://www.signalpeptide.de/index.php?sess=&m=myprotein&s=details&id=22229&listname= Signal Peptide Website] || - |
||
+ | |- |
||
+ | | [http://www.uniprot.org/uniprot/p47863 uniprot] || - |
||
+ | |- |
||
+ | | PolyPhobius || - |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 0px 0" | signal peptide sequence |
||
+ | | style="border-style: solid; border-width: 0 0 0px 0" | - |
||
+ | |} |
||
+ | </figtable> |
||
+ | |||
+ | <br style="clear:both;"> |
||
+ | '''Table 14''': Signal peptides prediction of P47863. |
||
+ | <br style="clear:both;"> |
||
+ | |||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
+ | File:sp_P47863_AQP4_RAT.png|Figure 18: signal peptides prediction of Aquaporin-4 (P47863) by using SignalP. |
||
+ | </gallery> |
||
=== P11279 === |
=== P11279 === |
||
+ | For the protein Lysosome-associated membrane glycoprotein 1 ([http://www.uniprot.org/uniprot/P11279 P11279]), as expected from its name, a long signal peptide was found with very high score: |
||
− | <figtable id="tab:signa_png"> |
||
+ | |||
− | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0" align="left" width="380px" |
||
+ | # Measure Position Value Cutoff signal peptide? |
||
− | |[[File:sp_P04062_GLCM_HUMAN.png|thumb|200px|left|signal peptides prediction of Glucosylceramidase (P04062) by using SignalP.]] |
||
+ | max. C 29 0.910 |
||
− | |[[File:sp_P02768_ALBU_HUMAN.png|thumb|200px|left|signal peptides prediction of Serum albumin (P02768) by using SignalP.]] |
||
+ | max. Y 29 0.939 |
||
− | |[[File:sp_P47863_AQP4_RAT.png|thumb|200px|left|signal peptides prediction of Aquaporin-4 (P47863) by using SignalP.]] |
||
+ | max. S 21 0.996 |
||
− | |[[File:sp_P11279_LAMP1_HUMAN.png|thumb|200px|left|signal peptides prediction of Lysosome-associated membrane glycoprotein 1 (P11279) by using SignalP.]] |
||
+ | mean S 1-28 0.962 |
||
+ | D 1-28 0.952 0.450 YES |
||
+ | Name=sp_P11279_LAMP1_HUMAN SP='YES' Cleavage site between pos. 28 and 29: ASA-AM D=0.952 D-cutoff=0.450 Networks=SignalP-noTM |
||
+ | |||
+ | <figtable id="tab:ss_sinal_P11279"> |
||
+ | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 1px 0" align="left" width="480px" |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | P11279 |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | position |
||
+ | |- |
||
+ | |- |
||
+ | | [http://www.cbs.dtu.dk/services/SignalP/ SignalP] || 1-28 |
||
+ | |- |
||
+ | | [http://www.signalpeptide.de/index.php?sess=&m=myprotein&s=details&id=17551&listname= Signal Peptide Website] || 1-28 |
||
+ | |- |
||
+ | | [http://www.uniprot.org/uniprot/p11279 uniprot] || 1-28 |
||
+ | |- |
||
+ | | PolyPhbius || 1-28 |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | |
||
+ | | style="border-style: solid; border-width: 0 0 1px 0" | |
||
+ | |- |
||
+ | | style="border-style: solid; border-width: 0 0 0px 0" | signal peptide sequence |
||
+ | | style="border-style: solid; border-width: 0 0 0px 0" | MAAPGSARRPLLLLLLLLLLGLMHCASA |
||
|} |
|} |
||
</figtable> |
</figtable> |
||
+ | |||
+ | <br style="clear:both;"> |
||
+ | '''Table 15''': Signal peptides prediction of P11279. |
||
+ | <br style="clear:both;"> |
||
+ | |||
+ | <gallery widths=250px heights=150px perrow=3 caption=""> |
||
+ | File:sp_P11279_LAMP1_HUMAN.png|Figure 19: signal peptides prediction of Lysosome-associated membrane glycoprotein 1 (P11279) by using SignalP. |
||
+ | File:p11279.png |Figure 20: signal peptides of Lysosome-associated membrane glycoprotein 1 (P11279) showing in [http://www.signalpeptide.de/index.php?sess=&m=myprotein&s=details&id=17551&listname= Signal Peptide Website]. |
||
+ | </gallery> |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
Line 557: | Line 659: | ||
[http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar GOPET] was used to predict GO terms for our protein Glucosylceramidase ([http://www.uniprot.org/uniprot/P04062 P04062]). The results was shown in <xr id="tab:ss_gopet"/> where only three GO terms were returned. Resetting the search options like lower the Confidence threshold did not change the prediction results. <br/> |
[http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar GOPET] was used to predict GO terms for our protein Glucosylceramidase ([http://www.uniprot.org/uniprot/P04062 P04062]). The results was shown in <xr id="tab:ss_gopet"/> where only three GO terms were returned. Resetting the search options like lower the Confidence threshold did not change the prediction results. <br/> |
||
− | All three GO terms belong to "Molecular Function Ontology" and are all involved in "hydrolase" or "glycosyl bonds". Since we have already known the function of our protein Glucosylceramidase ([http://www.uniprot.org/uniprot/P04062 P04062]) is to cleaves the glucosidic bonds of glucosylceramide and synthetic beta-glucosides. All these three predicted GO terms are correct. However, if compared with the GO terms list from uniprot or the prediction from ProtFun which is shown below, GOPET returned quite fewer results. It implies that GOPET prefers to report only the GO terms with high confidence.<br/> |
+ | All three GO terms belong to "Molecular Function Ontology" and are all involved in "hydrolase" or "glycosyl bonds". Since we have already known the function of our protein Glucosylceramidase ([http://www.uniprot.org/uniprot/P04062 P04062]) is to cleaves the glucosidic bonds of glucosylceramide and synthetic beta-glucosides. All these three predicted GO terms are correct. However, if compared with the GO terms list from uniprot([http://www.ebi.ac.uk/QuickGO/GProtein?ac=P04062 QuickGo]) or the prediction from ProtFun which is shown below, GOPET returned quite fewer results. It implies that GOPET prefers to report only the GO terms with high confidence.<br/> |
<figtable id="tab:ss_gopet"> |
<figtable id="tab:ss_gopet"> |
||
Line 565: | Line 667: | ||
! style="border-style: solid; border-width: 0 0 2px 0" align="left"| Confidence |
! style="border-style: solid; border-width: 0 0 2px 0" align="left"| Confidence |
||
! style="border-style: solid; border-width: 0 0 2px 0" align="left"| GO term |
! style="border-style: solid; border-width: 0 0 2px 0" align="left"| GO term |
||
− | ! style="border-style: solid; border-width: 0 0 2px 0" align="left"| |
+ | ! style="border-style: solid; border-width: 0 0 2px 0" align="left"| validate in [http://www.ebi.ac.uk/QuickGO/GProtein?ac=P04062 QuickGo] |
|- |
|- |
||
| GO:0016787 || F || 98% || hydrolase activity || yes |
| GO:0016787 || F || 98% || hydrolase activity || yes |
||
Line 703: | Line 805: | ||
<gallery widths=250px heights=150px perrow=3 caption=""> |
<gallery widths=250px heights=150px perrow=3 caption=""> |
||
− | File:pfam_p04062.png |One domain Glycoside hydrolase family 30 (green) was found by Pfam , the signal peptide (yellow) at the top (yellow). |
+ | File:pfam_p04062.png |Figure 21: One domain Glycoside hydrolase family 30 (green) was found by Pfam , the signal peptide (yellow) at the top (yellow). |
− | File:PDB_1nof_EBI.jpg |Structure of a xylanase from glycoside hydrolase family 5. |
+ | File:PDB_1nof_EBI.jpg |Figure 22: Structure of a xylanase from glycoside hydrolase family 5. |
</gallery> |
</gallery> |
||
Latest revision as of 02:45, 22 May 2012
Contents
Secondary structure
Knowing the secondary structure of a protein can shed light on its function since structure implies function. If the structure of a protein is known, secondary structure elements (helix, sheet, coiled) can be assigned to its residues depending on their affinity to form hydrogen bonds. DSSP <ref name="DSSP">Kabsch W, Sander C (1983). "Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features". Biopolymers</ref> is the most common method to perform such secondary structure assignments. If the structure of a protein is unknown, secondary structure elements be be predicted by tools like PSIPRED <ref name="PSIPRED">Liam J. McGuffin, Kevin Bryson, and David T. Jones (2000). "The PSIPRED protein structure prediction server". Bioinformatics</ref> or Reprof<ref name="Reprof">B Rost, C Sander (1993). "Prediction of protein secondary structure at better than 70% accuracy". J. Mol. Bio.</ref>. The aim of this task was to analyse the secondary structure of different proteins and the compare the secondary structure predictions of PSIPRED and Reprof with the DSSP secondary structure assignments. Following sequences were taken into account: <figtable id="tab:ss_sequences">
NAME | UniProtKB | PDB |
---|---|---|
Glucosylceramidase | P04062 | 1OGS |
Ribonuclease inhibitor | P10775 | 1DFJ |
Divalent-cation tolerance protein CutA | Q9X0E6 | 1KR4 |
Serine/threonine-protein phosphatase | Q08209 | 1AUI |
</figtable> Information about program calls and implementation details can be found in our protocol.
Predictions
For being able to better compare the different output formats, we mapped the secondary structure definitions of all three methods onto the three letters H (helix), E (sheet), and C (coiled) according to table <xr id="tab:ss_mapping"/>. Regions of the UniProt sequences which were not present in the PDB file as well as regions where no DSSP assignment was possible were ignored. <figtable id="tab:ss_mapping">
Method | H | E | C |
---|---|---|---|
DSSP | H,G,I | E,B | T,S,' ' |
PSIPRED | H | E | C |
Reprof | H | E | L |
</figtable>
P04062
Glycosylcermidase (P04062) is located the the membrane of lysosomes. It exhibits two domains which belong to the (1) glycosyl hydrolase domain fold and (2) the TIM beta/alpha-barrel fold. Both domains have hydrophobic beta sheets which anchor the protein in the membrane. <xr id="fig:ss_P04062"/> depicts the secondary structure elements of the corresponding crystal structure which coincide with the DSSP assignments. The following section shows the secondary structure annotations of the different methods: The PSIPRED predictions better coincide with the DSSP assigments than the Reprof predictions do. Reprof predicts sheets instead of helices in several regions. The residues of the beta-barell sheets (the tube in the middle of <xr id="fig:ss_P04062"/>) are marked by asterisks.
</figure>
40 ARPCIPKSFGYSSVVCVCNATYCDSFDPPTFPALGTFSRYESTRSGRRMELSMGPIQANHTGTGLLLTLQPEQKFQKVKG Reprof CCCCCCCCCCCEEEEEEECCEECCCCCCCCCCCCCEEEEEEECCCCCEEEEECCCEECCCCCCEEEEEECCCCEEEEEEC PSIPRED CCCCCCCCCCCCCEEEEECCHHCCCCCCCCCCCCCEEEEEEECCCCCCHHCCCCCCCCCCCCCCCEEEECCCCCCEEEEE DSSP CECCCEEECCCCCEEEEEECCCCCECCCCCCCCCCEEEEEEEECCCCCCEEEEEECECCCCCCCEEEEEEEEEEEEECCE 120 FGGAMTDAAALNILALSPPAQNLLLKSYFSEEGIGYNIIRVPMASCDFSIRTYTYADTPDDFQLHNFSLPEEDTKLKIPL Reprof CCCCCCHHHHHEEEEECCCCCCEEEEEEECCCCCEEEEEEECCCCCCEEEEEEEECCCCCCCEEEEECCCCCCCCCEEEE PSIPRED EEEEHHHHHHHHHHHCCHHHHHHHHHHHCCCCCCEEEEEEEEECCCCCCCCCCCCCCCCCCCCCCCCCCCCHHHHHHHHH DSSP EEEECCHHHHHHHCCCCHHHHHHHHHHHHCCCCCCCCEEEEEECCCCCCCCCCCCCCCCCCCCCCCCCCCHHHHCCHHHH 200 IHRALQLAQRPVSLLASPWTSPTWLKTNGAVNGKGSLKGQPGDIYHQTWARYFVKFLDAYAEHKLQFWAVTAENEPSAGL Reprof EEHHHHHCCCCCEEEECCCCCCCEEEECCCECCCEECCCCCCCCCCHHHHHHHHHHHHHCCCCEEEEEEEEECCCCCCCE PSIPRED HHHHHHHHCCCEEEEEEECCCCHHEEECCCCCCCCCCCCCCCHHHHHHHHHHHHHHHHHHHHCCCEEEEEEECCCCCCCC DSSP HHHHHHHCCCCCEEEEEECCCCHHHECCCCCCCCCEECCCCCCHHHHHHHHHHHHHHHHHHHCCCCCCEEECCCCCCHHH ****** *** 280 LSGYPFQCLGFTPEHQRDFIARDLGPTLANSTHHNVRLLMLDDQRLLLPHWAKVVLTDPEAAKYVHGIAVHWYLDFLAPA Reprof ECCCCEEEECCCCCCCCCEEEECCCCCCCCCCCCEEEEEEECCCCEECCCEEEEEECCCCCCEEEEEEEEEEEEEECCCC PSIPRED CCCCCCCCCCCCHHHHHHHHHHHHHHHHHHCCCCCEEEEEECCCCCCHHHHHHHHHCCHHHHHHCCEEEEECCCCCCCCH DSSP CCCCCCCCCECCHHHHHHHHHHCHHHHHHCCCCCCCEEEEEEEEHHHCCHHHHHHHCCHHHHCCCCEEEEEEECCCCCCH ******** ******* 360 KATLGETHRLFPNTMLFASEACVGSKFWEQSVRLGSWDRGMQYSHSIITNLLYHVVGWTDWNLALNPEGGPNWVRNFVDS Reprof CCECCCCCECCCCCEEEECCCCCCCEEEEEEEEECCCCCCCEEECEHHHHEEEEEEECCCCEEEECCCCCCEEEEECCCC PSIPRED HHHHHHHHHHCCCCEEEEECCCCCCCCCCCCCCCCCHHHHHHHHHHHHHHHHHHHHEEEEEEEEECCCCCCCCCCCCCCC DSSP HHHHHHHHHHCCCCEEEEEEEECCCCCCCCCCCCCCHHHHHHHHHHHHHHHHCCEEEEEEEECCECCCCCCCCCCCCCCC ******** ******** 440 PIIVDITKDTFYKQPMFYHLGHFSKFIPEGSQRVGLVASQKNDLDAVALMHPDGSAVVVVLNRSSKDVPLTIKDPAVGFL Reprof CEEEEECCCCCCCCCEEEECCCEEEECCCCCEEEEEEEECCCCCEEEEEECCCCCEEEEEEECCCCCCEEEECCCCEEEE PSIPRED CEEEECCCCEEEECHHHHHHHHHHHHCCCCCEEEEEECCCCCCEEEEEEECCCCCEEEEEEECCCCCEEEEEEECCCEEE DSSP CEEEEHHHCEEEECHHHHHHHHHHCCCCCCCEEEEEEECCCCCEEEEEEECCCCCEEEEEEECCCCCEEEEEEECCCEEE 520 ETISPGYSIHTYLWRRQ Reprof EEECCCCEEEEEEEECC PSIPRED EEECCCCEEEEEEEECC DSSP EEEECCCEEEEEEECCC |
<figure id="fig:ss_P04062"> |
P10775
<xr id="fig:ss_P10775"/> depicts the crystal structure 1DFJ which refers to P10775. It has two domains: d1dfji_ is a repeat domain consisting of altering alpha-helices and parallel beta-sheets. d1dfje_ contains long curved antiparallel beta-sheets and three alpha-helices. The alternating HHH and EEE regions in the following secondary structure annotations suit well with repetitive structure shown in <xr id="fig:ss_P10775"/>. Again, the PSIPRED predictions better match the DSSP assignments than Reprof.
</figure>
1 MNLDIHCEQLSDARWTELLPLLQQYEVVRLDDCGLTEEHCKDIGSALRANPSLTELCLRTNELGDAGVHLVLQGLQSPTC Reprof CCCCECHHCCCCCHHHHHHHHHHHCCEEEECCCCCCHHHHHHHHHHHHCCCCHHHHHHHHCCCCCCCHEEEHCCCCCCCC PSIPRED CEEECCCCCCCHHHHHHHHHHHCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCCEEECCCCCCCHHHHHHHHHHHHHCCC DSSP CECCCECCCCCHHHHHHHHHHHCCCCEEECECCCCCHHHHHHHHHHHHCCCCCCEEECCCCCCHHHHHHHHHHHHCCCCC 81 KIQKLSLQNCSLTEAGCGVLPSTLRSLPTLRELHLSDNPLGDAGLRLLCEGLLDPQCHLEKLQLEYCRLTAASCEPLASV Reprof EEEEECCCCCCCCHCCCCCCHHHHCHCHHHHHHCCCCCCCCHHHHHHHHHCCCCCHCCHHHHHHHHHHCCCCCCHHHHHH PSIPRED CCCEEECCCCCCCHHHHHHHHHHHHCCCCCCEEECCCCCCCHHHHHHHHHHHHCCCCCCCEEECCCCCCCHHHHHHHHHH DSSP CCCEEECCCCCCCCCHHHHHHHHHHHCCCCCEEECCCCCCHHHHHHHHHHHHHCCCCCCCEEECCCCCCCHHHHHHHHHH 161 LRATRALKELTVSNNDIGEAGARVLGQGLADSACQLETLRLENCGLTPANCKDLCGIVASQASLRELDLGSNGLGDAGIA Reprof HHHHHHHHHHCCCCCCHHHHHHHHHCCCCCCHHHHHHHHHHCCCCCCCCCHHHHHHHHHCHCCHHHCCCCCCCCCHHHHH PSIPRED HHHCCCCCEEECCCCCCCHHHHHHHHHHHHCCCCCCCEEECCCCCCCHHHHHHHHHHHHCCCCCCEEECCCCCCCHHHHH DSSP HHHCCCCCEEECCCCCCHHHHHHHHHHHHHCCCCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCCEEECCCCCCHHHHHH 241 ELCPGLLSPASRLKTLWLWECDITASGCRDLCRVLQAKETLKELSLAGNKLGDEGARLLCESLLQPGCQLESLWVKSCSL Reprof HHCCCCCCCHHHHCHHEEEHCCCCCHHHHHHHHHHHHHHHHHHHHHHCCCCCCHHHHHHHHHHCCCCCCHHHHHHHHCHH PSIPRED HHHHHHHHCCCCCCEEECCCCCCCHHHHHHHHHHHHCCCCCCEEECCCCCCCHHHHHHHHHHHHCCCCCCCEEECCCCCC DSSP HHHHHHCCCCCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCCEEECCCCCCHHHHHHHHHHHHHCCCCCCCEEECCCCCC 321 TAACCQHVSLMLTQNKHLLELQLSSNKLGDSGIQELCQALSQPGTTLRVLCLGDCEVTNSGCSSLASLLLANRSLRELDL Reprof HHHHHHHHHHHHHCCHHHHHHHCCCCCCCCHHHHHHHHHHCCCCCEEEEEEECCCCCCCCCHHHHHHHHHHHCCHHHHCC PSIPRED CHHHHHHHHHHHHHCCCCCEEECCCCCCCHHHHHHHHHHHCCCCCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCCEEEC DSSP EHHHHHHHHHHHHHCCCCCEEECCCCECHHHHHHHHHHHHHCCCCCCCEEECCCCCCEHHHHHHHHHHHHHCCCCCEEEC 401 SNNCVGDPGVLQLLGSLEQPGCALEQLVLYDTYWTEEVEDRLQALEGSKPGLRVIS Reprof CCCCCCCHHHHHHHCCCCCCCCHHHHHHHCCCCCCHHHHHHHHHHHCCCCCCCECC PSIPRED CCCCCCHHHHHHHHHHHHCCCCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCEECC DSSP CCCECCHHHHHHHHHHHCCCCCCCCEEECCCCCCCHHHHHHHHHHHHHCCCCEEEC |
<figure id="fig:ss_P10775"> |
Q9X0E6
<xr id="fig:ss_Q9X0E6"/> depicts the d1kr4a_ domain of 1KR4 which is made of three alpha-helices interrupted by beta-sheets. Reprof predicts too long helices.
</figure>
2 ILVYSTFPNEEKALEIGRKLLEKRLIACFNAFEIRSGYWWKGEIVQDKEWAAIFKTTEEKEKELYEELRKLHPYETPAIF Reprof EEEEECCCCHHHHHHHHHHHHHHHHHHHHCHCHHHCCCEEECEEECCHHHHHHHCCCHHHHHHHHHHHHHCCCCCCCHHE PSIPRED EEEEEECCCHHHHHHHHHHHHHCCCEEEEEEEEEEEEEEECCCEEEEEEEEEEECCCHHHHHHHHHHHHHHCCCCCCEEE DSSP EEEEEEECCHHHHHHHHHHHHHCCCCCEEEEEEEEEEEEECCEEEEEEEEEEEEEEEHHHHHHHHHHHHHHCCCCCCCEE 82 TLKVENVLTEYMNWLRESVL Reprof HHHHHHHHHHHHHHHHHHCC PSIPRED EEECCCCCHHHHHHHHHHCC DSSP EECCCCEEHHHHHHHHHHCC |
<figure id="fig:ss_Q9X0E6"> |
Q08209
Q08209 contains the domain d1auib_ and d1auia_ which are mainly assembled of alpha-helices (<xr id="fig:ss_Q08209"/>). PSIPRED predicts these alpha-helices considerably better than Reprof which suggests beta-sheets in some regions.
</figure>
14 TDRVVKAVPFPPSHRLTAKEVFDNDGKPRVDILKAHLMKEGRLEESVALRIITEGASILRQEKNLLDIDAPVTVCGDIHG Reprof CCCEEEEECCCCCCCEEEEEEECCCCCCEEEEEHHHECCCCCCCEEEEEEEEECCCCEEECCCCCCCCCCCEEEEECCCC PSIPRED CCCCCCCCCCCCCCCCCHHHCCCCCCCCCHHHHHHHHHCCCCCCHHHHHHHHHHHHHHHHHCCCEEEECCCEEEECCCCC DSSP CCCCCCCCCCCCCCCECHHHHECCCCCECHHHHHHHHHCCCCECHHHHHHHHHHHHHHHHCCCCEEEECCCEEEECCCCC 94 QFFDLMKLFEVGGSPANTRYLFLGDYVDRGYFSIECVLYLWALKILYPKTLFLLRGNHECRHLTEYFTFKQECKIKYSER Reprof HHHHHHEEEEECCCCCCCEEEEEEEECCCCEEEEEEEHHHHHHHCCCCCEEEEEECCCCCCEEEEEEEEEEEEEEEEECH PSIPRED HHHHHHHHHHHCCCCCCCCEEECCCCCCCCCCCHHHHHHHHHHHHCCCCCEEEECCCCHHHHHHCCCCHHHHHHHHCCHH DSSP CHHHHHHHHHHHCCCCCCCEEECCCCCCCCCCHHHHHHHHHHHHHHCCCCEEECCCCCCCHHHHHHCCHHHHHHHHCCHH 174 VYDACMDAFDCLPLAALMNQQFLCVHGGLSPEINTLDDIRKLDRFKEPPAYGPMCDILWSDPLEDFGNEKTQEHFTHNTV Reprof HHHHHHHHCCCCCHHHHHCCCEEEEECCCCCCCCCHHHHHHHHCCCCCCCCCCCEEEEECCCCCCCCCCCCCCECCCCCE PSIPRED HHHHHHHHHCCCHHHHHCCCCEEEECCCCCCCCCCHHHHHHCCCCCCCCCCCHHHHHHCCCCCCCCCCCCCCCCCCCCCC DSSP HHHHHHHHHCCCCCEEEECCCEEEECCCCCCCCCCHHHHHHCCCCCCCCCCCHHHHHHHCEECCCCCCCCCCCCEEECCC 254 RGCSYFYSYPAVCEFLQHNNLLSILRAHEAQDAGYRMYRKSQTTGFPSLITIFSAPNYLDVYNNKAAVLKYENNVMNIRQ Reprof CCEEEEECCCCEEEEHCCCCHHHHEHHHCCCCCCEEEEEECCCCCCCEEEEEEECCCEEEEECCCEEEEEECCCEEEEEE PSIPRED CCCCCCCCHHHHHHHHHHCCCCEEEEHHHHHHHHHHHHHCCCCCCCCCEEEEEECCCCCCCCCCEEEEEEECCCCCEEEE DSSP CCCCEEECHHHHHHHHHHCCCCEEEECCCCCCCCEEECCECCCCCCECEEEECCCCCHHHCCCCCEEEEEEECCEEEEEE 334 FNCSPHPYWLPNFMDVFTWSLPFVGEKVTEMLVNVLNICSSFEEAKGLDRINERMPPR Reprof ECCCCCCCCCCCCCEEEEEECCCCCHHHHHHHHHHHEECCEHHHHCCCCHHCCCCCCC PSIPRED ECCCCCCCCCCCCCCCCCCCHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHCCCCC DSSP ECCCCCCCCCHHHCCHHHHHHHHHHHHHHHHHHHHHCCCCCHHHHHHHHHHHHCCCCC |
<figure id="fig:ss_Q08209"> |
Prediction accuracy/precision
We compared the prediction performance of PSIPRED and Reprof via the Q3 score and the precision of the three secondary structure states H,E, and C. The Q3 score is identical to the accuracy, i.e. the number of correctly predicted states divided by the length of the protein. The precision of state X is the fraction of correct predictions of X, formally precision(X)=TP(X)/(TP(X)+FP(X)). <xr id="ss_acc"/> shows the results: PSIPRED clearly outperforms Reprof in all for cases. PSIPRED achieves an average accuracy of 87% which is significantly higher than 58% in case of Reprof. <figtable id="ss_acc">
Method | Q3 | Precision H | Precision E | Precision C |
---|---|---|---|---|
P04062 | ||||
PSIPRED | 0.831 | 0.830 | 0.872 | 0.810 |
Reprof | 0.553 | 1.000 | 0.455 | 0.592 |
P10775 | ||||
PSIPRED | 0.941 | 0.959 | 0.960 | 0.919 |
Reprof | 0.603 | 0.589 | 0.417 | 0.644 |
Q9X0E6 | ||||
PSIPRED | 0.890 | 1.000 | 0.895 | 0.720 |
Reprof | 0.580 | 0.562 | 0.917 | 0.458 |
Q08209 | ||||
PSIPRED | 0.833 | 0.842 | 0.902 | 0.812 |
Reprof | 0.579 | 0.762 | 0.293 | 0.743 |
</figtable>
Disorder
Disordered regions are regions with a varying three-dimensional structure. Nevertheless, such regions can be functionally highly important: regulation, signalling, and flexible ligand binding are only some examples. DisProt<ref name="DisProt">Vucetic S, Obradovic Z, Vacic V, et al. (2005). "DisProt: a database of protein disorder". Bioinformatics</ref> is a curated databases of proteins with experimentally determined disordered regions. IUPred<ref name="IUPred">Zsuzsanna Dosztányi, Veronika Csizmók, Péter Tompa and István Simon (2005). "The Pairwise Energy Content Estimated from Amino Acid Composition Discriminates between Folded and Intrinsically Unstructured Proteins". J. Mol. Biol.</ref> is a method for predicting disordered regions ab-initio, i.e. based solely on the protein sequence. We compared the predictions of IUPred with the annotations in the DisProt database for all four example proteins. IUPred was called to predict long, global disorder regions (confer the protocol for details). Residues involved in disordered regions were defined as those with a probability of at least 50%. These residues were compared to the DisProt annotations: either by the UniProt entry itself if available, or by a significant homolog (e-value < 1e-3) for which a DisProt entry existed. We measured the performance of IUPred via the precision (TP/(TP+FP)), sensitivity (TP/(TP+FN)), and specificity (TN/(TN+FP)).
P04062
Neither Glycosylceramidase (P04062) nor a homologous protein is annotated in DisProt. This might be due to lacking experimental data or, which is more likely, due to lacking disordered regions. The latter assumption is supported by the highly structured protein complex 1OGS (<xr id="fig:ss_P04062"/>). However, IUOred predicts some disordered residues with a probability >= 50%.
<figtable id="disorder_P04062">
Method | Disorder regions | |
---|---|---|
IUPred | 2, 3, 6, 90-93, 229-231, 235, 236 | |
DisProt | ||
Precision: 0% | Sensitivity: undef | Specificity: 98% |
</figtable>
P10775
P10775 is not annotated in DisProt itself, but there is a significant homolog (DP00554) with a disordered region from 31 to 50. This region is, however, is not covered by the pairwise alignment of P10775 and DP00554. Hence, IUpred correctly did not predict any disordered region (Specificity=100%).
Q9X0E6
There is neither an entry in DisProt which suggests a disordered region in Q9X0E6 and nor does IUPred predict such a region.
Q08209
Five disordered regions are annotated in DisProt for Q08209. These regions mainly cover the C-terminal end which exhibits several rather arbitrarily arranged alpha-helices (<xr id="fig:ss_Q08209"/>). All residues predicted by IUPred with a probability >= 50% are covered by DisProt annotations (precision 100%), but IUpred predicts only about half of all disordered regions (sensitivity 52%).
</figure>
<figtable id="disorder_Q08209">
</figtable> |
<figure id="fig:disorder_Q08209"> |
Transmembrane helices
PolyPhobius was used to predict transmembrane helices for our protein P04062 and other three proteins P35462, Q9YDF8 and P47863. The scripts which were used to do the prediction can be found here [protocal]. The prediction results from PolyPhobius were then compared with the membrane assignment of the structures for these proteins in OPM and/or PDBTM. For that purpose, a corresponding pdb structure for each protein was needed and was taken from uniprot(<xr id="tab:ss_sequences_trans"/>):
<figtable id="tab:ss_sequences_trans">
NAME | UniProtKB | PDB |
---|---|---|
Glucosylceramidase | P04062 | 1OGS |
D(3) dopamine receptor | P35462 | 3PBL |
DVoltage-gated potassium channel | Q9YDF8 | 1ORQ 1ORS |
Aquaporin-4 | P47863 | 2D57 |
</figtable>
Table 6: The proteins used to predict transmenbrane helix.
P04062
Our protein Glucosylceramidase is located in the lysosome, therefore contains no transmembrane regions. As expected, PolyPhobius has not reported any transmembrane regions. Instead, the most regions were predicted lying in the non-cytoplam and a signal peptide was found. Compared with that in uniprot, PolyPhobius has returned correct prediction (<xr id="tab:P04062_trans"/>):
<figtable id="tab:P04062_trans">
P04062 | Signal peptide | Others |
---|---|---|
PolyPhobius | 1-39 | 40-536 NON CYTOPLASMIC |
Uniprot | 1-39 | 40–536 Chain |
</figtable>
Table 7: Transmenbrane helix prediction results for P04062.
P35462
For the protein "dopamine receptor( P35462, )", only one pdb structure(3PBL) was found in uniprot. All three prediction tools PolyPhobius, OPM and PDBTM returned 7 regions of transmenbrane helix, which were confirmed by uniport. In <xr id="tab:P35462_trans"/> we can see that even the start and stop position of such regions showed very little deviation (maximally 4 bases). Figure 6 and 7 showed us the visualization of transmenbrane helix in OPM and PDBTM.
<figtable id="tab:P35462_trans">
P35462 | TRANSMEM 1 | TRANSMEM 2 | TRANSMEM 3 | TRANSMEM 4 | TRANSMEM 5 | TRANSMEM 6 | TRANSMEM 7 |
---|---|---|---|---|---|---|---|
PolyPhobius | 30-55 | 66-88 | 105-126 | 150-170 | 188-212 | 329-352 | 367-386 |
Uniprot | 33–55 | 66–88 | 105–126 | 150–170 | 188–212 | 330–351 | 367–388 |
OPM | 34-52 | 67-91 | 101-126 | 150–170 | 187-209 | 330-351 | 363-386 |
PDBTM | 35-52 | 68-84 | 109-123 | 152-166 | 191-206 | 334-347 | 368-382 |
</figtable>
Table 8: Transmenbrane helix prediction results for P35462.
Q9YDF8
For the protein KvAP("Voltage-gated potassium channel,Q9YDF8)", since no single pdb structure was found to cover the whole sequence, two pdb structures (1ORQ,1ORS) were used. Not like what happened to the previous protein, the prediction results varied among different prediction tools(<xr id="tab:Q9YDF8_trans"/>).
PolyPhobius returned the closest prediction as compared to that from uniprot. It reported one transmenbrane region less. By using two pdb structures, OPM has returned all the transmenbrane regions as those listed in uniprot, however the most of them showed quited big deviation of the start and stop positions. Although PolyPhobius reported one transmenbrane region less, the most showed much more accurate prediction of position. In PDBTM, even two pdb structures were used, still two regions were missed.
<figtable id="tab:Q9YDF8_trans">
Q9YDF8 | TRANSMEM 1 | TRANSMEM 2 | TRANSMEM 3 | TRANSMEM 4 | TRANSMEM 5 | TRANSMEM 6 | TRANSMEM 7 | TRANSMEM 8 |
---|---|---|---|---|---|---|---|---|
PolyPhobius | 42-60 | 68-88 | - | 108-129 | 137-157 | 163-184 | 196-213 | 224-244 |
Uniprot | 39–63 | 68–92 | 97–105(Intramembrane) | 109–125 | 129–145 | 160–184 | 196–208(Intramembrane) | 222 – 253 |
OPM_1ORQ | - | - | - | - | - | 153-172 | 183-195 | 207-225 |
OPM_1ORS | 25-46 | 55-78 | 86-97 | 100-107 | 117-148 | - | - | - |
PDBTM_1ORQ | 21-52 | 57-80 | - | - | 151-171 | - | - | 209-236 |
PDBTM_1ORS | 27-50 | 55-75 | 88-107 | 118-142 | - | - | - | - |
</figtable>
Table 10: Transmenbrane helix prediction results for Q9YDF8.
P47863
For the protein Aquaporin-4 ("Mercurial-insensitive water channel,P47863 )",the pdb structures 2D57 from X-ray with relative higher resolution was chosen. The prediction results did not vary too much among different prediction tools(<xr id="tab:P47863_trans"/>).
Not like PolyPhobius and that stated in uniprot, OPM reported two additional transmenbrane regions. PDBTM also returned two similar regions which were recognized as "Membrane loop". At transmembrane region 2, all three prediction tools showed an at least 5 bases leftward shift.
<figtable id="tab:P47863_trans">
P47863 | TRANSMEM_1 | TRANSMEM_2 | TRANSMEM_3 | TRANSMEM_4 | TRANSMEM_5 | TRANSMEM_6 | TRANSMEM_7 | TRANSMEM_8 |
---|---|---|---|---|---|---|---|---|
PolyPhobius | 34-58 | 70-91 | - | 115-136 | 156-177 | 188-208 | - | 231-252 |
Uniprot | 37–57 | 65–85 | - | 116–136 | 156–176 | 185–205 | - | 232–252 |
OPM | 34-56 | 70-88 | 98-107 | 112-136 | 156-178 | 189-203 | 214-223 | 231-252 |
PDBTM | 39-55 | 72-89 | 95-106 (Membrane loop) | 116-133 | 158-177 | 188-205 | 209-222 (Membrane loop) | 231-248 |
</figtable>
Table 12: Transmenbrane helix prediction results for P47863.
Signal peptides
The online server of SignalP with version 4.0 was used to predict signal peptides for our protein P04062 and other three proteins P02768, P47863 and P11279(<xr id="tab:ss_sequences_signa"/>): We also used Polyphobius and the information from uniprot and Signal Peptide Website to compare and validate the prediction results.
SignalP used three different scores, C, S and Y. Two additional scores, S-mean and the D-score, are reported in the output. The algorithm of SignalP employed two different neural networks, one for predicting the actual signal peptide and one for predicting the position of the signal peptidase cleavage site.
- The S-score is reported for every single amino acid position. If S-score is high, it suggests that the corresponding amino acid is part of a signal peptide. On the other hand, a low score indicate that the amino acid is part of a mature protein.
- The C-score is so called cleavage site score which should only be significantly high at the cleavage site.
- The Y-max is calculated by combining the C-score with the S-score and produces a better cleavage site prediction than the raw C-score alone.
- The S-mean is the average of the S-score for the length of the predicted signal peptide.
- The D-score is calculated as a weighted average of the S-mean and the Y-max scores.
Observing the above scores will help us to locate the signal peptides.
<figtable id="tab:ss_sequences_signa">
NAME | UniProtKB | PDB |
---|---|---|
Glucosylceramidase | P04062 | 1OGS |
Serum albumin | P02768 | 1AO6 |
Aquaporin-4 | P47863 | 2D57 |
Lysosome-associated membrane glycoprotein 1 | P11279 | - |
</figtable>
Table 11: The proteins used to predict signal peptides.
P04062
By using the standard option, SignalP did not return a positive result for our protein Glucosylceramidase P04062, even the s score and c score were quite high. After lowering the Cutoff value to 0.35, SignalP returned a finding of signal peptide which was validated by other tools and database(<xr id="tab:ss_sinal_p04062"/>).
# Measure Position Value Cutoff signal peptide? max. C 40 0.305 max. Y 40 0.396 max. S 36 0.684 mean S 1-39 0.323 D 1-39 0.367 0.350 YES Name=sp_P04062_GLCM_HUMAN SP='YES' Cleavage site between pos. 39 and 40: ASG-AR D=0.367 D-cutoff=0.350 Networks=SignalP-TM
<figtable id="tab:ss_sinal_p04062">
P04062 | position |
SignalP | 1-39 |
Signal Peptide Website | 1-39 |
uniprot | 1-39 |
PolyPhobius | 1-39 |
signal peptide sequence | MEFSSPSREECPKPLSRVSIMAGSLTGLLLLQAVSWASG |
</figtable>
Table 12: Signal peptides prediction of P04062.
Figure 15: signal peptides of Glucosylceramidase (P04062) showing in Signal Peptide Website.
P02768
For the protein Serum albumin (P02768), SignalP returned positive prediction which was validated by using other tools.
# Measure Position Value Cutoff signal peptide? max. C 19 0.710 max. Y 19 0.798 max. S 2 0.929 mean S 1-18 0.890 D 1-18 0.848 0.450 YES Name=sp_P02768_ALBU_HUMAN SP='YES' Cleavage site between pos. 18 and 19: AYS-RG D=0.848 D-cutoff=0.450 Networks=SignalP-noTM
<figtable id="tab:ss_sinal_p02768">
P02768 | position |
SignalP | 1-18 |
Signal Peptide Website | 1-18 |
uniprot | 1-18 |
PolyPhobius | 1-18 |
signal peptide sequence | MKWVTFISLLFLFSSAYS |
</figtable>
Table 13: Signal peptides prediction of P02768.
Figure 17: signal peptides of Serum albumin (P02768) showing in Signal Peptide Website.
P47863
For the protein Aquaporin-4 P47863, as expected from its alternative name "Mercurial-insensitive water channel", SignalP returned very low scores for which it was not meaningful to lower the cutoff value. All the other tools reported no finding either (<xr id="tab:ss_sinal_P47863"/>).
# Measure Position Value Cutoff signal peptide? max. C 41 0.209 max. Y 41 0.164 max. S 56 0.183 mean S 1-40 0.139 D 1-40 0.154 0.500 NO Name=sp_P47863_AQP4_RAT SP='NO' D=0.154 D-cutoff=0.500 Networks=SignalP-TM
<figtable id="tab:ss_sinal_P47863">
P47863 | position |
SignalP | - |
Signal Peptide Website | - |
uniprot | - |
PolyPhobius | - |
signal peptide sequence | - |
</figtable>
Table 14: Signal peptides prediction of P47863.
P11279
For the protein Lysosome-associated membrane glycoprotein 1 (P11279), as expected from its name, a long signal peptide was found with very high score:
# Measure Position Value Cutoff signal peptide? max. C 29 0.910 max. Y 29 0.939 max. S 21 0.996 mean S 1-28 0.962 D 1-28 0.952 0.450 YES Name=sp_P11279_LAMP1_HUMAN SP='YES' Cleavage site between pos. 28 and 29: ASA-AM D=0.952 D-cutoff=0.450 Networks=SignalP-noTM
<figtable id="tab:ss_sinal_P11279">
P11279 | position |
SignalP | 1-28 |
Signal Peptide Website | 1-28 |
uniprot | 1-28 |
PolyPhbius | 1-28 |
signal peptide sequence | MAAPGSARRPLLLLLLLLLLGLMHCASA |
</figtable>
Table 15: Signal peptides prediction of P11279.
Figure 20: signal peptides of Lysosome-associated membrane glycoprotein 1 (P11279) showing in Signal Peptide Website.
GO terms
GOPET
GOPET was used to predict GO terms for our protein Glucosylceramidase (P04062). The results was shown in <xr id="tab:ss_gopet"/> where only three GO terms were returned. Resetting the search options like lower the Confidence threshold did not change the prediction results.
All three GO terms belong to "Molecular Function Ontology" and are all involved in "hydrolase" or "glycosyl bonds". Since we have already known the function of our protein Glucosylceramidase (P04062) is to cleaves the glucosidic bonds of glucosylceramide and synthetic beta-glucosides. All these three predicted GO terms are correct. However, if compared with the GO terms list from uniprot(QuickGo) or the prediction from ProtFun which is shown below, GOPET returned quite fewer results. It implies that GOPET prefers to report only the GO terms with high confidence.
<figtable id="tab:ss_gopet">
GOid | Aspect | Confidence | GO term | validate in QuickGo |
---|---|---|---|---|
GO:0016787 | F | 98% | hydrolase activity | yes |
GO:0004348 | F | 97% | glucosylceramidase activity | yes |
GO:0016798 | F | 97% | hydrolase activity acting on glycosyl bonds | yes |
</figtable>
Table 16: GO terms prediction results for P04062 in GOPET.
ProtFun
ProtFun2.0 was also used to predict GO terms for our protein Glucosylceramidase (P04062). The results were listed in <xr id="tab:ss_protfun"/>. In the table, for each predicted terms, there are two scores. The first one is the estimated probability that the entry belongs to the class.It is influenced by the prior probability of that class. The second number represents the odds that the sequence belongs to that class/category. It is independent of the prior probability.
Compared with GOPET which reported only GO identifier, ProtFun returned much more information about the input protein:
At first, the input protein was classified into 12 different functional categories based on the scheme developed by Monica Riley for E. coli in 1993<ref name="protfun_scheme">Monica Riley (1993). "Functions of the gene products of Escherichia coli.". Microbiol Rev[1]</ref>. As we can see in the table under "Functional category" section where all the functional categories were listed with their probabilities and odds. For our protein Glucosylceramidase (P04062), "Cell_envelope" in green was considered by ProrFun as the most possible functional category that the protein should belong to.
At the second part, the protein was classified as enzyme or non-enzyme where our protein was correctly predicted as Enzyme.
At the third part,the protein was further classified into 6 different enzyme classes and possible EC numbers were also given. We can see that in that section, ProtFun did not return a significant prediction. No one class from them had a significant higher probability and odds than the others. Even though, Hydrolase (EC 3.-.-.-) with the highest probability 0.272 gave us a good suggestion since we have already known that our protein has EC=3.2.1.45.
At the last part, we can see different Gene Ontology categories with the prediction scores. The Gene Ontology categories listed here seemed different with that from the Gene Ontology. For our protein, ProtFun returned that its most possible Gene Ontology is "Immune_response" which stayed in question. Since there are no further information available about that in ProtFun, it was not easy to judge that prediction.
<figtable id="tab:ss_protfun">
Functional category | Prob | Odds |
Amino_acid_biosynthesis | 0.035 | 1.593 |
Biosynthesis_of_cofactors | 0.182 | 2.528 |
Cell_envelope | 0.504 | 8.262 |
Cellular_processes | 0.032 | 0.438 |
Central_intermediary_metabolism | 0.382 | 6.063 |
Biosynthesis_of_cofactors | 0.182 | 2.528 |
Energy_metabolism | 0.067 | 0.740 |
Fatty_acid_metabolism | 0.027 | 2.088 |
Purines_and_pyrimidines | 0.538 | 2.213 |
Regulatory_functions | 0.031 | 0.191 |
Replication_and_transcription | 0.126 | 0.471 |
Translation | 0.082 | 1.863 |
Transport_and_binding | 0.560 | 1.365 |
Enzyme/nonenzyme | Prob | Odds |
Enzyme | 0.773 | 2.698 |
Nonenzyme | 0.227 | 0.318 |
Enzyme class | Prob | Odds |
Oxidoreductase (EC 1.-.-.-) | 0.083 | 0.399 |
Transferase (EC 2.-.-.-) | 0.228 | 0.660 |
Hydrolase (EC 3.-.-.-) | 0.272 | 0.859 |
Lyase (EC 4.-.-.-) | 0.045 | 0.961 |
Isomerase (EC 5.-.-.-) | 0.011 | 0.345 |
Ligase (EC 6.-.-.-) | 0.017 | 0.332 |
Gene Ontology category | Prob | Odds |
Signal_transducer | 0.054 | 0.251 |
Receptor | 0.027 | 0.158 |
Hormone | 0.001 | 0.206 |
Structural_protein | 0.002 | 0.087 |
Transporter | 0.024 | 0.222 |
Ion_channel | 0.018 | 0.307 |
Voltage-gated_ion_channel | 0.004 | 0.195 |
Cation_channel | 0.012 | 0.268 |
Transcription | 0.070 | 0.550 |
Transcription_regulation | 0.030 | 0.237 |
Stress_response | 0.085 | 0.962 |
Immune_response | 0.153 | 1.804 |
Growth_factor | 0.005 | 0.376 |
Metal_ion_transport | 0.009 | 0.020 |
</figtable>
Table 17: GO terms prediction results for P04062 in Protfun2.0.
Pfam
At last, Pfam was used to find pfam family for our protein Glucosylceramidase (P04062). As we can see in the following figure, one domain was found by Pfam: Glycoside hydrolase family 30.
It is a group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycoside hydrolases, based on sequence similarity, has led to the definition of >100 different families. One can find this system on the CAZy web site.
One of the famous enzyme under this group is the mammalian Glucosylceramidase (P04062) which is the protein we are working with. It cleaves the glucosidic bonds of glucosylceramide and synthetic beta-glucosides.Any one of over 50 different mutations in the gene of glucocerebrosidase have been found to affect activity of this hydrolase, producing variants of Gaucher disease, the most prevalent lysosomal storage disease.
References
<references/>