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. And with both together we would classify it as damaging, which is correct.
Mutation 7
We predicted the seventh mutation as non-neutral in sequence-based mutation analysis, which is correct, and as neutral in structure-based mutation analysis. The comparisons of the substitution matrices etc. indicated a neutral mutation but all prediction tools in sequence-based analysis classified the mutation as damaging, which was the reason why we decided to predict it as non-neutral. But it was hard to choose the right prediction. In structure-based mutation analysis all results lead to the assumption that the mutation is harmless. As it is the most common mutation in Gaucher Disease it is really interesting that we have so many problems to classify it correctly. Also with both methods we would tend to classify it as neutral. To get a better insight into the mutation's effect, we made a Molecular Dynamics analysis. Maybe it helps to understand the damaging effect.
Mutation 8
We predicted mutation eight as damaging in the sequence-based mutation analysis and in the structure-based mutation analysis which is correct. It is listed in the HGMD and associated to Gaucher Disease 2. In both cases it was easy to decide to classify the mutation as damaging.
Mutation 9
The ninth mutation is damaging, which we also predicted in the sequence-based mutation analysis. But in structure-based mutation analysis we classified it as neutral. In sequence-based analysis we had many results that indicated a damaging mutation. The amino acid properties, the PSSM and also the prediction tools gave us enough reason to classify it like that. In structure-based analysis it was hard to decide. There were some reasons for classifying it as damaging and some to classify it as harmless. It seems that we did not judges the importance of the different results correctly. On the other hand there are so many signs in sequence-based analysis, that this should be enough to classify the mutation as damaging.
Mutation 10
For mutation ten we also have different results. In sequence-based analysis we predicted it as damaging whereas in structure-based analysis as harmless, which is correct. In sequence-based analysis all results lead us to the putative clear prediction, that the mutation is damaging. But that is wrong. In structure-based analysis we found only a little higher energy but no other signs for a damaging mutation. So we classified it as neutral. All in all we also would classify it as damaging with the results of the structure-based analysis because the sequence-based analysis seems to be clear. But that's wrong. For this mutation we also made a Molecular Dynamics anaysis as we have such different results in sequence-based and structure-based analysis.
Summary
sequence-based mutation analysis | structure-based mutation analysis | ||
correctly predicted | wrong predicted | correctly predicted | wrong predicted |
7 | 3 | 6 | 4 |