Difference between revisions of "Rs61747114"
(→Gromacs Energy Comparison) |
(→SNAP Prediction) |
||
(19 intermediate revisions by one other user not shown) | |||
Line 20: | Line 20: | ||
== Sequence-based Mutation Analysis == |
== Sequence-based Mutation Analysis == |
||
− | === Pysicochemical |
+ | === Pysicochemical Properties === |
First of all, we explored the amino acid properties and compared them for the original and the mutated amino acid. Therefore we created the possible effect that the mutation could have on the protein. |
First of all, we explored the amino acid properties and compared them for the original and the mutated amino acid. Therefore we created the possible effect that the mutation could have on the protein. |
||
Line 31: | Line 31: | ||
|aliphatic, hydrophobic, neutral |
|aliphatic, hydrophobic, neutral |
||
|aromatic, hydrophobic, neutral |
|aromatic, hydrophobic, neutral |
||
− | |Leu is an aliphatic amino acid, |
+ | |Leu is an aliphatic amino acid, whereas Phe is an aromatic amino acid. This means, that Phe has an aromatic ring in its structure. But both amino acids are relatively big and so it is possible, that the exchange of this amino acids does not change the structure of the protein that much. Therefore, we suggest it is possible, that the protein works. |
|- |
|- |
||
|} |
|} |
||
Line 39: | Line 39: | ||
---- |
---- |
||
− | === |
+ | === Visualization of the Mutation === |
− | In the next step, we created the visualization of the |
+ | In the next step, we created the visualization of the mutation with PyMol. Therefore we created a picture for the original amino acid (Figure 1), for the new mutated amino acid (Figure 2) and finally for both together in one picture whereas the mutation is white colored (Figure 3). The following pictures display that the mutated amino acid Phenylalanine looks very different to Leucine. Leucine forks at the end of the rest whereas Phenylalanine has an huge aromatic ring. The only thing that agrees is the orientation and the length of the residues before the fork or the aromatic ring. All in all the difference probably prevails the agreement and therefore the mutation will probably cause changes in protein structure and function. |
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
− | |picture original |
+ | |picture original amino acid |
− | |picture mutated |
+ | |picture mutated amino acid |
|combined picture |
|combined picture |
||
|- |
|- |
||
− | |[[Image:L248.png|thumb|150px|Amino acid Leucine]] |
+ | |[[Image:L248.png|thumb|150px|Figure 1: Amino acid Leucine]] |
− | |[[Image:248W.png|thumb|150px|Amino acid Phenylalanine]] |
+ | |[[Image:248W.png|thumb|150px|Figure 2: Amino acid Phenylalanine]] |
− | |[[Image:L248W.png|thumb|150px|Picture which visualize the mutation]] |
+ | |[[Image:L248W.png|thumb|150px|Figure 3: Picture which visualize the mutation]] |
|- |
|- |
||
|} |
|} |
||
Line 58: | Line 58: | ||
---- |
---- |
||
− | === |
+ | === Substitution Matrices Values === |
− | Afterwards, we looked at the values of the substitution matrices PAM1, PAM250 and BLOSSUM62. Therefore we looked detailed at the three values: the value for |
+ | Afterwards, we looked at the values of the substitution matrices PAM1, PAM250 and BLOSSUM62. Therefore we looked detailed at the three values: the value for this amino acid substitution, the most frequent value for the substitution of the examined amino acid and the rarest substitution. |
− | In this case, the substitution of Leucine to Phenylalanine has average values that are almost as far away from the rarest value as |
+ | In this case, the substitution of Leucine to Phenylalanine has average values that are almost as far away from the rarest value as from the most frequent value for all three substitution matrices. Only BLOSOUM62 is a little bit closer to the most frequent substitution. Therefore, the values from both PAMs are not really significant and therefore we are not able to determine effects on the protein for these two matrices. Otherwise, according to BLOSOUM62 a mutation at this position will probably not cause structural changes which can affect functional changes. |
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
Line 69: | Line 69: | ||
|colspan="3" | BLOSOUM 62 |
|colspan="3" | BLOSOUM 62 |
||
|- |
|- |
||
− | |value |
+ | |value amino acid |
|most frequent substitution |
|most frequent substitution |
||
|rarest substitution |
|rarest substitution |
||
− | |value |
+ | |value amino acid |
|most frequent substitution |
|most frequent substitution |
||
|rarest substitution |
|rarest substitution |
||
− | |value |
+ | |value amino acid |
|most frequent substitution |
|most frequent substitution |
||
|rarest substitution |
|rarest substitution |
||
Line 105: | Line 105: | ||
|colspan="3" | PSSM |
|colspan="3" | PSSM |
||
|- |
|- |
||
− | |value |
+ | |value amino acid |
|most frequent substitution |
|most frequent substitution |
||
|rarest substitution |
|rarest substitution |
||
Line 121: | Line 121: | ||
=== Conservation Analysis with Multiple Alignments === |
=== Conservation Analysis with Multiple Alignments === |
||
− | As a next step we created a multiple alignment which contains the HEXA sequence and 9 other mammalian homologous sequences from uniprot. Afterwards we looked at the position of the different mutations and looked at the conservation level |
+ | As a next step we created a multiple alignment which contains the HEXA sequence and 9 other mammalian homologous sequences from [[http://www.uniprot.org UniProt]]. Afterwards we looked at the position of the different mutations and looked at the conservation level at this position, which can be seen in Figure 4. The regarded mutation is presented by the second colored column. Here we can see, that all the other mammalians have also the amino acid Leucine on this position. Therefore, the amino acid at this position is highly conserved and a mutation there will cause probably huge structural and functional changes for the protein. |
− | [[Image:mut_5.png|thumb|center|600px|Mutation in the multiple alignment]] |
+ | [[Image:mut_5.png|thumb|center|600px|Figure 4: Mutation in the multiple alignment]] |
Line 131: | Line 131: | ||
=== Secondary Structure Mutation Analysis === |
=== Secondary Structure Mutation Analysis === |
||
− | As a next step we compared the different results of the secondary structure prediction tools JPred and PsiPred. Afterwards we can examine in which secondary structure element and where therein the mutation takes place. This can give an overview of how |
+ | As a next step we compared the different results of the secondary structure prediction tools JPred and PsiPred. Afterwards we can examine in which secondary structure element and where therein the mutation takes place. This can give an overview of how drastic the mutation can be. In this case both tools agree and predict at the position of the mutation a coil. This has as result, that the mutation at this position would not destroy or split a secondary structure element. It will probably only changes the coil between two secondary structure elements, but this can sometimes also cause a change of the the following secondary structure. We think that a drastic change of the protein structure and its function is unlikely because the mutation does not affect a secondary structure element. The change of the coil will probably only take places between two secondary structure elements which will probably not change. |
JPred: |
JPred: |
||
Line 140: | Line 140: | ||
''' Comparison with the real Structure: ''' |
''' Comparison with the real Structure: ''' |
||
− | Afterwards we also visualize the position of the |
+ | Afterwards we also visualize the position of the mutation (red) in the real 3D-structure of PDB and compare it with the predicted secondary structure (Figure 5 and Figure 6). The visualization can therefore like above the predicted secondary structure display if the mutation is in a secondary structure element or in some other regions. |
− | Here in this case the mutation position disagrees with the position of the predicted secondary structure and is within an alpha-helix. This means a mutation will probably destroy or split the alpha helix which affects |
+ | Here in this case the mutation position disagrees with the position of the predicted secondary structure and is within an alpha-helix. This means a mutation will probably destroy or split the alpha helix which affects drastic structural changes on the protein. We think that a structural change is very likely, because it is within a secondary structure element and will therefore cause extreme changes. |
{| |
{| |
||
− | | [[Image:248_mut.png|thumb|250px|Mutation at position 248]] |
+ | | [[Image:248_mut.png|thumb|250px|Figure 5: Mutation at position 248]] |
− | | [[Image:248_mut_detail.png|thumb|250px|Mutation at position 248 - detailed view]] |
+ | | [[Image:248_mut_detail.png|thumb|250px|Figure 6: Mutation at position 248 - detailed view]] |
|} |
|} |
||
Line 155: | Line 155: | ||
=== SNAP Prediction === |
=== SNAP Prediction === |
||
− | Next, we looked at the result of the SNAP prediction. For this prediction we took the amino acid of the certain position and checked every possible amino acid mutation. Afterwards we extract the result for Phenylalanine which is the real mutation in this case. SNAP has a result that the exchange from Leucine to Phenylalanine at this position is neutral with a relative high accuracy. This means that this certain mutation on this position cause very likely no structural and functional changes of the protein. |
+ | Next, we looked at the result of the SNAP prediction. For this prediction we took the amino acid of the certain position and checked every possible amino acid mutation. Afterwards we extract the result for Phenylalanine which is the real mutation in this case. SNAP has as a result that the exchange from Leucine to Phenylalanine at this position is neutral with a relative high accuracy. This means that this certain mutation on this position cause very likely no structural and functional changes of the protein. |
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
Line 178: | Line 178: | ||
=== SIFT Prediction === |
=== SIFT Prediction === |
||
− | Next, we used SIFT Prediction which displays if a mutation is neutral or not. Therefore, it first shows a row which contains a score for the particular |
+ | Next, we used SIFT Prediction which displays if a mutation is neutral or not. Therefore, it first shows a row which contains a score for the particular mutation position to a certain amino acid. The amino acid which are not tolerated at this position are colored red. Besides, it also constructs a table which lists the amino acids that are predicted as tolerated and not-tolerated (compare Figure 8). |
− | In this case, there are seven substitutions that are tolerated: Tyrosine, Isoleucine, Aspartic acid, Valine, Phenylalanine, Methionine and Leucine. The substitution to |
+ | In this case, there are seven substitutions that are tolerated: Tyrosine, Isoleucine, Aspartic acid, Valine, Phenylalanine, Methionine and Leucine. The substitution to Phenylalanine is tolerated at this position. This means that this mutation at this position is probably neutral and will not cause any structural and function changes of the protein. |
SIFT Matrix:<br> |
SIFT Matrix:<br> |
||
Line 186: | Line 186: | ||
{| |
{| |
||
− | | [[Image:sift_legend.png|center]] |
+ | | [[Image:sift_legend.png|center|800px|thumb|Figure 7: Legend]] |
|- |
|- |
||
− | | [[Image:248_sift.png.png|center |
+ | | [[Image:248_sift.png.png|center|thumb|800px|Figure 8: SIFT Table<br> |
+ | Threshold for intolerance is 0.05.<BR>Amino acid color code: non-polar, <font color=green>uncharged polar</font>, <font color=red>basic</font>, <font color=blue>acidic</font>. <BR>Capital letters indicate amino acids appearing in the alignment, lower case letters result from prediction.]] |
||
|} |
|} |
||
− | SIFT Table<br> |
||
− | Threshold for intolerance is 0.05.<BR>Amino acid color code: nonpolar, <font color=green>uncharged polar</font>, <font color=red>basic</font>, <font color=blue>acidic</font>. <BR>Capital letters indicate amino acids appearing in the alignment, lower case letters result from prediction. |
||
<br> |
<br> |
||
{| class="wikitable centered" |
{| class="wikitable centered" |
||
Line 206: | Line 205: | ||
=== PolyPhen2 Prediction === |
=== PolyPhen2 Prediction === |
||
− | Finally, we also regarded the PolyPhen2 prediction for this |
+ | Finally, we also regarded the PolyPhen2 prediction for this mutation. This prediction visualizes how strongly damaging the mutation probably will be. Therefore it gives the result for two possible cases: HumDiv and HumVar. HumDiv is the preferred model for evaluation rare allelse, dense mapping of regions identified by genome-wide association studies and analysis of neutral selection. In contrast, HumVar is the preferred model for diagnostic of Mendelian diseases which require distinguishing mutations with drastic effects from all remaining human variations including abundant mildly deleterious alleles. We decided to look at both possible models, which agreed the most cases. |
− | In this case, the models disagree. The first model HumDic predicts that the mutation is benign whereas the second model HumVar predicts that the mutation is possibly damaging. When we look at the scores we can see that the difference is with 0.213 not so high, which displays that they are moving in the same direction. Looking at both together, there is probably no huge structural change but with a small probability there can possibly be a damage of the protein. We assume with the help of |
+ | In this case, the models disagree. The first model HumDic predicts that the mutation is benign (which could be seen on figure 9) whereas the second model HumVar predicts that the mutation is possibly damaging (compare Figure 10). When we look at the scores we can see that the difference is with 0.213 not so high, which displays that they are moving in the same direction. Looking at both together, there is probably no huge structural change but with a small probability there can possibly be a damage of the protein. We assume with the help of these two results that the mutation is probably neutral and will cause no damage. |
{| |
{| |
||
− | | [[Image:mut_5_humdiv.png|thumb|450px|HumDiv prediction]] |
+ | | [[Image:mut_5_humdiv.png|thumb|450px|Figure 9: HumDiv prediction]] |
− | | [[Image:mut_5_humvar.png|thumb|450px|HumVar prediction]] |
+ | | [[Image:mut_5_humvar.png|thumb|450px|Figure 10: HumVar prediction]] |
|} |
|} |
||
Line 223: | Line 222: | ||
=== Mapping onto Crystal Structure === |
=== Mapping onto Crystal Structure === |
||
− | [[Image:mut248_active.png|thumb|center|400px|Visualization of the mutation and important functional sites |
+ | [[Image:mut248_active.png|thumb|center|400px|Figure 11: Visualization of the mutation and important functional sites<br>Color declaration: <br> |
+ | * <font color=red>red</font>: position of mutation<br> |
||
+ | * <font color=green>green</font>: position of active side<br> |
||
+ | * <font color=yellow>yellow</font>: position of glycolysation<br> |
||
+ | * <font color=cyan>cyan</font>: position of Cysteine]] |
||
− | Color declaration: |
||
− | * <font color=red>red</font>: position of mutation |
||
− | * <font color=green>green</font>: position of active side |
||
− | * <font color=yellow>yellow</font>: position of glycolysation |
||
− | * <font color=cyan>cyan</font>: position of Cystein |
||
+ | First of all, we colored the important residues and also the mutated residue in the crystal structure, to see if the mutation is near of far away from the functional residues. As you can see on Figure 11, the mutation is located within a helix but far away from the functional residues. The mutation is located at the end or begin of a helix. Therefore, of course it is possible that a mutation of this residue causes changes in the structure of the protein. But it is at the side of a helix and therefore the damage should not be that dramatically. |
||
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
Line 236: | Line 235: | ||
=== SCWRL Prediction === |
=== SCWRL Prediction === |
||
+ | |||
+ | Next, we analysed the mutation in more detail. Therefore we looked directly to the structure and orientation of the residue. |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
− | |picture original |
+ | |picture original amino acid |
− | |picture mutated |
+ | |picture mutated amino acid |
|combined picture |
|combined picture |
||
|- |
|- |
||
− | |[[Image:mut248_org_aa.png|thumb|150px|Amino acid Leucine]] |
+ | |[[Image:mut248_org_aa.png|thumb|150px|Figure 12: Amino acid Leucine]] |
− | |[[Image:mut248_aa.png|thumb|150px|Amino acid Phenylalanine]] |
+ | |[[Image:mut248_aa.png|thumb|150px|Figure 13: Amino acid Phenylalanine]] |
− | |[[Image:mut248_both.png|thumb|150px|Picture which visualize the mutation]] |
+ | |[[Image:mut248_both.png|thumb|150px|Figure 14: Picture which visualize the mutation]] |
|- |
|- |
||
|} |
|} |
||
+ | |||
+ | On the pictures we can see, that the amino acids themselves are very different (Figure 12 and Figure 13), but also the location and orientation of the amino acids differ dramatically. Furthermore the two amino acids have a totally different structure (Figure 14) and it is very likely that the properties of these two amino acids differ extremely. |
||
+ | Therefore it is very likely, that this mutation destroys the helix structure. |
||
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
Line 252: | Line 256: | ||
=== FoldX Energy Comparison === |
=== FoldX Energy Comparison === |
||
+ | |||
+ | One important point in the analysis of mutations is to look at the energy of the protein with the original and with the mutated amino acid. Often the energy increases dramatically with the mutated amino acid. This means, that the protein becomes very unstable and therefore, it is often possible that the protein can not bind its ligands any longer. Otherwise, it is also possible, that the protein with the mutated amino acid has a lower energy than the original protein. This means, that the protein is too rigid and loses its flexibility. Than it is also possible, that the protein can not bind the ligands any longer. |
||
+ | |||
+ | Therefore, we compared the energy of our protein with different methods. Here we want to present the result of FoldX. |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
Line 263: | Line 271: | ||
|- |
|- |
||
|} |
|} |
||
+ | |||
+ | |||
+ | The energy of the mutated structure is a little bit higher than the energy of the original structure, but the difference is not that much. |
||
+ | To have the possibility to compare the energy values with values from other programs, we calculated a ratio between these two values. |
||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |Ratio of the original structure |
||
+ | |Ratio of the mutated structure |
||
+ | |Difference |
||
+ | |- |
||
+ | |100 |
||
+ | |99.85 |
||
+ | |0.15 |
||
+ | |- |
||
+ | |} |
||
+ | |||
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
Line 269: | Line 293: | ||
=== Minimise Energy Comparison === |
=== Minimise Energy Comparison === |
||
+ | Next we use the minimise energy tool to compare the energy values of the two different structures. |
||
+ | |||
+ | ''' Comparing Energy: ''' |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
Line 279: | Line 306: | ||
|} |
|} |
||
+ | As before, in this case the mutated structure has a higher energy than the original structure. But again, the difference between the two energy values is not very high. To have the possibility to compare the energy values of the different tools, we calculate again the ratio between these two values. |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
− | | |
+ | |Ratio of the original structure |
− | | |
+ | |Ratio of the mutated structure |
+ | |Difference |
||
+ | |- |
||
+ | |100 |
||
+ | |99.96 |
||
+ | |0.04 |
||
+ | |- |
||
+ | |} |
||
+ | |||
+ | As before, the ratio between these two energy values is very low. |
||
+ | |||
+ | ''' Comparing Structure: ''' |
||
+ | |||
+ | This tool also gives as output a PDB file with the position of the original and the mutated amino acid. |
||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |picture original amino acid |
||
+ | |picture mutated amino acid |
||
|combined picture |
|combined picture |
||
|- |
|- |
||
− | |[[Image:mmut248_org.png|thumb|150px|Amino acid Leucine]] |
+ | |[[Image:mmut248_org.png|thumb|150px|Figure 15: Amino acid Leucine]] |
− | |[[Image:mmut248_mut.png|thumb|150px|Amino acid Phenylalanine]] |
+ | |[[Image:mmut248_mut.png|thumb|150px|Figure 16: Amino acid Phenylalanine]] |
− | |[[Image:mmut248_both.png|thumb|150px|Picture which visualize the mutation]] |
+ | |[[Image:mmut248_both.png|thumb|150px|Figure 17: Picture which visualize the mutation]] |
|- |
|- |
||
|} |
|} |
||
+ | |||
+ | Again we can see that these two amino acids (Figure 15, Figure 16 and Figure 17) differ extremely in location and orientation within the protein. The result is the same as we got wit SCWRL (Figure 12, Figure 13 and Figure 14). |
||
+ | |||
+ | ''' Visualization of H-bonds and Clashes: ''' |
||
+ | |||
+ | To get more insight in the effects of the mutated amino acid on the structure, we also analysed the H-bonds and clashes of the Asparagine residue. |
||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |H-bonds of the original structure |
||
+ | |H-bonds near the mutation |
||
+ | |Clashes of the mutation |
||
+ | |- |
||
+ | |[[Image:orig248.png|thumb|150px|Figure 18: H-bonds of the original amino acid (colored in magenta)]] |
||
+ | |[[Image:hbond248.png|thumb|150px|Figure 19: H-bonds of the mutated amino acids (colored in red)]] |
||
+ | |[[Image:clash248.png|thumb|150px|Figure 20: Possible clashes]] |
||
+ | |- |
||
+ | |} |
||
+ | |||
+ | If we compare the H-Bonds of the original (Figure 18) and the mutated amino acids (Figure 19), we can see that there are no H-Bonds between the residues and the rest of the protein. Therefore the damage of the protein is not caused by missing H-Bonds. |
||
+ | Furthermore, we can see that there are no clashes in the protein if the amino acid was mutated (Figure 20). |
||
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
− | --- |
+ | ---- |
=== Gromacs Energy Comparison === |
=== Gromacs Energy Comparison === |
||
+ | |||
+ | The last tool we used for the energy analysis was Gromacs. |
||
+ | |||
+ | ''' Comparing Energy: ''' |
||
+ | |||
+ | To analyse the energy values calculated by Gromacs, we used the AMBER99SB-ILDN force field. |
||
+ | |||
+ | Here are the values of the original structure: |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
Line 305: | Line 378: | ||
|- |
|- |
||
|Bond |
|Bond |
||
+ | |1091.57 |
||
− | | |
||
− | | |
+ | |270 |
+ | | -nan |
||
− | | |
||
+ | | -1622.75 |
||
− | | |
||
|- |
|- |
||
|Angle |
|Angle |
||
+ | |3326.81 |
||
− | | |
||
− | | |
+ | |62 |
+ | | -nan |
||
− | | |
||
+ | |404.076 |
||
− | | |
||
|- |
|- |
||
|Potential |
|Potential |
||
+ | | -61304.1 |
||
− | | |
||
− | | |
+ | |960 |
+ | | -nan |
||
− | | |
||
+ | | -6402.44 |
||
− | | |
||
|- |
|- |
||
|} |
|} |
||
+ | |||
− | <br><br> |
||
+ | Here you can see the values which gromacs calculated for the structure with the mutated amino acid: |
||
+ | |||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |Energy |
||
+ | |Average |
||
+ | |Err.Est. |
||
+ | |RMSD |
||
+ | |Tot-Drift |
||
+ | |- |
||
+ | |Bond |
||
+ | |920.202 |
||
+ | |180 |
||
+ | |2165.56 |
||
+ | | -1059.94 |
||
+ | |- |
||
+ | |Angle |
||
+ | |3202.18 |
||
+ | |61 |
||
+ | |233.689 |
||
+ | |406.391 |
||
+ | |- |
||
+ | |Potential |
||
+ | | -48802.5 |
||
+ | |830 |
||
+ | |3227.37 |
||
+ | | -5533.99 |
||
+ | |- |
||
+ | |} |
||
+ | |||
+ | One difference between gromacs and the other tools we used is, that gromacs also calculated the energy for the bonds and the angles. To compare the energies between the different tools we only consider the potential energy in our analysis, because the potential energy is the energy of the complete protein. Therefore, we calculated the ratio between the energies only for the potential energy. |
||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |Ratio of the original structure |
||
+ | |Ratio of the mutated structure |
||
+ | |Difference |
||
+ | |- |
||
+ | |100 |
||
+ | |79.60 |
||
+ | |20.40 |
||
+ | |- |
||
+ | |} |
||
+ | |||
+ | In comparison with the other two methods we used for the energy analysis the difference between the gromacs energy values is much higher. |
||
+ | |||
+ | ''' Comparing Structure: ''' |
||
{| border="1" style="text-align:center; border-spacing:0;" |
{| border="1" style="text-align:center; border-spacing:0;" |
||
− | |picture original |
+ | |picture original amino acid |
− | |picture mutated |
+ | |picture mutated amino acid |
|combined picture |
|combined picture |
||
|- |
|- |
||
− | |[[Image:gro_mut248_org.png|thumb|150px|Amino acid Leucine]] |
+ | |[[Image:gro_mut248_org.png|thumb|150px|Figure 21: Amino acid Leucine]] |
− | |[[Image:gro_mut248_mut.png|thumb|150px|Amino acid Phenylalanine]] |
+ | |[[Image:gro_mut248_mut.png|thumb|150px|Figure 22: Amino acid Phenylalanine]] |
− | |[[Image:gro_mut248_both.png|thumb|150px|Picture which visualize the mutation]] |
+ | |[[Image:gro_mut248_both.png|thumb|150px|Figure 23: Picture which visualize the mutation]] |
|- |
|- |
||
|} |
|} |
||
+ | As we saw before, the two amino acids differ extremely. |
||
+ | ''' Visualization of H-bonds and Clashes: ''' |
||
− | Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
+ | Gromacs also offers the user some pictures with the mutated amino acid. |
||
− | --- |
||
+ | |||
+ | {| border="1" style="text-align:center; border-spacing:0;" |
||
+ | |H-bonds of the original structure |
||
+ | |H-bonds near the mutation |
||
+ | |Clashes of the mutation |
||
+ | |- |
||
+ | |[[Image:orig248.png|thumb|150px|Figure 24: H-bonds of the original amino acid (colored in magenta)]] |
||
+ | |[[Image:gro_hbond248.png|thumb|150px|Figure 25: H-bonds of the mutated amino acid (colored in red)]] |
||
+ | |[[Image:gro_clash248.png|thumb|150px|Figure 26: Possible clashes]] |
||
+ | |- |
||
+ | |} |
||
+ | |||
+ | If we look at the pictures (Figure 24, Figure 25) we can see, that there are no H-bonds between the amino acids (the original and mutated) and the protein. Therefore, a damage of the protein can not be explained by missing or additional H-Bonds. |
||
+ | Furthermore, it is not possible to find any clashes between the mutated amino acid and the rest of the protein, what can be seen in Figure 26. |
||
+ | |||
+ | Back to [[http://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Structure-based_mutation_analysis_HEXA Structure-based mutation analysis]] |
||
+ | ---- |
Latest revision as of 22:38, 31 August 2011
Contents
General Information
SNP-id | rs61747114 |
Codon | 248 |
Mutation Codon | Leu -> Phe |
Mutation Triplet | CTT -> TTT |
Back to [Sequence-based mutation analysis]
Sequence-based Mutation Analysis
Pysicochemical Properties
First of all, we explored the amino acid properties and compared them for the original and the mutated amino acid. Therefore we created the possible effect that the mutation could have on the protein.
Leu | Phe | consequences |
aliphatic, hydrophobic, neutral | aromatic, hydrophobic, neutral | Leu is an aliphatic amino acid, whereas Phe is an aromatic amino acid. This means, that Phe has an aromatic ring in its structure. But both amino acids are relatively big and so it is possible, that the exchange of this amino acids does not change the structure of the protein that much. Therefore, we suggest it is possible, that the protein works. |
Back to [Sequence-based mutation analysis]
Visualization of the Mutation
In the next step, we created the visualization of the mutation with PyMol. Therefore we created a picture for the original amino acid (Figure 1), for the new mutated amino acid (Figure 2) and finally for both together in one picture whereas the mutation is white colored (Figure 3). The following pictures display that the mutated amino acid Phenylalanine looks very different to Leucine. Leucine forks at the end of the rest whereas Phenylalanine has an huge aromatic ring. The only thing that agrees is the orientation and the length of the residues before the fork or the aromatic ring. All in all the difference probably prevails the agreement and therefore the mutation will probably cause changes in protein structure and function.
picture original amino acid | picture mutated amino acid | combined picture |
Back to [Sequence-based mutation analysis]
Substitution Matrices Values
Afterwards, we looked at the values of the substitution matrices PAM1, PAM250 and BLOSSUM62. Therefore we looked detailed at the three values: the value for this amino acid substitution, the most frequent value for the substitution of the examined amino acid and the rarest substitution.
In this case, the substitution of Leucine to Phenylalanine has average values that are almost as far away from the rarest value as from the most frequent value for all three substitution matrices. Only BLOSOUM62 is a little bit closer to the most frequent substitution. Therefore, the values from both PAMs are not really significant and therefore we are not able to determine effects on the protein for these two matrices. Otherwise, according to BLOSOUM62 a mutation at this position will probably not cause structural changes which can affect functional changes.
PAM 1 | Pam 250 | BLOSOUM 62 | ||||||
value amino acid | most frequent substitution | rarest substitution | value amino acid | most frequent substitution | rarest substitution | value amino acid | most frequent substitution | rarest substitution |
13 | 45 (Met) | 0 (Asp, Cys) | 13 | 20 (Met) | 2 (Cys) | 0 | 2 (Ile, Met) | -4 (Asp, Gly) |
Back to [Sequence-based mutation analysis]
PSSM Analysis
Besides, we looked additional at the position specific scoring matrix (PSSM) for our sequence. In contrast to PAM and BLOSOUM, the PSSM contains a specific substitution rate for each position in the sequence. Therefore, the PSSM is more position specific than PAM or BLOSOUM. We extracted the substitution value for the underlying mutation, the value for the most frequent substitution and the rarest substitution.
In this case the substitution rate for Leucine to Phenylalanine at this position is very low and near the value for the rarest substitution. This means this substitution at this position is likely very uncommon which indicates that this substitution has bad effects as a consequence. Therefore, we concluded that this mutation will probably cause protein structure changes as well as functional changes.
PSSM | ||
value amino acid | most frequent substitution | rarest substitution |
-3 | 3 | -5 |
Back to [Sequence-based mutation analysis]
Conservation Analysis with Multiple Alignments
As a next step we created a multiple alignment which contains the HEXA sequence and 9 other mammalian homologous sequences from [UniProt]. Afterwards we looked at the position of the different mutations and looked at the conservation level at this position, which can be seen in Figure 4. The regarded mutation is presented by the second colored column. Here we can see, that all the other mammalians have also the amino acid Leucine on this position. Therefore, the amino acid at this position is highly conserved and a mutation there will cause probably huge structural and functional changes for the protein.
Back to [Sequence-based mutation analysis]
Secondary Structure Mutation Analysis
As a next step we compared the different results of the secondary structure prediction tools JPred and PsiPred. Afterwards we can examine in which secondary structure element and where therein the mutation takes place. This can give an overview of how drastic the mutation can be. In this case both tools agree and predict at the position of the mutation a coil. This has as result, that the mutation at this position would not destroy or split a secondary structure element. It will probably only changes the coil between two secondary structure elements, but this can sometimes also cause a change of the the following secondary structure. We think that a drastic change of the protein structure and its function is unlikely because the mutation does not affect a secondary structure element. The change of the coil will probably only take places between two secondary structure elements which will probably not change.
JPred: ...HHHHHHHHCCCEEEECCCCCHHHHHHHHHCCCCCCCCCCCCCCCCCCCCCCCCCC... PsiPred: ...HHHHHHHHCCCEEEECCCCCHHHHHHHHCCCCCCCCCCCCCCCCCCCCCCCCCCH...
Comparison with the real Structure:
Afterwards we also visualize the position of the mutation (red) in the real 3D-structure of PDB and compare it with the predicted secondary structure (Figure 5 and Figure 6). The visualization can therefore like above the predicted secondary structure display if the mutation is in a secondary structure element or in some other regions.
Here in this case the mutation position disagrees with the position of the predicted secondary structure and is within an alpha-helix. This means a mutation will probably destroy or split the alpha helix which affects drastic structural changes on the protein. We think that a structural change is very likely, because it is within a secondary structure element and will therefore cause extreme changes.
Back to [Sequence-based mutation analysis]
SNAP Prediction
Next, we looked at the result of the SNAP prediction. For this prediction we took the amino acid of the certain position and checked every possible amino acid mutation. Afterwards we extract the result for Phenylalanine which is the real mutation in this case. SNAP has as a result that the exchange from Leucine to Phenylalanine at this position is neutral with a relative high accuracy. This means that this certain mutation on this position cause very likely no structural and functional changes of the protein.
Substitution | Prediction | Reliability Index | Expected Accuracy |
F | Neutral | 3 | 78% |
A detailed list of all possible substitutions can be found [here]
Back to [Sequence-based mutation analysis]
SIFT Prediction
Next, we used SIFT Prediction which displays if a mutation is neutral or not. Therefore, it first shows a row which contains a score for the particular mutation position to a certain amino acid. The amino acid which are not tolerated at this position are colored red. Besides, it also constructs a table which lists the amino acids that are predicted as tolerated and not-tolerated (compare Figure 8).
In this case, there are seven substitutions that are tolerated: Tyrosine, Isoleucine, Aspartic acid, Valine, Phenylalanine, Methionine and Leucine. The substitution to Phenylalanine is tolerated at this position. This means that this mutation at this position is probably neutral and will not cause any structural and function changes of the protein.
SIFT Matrix:
Each entry contains the score at a particular position (row) for an amino acid substitution (column). Substitutions predicted to be intolerant are highlighted in red.
|
Back to [Sequence-based mutation analysis]
PolyPhen2 Prediction
Finally, we also regarded the PolyPhen2 prediction for this mutation. This prediction visualizes how strongly damaging the mutation probably will be. Therefore it gives the result for two possible cases: HumDiv and HumVar. HumDiv is the preferred model for evaluation rare allelse, dense mapping of regions identified by genome-wide association studies and analysis of neutral selection. In contrast, HumVar is the preferred model for diagnostic of Mendelian diseases which require distinguishing mutations with drastic effects from all remaining human variations including abundant mildly deleterious alleles. We decided to look at both possible models, which agreed the most cases.
In this case, the models disagree. The first model HumDic predicts that the mutation is benign (which could be seen on figure 9) whereas the second model HumVar predicts that the mutation is possibly damaging (compare Figure 10). When we look at the scores we can see that the difference is with 0.213 not so high, which displays that they are moving in the same direction. Looking at both together, there is probably no huge structural change but with a small probability there can possibly be a damage of the protein. We assume with the help of these two results that the mutation is probably neutral and will cause no damage.
Back to [Sequence-based mutation analysis]
Structure-based Mutation Analysis
Mapping onto Crystal Structure
First of all, we colored the important residues and also the mutated residue in the crystal structure, to see if the mutation is near of far away from the functional residues. As you can see on Figure 11, the mutation is located within a helix but far away from the functional residues. The mutation is located at the end or begin of a helix. Therefore, of course it is possible that a mutation of this residue causes changes in the structure of the protein. But it is at the side of a helix and therefore the damage should not be that dramatically.
Back to [Structure-based mutation analysis]
SCWRL Prediction
Next, we analysed the mutation in more detail. Therefore we looked directly to the structure and orientation of the residue.
picture original amino acid | picture mutated amino acid | combined picture |
On the pictures we can see, that the amino acids themselves are very different (Figure 12 and Figure 13), but also the location and orientation of the amino acids differ dramatically. Furthermore the two amino acids have a totally different structure (Figure 14) and it is very likely that the properties of these two amino acids differ extremely. Therefore it is very likely, that this mutation destroys the helix structure.
Back to [Structure-based mutation analysis]
FoldX Energy Comparison
One important point in the analysis of mutations is to look at the energy of the protein with the original and with the mutated amino acid. Often the energy increases dramatically with the mutated amino acid. This means, that the protein becomes very unstable and therefore, it is often possible that the protein can not bind its ligands any longer. Otherwise, it is also possible, that the protein with the mutated amino acid has a lower energy than the original protein. This means, that the protein is too rigid and loses its flexibility. Than it is also possible, that the protein can not bind the ligands any longer.
Therefore, we compared the energy of our protein with different methods. Here we want to present the result of FoldX.
Original total energy | Total energy for the mutated protein | Strongest energy changes within the mutated protein |
-154.17 | -153.78 | - |
The energy of the mutated structure is a little bit higher than the energy of the original structure, but the difference is not that much.
To have the possibility to compare the energy values with values from other programs, we calculated a ratio between these two values.
Ratio of the original structure | Ratio of the mutated structure | Difference |
100 | 99.85 | 0.15 |
Back to [Structure-based mutation analysis]
Minimise Energy Comparison
Next we use the minimise energy tool to compare the energy values of the two different structures.
Comparing Energy:
Original total energy | Total energy for the mutated protein |
-9610.467157 | -9606.588566 |
As before, in this case the mutated structure has a higher energy than the original structure. But again, the difference between the two energy values is not very high. To have the possibility to compare the energy values of the different tools, we calculate again the ratio between these two values.
Ratio of the original structure | Ratio of the mutated structure | Difference |
100 | 99.96 | 0.04 |
As before, the ratio between these two energy values is very low.
Comparing Structure:
This tool also gives as output a PDB file with the position of the original and the mutated amino acid.
picture original amino acid | picture mutated amino acid | combined picture |
Again we can see that these two amino acids (Figure 15, Figure 16 and Figure 17) differ extremely in location and orientation within the protein. The result is the same as we got wit SCWRL (Figure 12, Figure 13 and Figure 14).
Visualization of H-bonds and Clashes:
To get more insight in the effects of the mutated amino acid on the structure, we also analysed the H-bonds and clashes of the Asparagine residue.
H-bonds of the original structure | H-bonds near the mutation | Clashes of the mutation |
If we compare the H-Bonds of the original (Figure 18) and the mutated amino acids (Figure 19), we can see that there are no H-Bonds between the residues and the rest of the protein. Therefore the damage of the protein is not caused by missing H-Bonds. Furthermore, we can see that there are no clashes in the protein if the amino acid was mutated (Figure 20).
Back to [Structure-based mutation analysis]
Gromacs Energy Comparison
The last tool we used for the energy analysis was Gromacs.
Comparing Energy:
To analyse the energy values calculated by Gromacs, we used the AMBER99SB-ILDN force field.
Here are the values of the original structure:
Energy | Average | Err.Est. | RMSD | Tot-Drift |
Bond | 1091.57 | 270 | -nan | -1622.75 |
Angle | 3326.81 | 62 | -nan | 404.076 |
Potential | -61304.1 | 960 | -nan | -6402.44 |
Here you can see the values which gromacs calculated for the structure with the mutated amino acid:
Energy | Average | Err.Est. | RMSD | Tot-Drift |
Bond | 920.202 | 180 | 2165.56 | -1059.94 |
Angle | 3202.18 | 61 | 233.689 | 406.391 |
Potential | -48802.5 | 830 | 3227.37 | -5533.99 |
One difference between gromacs and the other tools we used is, that gromacs also calculated the energy for the bonds and the angles. To compare the energies between the different tools we only consider the potential energy in our analysis, because the potential energy is the energy of the complete protein. Therefore, we calculated the ratio between the energies only for the potential energy.
Ratio of the original structure | Ratio of the mutated structure | Difference |
100 | 79.60 | 20.40 |
In comparison with the other two methods we used for the energy analysis the difference between the gromacs energy values is much higher.
Comparing Structure:
picture original amino acid | picture mutated amino acid | combined picture |
As we saw before, the two amino acids differ extremely.
Visualization of H-bonds and Clashes:
Gromacs also offers the user some pictures with the mutated amino acid.
H-bonds of the original structure | H-bonds near the mutation | Clashes of the mutation |
If we look at the pictures (Figure 24, Figure 25) we can see, that there are no H-bonds between the amino acids (the original and mutated) and the protein. Therefore, a damage of the protein can not be explained by missing or additional H-Bonds. Furthermore, it is not possible to find any clashes between the mutated amino acid and the rest of the protein, what can be seen in Figure 26.
Back to [Structure-based mutation analysis]