The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSIBLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

/pdf/gapped-blast-and-psi-blast-a-new-generation-of-protein-2ybl67x5yx.pdf

Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

/pdf/clustal-w-improving-the-sensitivity-of-progressive-multiple-2na2s9sq2h.pdf

Clustal w: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

CLUSTAL X is a new windows interface for the widely-used progressive multiple sequence alignment program CLUSTAL W. The new system is easy to use, providing an integrated system for performing multiple sequence and profile alignments and analysing the results. CLUSTAL X displays the sequence alignment in a window on the screen. A versatile sequence colouring scheme allows the user to highlight conserved features in the alignment. Pull-down menus provide all the options required for traditional multiple sequence and profile alignment. New features include: the ability to cut-and-paste sequences to change the order of the alignment, selection of a subset of the sequences to be realigned, and selection of a sub-range of the alignment to be realigned and inserted back into the original alignment. Alignment quality analysis can be performed and low-scoring segments or exceptional residues can be highlighted. Quality analysis and realignment of selected residue ranges provide the user with a powerful tool to improve and refine difficult alignments and to trap errors in input sequences. CLUSTAL X has been compiled on SUN Solaris, IRIX5.3 on Silicon Graphics, Digital UNIX on DECstations, Microsoft Windows (32 bit) for PCs, Linux ELF for x86 PCs, and Macintosh PowerMac.

/pdf/the-clustal-x-windows-interface-flexible-strategies-for-3qez0eyvqt.pdf

The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.

We describe MUSCLE, a new computer program for creating multiple alignments of protein sequences. Elements of the algorithm include fast distance estimation using kmer counting, progressive alignment using a new profile function we call the logexpectation score, and refinement using treedependent restricted partitioning. The speed and accuracy of MUSCLE are compared with T-Coffee, MAFFT and CLUSTALW on four test sets of reference alignments: BAliBASE, SABmark, SMART and a new benchmark, PREFAB. MUSCLE achieves the highest, or joint highest, rank in accuracy on each of these sets. Without refinement, MUSCLE achieves average accuracy statistically indistinguishable from T-Coffee and MAFFT, and is the fastest of the tested methods for large numbers of sequences, aligning 5000 sequences of average length 350 in 7 min on a current desktop computer. The MUSCLE program, source code and PREFAB test data are freely available at http://www.drive5. com/muscle.

MUSCLE: multiple sequence alignment with high accuracy and high throughput

The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large-scale molecular assemblies such as viral coats, and Collaboratory, which allows researchers to share a Chimera session interactively despite being at separate locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and associated structures; ViewDock, for screening docked ligand orientations; Movie, for replaying molecular dynamics trajectories; and Volume Viewer, for display and analysis of volumetric data. A discussion of the usage of Chimera in real-world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.

/pdf/ucsf-chimera-a-visualization-system-for-exploratory-research-59xu54fhxg.pdf

UCSF Chimera--a visualization system for exploratory research and analysis.

Lorist6, a cosmid vector with BamHI, NotI, ScaI and HindIII cloning sites and altered neomycin phosphotransferase gene expression

Cyclins, TFIIB and RB play major roles in cell cycle and/or gene regulation. Earlier work has suggested common ancestry for the TFIIB repeats and RB pocket B which share 20% sequence identity. We now report that database searches with profiles based on a multiple alignment of cyclin core regions (the 'cyclin box') detect the TFIIB repeats with equivalent scores to divergent cyclins. Several features of the sequences support the notion of common ancestry: e.g. cyclins A/B, C and D share approximately 20-30% identity but each have approximately 15-20% identity with vertebrate TFIIB, showing that conserved cyclin features underlie the match. These results suggest the presence of a domain superfamily, which we term the TR domain, in nuclear regulatory proteins belonging to the TFIIB, cyclin and RB families, that has been duplicated many times during eukaryotic evolution. The TR domain appears to function in protein-protein interactions.

/pdf/evidence-for-a-protein-domain-superfamily-shared-by-the-4jv5aqqv4m.pdf

Evidence for a protein domain superfamily shared by the cyclins, TFIIB and RB/p107.

/pdf/dna-binding-domain-ancestry-44n7lx5jj3.pdf

DNA-binding domain ancestry.

Formin was originally isolated as the gene affected by the murine limb deformity (ld) mutations, which disrupt the epithelial-mesenchymal interactions regulating patterning of the vertebrate limb autopod. More recently, a rapidly growing number of genes with similarity to formin have been isolated from many different species including fungi and plants. Genetic and biochemical analysis shows that formin family members function in cellular processes regulating either cytokinesis and/or cell polarisation. Another common feature among formin family members is their requirement in morphogenetic processes such as budding and conjugation of yeast, establishment of Drosophila oocyte polarity and vertebrate limb pattern formation. Vertebrate formins are predominantly nuclear proteins which control polarising activity in limb buds through establishment of the SHH/FGF-4 feedback loop. Formin acts in the limb bud mesenchyme to induce apical ectodermal ridge (AER) differentiation and FGF-4 expression in the posterior AER compartment. Finally, disruption of the epithelial-mesenchymal interactions controlling induction of metanephric kidneys in ld mutant embryos indicates that formin might function more generally in transduction of morphogenetic signals during embryonic pattern formation.

Formin defines a large family of morphoregulatory genes and functions in establishment of the polarising region

The solution structure of an active synthetic peptide containing both the leucine zipper and the adjacent basic domain of the yeast transcription factor GCN4 (residues 220-280) was determined by NMR. The two domains show structurally distinct behaviours. In the absence of DNA, the basic domain is, although very flexible, structured and fluctuating around a helical conformation. The leucine zipper region forms a long, uninterrupted helix. From a suitable set of NMR distances the three-dimensional structure of the leucine zipper monomeric sub-domain was calculated by distance geometry algorithms. The structure of the symmetrical parallel dimer was obtained by model building using the NMR information. A smaller peptide with the sequence of the isolated basic region (residues 1-35 of the 61 residue peptide) was also synthesized. Circular dichroism studies showed 30-40% helicity. A flexible helix spans the region between residues 8 and 21. The comparison of our results with suggested models is discussed in detail.

Toby J. Gibson

Papers

Lorist6, a cosmid vector with BamHI, NotI, ScaI and HindIII cloning sites and altered neomycin phosphotransferase gene expression

Evidence for a protein domain superfamily shared by the cyclins, TFIIB and RB/p107.

DNA-binding domain ancestry.

Formin defines a large family of morphoregulatory genes and functions in establishment of the polarising region

Solution structure of the DNA-binding domain of the yeast transcriptional activator protein GCN4