TY - JOUR
T1 - Revisiting gap locations in amino acid sequence alignments and a proposal for a method to improve them by introducing solvent accessibility
AU - Hijikata, Atsushi
AU - Yura, Kei
AU - Noguti, Tosiyuki
AU - Go, Mitiko
PY - 2011/6
Y1 - 2011/6
N2 - In comparative modeling, the quality of amino acid sequence alignment still constitutes a major bottleneck in the generation of high quality models of protein three-dimensional (3D) structures. Substantial efforts have been made to improve alignment quality by revising the substitution matrix, introducing multiple sequences, replacing dynamic programming with hidden Markov models, and incorporating 3D structure information. Improvements in the gap penalty have not been a major focus, however, following the development of the affine gap penalty and of the secondary structure dependent gap penalty. We revisited the correlation between protein 3D structure and gap location in a large protein 3D structure data set, and found that the frequency of gap locations approximated to an exponential function of the solvent accessibility of the inserted residues. The nonlinearity of the gap frequency as a function of accessibility corresponded well to the relationship between residue mutation pattern and residue accessibility. By introducing this relationship into the gap penalty calculation for pairwise alignment between template and target amino acid sequences, we were able to obtain a sequence alignment much closer to the structural alignment. The quality of the alignments was substantially improved on a pair of sequences with identity in the "twilight zone" between 20 and 40%. The relocation of gaps by our new method made a significant improvement in comparative modeling, exemplified here by the Bacillus subtilis yitF protein. The method was implemented in a computer program, ALAdeGAP (ALignment with Accessibility dependent GAp Penalty), which is available at.
AB - In comparative modeling, the quality of amino acid sequence alignment still constitutes a major bottleneck in the generation of high quality models of protein three-dimensional (3D) structures. Substantial efforts have been made to improve alignment quality by revising the substitution matrix, introducing multiple sequences, replacing dynamic programming with hidden Markov models, and incorporating 3D structure information. Improvements in the gap penalty have not been a major focus, however, following the development of the affine gap penalty and of the secondary structure dependent gap penalty. We revisited the correlation between protein 3D structure and gap location in a large protein 3D structure data set, and found that the frequency of gap locations approximated to an exponential function of the solvent accessibility of the inserted residues. The nonlinearity of the gap frequency as a function of accessibility corresponded well to the relationship between residue mutation pattern and residue accessibility. By introducing this relationship into the gap penalty calculation for pairwise alignment between template and target amino acid sequences, we were able to obtain a sequence alignment much closer to the structural alignment. The quality of the alignments was substantially improved on a pair of sequences with identity in the "twilight zone" between 20 and 40%. The relocation of gaps by our new method made a significant improvement in comparative modeling, exemplified here by the Bacillus subtilis yitF protein. The method was implemented in a computer program, ALAdeGAP (ALignment with Accessibility dependent GAp Penalty), which is available at.
KW - ALAdeGAP
KW - Amino acid sequence alignment
KW - Comparative modeling
KW - Position dependent gap penalty
KW - Solvent accessibility
UR - http://www.scopus.com/inward/record.url?scp=79955705563&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=79955705563&partnerID=8YFLogxK
U2 - 10.1002/prot.23011
DO - 10.1002/prot.23011
M3 - Article
C2 - 21465562
AN - SCOPUS:79955705563
VL - 79
SP - 1868
EP - 1877
JO - Proteins: Structure, Function and Genetics
JF - Proteins: Structure, Function and Genetics
SN - 0887-3585
IS - 6
ER -