SFESA: a web server for pairwise alignment refinement by secondary structure shifts
BMC Bioinformatics volume 16, Article number: 282 (2015)
Protein sequence alignment is essential for a variety of tasks such as homology modeling and active site prediction. Alignment errors remain the main cause of low-quality structure models. A bioinformatics tool to refine alignments is needed to make protein alignments more accurate.
We developed the SFESA web server to refine pairwise protein sequence alignments. Compared to the previous version of SFESA, which required a set of 3D coordinates for a protein, the new server will search a sequence database for the closest homolog with an available 3D structure to be used as a template. For each alignment block defined by secondary structure elements in the template, SFESA evaluates alignment variants generated by local shifts and selects the best-scoring alignment variant. A scoring function that combines the sequence score of profile-profile comparison and the structure score of template-derived contact energy is used for evaluation of alignments. PROMALS pairwise alignments refined by SFESA are more accurate than those produced by current advanced alignment methods such as HHpred and CNFpred. In addition, SFESA also improves alignments generated by other software.
SFESA is a web-based tool for alignment refinement, designed for researchers to compute, refine, and evaluate pairwise alignments with a combined sequence and structure scoring of alignment blocks. To our knowledge, the SFESA web server is the only tool that refines alignments by evaluating local shifts of secondary structure elements. The SFESA web server is available at http://prodata.swmed.edu/sfesa.
Homology modeling that constructs a structural model of a “query” protein based on its similarity to a homologous protein with known 3-dimensional structure (the “template”) remains the most reliable method of structure prediction. In most homology modeling methods, an essential step requires the input or construction of a pairwise sequence alignment between the query and the template, from which structurally equivalent residue pairs are deduced. Pairwise alignment is also the foundation for most multiple sequence alignment (MSA) methods. For example, the progressive method for MSA construction assembles a multiple sequence alignment by a series of pairwise alignments of sequences or pre-aligned groups .
Early methods of pairwise protein alignments apply dynamic programming algorithms that rely on general substitution matrices of amino acid residues and pre-defined gap penalties [2, 3]. Heuristic pairwise alignment tools such as BLAST  excel in speed and are suitable for sequence database searches. Numerical sequence profiles have been designed to incorporate information of homologous proteins to help aligning divergent sequences. PSI-BLAST  and HMMER  are examples of sequence-profile comparison methods that are generally more accurate than methods of sequence-sequence comparison. The subsequent development of profile-profile comparison methods [7–10] further enhanced alignment quality and the ability to detect homologous relationships. In addition to amino acid sequence profiles, predicted structural information, e.g., secondary structure and solvent accessibility, was also included in various alignment methods [11–13]. Three-dimensional structure information has been used in alignment construction methods in various ways, such as those based on structure-dependent profiles [14, 15] and a Monte Carlo-based alignment method that samples a set of moves of gapless alignment stretches and scores based on a template contact map .
Despite continuous method development in the alignment field, obtaining high-quality alignments for distantly related proteins remains a challenge. Alignment errors are still the main cause for the low quality of models built by homology. One common type of alignment error is the local misalignment, often by only a few residues, of secondary structure elements (α-helices and β-strands). Such errors often reflect the periodic nature of regular secondary structures. For example, many α-helices can be shifted by three or four residues while still maintaining a similar pattern of hydrophobic residues and polar residues. Therefore, one possible direction for refining an alignment lies in the generation of alignment variants by locally shifting secondary structure elements and evaluating the sequence and structure fitness of these alignment variants to determine which one is more likely to be correct.
Here we describe the SFESA web server, which refines pairwise protein alignments by evaluating alignment variants resulting from locally shifting secondary structure elements. The SFESA web server enables researchers to compute, refine, and evaluate pairwise alignments with a combined sequence and structure scoring of alignment blocks. The previous version of SFESA required the upload of a predefined template structure. In contrast, the new web server allows for a template to be specified by its PDB and chain identifiers. Furthermore, if no structure is provided, the SFESA server will search the database of sequences with experimentally determined 3D structures for the closest template, and this will then be used in the alignment refinement. The server facilitates further analysis of alignments at the level of secondary structure, providing detailed results of sequence and structure scores for local shifts of secondary structure elements. To our knowledge, the SFESA web server is the only online tool that refines alignments by evaluating local shifts of secondary structure elements.
Overview of the SFESA alignment refinement method and procedure
Recently we developed SFESA , a method that refines pairwise protein sequence alignment by evaluating alignment variants generated from local shifts of secondary structure elements. SFESA first delineates alignment blocks from a starting pairwise protein alignment. Each alignment block corresponds to a regular secondary structural element (α-helix or β-strand as delineated by PALSSE ) in the template and the corresponding aligned region in the query. For each alignment block, SFESA generates a set of alignment variants by locally shifting query residues relative to template residues. Then, both a profile-based sequence score and a contact-based structure score of the aligned residue pairs in the original alignment block and the alignment variants are calculated. We have shown that the best-scoring alignment variant has the highest probability of being correct, e.g., showing the best agreement with the structure-based alignment.
SFESA uses two local shifting strategies to generate alignment variants with different treatments of gaps in the original alignment block. In the first strategy, up to 8 alignment variants are generated by shifting query residues up to four positions left or right relative to the template while maintaining the gap pattern in the original alignment block. However, we observed that gaps rarely occur in the middle of secondary structure elements in structure-based alignments. Therefore, in the second strategy, SFESA preprocesses the gap pattern in the original alignment block by eliminating gaps in the middle of the secondary structure elements. To achieve this, residues of an alignment block in both the query and template are shifted all the way to the left or right while all gaps are placed on the opposite side. Two preprocessed alignment blocks are generated: one by shifting residues to the left and filling the right side with gaps and the other by shifting residues to the right and filling the left side with gaps. Each of these two alignment variants is then used as a starting point to generate 8 additional alignment variants by ±4 shifts while keeping the modified gap patterns. This procedure gives rise to up to 18 (1 + 8 + 1 + 8) unique alignment variants (for details, see ).
For the sequence score, we use the profile-profile COMPASS score . Sequence profiles are generated from PSI-BLAST multiple sequence alignments . For the structure score, we define residue contacts based on the structure of the template. A residue contact is defined as a residue pair within a distance cutoff. In the template of an alignment, the residue contacts can be identified using the known structure of the template. We then evaluate the contact energy of corresponding contact residue pairs in the query that are inferred from query-template alignment. For example, if residue i in the template makes contact with residues j, k, and m in the template structure (i.e., contact pairs are (i, j), (i, k), and (i, m)), and the corresponding aligned residues for i, j, k, and m in the query are i’ , j’ , k’ , and m’ , respectively, then the inferred contact pairs in the query are (i’ , j’), (i’ , k’), and (i’, m’). The structure score for the aligned residue pair i and i’ is CE(i’ , j’) + CE(i’ , k’) + CE(i’ , m’), reflecting the structural fitness of the inferred query contact residue pairs. Here, CE is a matrix of the contact energy for residue pairs. We used two contact energy matrices: one is derived by Miyazawa and Jernigan  with contacts defined as residue pairs with side chain centers less than 6.5 Å, and the other is developed by us to best discriminate correct alignment variants from incorrect alignment variants (for details, see ). Regarding our derived contact matrix, the cutoff for contact definition is 6.5 Å between any side chain atoms of two residues.
In practice, the SFESA method uses a two-filter strategy to compare the scores of the original alignment block and the alignment variants and determines whether the original alignment block should be kept or changed to one of the alignment variants. The first filter checks if there are any alignment variants with a higher combined score I (Scomb_I, a linear combination of sequence score and structure score) than the original alignment block. If none of the alignment variants has a Scomb_I higher than the original alignment block, SFESA rejects all the alignment variants and keeps the original alignment block. Otherwise, the alignment variant with the highest Scomb_I is selected and passed to the second filter. In the second filter, SFESA uses combined score II (Scomb_II, a linear combination of sequence score and structure score) or an SVM score (SSVM) to compare the selected alignment variant and the original alignment block. If the selected alignment variant still has a higher Scomb_II or SSVM, SFESA will accept this alignment variant. Otherwise, SFESA keeps the original alignment block. The weights of the sequence score and structure score in Scomb_I and Scomb_II are optimized separately. SSVM is a score reported by a support vector machine (SVM) that was trained to differentiate correct alignment variants from incorrect alignment variants by using a number of features including a COMPASS-based sequence score , a contact-based structure score, a solvent accessibility score and a secondary structure score (for details, see ). The solvent accessibility score is based on a three-by-three relative solvent accessibility substitution matrix derived from FAST  structural alignments of SCOP  domains. Similarly, the secondary structure score is based on a three-by-three secondary structure substitution matrix derived from FAST  structural alignments of SCOP  domains (for details, see ). The secondary structure is predicted by PSIPRED  for the query; the secondary structure information in DSSP  is used for the template. For each alignment block, starting from the N-terminus and proceeding to the C-terminus, SFESA decides whether to keep the original alignment block or to accept one of the alignment variants.
The SFESA web server
The SFESA web server is a tool for constructing, refining, and evaluating pairwise protein alignments (Fig. 1). The workflow of the server is shown in Fig. 1. Compared to the previously reported version of SFESA , in which a user must provide a structure for the template sequence, the updated server will search against our inhouse protein structure database to find the closest (to either sequence) homolog with available 3D structure to improve the alignment.
Users can input or upload sequences for the query and template either as a pairwise alignment or as two unaligned sequences in FASTA format. If two unaligned sequences are provided, the server uses PROMALS  to automatically construct a pairwise alignment. Input of a 3-dimensional structure (in pdb format) with high sequence similarity to the template is optional but recommended. A user can input a PDB identifier and a chain identifier, instead of a coordinate set, to directly use the structure from the RCSB PDB database . If no structure is provided for the template, the server uses BLAST  to automatically search for homologs for either the query or template in a database of representative spatial structures and selects the best hit as the homologous structure, used as a template for structure score calculation.
Four SFESA alignment refinement modes are available in the web server: SFESA (O) uses up to 8 variants generated by ±4 shifts that keep the gap patterns of the original alignment block and the Miyazawa-Jernigan (MJ)  contact matrix for structure score calculation; SFESA (O + G) uses up to 18 variants by considering gap shifts and the MJ contact matrix; and SFESA (O + G + M) (default) uses a newly derived contact matrix in addition to gap processing; SFESA (O + G + M + S) differs from SFESA (O + G + M) in that an SVM-derived score is used in the second filtering step instead of Scomb_II.
Several parameters are provided. One parameter is the sequence identity threshold between the template sequence and its homolog with a known 3D structure. SFESA refinement is applied only when the sequence identity between the template and its structure homolog is higher than the threshold (default = 0.5). Another parameter is the maximal number of residue positions to shift (default = 4, i.e., shifts are applied from −4 to +4 positions). Increasing this parameter generates more alignment variants, but also increases the probability that a wrong variant is accepted. The third parameter is the threshold for the fraction of non-gapped residue pairs above which an alignment block is used in the refinement process (default = 0.5). We also provide parameters for running and processing PSI-BLAST  results to generate the sequence profile used for the sequence score calculation, such as the number of iterations, the e-value inclusion cutoff, and a sequence identity cutoff to remove divergent hits.
The output page of the SFESA web server includes the starting alignment (the input alignment or in the case of the input of unaligned sequences, the automatically generated PROMALS alignment), the refined alignment, and the refinement details for each evaluated alignment block. Figure 2a, b, and c shows one example of the output page. The first part of the output page (Fig. 2a) contains the starting alignment and the refined alignment with colored alignment blocks. PSIPRED  predicts secondary structure elements of the query and secondary structure elements of the template are based on PALSSE  and DSSP . These predicted elements are shown above the query sequence and below the template sequence, respectively. Evaluated alignment blocks are depicted in red and orange for α-helices and blue and dark green for β-strands to distinguish them. In the SFESA-refined alignment, the modified alignment blocks are marked with underscores.
The second part of the output page (Fig. 2b) is a table summarizing the refinement results of the evaluated alignment blocks, numbered from the N-terminus to the C-terminus. Each row in the table provides the element start and end position numbers in the template, the element secondary structure type, the original alignment block, the shift result, and the refined alignment block. The shift result column shows Gap Mode and Shift Number. Gap Mode can be “Left” (gap pattern preprocessed by moving residues all the way to the left), “Right” (gap pattern preprocessed by moving residues all the way to the right), or “Original” (no gap preprocessing). Shift Number (in brackets) is the number of positions the residues in the query are shifted by, relative to the template. The “+” and “−” signs in Shift Number denote that the query residues in the alignment block are shifted towards the C-terminus or the N-terminus, respectively. If no alignment variant was accepted for an alignment block (i.e., the original refinement retained), “No shift” is shown in the shift result column and “-” is shown in the column of Refined Alignment Block. The third part of the output page contains tables with scoring details for the alignment variants. A table is provided for each alignment block evaluated by SFESA and presents each alignment variant and its sequence score, structure score, and combined scores I and II. Figure 2c provides an example for alignment block number 4. The residues in the original alignment block are colored blue and pink for the query and the template, respectively. The scoring details for alignment variants may help users manually evaluate and select alternative alignments.
Test on the SABmark benchmark
We tested SFESA on the SABmark benchmark database of alignments . The SABmark database has two alignment sets with different levels of sequence similarity and thus different levels of difficulty for alignment. The Twilight Zone set was created by selecting structurally similar domains at the SCOP  fold level and contains sequences with very low to low similarity. The Superfamilies set was created by selecting structurally similar domains at the SCOP superfamily level and contains sequences with low to intermediate similarity. We tested the ability of SFESA to refine alignments generated by several alignment methods, such as PROMALS , HHpred , and CNFpred . Here HHpred was used in the global alignment mode because its local alignment mode often results in short alignments and shows lower alignment accuracy than global alignments.
We used the reference-dependent Q-score as the assessment. The Q-score is the fraction of correctly aligned residue pairs in a test alignment among all aligned residue pairs in a reference alignment. In this paper, the range of Q-score values is from 0 to 100 (e.g. 100 means 100 % agreement with reference).
SFESA can improve the PROMALS Q-score from 46 to 48 for the Twilight Zone set and from 71 to 72 for the Superfamilies set (Table 1). PROMALS-based SFESA outperforms other advanced alignment methods, such as HHpred and CNFpred. In practice, SFESA (O + G + M) and SFESA (O + G + M + S) produced similar results that are on average better than SFESA (O) and SFESA (O + G). Furthermore, SFESA also improves alignments generated by other methods (Table 1), including HHpred and CNFpred.
An example of an alignment improved by the SFESA server
In the example shown in Fig. 2, the input consisted of two SCOP domains, d1ja1a3 (query) and d2piaa2 (template), and the 3D structure of the template. SFESA used PROMALS to obtain the starting alignment, which was refined to generate the refined alignment with the default option SFESA (O + G + M). Out of the seven alignment blocks evaluated by SFESA, five alignment blocks were kept without shifts and two alignment blocks were modified according to SFESA refinement scores (Fig. 2b). Both of these modified alignment blocks are in better agreement with the Dali structural alignment  of the query and the template compared to the original alignment blocks. We generated structure models for the query based on the starting alignment (Fig. 2d, left panel) and the refined alignment (Fig. 2d, right panel). Both models (in light grey and red ribbons) were superimposed upon the real structure of the query (in dark grey and green ribbons). The GDT-TS scores  for models generated from the starting alignment and the refined alignment are 57.7 and 67.1, respectively. The query secondary structure element in the fourth evaluated alignment block is highlighted in both structure superpositions (green for the real structure and red for the model). This element, misaligned by two residues in the starting alignment (Fig. 2d, left panel), has been corrected in the refined alignment (Fig. 2d, right panel). As a result, the RMSD for this secondary structure element between the model and the real structure improved from 5.3 Å for the model generated by the starting alignment to 2.0 Å for the model generated by the refined alignment.
Despite many significant research efforts, it is still challenging to correctly align weakly similar but homologous protein sequences. Alignment errors remain the main reason for the poor quality of homology-based models. Refining the alignments generated by automatic methods is a promising approach for increasing alignment quality. We found that secondary structure elements are often misaligned by only a few residues and that more accurate solutions can be identified within a limited set of local shifts of secondary structure elements. Therefore, we developed the SFESA method in order to refine alignments by evaluating the alignment variants generated by local shifts of template-defined secondary structures.
In the SFESA scoring system, both a profile-based sequence score and a novel contact-based structure score of the aligned residue pairs in the original alignment block and the alignment variants are calculated. Thus, an insufficient number of contacts can limit the quality of the alignment refinement. We found that structure scoring works well when there are sufficient contacts in the template as well as sufficient corresponding aligned residues in the query . However, if a secondary structure element is involved in too few contacts (e.g. exposed edge β-strands), these contacts are insufficient to define a complete structural environment. SFESA is less effective in these cases. This observation suggests that dedicated efforts on misaligned blocks with insufficient contacts are required to improve alignments further.
SFESA is a web-based tool to compute, refine, and evaluate pairwise alignments with a combined sequence and structure scoring of alignment blocks. Taking a pairwise alignment as input, the SFESA web server searches against an in-house database of protein spatial structures to find the closest homolog of either sequence. It then refines the pairwise alignment by combining the sequence profile similarity and residue-residue contact information that were obtained from the homolog with the structure. Finally, it facilitates further analysis of the alignment results at the level of secondary structure, providing details about scoring for all shifts of secondary structure elements.
Availability and requirements
Project name: SFESA, a pairwise alignment refinement tool.
Project home page: http://prodata.swmed.edu/sfesa.
Operating system(s): Platform independent.
Programming language: Perl.
Other requirements: Internet connection, a browser.
Any restrictions to use by non-academics: None.
Feng DF, Doolittle RF. Progressive sequence alignment as a prerequisite to correct phylogenetic trees. J Mol Evol. 1987;25(4):351–60.
Needleman SB, Wunsch CD. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970;48(3):443–53.
Smith TF, Waterman MS. Identification of common molecular subsequences. J Mol Biol. 1981;147(1):195–7.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. doi:10.1016/S0022-2836(05)80360-2.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14(9):755–63.
Sadreyev R, Grishin N. COMPASS: a tool for comparison of multiple protein alignments with assessment of statistical significance. J Mol Biol. 2003;326(1):317–36.
Yona G, Levitt M. Within the twilight zone: a sensitive profile-profile comparison tool based on information theory. J Mol Biol. 2002;315(5):1257–75. doi:10.1006/jmbi.2001.5293.
Rychlewski L, Jaroszewski L, Li W, Godzik A. Comparison of sequence profiles. Strategies for structural predictions using sequence information. Protein Sci. 2000;9(2):232–41. doi:10.1110/ps.9.2.232.
Gniewek P, Kolinski A, Gront D. Optimization of profile-to-profile alignment parameters for one-dimensional threading. J Comput Biol. 2012;19(7):879–86. doi:10.1089/cmb.2011.0307.
Pei J, Grishin NV. PROMALS: towards accurate multiple sequence alignments of distantly related proteins. Bioinformatics. 2007;23(7):802–8. doi:10.1093/bioinformatics/btm017.
Soding J. Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005;21(7):951–60. doi:10.1093/bioinformatics/bti125.
Ma J, Peng J, Wang S, Xu J. A conditional neural fields model for protein threading. Bioinformatics. 2012;28(12):i59–66. doi:10.1093/bioinformatics/bts213.
Bowie JU, Luthy R, Eisenberg D. A method to identify protein sequences that fold into a known three-dimensional structure. Science. 1991;253(5016):164–70.
Kelley LA, MacCallum RM, Sternberg MJ. Enhanced genome annotation using structural profiles in the program 3D-PSSM. J Mol Biol. 2000;299(2):499–520. doi:10.1006/jmbi.2000.3741.
Gniewek P, Kolinski A, Kloczkowski A, Gront D. BioShell-Threading: versatile Monte Carlo package for protein 3D threading. BMC Bioinformatics. 2014;15:22. doi:10.1186/1471-2105-15-22.
Tong J, Pei J, Otwinowski Z, Grishin NV. Refinement by shifting secondary structure elements improves sequence alignments. Proteins. 2015;83(3):411–27. doi:10.1002/prot.24746.
Majumdar I, Krishna SS, Grishin NV. PALSSE: a program to delineate linear secondary structural elements from protein structures. BMC Bioinformatics. 2005;6:202. doi:10.1186/1471-2105-6-202.
Miyazawa S, Jernigan RL. An empirical energy potential with a reference state for protein fold and sequence recognition. Proteins. 1999;36(3):357–69. doi:10.1002/(SICI)1097-0134(19990815).
Zhu J, Weng Z. FAST: a novel protein structure alignment algorithm. Proteins. 2005;58(3):618–27. doi:10.1002/prot.20331.
Andreeva A, Howorth D, Chandonia JM, Brenner SE, Hubbard TJ, Chothia C, et al. Data growth and its impact on the SCOP database: new developments. Nucleic Acids Res. 2008;36(Database issue):D419–25.
Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292(2):195–202. doi:10.1006/jmbi.1999.3091.
Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983;22(12):2577–637. doi:10.1002/bip.360221211.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42.
Van Walle I, Lasters I, Wyns L. SABmark--a benchmark for sequence alignment that covers the entire known fold space. Bioinformatics. 2005;21(7):1267–8. doi:10.1093/bioinformatics/bth493.
Holm L, Sander C. Mapping the protein universe. Science. 1996;273(5275):595–603.
Zemla A. LGA: A method for finding 3D similarities in protein structures. Nucleic Acids Res. 2003;31(13):3370–4.
We thank Raquel Bromberg for critical reading of the manuscript. We also would like to thank Ming Tang for providing technical support with setting up the server.
This work was supported by the National Institutes of Health (GM094575 to NVG) and the Welch Foundation (I-1505 to NVG).
The authors declare that they have no competing interests.
JT, JP, and NVG designed and developed the algorithms and workflow. JT and JP designed the web server. JT performed the analysis and drafted the manuscript. NVG conceptualized and oversaw the development of the tool. All authors contributed in discussions and approved the final draft of the manuscript.
About this article
Cite this article
Tong, J., Pei, J. & Grishin, N.V. SFESA: a web server for pairwise alignment refinement by secondary structure shifts. BMC Bioinformatics 16, 282 (2015). https://doi.org/10.1186/s12859-015-0711-0