Abstract:We have tested CLUSTAL W in a wide variety of situations, and it is capable of handling some very difficult protein alignment problems. If the data set consists of enough closely related sequences so that the first alignments are accurate, then CLUSTAL W will usually find an alignment that is very close to ideal. Problems can still occur if the data set includes sequences of greatly different lengths or if some sequences include long regions that are impossible to align with the rest of the data set. Trying to balance the need for long insertions and deletions in some alignments with the need to avoid them in others is still a problem. The default values for our parameters were tested empirically using test cases of sets of globular proteins where some information as to the correct alignment was available. The parameter values may not be very appropriate with nonglobular proteins. We have argued that using one weight matrix and two gap penalties is too simplistic to be of general use in the most difficult cases. We have replaced these parameters with a large number of new parameters designed primarily to help encourage gaps in loop regions. Although these new parameters are largely heuristic in nature, they perform surprisingly well and are simple to implement. The underlying speed of the progressive alignment approach is not adversely affected. The disadvantage is that the parameter space is now huge; the number of possible combinations of parameters is more than can easily be examined by hand. We justify this by asking the user to treat CLUSTAL W as a data exploration tool rather than as a definitive analysis method. It is not sensible to automatically derive multiple alignments and to trust particular algorithms as being capable of always getting the correct answer. One must examine the alignments closely, especially in conjunction with the underlying phylogenetic tree (or estimate of it) and try varying some of the parameters. Outliers (sequences that have no close relatives) should be aligned carefully, as should fragments of sequences. The program will automatically delay the alignment of any sequences that are less than 40% identical to any others until all other sequences are aligned, but this can be set from a menu by the user. It may be useful to build up an alignment of closely related sequences first and to then add in the more distant relatives one at a time or in batches, using the profile alignments and weighting scheme described earlier and perhaps using a variety of parameter settings. We give one example using SH2 domains. SH2 domains are widespread in eukaryotic signalling proteins where they function in the recognition of phosphotyrosine-containing peptides. In the chapter by Bork and Gibson ([11], this volume), Blast and pattern/profile searches were used to extract the set of known SH2 domains and to search for new members. (Profiles used in database searches are conceptually very similar to the profiles used in CLUSTAL W: see the chapters [11] and [13] for profile search methods.) The profile searches detected SH2 domains in the JAK family of protein tyrosine kinases, which were thought not to contain SH2 domains. Although the JAK family SH2 domains are rather divergent, they have the necessary core structural residues as well as the critical positively charged residue that binds phosphotyrosine, leaving no doubt that they are bona fide SH2 domains. The five new JAK family SH2 domains were added sequentially to the existing alignment of 65 SH2 domains using the CLUSTAL W profile alignment option. Figure 6 shows part of the resulting alignment. Despite their divergent sequences, the new SH2 domains have been aligned nearly perfectly with the old set. No insertions were placed in the original SH2 domains. In this example, the profile alignment procedure has produced better results than a one-step full alignment of all 70 SH2 domains, and in considerably less time. (ABSTRACT TRUNCATED)

ClustalXeed: a GUI-based grid computation version for high performance and terabyte size multiple sequence alignment

Gene Sequence Alignment on a Public Computing Platform

HAlign-II: Efficient Ultra-Large Multiple Sequence Alignment and Phylogenetic Tree Reconstruction with Distributed and Parallel Computing

Parallel Algorithms for Large-Scale Biological Sequence Alignment on Xeon-Phi Based Clusters

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

Fast, Scalable Generation of High-Quality Protein Multiple Sequence Alignments Using Clustal Omega

Cluster-Distribute-Align-Merge: A General Algorithm to Speed Up Multiple Sequence Alignment on Multi-Core Computers

MAGUS: Multiple sequence Alignment using Graph clUStering

CUDA ClustalW: An efficient parallel algorithm for progressive multiple sequence alignment on Multi-GPUs

Kalign – an accurate and fast multiple sequence alignment algorithm

GPU-ClustalW: using graphics hardware to accelerate multiple sequence alignment

Sample-Align-D: A High Performance Multiple Sequence Alignment System using Phylogenetic Sampling and Domain Decomposition

Load balancing and parallel multiple sequence alignment with tree accumulation

Memory Efficient Pair-Wise Genome Alignment Algorithm - A Small-Scale Application With Grid Potential

ClustalY: Speed Up the Guide Tree Building for ClustalW

Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences

Using CLUSTAL for multiple sequence alignments

A data parallel strategy for aligning multiple biological sequences on multi-core computers

HAlign 4: A New Strategy for Rapidly Aligning Millions of Sequences

CD-HIT: accelerated for clustering the next-generation sequencing data

Distributed Sequence Alignment Applications for the Public Computing Architecture