Difference between revisions of "Sequence-based mutation analysis (Phenylketonuria)"

From Bioinformatikpedia
(→‎Analyze SNPs)
(→‎Analyze SNPs)
Line 36: Line 36:
 
Looking at the substitution matrices you can see that the substitutions analyzed in this task never have the worst values.<br>
  
Looking at the substitution matrices you can see that the substitutions analyzed in this task never have the worst values.<br>
 
As 2pah begins at position 118 for amino acids 20 and 103 the substitutions cannot be shown with pymol.
  
As 2pah begins at position 118 for amino acids 20 and 103 the substitutions cannot be shown with pymol.
  +
<figure id ="all_mutations"> [[File:pah_all_mutations.png|'''<caption>''' All eight substitutions shown in 2pah (green), where the original amino acids are shown in yellow and the mutated one in purple.</caption>]]</figure>
−
===Mammalian Homologous Sequences===
  
−
We used a normal BLAST search on the mammalian database of Uniprot and filter the [http://www.uniprot.org/blast/uniprot/2013062760F4KDGT75 results] per hand for double entries. Altogether we found 22 homologues sequences (<xr id="IDs"/>).
  
−
<figtable id="IDs">
  
−
{| border="1" style="text-align:center;" cellpadding="5" cellspacing="0" align="center"
  
−
|-
  
−
! colspan="8" style="background:#32CD32;" | UniProt-IDs homologue to P00439
  
−
|-
  
−
| [http://www.uniprot.org/uniprot/H2Q6R0 H2Q6R0] || [http://www.uniprot.org/uniprot/G3S964 G3S964] || [http://www.uniprot.org/uniprot/G1R3M2 G1R3M2] || [http://www.uniprot.org/uniprot/G7PJC2 G7PJC2] || [http://www.uniprot.org/uniprot/F7HMW9 F7HMW9] || [http://www.uniprot.org/uniprot/F7I717 F7I717] || [http://www.uniprot.org/uniprot/F7BKF9 F7BKF9]
 
−
|-
  
−
|[http://www.uniprot.org/uniprot/F6XY00 <span style="color:green">F6XY00</span>] || [http://www.uniprot.org/uniprot/E2R366 E2R366] || [http://www.uniprot.org/uniprot/G1T8B6 G1T8B6] || [http://www.uniprot.org/uniprot/H2NIF5 H2NIF5] || [http://www.uniprot.org/uniprot/M3YKN3 M3YKN3] || [http://www.uniprot.org/uniprot/M9P0Q7 M9P0Q7] || [http://www.uniprot.org/uniprot/H0WTI6 H0WTI6]
  
−
|-
  
−
|[http://www.uniprot.org/uniprot/M3W9R1 M3W9R1] || [http://www.uniprot.org/uniprot/G3TIW0 G3TIW0] || [http://www.uniprot.org/uniprot/G1LIM6 G1LIM6] || [http://www.uniprot.org/uniprot/Q2KIH7 Q2KIH7] || [http://www.uniprot.org/uniprot/G1P4I7 G1P4I7] || [http://www.uniprot.org/uniprot/P16331 P16331] || [http://www.uniprot.org/uniprot/M9P0Y5 M9P0Y5]
  
−
|}
  
−
<center><small>'''<caption>''' UniProt IDs of the 22 mammalian homologues sequences of the PAH protein. The only ID that has a sequence that is not 100% identical to P00439 is marked in green.</caption></small></center>
  
−
</figtable>
  
−  
−
For the creation of the multiple alignment we used clustalw. Nevertheless, only F6XY00 has a slightly different sequence. Like E2R366 F6XY00 is a homologue sequence of the human PAH in canis familiaris. On the positions of our substitutions the sequences are conserved.
  
−  
 
===Ala-259-Val===
  
===Ala-259-Val===
 
<figure id="A259V">
  
<figure id="A259V">
Line 324: Line 307:
 
|colspan="2"|'''&rArr;''' prediction: '''non-neutral'''
  
|colspan="2"|'''&rArr;''' prediction: '''non-neutral'''
 
|}
  
|}
  +
===Mammalian Homologous Sequences===
  +
We used a normal BLAST search on the mammalian database of Uniprot and filter the [http://www.uniprot.org/blast/uniprot/2013062760F4KDGT75 results] per hand for double entries. Altogether we found 22 homologues sequences (<xr id="IDs"/>).
  +
<figtable id="IDs">
  +
{| border="1" style="text-align:center;" cellpadding="5" cellspacing="0" align="center"
  +
|-
  +
! colspan="8" style="background:#32CD32;" | UniProt-IDs homologue to P00439
  +
|-
  +
| [http://www.uniprot.org/uniprot/H2Q6R0 H2Q6R0] || [http://www.uniprot.org/uniprot/G3S964 G3S964] || [http://www.uniprot.org/uniprot/G1R3M2 G1R3M2] || [http://www.uniprot.org/uniprot/G7PJC2 G7PJC2] || [http://www.uniprot.org/uniprot/F7HMW9 F7HMW9] || [http://www.uniprot.org/uniprot/F7I717 F7I717] || [http://www.uniprot.org/uniprot/F7BKF9 F7BKF9]
  +
|-
  +
|[http://www.uniprot.org/uniprot/F6XY00 <span style="color:green">F6XY00</span>] || [http://www.uniprot.org/uniprot/E2R366 E2R366] || [http://www.uniprot.org/uniprot/G1T8B6 G1T8B6] || [http://www.uniprot.org/uniprot/H2NIF5 H2NIF5] || [http://www.uniprot.org/uniprot/M3YKN3 M3YKN3] || [http://www.uniprot.org/uniprot/M9P0Q7 M9P0Q7] || [http://www.uniprot.org/uniprot/H0WTI6 H0WTI6]
  +
|-
  +
|[http://www.uniprot.org/uniprot/M3W9R1 M3W9R1] || [http://www.uniprot.org/uniprot/G3TIW0 G3TIW0] || [http://www.uniprot.org/uniprot/G1LIM6 G1LIM6] || [http://www.uniprot.org/uniprot/Q2KIH7 Q2KIH7] || [http://www.uniprot.org/uniprot/G1P4I7 G1P4I7] || [http://www.uniprot.org/uniprot/P16331 P16331] || [http://www.uniprot.org/uniprot/M9P0Y5 M9P0Y5]
  +
|}
  +
<center><small>'''<caption>''' UniProt IDs of the 22 mammalian homologues sequences of the PAH protein. The only ID that has a sequence that is not 100% identical to P00439 is marked in green.</caption></small></center>
  +
</figtable>
  +
  +
For the creation of the multiple alignment we used clustalw. Nevertheless, only F6XY00 has a slightly different sequence. Like E2R366 F6XY00 is a homologue sequence of the human PAH in canis familiaris. On the positions of our substitutions the sequences are conserved.
   
 
== Prediction ==
  
== Prediction ==

Revision as of 09:20, 3 July 2013

Summary

...

Mutation dataset

Lab journal

SNPs
AA - three letter AA - one letter Nucleotides
Ala259Val A259V C776T
Arg123Ile R123I G368T
Gln20Leu Q20L A59T
Gln172His Q172H G516T
Gly103Ser G103S G307A
Ile421Thr I421T T1262C
Lys341Thr K341T A1022C
Phe392Ser F392S T1175C
Pro416Gln P416Q C1247A
Thr266Ala T266A A796C

Analyze SNPs

The complete results of the prediction tools can be found in the Lab journal.
Looking at the substitution matrices you can see that the substitutions analyzed in this task never have the worst values.
As 2pah begins at position 118 for amino acids 20 and 103 the substitutions cannot be shown with pymol. <figure id ="all_mutations"> All eight substitutions shown in 2pah (green), where the original amino acids are shown in yellow and the mutated one in purple.</figure>

Ala-259-Val

<figure id="A259V">

Mutation of alanine (yellow) to valine (purple) at position 259 of 2pah (green).

</figure>

Properties alanine is small, non-polar, neutral and hydrophobic
valine is aliphatic, small, non-polar, neutral and hydrophobic
Structure HELIX
Substitution Matrices middle ranged value → neutral substitution
PSSM alanine seems to be highly conserved with a value of 70. Valine has a value of 1.
SIFT TOLERATED with a score of 0.16
PolyPhen2 probably damaging with a score of 1.000
SNAP non-neutral, RI = 2
MutationTaster disease causing
⇒ prediction: neutral/non-neutral???


Arg-123-Ile

<figure id="R123I">

Mutation of arginine (yellow) to isoleucine (purple) at position 123 of 2pah (green).

</figure>

Properties arginine is positively charged, polar and hydrophilic
isoleucine is aliphatic, neutral, non-polar and hydrophobic
Substitution Matrices low value → bad substitution
Structure LOOP
PSSM arginine seems to be highly conserved with a value of 49. Isoleucine has a value of 7.
SIFT AFFECT PROTEIN FUNCTION with a score of 0.00
PolyPhen2 possibly damaging with a score of 0.807 / 0.582
SNAP neutral, RI = 0
MutationTaster disease causing
⇒ prediction: non-neutral

Gln-20-Leu

Properties glutamine is neutral, polar and hydrophilic
leucine is aliphatic, neutral, non-polar and hydrophilic
Substitution Matrices low value → bad substitution
Structure LOOP
PSSM seems to be not conserved as the value for aspartic acid (17) is higher than for glutamine (14).

leucine has a value of 6.

SIFT TOLERATED with a score of 0.90
PolyPhen2 benign with a score of 0.000
SNAP neutral, RI = 0
MutationTaster disease causing
⇒ prediction: neutral

Gln-172-His

<figure id="Q172H">

Mutation of glutamine (yellow) to histidine (purple) at position 172 of 2pah (green).

</figure>

Properties glutamine is neutral, polar and hydrophilic
histidine is positively charged, polar and hydrophobic
Substitution Matrices middle ranged → neutral substitution
Structure LOOP
PSSM ...
...
SIFT AFFECT PROTEIN FUNCTION with a score of 0.03.
PolyPhen2 possibly damaging with a score of 0.705 / benign with a score of 0.170
SNAP neutral, RI = 4
MutationTaster disease causing
⇒ prediction: neutral

Gly-103-Ser

Properties glycine is small, neutral, non-polar and hydrophilic
serine is small, neutral, polar and hydrophilic
Substitution Matrices middle ranged value → neutral substitution
Structure HELIX
PSSM glycine seems to be completely unconserved as every amino acid has the same value 5.
SIFT TOLERATED with a score of 0.05
PolyPhen2 benign with a score of 0.003 / 0.006
SNAP neutral, RI = 6
MutationTaster disease causing
⇒ prediction: neutral

Ile-421-Thr

<figure id="I421T">

Mutation of isoleucine (yellow) to threonine (purple) at position 421 of 2pah (green).

</figure>

Properties isoleucine is aliphatic, neutral, non-polar and hydrophobic
threonine is small, neutral, polar and hydrophobic
Substitution Matrices low value for Blosum62, middle ranged value for PAM1/250 → bad or neutral substitution
Structure STRAND
PSSM isoleucine with a value of 32 seems to be not conserved as valine has a value of 41.

threonine has a value of 0.

SIFT AFFECT PROTEIN FUNCTION with a score of 0.00.
PolyPhen2 possibly damaging with a score of 0.667 / probably damaging with a score of 0.913
SNAP non-neutral, RI = 2
MutationTaster disease causing
⇒ prediction: non-neutral

Lys-341-Thr

<figure id="K341T">

Mutation of lysine (yellow) to threonine (purple) at position 341 of 2pah (green).

</figure>

Properties lysine is positively charged, polar and hydrophobic
threonine is small, neutral, polar and hydrophobic
Substitution Matrices low value for Blosum62, middle ranged value for PAM1/250 → bad or neutral substitution
Structure STRAND
PSSM lysine seems to be conserved with a value of 39.

threonine has a value of 1.

SIFT AFFECT PROTEIN FUNCTION with a score of 0.00.
PolyPhen2 probably damaging with a score of 1.000 / 0.996
SNAP non-neutral, RI = 3
MutationTaster disease causing
⇒ prediction: non-neutral

Phe-392-Ser

<figure id="F392S">

Mutation of phenylalanine (yellow) to serine (purple) at position 392 of 2pah (green).

</figure>

Properties phenylalanine is aromatic, neutral, non-polar and hydrophobic
serine is small, neutral, polar and hydrophilic
Substitution Matrices low value → bad substitution
Structure HELIX
PSSM phenylalanine seems to be highly conserved with a value of 49. Serine has a value of 0.
SIFT AFFECT PROTEIN FUNCTION with a score of 0.00.
PolyPhen2 probably damaging with a score of 1.000
SNAP non-neutral, RI = 1
MutationTaster disease causing
⇒ prediction: non-neutral

Pro-416-Gln

<figure id="P416Q">

Mutation of proline (yellow) to glutamine (purple) at position 416 of 2pah (green).

</figure>

Properties proline is small, neutral, non-polar and hydrophilic
glutamine is neutral, polar and hydrophilic
Substitution Matrices low value → bad substitution
Structure LOOP
PSSM proline seems to be highly conserved with a value of 66. Glutamine has a value of 0.
SIFT AFFECT PROTEIN FUNCTION with a score of 0.00.
PolyPhen2 probably damaging with a score of 0.996 / 0.985
SNAP non-neutral, RI = 1
MutationTaster disease causing
⇒ prediction: non-neutral

Thr-266-Ala

<figure id="T266A">

Mutation of threonine (yellow) to alanine (purple) at position 266 of 2pah (green).

</figure>

Properties threonine is small, neutral, polar and hydrophobic
alanine is small, non-polar, neutral and hydrophobic
Substitution Matrices middle ranged value &rarr neutral substitution
Structure HELIX
PSSM threonine seems to be conserved with a value of 43. Alanine has a value of 29.
SIFT AFFECT PROTEIN FUNCTION with a score of 0.00.
PolyPhen2 probably damaging with a score of 1.000
SNAP neutral, RI = 2
MutationTaster disease causing
⇒ prediction: non-neutral

Mammalian Homologous Sequences

We used a normal BLAST search on the mammalian database of Uniprot and filter the results per hand for double entries. Altogether we found 22 homologues sequences (<xr id="IDs"/>). <figtable id="IDs">

UniProt-IDs homologue to P00439
H2Q6R0 G3S964 G1R3M2 G7PJC2 F7HMW9 F7I717 F7BKF9
F6XY00 E2R366 G1T8B6 H2NIF5 M3YKN3 M9P0Q7 H0WTI6
M3W9R1 G3TIW0 G1LIM6 Q2KIH7 G1P4I7 P16331 M9P0Y5
UniProt IDs of the 22 mammalian homologues sequences of the PAH protein. The only ID that has a sequence that is not 100% identical to P00439 is marked in green.

</figtable>

For the creation of the multiple alignment we used clustalw. Nevertheless, only F6XY00 has a slightly different sequence. Like E2R366 F6XY00 is a homologue sequence of the human PAH in canis familiaris. On the positions of our substitutions the sequences are conserved.

Prediction

SNP-Prediction
SNP Prediction Validation
Ala259Val neutral/non-neutral??? HGMD
Arg123Ile non-neutral ? dbSNP
Gln20Leu neutral wrong HGMD
Gln172His neutral ? dbSNP
Gly103Ser neutral wrong HGMD
Ile421Thr non-neutral ? dbSNP
Lys341Thr non-neutral correct HGMD
Phe392Ser non-neutral ? dnSNP
Pro416Gln non-neutral correct HGMD
Thr266Ala non-neutral correct dbSNP

References

<references/>

Retrieved from "https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Sequence-based_mutation_analysis_(Phenylketonuria)&oldid=33325"