02070nas a2200205   4500008004100000245006600041210006500107300001100172490000700183520128000190653004201470653002501512653010401537653004001641100002401681700002401705700001601729700001301745856010601758       2006                        eng d00aVariable gap penalty for protein sequence-structure alignment0 aVariable gap penalty for protein sequencestructure alignment  a129-330 v193 aThe penalty for inserting gaps into an alignment between two protein sequences is a major determinant of the alignment accuracy. Here, we present an algorithm for finding a globally optimal alignment by dynamic programming that can use a variable gap penalty (VGP) function of any form. We also describe a specific function that depends on the structural context of an insertion or deletion. It penalizes gaps that are introduced within regions of regular secondary structure, buried regions, straight segments and also between two spatially distant residues. The parameters of the penalty function were optimized on a set of 240 sequence pairs of known structure, spanning the sequence identity range of 20-40%. We then tested the algorithm on another set of 238 sequence pairs of known structures. The use of the VGP function increases the number of correctly aligned residues from 81.0 to 84.5% in comparison with the optimized affine gap penalty function; this difference is statistically significant according to Student’s t-test. We estimate that the new algorithm allows us to produce comparative models with an additional approximately 7 million accurately modeled residues in the approximately 1.1 million proteins that are detectably related to a known structure.10aAlgorithms Amino Acid Sequence Models10aAmino Acid *Software10aMolecular Molecular Sequence Data Proteins/*chemistry Sequence Alignment/*methods Sequence Analysis10aProtein/*methods *Sequence Homology1 aMadhusudhan, M., S.1 aMarti-Renom, M., A.1 aSanchez, R.1 aSali, A.  uhttp://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=16423846