Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Highly Accessed Research article

The meaning of alignment: lessons from structural diversity

Walter Pirovano, K Anton Feenstra and Jaap Heringa*

Author Affiliations

Centre for Integrative Bioinformatics VU (IBIVU), VU University Amsterdam, De Boelelaan 1081A, 1081HV Amsterdam, the Netherlands

For all author emails, please log on.

BMC Bioinformatics 2008, 9:556  doi:10.1186/1471-2105-9-556


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/9/556


Received:20 August 2008
Accepted:23 December 2008
Published:23 December 2008

© 2008 Pirovano et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Protein structural alignment provides a fundamental basis for deriving principles of functional and evolutionary relationships. It is routinely used for structural classification and functional characterization of proteins and for the construction of sequence alignment benchmarks. However, the available techniques do not fully consider the implications of protein structural diversity and typically generate a single alignment between sequences.

Results

We have taken alternative protein crystal structures and generated simulation snapshots to explicitly investigate the impact of structural changes on the alignments. We show that structural diversity has a significant effect on structural alignment. Moreover, we observe alignment inconsistencies even for modest spatial divergence, implying that the biological interpretation of alignments is less straightforward than commonly assumed. A salient example is the GroES 'mobile loop' where sub-Ångstrom variations give rise to contradictory sequence alignments.

Conclusion

A comprehensive treatment of ambiguous alignment regions is crucial for further development of structural alignment applications and for the representation of alignments in general. For this purpose we have developed an on-line database containing our data and new ways of visualizing alignment inconsistencies, which can be found at http://www.ibi.vu.nl/databases/stralivari webcite.

Background

Sequence comparison has become a major tool for biological research in the post-genomic era, forming the basis for functional annotation, classification, and analysis of evolutionary relationships. At the residue level, however, the relation between sequence, structure and function can often be obscure, and examples abound of proteins with a clear functional and homologous relationship but sharing negligible similarity at the sequence level.

Structural alignment therefore is the method of choice for reliable homology assessment and derived features like functional classification and phylogeny. This importance is reflected in the number of tools available for structural alignment, such as DALI [1], SSAP [2], STRUCTAL [3], MAMMOTH [4], CE [5] and COMPARER [6] (for recent reviews on the topic, see Kolodny et al. [7] and Mayr et al. [8]). Databases for functional classification such as CATH [9], FSSP [10] and PASS2 [11] each derive directly from the use of one or more of these methods, whereas for SCOP expert input in the structural classification is deemed critical [12]. Structural alignments are also routinely used for benchmarking sequence alignment methods. A number of databases have been developed for this purpose, among which BAliBASE [13], HOMSTRAD [14] and SABmark [15] are widely used. These databases often rely on expert knowledge and include a notion of 'core blocks', i.e. where alignment ambiguity does not occur and hence can be trusted. The general problem of uncertainty in sequence alignment has recently been discussed by Wong et al. [16]. Due to the complexity of interpreting non-trivial alignment regions, these are often omitted in large-scale evolutionary analyses, even though there is ample evidence for their fundamental importance [16,17]. An approach to pinpointing alignment ambiguity is the generation of ensembles of suboptimal alignments [18], but computational demands remain prohibitive for genome wide studies.

Recent structural alignment methods have started to place emphasis on dealing with structural flexibility, such as FATCAT [19], MultiProt [20], MATT [21] and RAPIDO [22]. This may increase the consistency of alignments produced by each of these methods, but does not address the intrinsic ambiguity arising from structural divergence. The fundamental issue is whether a one-to-one equivalence exists between residues from different proteins that could be expressed as one definite alignment between sequences [18]. This is illustrated in Figure 1, where we show that a single insertion can lead to ambiguity in the functional correspondence between most residues in the loop.

thumbnailFigure 1. Dealing with structural flexibility: a single insertion (5', left) can lead to ambiguity in the pairwise residue alignment between the loops (right). Therefore, a simple one-to-one functional equivalence between residues from different proteins may not exist.

To further elucidate the effect of structural diversity on structural alignment, we prepared two distinct comprehensive sets of alternative structures for proteins from the HOMSTRAD database of homologous protein families. The first set comprises proteins for which alternative crystal structures are available. The other set is derived from molecular dynamics simulations to explore a more extensive spectrum of possible structures. An overview of our analysis procedure is outlined in Figure 2.

thumbnailFigure 2. Overview of the approach. SP scores are calculated to describe the differences at the sequence level between the reference and alternative structural alignments. In addition each alternative structure (either obtained with molecular simulation or from the PDB) is fit onto the reference structure and root mean square deviations (RMSDs) are calculated.

Our main results show that in many cases structural variation strongly affects structural alignments, even for highly similar sequences. Moreover, the derived alignment appears to be highly sensitive to even small conformational changes of the proteins. The uncertainty in pairing up structural equivalent residues makes it difficult to determine which alignment alternative would describe most closely the functional relationship between the proteins. To address this issue, we show how alternative alignment visualizations may be used to exploit the information contained within variable alignment regions.

Results and discussion

Structural diversity and alignment stability

The relation between the variation of the alternative structures (RMSD) and the corresponding alignment similarity (SP score) is shown in Figure 3 (bottom panel). It is clear that the structural variation between crystal structures (in orange) is much smaller (up to 3–4 Å RMSD) than that of the simulation snapshots (in blue; up to 10 Å RMSD). A crucial aspect is that even for small (<1 Å RMSD) and modest (1–3 Å RMSD) structural differences, alignments can easily vary up to 20% and sometimes as much as 40% or more in their SP score. On the other hand, a considerable number of alignments appear robust to larger (up to 6 Å RMSD) and even extreme (up to 10 Å RMSD) structural variations. Additionally, for the crystal structures, the sequence similarity has no effect on the variation in structural alignments (Figure 3, top panel). For the simulation snapshots, however, there seems to be a slight but distinct tendency for more similar sequences to have less variation in structural alignments, but this can be mainly attributed to the larger variations (>3 Å RMSD) in structure that arise from the simulations (see additional file 1). As an alternative for RMSD measurements we also tested the rho-score [23], a protein size-independent measure, which resulted in the same trend (data not shown).

Additional file 1. Figure S1: Combined effects of structural variation and sequence variation on the alignment.

Format: PDF Size: 4KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 3. Effects of structure and sequence variation on the alignment. The bottom panel shows structural difference (measured by the RMSD) versus alignment similarity, measured by the SP score, which is defined as the fraction of aligned reference residues pairs that are reproduced in the query alignment. The top panel shows distributions of SP scores for alignments sharing less and more than 50% sequence identity. Orange (lighter) refers to alternative crystal structures while blue (darker) refers to alternative structures obtained from molecular simulations.

A quite interesting example of the impact that small structural variations can have on the structural alignment is found in the GroES so-called 'mobile loop', which is the main region for interaction with GroEL and therefore is a crucial component of the GroEL/ES chaperonin machinery [24]. The structural variations for this loop in E. coli GroES (Figure 4C, shown in blue) are almost negligible (whole protein Cα RMSDs 0.42 ± 0.13 Å). It is therefore surprising that the corresponding DALI sequence alignments with M. tuberculosis GroES show remarkable variation in this region (Figure 4A). To pinpoint the source of this variation, we also used three other structural alignment programs: CE, MATT and FATCAT. The latter two explicitly take structural flexibility into account and this leads to more consistent alignments in the variable loop (alignment positions 20–69, Figure 4A). On the other hand, two regions (84–89 and 107–109, Figure 4A) are aligned consistently by DALI but show inconsistencies when aligned by CE and the two flexibility-aware methods. Strikingly, there is no overall consistency between the four methods, which is in line with several other studies where several structural alignment methods are compared [7,8,18]. It should be stressed that the focus of this paper is not on comparing the performance of the various methods but rather on the effects of structural diversity. A comprehensive overview of the GroES variability is given by the alignment matrix and the consistency plots (Figure 4B). The alignment matrix scores the occurrence of aligned residue pairs over all alignments, similar to the dot-plot [25,26]. Consistency plots show for each residue the standard deviation from the alignment position of the consensus pair. The alignment matrix and associated consistency plots allow a detailed visualization of the variability while enabling easy interpretation of the ensemble of alternative alignments.

thumbnailFigure 4. An example of the impact of tiny structural variations in the GroES 'mobile loop' that lead to quite dramatic variations in the alignment. A) The 'master-slave' alignments with 1p3h-N as master, variable regions marked at the top, and secondary structure with the mobile loop shown in red at the bottom. B) Alignment matrix with consistency plots along both axes give an overview of variability in each of the alignments from A). C) The different GroEs structures with 1p3h-N in red and 1aon-O and alternatives in blue. Alignment image created using JalView [33] with 'Zappo' colouring; secondary structure assignment according to Xu, et al. [24]. Protein structure rendered using SwissPDBViewer [34] and PovRay [35].

Although alignment uncertainty has been shown to have a great impact on large scale sequence analysis [16,17], the relation with structural variation has not been widely explored [27]. This is remarkable given that structural alignments are generally employed to benchmark sequence alignment methods. We demonstrate that in many cases structural alignments can vary dramatically even for small structural changes. Trends observed in the set of crystal structures corroborate those observed in the set of simulation snapshots, albeit alignment differences in the latter set are more pronounced due to larger structural variations.

A depository for alignment variability

It is questionable whether a single reference alignment captures the full width of naturally occurring sequence variability [18]. Yet, current visualization and alignment methods are not designed to take variable regions into account, and they are typically ignored in sequence alignment benchmark protocols. Since variable regions are often important structurally and/or functionally, new approaches for visualization, alignment and benchmarking are desirable.

To this end we have constructed a database of 'flexible' reference alignments. This database is available online http://www.ibi.vu.nl/databases/stralivari webcite and contains all structures and alignments used in this study. For each alignment in our database, variation is visualized using alignment matrices and consistency plots as shown in Figure 4B. In addition the database contains the ensemble 'master-slave' alignments as shown in Figure 4A. This pinpoints alignment regions that are affected by variability.

Conclusion

Structural variation, as presented here by alternative crystal structures and molecular dynamics simulations, has a profound effect on structural alignment. The sensitivity to structural variation is a bottleneck for the effective application of structural alignment approaches. This undermines the current basis of all sequence alignment methodologies and is an underestimated problem for the homology assessment used in structural and functional classification. The GroES 'mobile loop' example demonstrates how functionally essential protein regions can coincide with variable structural alignment segments. Our database should therefore be useful for alignment verification and delineation of functionally important protein regions.

Methods

The HOMSTRAD database of homologous structure alignments [14] was used as a source to select homologous proteins with known structure. HOMSTRAD families containing two homologous proteins (A and B in Figure 2) were selected. The corresponding structures were retrieved from the PDB [28] and taken as reference. For each reference structure, after equilibration, molecular dynamics simulations were performed for up to 10 ns, and snapshot structures were stored every 1 ns. Standard solvated conditions in the Gromos 43a1 forcefield [29] and the Gromacs simulation package [30] were used (details summarized in additional file 2). In addition, for each reference structure, we retrieved all alternative PDB structures with 100% sequence identity. In the subsequent analysis only the residues corresponding to the HOMSTRAD sequences were used.

Additional file 2. Table S1: Molecular Dynamics simulation set-up.

Format: PDF Size: 73KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

From each pair of reference HOMSTRAD structures, we constructed reference alignments with the widely used structural alignment tool DALI [1]. We also used DALI to create pairwise alignments between each reference structure and the alternatives of the other reference structure (PDB and snapshots). The sequence differences between the alignments were calculated using Sum-of-Pairs (SP) scoring implemented in the BAliBASE alignment comparison tool [13]. SP scores range from 0 (non-identical) to 1 (identical sequence alignments). Finally we calculated the root mean square deviation (RMSD) between the Cα atoms of the alternative structures and their reference structure using the McLachlan algorithm [31] as implemented in the program ProFit version 2.5.3 [32].

Our final database consists of 496 proteins (divided over 341 families) for which 3309 snapshot structures could be made and 565 proteins (divided over 395 families) for which we found in total 2998 alternative crystal structures with redundant sequences. A full list of all aligned structures and relevant details is provided in additional file 3.

Additional file 3. Table S2: Details of aligned Crystal Structures (a) and Simulation Snapshots (b).

Format: PDF Size: 95KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Authors' contributions

All authors designed the research, analyzed the results and wrote the paper. WP and KAF performed the research. All authors read and approved the final manuscript.

Acknowledgements

We like to thank Sander W. Timmer and Anneke van der Reijden for development of the data-analysis scripts and Bernd W. Brandt for the set of redundant protein structures. Financial support was provided by the Netherlands Bioinformatics Centre, BioRange Bioinformatics research programmes SP 3.2.2 and SP 2.3.1.

References

  1. Holm L, Park J: DaliLite workbench for protein structure comparison.

    Bioinformatics 2000, 16(6):566-567. PubMed Abstract | Publisher Full Text OpenURL

  2. Taylor WR, Orengo CA: Protein structure alignment.

    J Mol Biol 1989, 208(1):1-22. PubMed Abstract | Publisher Full Text OpenURL

  3. Gerstein M, Levitt M: Comprehensive assessment of automatic structural alignment against a manual standard, the scop classification of proteins.

    Protein Sci 1998, 7(2):445-456. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Lupyan D, Leo-Macias A, Ortiz AR: A new progressive-iterative algorithm for multiple structure alignment.

    Bioinformatics 2005, 21(15):3255-3263. PubMed Abstract | Publisher Full Text OpenURL

  5. Shindyalov IN, Bourne PE: Protein structure alignment by incremental combinatorial extension (CE) of the optimal path.

    Protein Eng 1998, 11(9):739-747. PubMed Abstract | Publisher Full Text OpenURL

  6. Sali A, Blundell TL: Definition of general topological equivalence in protein structures. A procedure involving comparison of properties and relationships through simulated annealing and dynamic programming.

    J Mol Biol 1990, 212(2):403-428. PubMed Abstract | Publisher Full Text OpenURL

  7. Kolodny R, Koehl P, Levitt M: Comprehensive evaluation of protein structure alignment methods: scoring by geometric measures.

    J Mol Biol 2005, 346(4):1173-1188. PubMed Abstract | Publisher Full Text OpenURL

  8. Mayr G, Domingues FS, Lackner P: Comparative analysis of protein structure alignments.

    BMC Struct Biol 2007, 7:50. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  9. Orengo CA, Michie AD, Jones S, Jones DT, Swindells MB, Thornton JM: CATH – a hierarchic classification of protein domain structures.

    Structure 1997, 5(8):1093-1108. PubMed Abstract | Publisher Full Text OpenURL

  10. Holm L, Ouzounis C, Sander C, Tuparev G, Vriend G: A database of protein structure families with common folding motifs.

    Protein Sci 1992, 1(12):1691-1698. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Bhaduri A, Pugalenthi G, Sowdhamini R: PASS2: an automated database of protein alignments organised as structural superfamilies.

    BMC Bioinformatics 2004, 5:35. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  12. Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures.

    J Mol Biol 1995, 247(4):536-540. PubMed Abstract | Publisher Full Text OpenURL

  13. Thompson JD, Plewniak F, Poch O: BAliBASE: a benchmark alignment database for the evaluation of multiple alignment programs.

    Bioinformatics 1999, 15(1):87-88. PubMed Abstract | Publisher Full Text OpenURL

  14. Mizuguchi K, Deane CM, Blundell TL, Overington JP: HOMSTRAD: a database of protein structure alignments for homologous families.

    Protein Sci 1998, 7(11):2469-2471. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. 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-1268. PubMed Abstract | Publisher Full Text OpenURL

  16. Wong KM, Suchard MA, Huelsenbeck JP: Alignment uncertainty and genomic analysis.

    Science 2008, 319(5862):473-476. PubMed Abstract | Publisher Full Text OpenURL

  17. Rokas A: Genomics. Lining up to avoid bias.

    Science 2008, 319(5862):416-417. PubMed Abstract | Publisher Full Text OpenURL

  18. Godzik A: The structural alignment between two proteins: is there a unique answer?

    Protein Sci 1996, 5(7):1325-1338. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Ye Y, Godzik A: Flexible structure alignment by chaining aligned fragment pairs allowing twists.

    Bioinformatics 2003, 19(Suppl 2):ii246-255. PubMed Abstract | Publisher Full Text OpenURL

  20. Shatsky M, Nussinov R, Wolfson HJ: A method for simultaneous alignment of multiple protein structures.

    Proteins 2004, 56(1):143-156. PubMed Abstract | Publisher Full Text OpenURL

  21. Menke M, Berger B, Cowen L: Matt: local flexibility aids protein multiple structure alignment.

    PLoS Comput Biol 2008, 4(1):e10. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Mosca R, Schneider TR: RAPIDO: a web server for the alignment of protein structures in the presence of conformational changes.

    Nucleic Acids Res 2008, (36 Web Server):W42-46. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Maiorov VN, Crippen GM: Size-independent comparison of protein three-dimensional structures.

    Proteins 1995, 22(3):273-283. PubMed Abstract | Publisher Full Text OpenURL

  24. Xu Z, Horwich AL, Sigler PB: The crystal structure of the asymmetric GroEL-GroES-(ADP)7 chaperonin complex.

    Nature 1997, 388(6644):741-750. PubMed Abstract | Publisher Full Text OpenURL

  25. Maizel JV Jr, Lenk RP: Enhanced graphic matrix analysis of nucleic acid and protein sequences.

    Proc Natl Acad Sci USA 1981, 78(12):7665-7669. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Zuker M: Suboptimal sequence alignment in molecular biology. Alignment with error analysis.

    J Mol Biol 1991, 221(2):403-420. PubMed Abstract | Publisher Full Text OpenURL

  27. Notredame C: Recent evolutions of multiple sequence alignment algorithms.

    PLoS Comput Biol 2007, 3(8):e123. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank.

    Nucleic Acids Res 2000, 28(1):235-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Hunenberger PH, Mark AE, van Gunsteren WF: Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations.

    J Mol Biol 1995, 252(4):492-503. PubMed Abstract | Publisher Full Text OpenURL

  30. Lindahl E, Hess B, Spoel D: GROMACS 3.0: a package for molecular simulation and trajectory analysis.

    J Mol Mod 2001, 7(8):306-317. OpenURL

  31. McLachlan A: Rapid comparison of protein structures.

    Acta Cryst 1982, A38:871-873. OpenURL

  32. ProFit [http://www.bioinf.org.uk/software/profit] webcite

  33. Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java alignment editor.

    Bioinformatics 2004, 20(3):426-427. PubMed Abstract | Publisher Full Text OpenURL

  34. Kaplan W, Littlejohn TG: Swiss-PDB Viewer (Deep View).

    Brief Bioinform 2001, 2(2):195-197. PubMed Abstract | Publisher Full Text OpenURL

  35. Persistence of Vision (TM) Raytracer [http://www.povray.org] webcite