Sequence and structure based mutation analysis of GBA
Contents
Introduction
In this section we want to combine the results of sequence- and structure-based mutation analysis. Therefore we use the results of task 6 and task 7.
Sequence-based mutation analysis
The following table summarizes our results of the sequence-based mutation analysis.
Mutation | Amino-Acid Properties | Substitution Matrices | PSSM | Conservation | Secondary Structure | SNAP | SIFT | PolyPhen-2 | |||
BLOSUM62 | PAM1 | PAM250 | HumDiv | HumVar | |||||||
1 | non-neutral | neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | neutral | neutral | non-neutral | non-neutral |
2 | non-neutral | neutral | neutral | neutral | non-neutral | non-neutral | neutral | non-neutral | neutral | non-neutral | non-neutral |
3 | neutral | neutral | neutral | neutral | non-neutral | neutral | neutral | neutral | neutral | neutral | neutral |
4 | non-neutral | neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral |
5 | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | neutral | non-neutral | non-neutral | non-neutral | non-neutral |
6 | neutral | neutral | neutral | neutral | neutral | non-neutral | neutral | neutral | neutral | neutral | neutral |
7 | neutral | neutral | neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | neutral |
8 | non-neutral | neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral |
9 | non-neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | neutral | non-neutral | non-neutral |
10 | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral | non-neutral |
Structure-based mutation analysis
The following table summarizes our results of the structure-based mutation analysis.
Mutation | SCWRL | Minimise | Gromacs | FoldX Energy | ||||||
Polar Interactions |
Clashes Holes |
Energy | Polar Interactions |
Clashes Holes |
Energy | Polar Interactions |
Clashes Holes |
Energy | ||
1 | neutral | neutral | neutral | neutral | neutral | neutral | neutral | non-neutral | neutral | neutral |
2 | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral |
3 | non-neutral | non-neutral | neutral | non-neutral | non-neutral | neutral | non-neutral | neutral | neutral | neutral |
4 | non-neutral | neutral | neutral | non-neutral | neutral | neutral | non-neutral | non-neutral | non-neutral | neutral |
5 | neutral | neutral | non-neutral | neutral | neutral | non-neutral | neutral | neutral | non-neutral | non-neutral |
6 | non-neutral | non-neutral | neutral | non-neutral | non-neutral | neutral | non-neutral | non-neutral | non-neutral | neutral |
7 | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral | neutral |
8 | non-neutral | non-neutral | neutral | non-neutral | neutral | non-neutral | non-neutral | neutral | non-neutral | neutral |
9 | neutral | non-neutral | neutral | neutral | neutral | non-neutral | neutral | non-neutral | non-neutral | neutral |
10 | neutral | non-neutral | neutral | neutral | non-neutral | neutral | neutral | neutral | non-neutral | neutral |
Predictions of sequence-based and structure-based mutation analysis
The following table shows the predictions we made with sequence-based and structure-based analysis. It is also shown, if we agree in both methods and if the mutation is listed in HGMD and therefore damaging.
Mutation | Sequence-based mutation analysis | Structure-based mutation analysis | in HGMD? | prediction |
1 | neutral | neutral | yes | wrong |
2 | non-neutral | neutral | yes | partly correct |
3 | neutral | neutral | no | correct |
4 | non-neutral | non-neutral | yes | correct |
5 | non-neutral | non-neutral | yes | correct |
6 | neutral | non-neutral | yes | partly correct |
7 | non-neutral | neutral | yes | partly correct |
8 | non-neutral | non-neutral | yes | correct |
9 | non-neutral | neutral | yes | partly correct |
10 | non-neutral | neutral | no | partly correct |
Discussion
Mutation 1
The first mutation is the only one we predicted totally wrong. In sequence- and in structure-based analysis we predicted it as neutral, but as it is listed in HGMD it is damaging. In sequence-based analysis the amino-acid properties, the PSSM, the conservation, the secondary structure and the prediction of Polyphen-2 indicated that the mutation would be damaging. So it was not easy to decide if we classify the mutation as neutral or damaging. But the affected amino acid is at the exterior of the protein and there were also many results that lead us to the prediction that the mutation is neutral. In structure-based mutation analysis almost all results lead us to the prediction that the mutation is harmless. Only one Gromacs showed a different surface. It is interesting, that we did not find any significant changes in structure-based mutation analysis as we would expect after the results in sequence-based mutation analysis. All in all there were too little signs for a damaging mutation. For this mutation our prediction was totally wrong. We failed in both methods. It would be interesting to use more methods to see if we would be able find the reason for the damaging effect.
Mutation 2
We predicted the second mutation partly correct. It is listed in HGMD and so damaging, as we also conluded in sequence-based analysis. In structure-based analysis we classified the mutation as harmless. In sequence-based analysis it was hard to classify. We had different results and so were not sure which ones are the most important. But mainly because of the change from an acidic to a neutral amino acid we predicted the mutation as damaging. In structure-based analysis all results indicated a neutral substitution. So we classified the mutation as harmless, which was wrong. The effect must be directly at the amino acid with no structural changes. Maybe the binding differs somehow or the loss of the acidic character is damaging.
Mutation 3
We predicted the third mutation correct in sequence- and structure-based mutation analysis. Although it was hard to decide in structure-based analysis we predicted it correctly as neutral.
Mutation 4
Also the fourth mutation was predicted correclty. In both tasks we classified it as damaging, which is correct. In sequence- and structure-based mutation analysis it was easy to decide.
Mutation 5
The fifth mutation we also predicted correctly as damaging. In sequence-based analysis it was absolutely clear to predict it as damaging and in structure-based analysis we also were sure to classify it correctly because of the high difference in energy comparison.
Mutation 6
For the sixth mutation we had to different results. In sequence-based mutation analysis we predicted it as neutral, which was wrong, in structure-based mutation analysis we predicted it as damaging, which is correct as it is listed in HGMD. In sequence-based analysis only the conservation indicated a damaging mutation but all the other methods lead us to the conclusion, that the mutation is neutral. But in structure-based analysis the polar interactions changed, the surface changed and also the energy was too high. So we decided to classify it as damaging, which is right. This is a good example to see that you have consider all possible effects of a mutation. Although the sequence-based analysis does not show a hint, that the mutation could be damaging it is very clear if you consider structure-based analysis.
Mutation 7
MD
Mutation 8
Mutation 9
Mutation 10
MD