当前位置: 首页 > 期刊 > 《核酸研究》 > 2006年第We期 > 正文
编号:11367377
OPAAS: a web server for optimal, permuted, and other alternative align
http://www.100md.com 《核酸研究医学期刊》
     Institute of Biomedical Sciences, Academia Sinica Taipei 11529, Taiwan

    *To whom correspondence should be addressed. Tel: +886 2 2789 9033; Fax: +886 2 2788 7641; Email: mjhwang@ibms.sinica.edu.tw

    ABSTRACT

    The large number of experimentally determined protein 3D structures is a rich resource for studying protein function and evolution, and protein structure comparison (PSC) is a key method for such studies. When comparing two protein structures, almost all currently available PSC servers report a single and sequential (i.e. topological) alignment, whereas the existence of good alternative alignments, including those involving permutations (i.e. non-sequential or non-topological alignments), is well known. We have recently developed a novel PSC method that can detect alternative alignments of statistical significance (alignment similarity P-value <10–5), including structural permutations at all levels of complexity. OPAAS, the server of this PSC method freely accessible at our website (http://opaas.ibms.sinica.edu.tw), provides an easy-to-read hierarchical layout of output to display detailed information on all of the significant alternative alignments detected. Because these alternative alignments can offer a more complete picture on the structural, evolutionary and functional relationship between two proteins, OPAAS can be used in structural bioinformatics research to gain additional insight that is not readily provided by existing PSC servers.

    INTRODUCTION

    Protein structure comparison (PSC) has been a staple method for obtaining information about a protein when its 3D structure is determined experimentally or predicted computationally. It is therefore not surprising that the development of new PSC algorithms has been continuing for more than two decades with no sign of ceasing (1–6). These efforts are needed not only to meet new scientific challenges but also to benefit maximally from the large number of new structures now pouring in from structural genomics projects (7,8). To these ends, a number of laboratories have created PSC servers in recent years to provide information beyond the basic PSC operations, including, e.g. those that do flexible alignment (9,10), those that discover recurring substructures or motifs (11,12), those that perform multiple structure alignment (13) and those that focus on fast structure feature extraction (14–16).

    Here we offer a new PSC server with the functionality to report statistically significant alternative alignments (17,18) and structural permutations (19,20) at all levels of complexity. Our method, named OPAAS, which has been detailed elsewhere (21,22), deduces the probabilities of aligning every possible pair of secondary structure elements (SSEs) between two protein structures prior to the search for a solution of their alignment. This deduction allows the ready identification of most, though not all, statistically significant alignment solutions, many of which being distinct alternatives to the ‘optimal’ solution, the target of conventional PSC operations. As we reported previously from a study of all-against-all database comparisons (22), about half of the alternative alignments were detectable only when permutation, i.e. non-topological alignment, was allowed. Moreover, many of the permuted alignments exhibited a permutation complexity higher than that of circular permutation, meaning that more than two separable regions of the protein structure could be aligned non-sequentially. To quantitatively measure the level of permutation complexity for all the alignments, we devised a permutation index (PI) as follows:

    where Si is the size (number of aligned amino acid residues) of the aligned region i and n is the total number of aligned regions. A region is an independently, and, within the region itself, topologically aligned part of an alignment. That is, within a region, all the aligned residues are ordered sequentially, which may or may not be interrupted by gaps, but these regions, if there are more than one, are aligned non-sequentially. It follows that an alignment without any permutation will have just one region, and will have, by definition, a PI value of 1.0. Also by definition, a circular permutation, which involves swapping two regions in a non-topological alignment (19,20), will receive a PI value >1.0 but not >2.0. PI hence furthermore let us know how much the sizes of the separately aligned regions differ. For example, given two permuted protein pairs having PI 3.0 and 2.5, respectively, we will know that they both have three aligned regions, but the sizes of the three regions are equal for the former and vary significantly for the latter.

    Both permuted and non-permuted alternative alignments are reported by the OPAAS server in a fashion that is easy for a non-specialist user to grasp the main significance of the comparison as one would with the ‘optimal’ alignment featured by other PSC servers. This is aided by the server's user-friendly interfaces described below, which use intuitive viewing directions, informative tables that can be sorted by different parameters, cascading information windows, and a structured user guide with examples.

    OPAAS WEB SERVER LAYOUT

    At the portal of the OPAAS web server (Figure 1a) lay two main structure comparison functions, ‘1 against SCOP90 dataset’ and ‘2 chains alignment’, and a Help webpage for a structured OPAAS user guide, which can be viewed on-line (http://opaas.ibms.sinica.edu.tw/help/opaas.html) as well as interactively in different contexts (see below).

    Figure 1 Snapshots of the OPAAS web server. (a) The portal page of the server. Two main applications can be chosen here. (b) The input page for the function of one-against-all search in the SCOP90 dataset. Three different input preferences and optional parameter settings are allowed in the left and right panel, respectively. (c) This table of matched results is retrieved from a database containing pre-computed all-against-all comparison results, which is available only with the one-against-all function. (d) The input page for the function of two chain alignment. (e) Output of a two chain alignment. The information of all alignment solutions (‘optimal’ and alternatives) is hierarchically laid out in three frames.

    One-against-all search on SCOP90

    The one-against-all on SCOP90 function is designed to find structural neighbors of a protein of interest in the structure classification of proteins (SCOP) (23) database. One of the following three input options (left panel in Figure 1b and user guide 3.1.1) is available for the search: a structure domain already in SCOP90 (SCOP version 1.55, <90% sequence identity non-redundant set), a structure in the current Protein Data Bank (PDB; the server updates its local PDB weekly) (24) with a specified chain (PDB ID option), or a structure in PDB format uploaded by user (User's structure option). All of the three inputs can be accompanied with optional parameter settings for customized output. The parameters that could be changed from default include minimum rough Z-score , minimum refined Z-score, maximum root mean square deviation (RMSD) of C superposition between two aligned protein chains, minimum aligned sequence identity, maximum number of shown matches, and sorting options (right panel in Figure 1b and user guide 3.1.2). Unlike the first input option, for which pre-computed results will be retrieved, the last two input options require entry of user's e-mail address (user guide 3.1.3) because the OPAAS server needs time for computation and will return the result via e-mail in, typically, minutes to hours, depending mainly on the stringency of the selected parameters.

    For the first input option, a table of matched results retrieved from a pre-computed database will be displayed directly on the web page (Figure 1c and user guide 3.2). In this table, Z-score, #a.a. aligned (number of aligned amino acids), and RMSD reported are those for the optimal alignment solution only, #alt. sol. is the number of alternative alignment solutions that satisfy the parameters specified by the user, and PI indicates the level of permutation complexity as described above. User can sort the table by the field selected (user guide 3.2.2) and download this table either as a plain text file or a comma separated values format file (user guide 3.2.3). User also can click on the matched entry to view details of each alignment solution, which will have exactly the same output as that from using the function of ‘2 chains alignment’ (user guide 3.2.4; see below). For the output of the last two input options, an e-mail of search result with the subject title ‘OPAAS result’ will be sent back to the e-mail address supplied (user guide 3.3). The e-mail contains a table of search result like that of the first input option described above, but to interactively view details of individual alignment solution, the second function, described below, needs to be separately invoked.

    Alignment of two protein chains

    For the function of ‘2 chains alignment’ (Figure 1d and user guide 4.1), the structures of the two protein chains could either be selected from PDB (i.e. entering PDB ID) or uploaded by user in PDB format. Display of the comparison results, which can be expected to follow immediately upon submission of the request, is split into three frames to show ‘basic information of the comparison’, ‘information of each solution’, and ‘detail alignment of selected region’ (Figure 1e and user guide 4.2).

    In the frame of ‘basic information of the comparison’, one of two tables gives some basic information about the two chains compared: size (number of amino acid residues), and number of SSE (user guide 4.2.1). User can click a hyperlink on the name of the two chains placed above the table to learn more about the compared proteins from the PDB website (user guide 4.2.2). The other table shows information of all the solutions (both the ‘optimal’ and alternative alignments, if any) of this comparison including #a.a. aligned (number of aligned residues), RMSD, refined Z-score, PI and region (number of permutedly aligned regions). Clicking on the solution number shows details of that solution in another two frames (user guide 4.2.3). The frame of ‘information of each solution’ shows the alignment of the selected solution graphically (user guide 4.2.4), both in a schematic representation and in a 3D superposition supported by chime plug-in (user guide 4.2.5). A file containing Cartesian coordinates of this alignment solution in standard PDB format can be exported. Different colors of the boxes in the diagram of the schematic representation, as well as of the traces in the 3D superposition, refer to different aligned regions. Clicking on the region box will show ‘detail alignment of the selected region’ in the third frame. The sequence alignment of this region can also be downloaded in MFA (Multi-FASTA Alignment) format by clicking ‘export alignment’ at the top of this frame.

    DISCUSSION

    The best way to compare two protein structures often depends on the question being asked (6), so having a server like OPAAS that can simultaneously analyze solutions beyond the ‘optimal’ alignment is useful. Although most of the published PSC algorithms can be modified to offer similar capability, to our knowledge, only two PSC servers give user the option to see alternative alignments: Prosup (25) and SARF2 (26), but Prosup is limited to topological alignments and neither offers one-against-all database searching service. Moreover, with an intuitive hierarchical layout of the comparison results and optional parameter settings to view most significant alignments (e.g. with similarity P-value set at 10–5, a typical comparison usually resulted in <5 such solutions), an informative summary that could lead to unexpected insight from unexpected alternative alignments is effectively produced. The main limitation of OPAAS, at its current version, arose from a compromise to trade for computational efficiency, which dictates that a structure must possess at least three SSEs to be compared (21,22); elimination of this limitation is in progress. Significantly, our server allows database search with an efficiency comparable that of the popular CE server (27), despite ours being run on a personal computer (Pentium IV) and being asked to find alternative alignments. The source code of OPAAS is also available at the server for free download for standalone computations and for incorporation of structure database other than SCOP.

    ACKNOWLEDGEMENTS

    Funding to pay the Open Access publication charges for this article was provided by the Academia Sinica and National Science Council of Taiwan.

    REFERENCES

    Orengo, C. (1994) Classification of protein folds Curr. Opin. Struct. Biol, . 4, 429–440 .

    Brown, N.P., Orengo, C.A., Taylor, W.R. (1996) A protein structure comparison methodology Comput. Chem, . 20, 359–380 .

    Gibrat, J.-F., Madej, T., Bryant, S.H. (1996) Surprising similarities in structure comparison Curr. Opin. Struct. Biol, . 6, 377–385 .

    Holm, L. and Sander, C. (1996) Mapping the protein universe Science, 273, 595–602 .

    Koehl, P. (2001) Protein structure similarities Curr. Opin. Struct. Biol, . 11, 348–353 .

    Sierk, M.L. and Kleywegt, G.J. (2004) Deja vu all over again: finding and analyzing protein structure similarities Structure, 12, 2103–2111 .

    Todd, A.E., Marsden, R.L., Thornton, J.M., Orengo, C.A. (2005) Progress of structural genomics initiatives: an analysis of solved target structures J. Mol. Biol, . 348, 1235–1260 .

    Chandonia, J.M. and Brenner, S.E. (2006) The impact of structural genomics: expectations and outcomes Science, 311, 347–351 .

    Shatsky, M., Nussinov, R., Wolfson, H.J. (2002) Flexible protein alignment and hinge detection Proteins, 48, 242–256 .

    Ye, Y. and Godzik, A. (2004) FATCAT: a web server for flexible structure comparison and structure similarity searching Nucleic Acids Res, . 32, W582–W585 .

    Chang, D.T., Chen, C.Y., Chung, W.C., Oyang, Y.J., Juan, H.F., Huang, H.C. (2004) ProteMiner-SSM: a web server for efficient analysis of similar protein tertiary substructures Nucleic Acids Res, . 32, W76–W82 .

    Shapiro, J. and Brutlag, D. (2004) FoldMiner and LOCK 2: protein structure comparison and motif discovery on the web Nucleic Acids Res, . 32, W536–W541 .

    Guda, C., Lu, S., Scheeff, E.D., Bourne, P.E., Shindyalov, I.N. (2004) CE-MC: a multiple protein structure alignment server Nucleic Acids Res, . 32, W100–W103 .

    Shyu, C.R., Chi, P.H., Scott, G., Xu, D. (2004) ProteinDBS: a real-time retrieval system for protein structure comparison Nucleic Acids Res, . 32, W572–W575 .

    Maiti, R., Van Domselaar, G.H., Zhang, H., Wishart, D.S. (2004) SuperPose: a simple server for sophisticated structural superposition Nucleic Acids Res, . 32, W590–W594 .

    Vlahovicek, K., Pintar, A., Parthasarathi, L., Carugo, O., Pongor, S. (2005) CX, DPX and PRIDE: WWW servers for the analysis and comparison of protein 3D structures Nucleic Acids Res, . 33, W252–W254 .

    Feng, Z.K. and Sippl, M.J. (1996) Optimum superimposition of protein structures: ambiguities and implications Fold Des, . 1, 123–132 .

    Godzik, A. (1996) The structural alignment between two proteins: is there a unique answer? Protein Sci, . 5, 1325–1338 .

    Cunningham, B.A., Hemperly, J.J., Hopp, T.P., Edelman, G.M. (1979) Favin versus concanavalin A: circularly permuted amino acid sequences Proc. Natl Acad. Sci, . 76, 3218–3222 .

    Heinemann, U. and Hahn, M. (1995) Circular permutation of polypeptide chains: implications for protein folding and stability Prog. Biophys. Mol. Biol, . 64, 121–143 .

    Shih, E.S.C. and Hwang, M.J. (2003) Protein structure comparison by probability-based matching of secondary structure elements Bioinformatics, . 19, 735–741 .

    Shih, E.SC. and Hwang, M.J. (2004) Alternative alignments from comparison of protein structures Proteins, 56, 519–527 .

    Murzin, A.G., Brenner, S.E., Hubbard, T., Chothia, C. (1995) SCOP: a structural classification of proteins database for the investigation of sequences and structures J. Mol. Biol, . 247, 536–540 .

    Berman, H.M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T.N., Weissig, H., Shindyalov, I.N., Bourne, P.E. (2000) The Protein Data Bank Nucleic Acids Res, . 28, 235–242 .

    Lackner, P., Koppensteiner, W.A., Sippl, M.J., Domingues, F.S. (2000) ProSup: a refined tool for protein structure alignment Protein Eng, . 13, 745–752 .

    Alexandrov, N.N. and Fischer, D. (1996) Analysis of topological and nontopological structural similarities in the PDB: new examples with old structures Proteins, 25, 354–365 .

    Shindyalov, I.N. and Bourne, P.E. (1998) Protein structure alignment by incremental combinatorial extension (CE) of the optimal path Protein Eng, . 11, 739–747 .(Edward S. C. Shih, Ruei-chi R. Gan and M)