Biomedical Engineering Reference
In-Depth Information
H -2 0 1 -1 -3 0 0 -2 8
I -1 -3 -3 -3 -1 -3 -3 -4 -3 4
L -1 -2 -3 -4 -1 -2 -3 -4 -3 2 4
K -1 2 0 -1 -3 1 1 -2 -1 -3 -2 5
M -1 -1 -2 -3 -1 0 -2 -3 -2 1 2 -1 5
F -2 -3 -3 -3 -2 -3 -3 -3 -1 0 0 -3 0 6
P -1 -2 -2 -1 -3 -1 -1 -2 -2 -3 -3 -1 -2 -4 7
S 1 -1 1 0 -1 0 0 0 -1 -2 -2 0 -1 -2 -1 4
T 0 -1 0 -1 -1 -1 -1 -2 -2 -1 -1 -1 -1 -2 -1 1 5
W -3 -3 -4 -4 -2 -2 -3 -2 -2 -3 -2 -3 -1 1 -4 -3 -2 11
Y -2 -2 -2 -3 -2 -1 -2 -3 2 -1 -1 -2 -1 3 -3 -2 -2 2 7
V 0 -3 -3 -3 -1 -2 -2 -3 -3 3 1 -2 1 -1 -2 -2 0 -3 -1 4
B -2 -1 3 4 -3 0 1 -1 0 -3 -4 0 -3 -3 -2 0 -1 -4 -3 -3 4
Z -1 0 0 1 -3 3 4 -2 0 -3 -3 1 -1 -3 -1 0 -1 -3 -2 -2 1 4
X 0 -1 -1 -1 -2 -1 -1 -1 -1 -1 -1 -1 -1 -1 -2 0 0 -2 -1 -1 -1 -1 -1
As noted in the earlier discussion of gap penalties, arbitrarily selecting opening and extension costs
so that the output "looks nice" from a mathematical perspective likely has no relevance to the actual
biology of the protein under study. It's commonly assumed that a better approach is to assign gap
extension and opening costs relative to the substitution matrix used for a given protein. If the gap
penalty figures are too high relative to the matrix scores, the gap penalty figures will override the
matrix scores, and gaps will never appear in the sequence alignment. Conversely, if gap penalty
figures are too low relative to the matrix scores, gaps will be used wherever possible in order to align
the sequences. That is, simply because a substitution matrix is used doesn't guarantee biologically
relevant results. The matrices and related calculations must be used appropriately, and in
consideration of the underlying biology.
