I am interested on how profile hidden markov model is used for protein analysis in bioinformatics. Can you guys explain why phmm is more preferable method comparing to simple MSA(Multiple Sequence Alignment) scoring or simple MSA probability scoring? For example we have three MSA sequences as based model and a target to be scored.
A B C B C A
A B B C C A
A B C B A A
A B B B C A
By MSA scoring we can count score for target sequence comparing to a profile MSA. For example above target sequence protein column will counted as match(A), match(B), mismatch(B-majorityMSA(C)), match(B-majorityMSA(B)), match(C-majorityMSA(C)), match(A). Just give score for match and mismatch and we can get total score for a target sequence comparing to a profile MSA sequences.
By MSA probability scoring we can count probability of each protein in each column so for example above
P(target_seq similar to MSA profile) = P(A in MSA sequences)*P(B in MSA sequences)*.... *P(A in MSA sequences) = 1*1*(1/3)*(2/3)*(2/3)*1 = 4/27
Isn't MSA scoring like above examples powerful enough?
Thanks a lot !!