|Home | About | Journals | Submit | Contact Us | Français|
Summary: Accurate alignment of protein sequences and/or structures is crucial for many biological analyses, including functional annotation of proteins, classifying protein sequences into families, and comparative protein structure modeling. Described here is a web interface to SALIGN, the versatile protein multiple sequence/structure alignment module of MODELLER. The web server automatically determines the best alignment procedure based on the inputs, while allowing the user to override default parameter values. Multiple alignments are guided by a dendrogram computed from a matrix of all pairwise alignment scores. When aligning sequences to structures, SALIGN uses structural environment information to place gaps optimally. If two multiple sequence alignments of related proteins are input to the server, a profile–profile alignment is performed. All features of the server have been previously optimized for accuracy, especially in the contexts of comparative modeling and identification of interacting protein partners.
Availability: The SALIGN web server is freely accessible to the academic community at http://salilab.org/salign. SALIGN is a module of the MODELLER software, also freely available to academic users (http://salilab.org/modeller).
As part of our effort to better understand cellular processes, we annotate the functions of proteins, classify them into families, determine or model their 3D structures, categorize these structures, etc. To achieve these aims, robust and accurate methods for aligning protein sequences and structures with one another are essential. Alignment methods can be divided into three categories based on the information used to align proteins: sequence–sequence, sequence–structure, and structure–structure alignments. The SALIGN web server is designed to provide a user-friendly interface to producing robust and accurate solutions for the alignment problems in these different categories. Although several programs and servers (e.g. Di Tommaso et al., 2011; Shatsky et al., 2004) align protein sequences or protein structures with one another, SALIGN can align sequences, structures, or a combination of sequences and structures.
The multi-purpose alignment module of MODELLER (Sali and Blundell, 1993), SALIGN (Madhusudhan, et al., 2006; 2009; Marti-Renom et al., 2004; Pieper et al., 2011), can produce alignments of two or more protein sequences, based on sequence and/or structure information. Given a set of input protein sequences and/or structures, all possible pairwise alignments are calculated using dynamic programming.
SALIGN provides two methods for constructing multiple alignments from pairwise alignments, ‘tree alignment’ and ‘progressive alignment’ (Madhusudhan et al., 2009). The tree algorithm first creates a dendrogram of the input proteins from a matrix of all pairwise alignment scores. This guides the multiple alignment by serially aligning, in pairwise fashion, the closest linked branches to each other. The progressive alignment algorithm is computationally less expensive, where two arbitrary sequences are first aligned to one another, followed by the addition of a third; in n-1 steps, a multiple alignment of n sequences is created.
If two pre-aligned blocks of sequences are to be aligned, the profile–profile alignment method is used (Marti-Renom et al., 2004). To align a block of sequences to a block of structures, the ALIGN2D algorithm is used, where affine gap penalties are replaced by an environment-dependent gap penalty function (Madhusudhan et al., 2006).
SALIGN's default settings suffice for many applications. It has been fine tuned and extensively tested for alignment accuracy (Davis et al., 2006; Madhusudhan et al., 2006; 2009; Marti-Renom et al., 2004; 2007; Pieper et al., 2011). Nevertheless, the interface allows the user to manipulate many options if so desired.
The main interface to the server allows the user to input protein structures and/or sequences. Structures can be uploaded in protein data bank (PDB) format or specified by their four-letter PDB codes. Sequences can be uploaded as files in the FASTA or MODELLER Protein Information Resource (PIR) format, or pasted into the sequence window. Alignment files in the MODELLER PIR format can also be used to specify the input structures; the server searches the PDB for the structures denoted in the alignment file. In the event no structure is found matching an alignment entry, it is treated as a sequence.
Given the input, the server decides upon the optimal alignment protocol. If structures were input, the user has the option of specifying a structure segment (e.g. a single chain or domain). Further, on the advanced view page, the user can override the default alignment protocol and parameter values. The parameters that can be tweaked depend on the selected alignment protocol. For instance, if a structure–sequence alignment is chosen, parameters related to the environment-dependent gap penalty function can be changed. An explanation for each of the tunable parameters is provided in the help pages or via links to the MODELLER manual.
After the successful completion of an alignment task, the user is e-mailed a link to an archive of output files. The archive contains the final alignment file, superimposed coordinates if structures were aligned, a dendrogram file if a tree was constructed, and the MODELLER input and log files. The input file(s) can be used with a standalone version of MODELLER (version 9 or higher). The log file contains information about the calculation, including root mean square deviation values, numbers of equivalent positions, and any warnings and errors. If structures have been aligned, the output page also includes a link that displays the aligned structures in the molecular graphics viewer Chimera (Pettersen et al., 2004). This functionality provides an instant visualization of the alignment (Fig. 1).
We thank Eric Pettersen for helpful discussions.
Funding: National Institutes of Health (NIH), U.S.A [U54 GM094662 and R01 GM54762 to A.S.]
Conflict of Interest: none declared.