Abstract:The FASTA program can search the NBRF protein sequence library (2.5 million residues) in less than 20 min on an IBM-PC microcomputer and unambiguously detect proteins that shared a common ancestor billions of years in the past. FASTA is both fast and selective because it initially considers only amino acid identities. Its sensitivity is increased not only by using the PAM250 matrix to score and rescore regions with large numbers of identities but also by joining initial regions. The results of searches with FASTA compare favorably with results using NWS-based programs that are 100 times slower. FASTA is slightly less sensitive but considerably more selective. It is not clear that NWS-based programs would be more successful in finding distantly related members of the G-protein-coupled receptor family. The joining step by FASTA to calculate the initn score is especially useful for sequences that share regions of sequence similarity that are separated by variable-length loops. FASTP and FASTA were designed to identify protein sequences that have descended from a common ancestor, and they have proved very useful for this task. In many cases, a FASTA sequence search will result in a list of high scoring library sequences that are homologous to the query sequence, or the search will result in a list of sequences with similarity scores that cannot be distinguished from the bulk of the library. In either case, the question of whether there are sequences in the library that are clearly related to the query sequence has been answered unambiguously. Unfortunately, the results often will not be so clear-cut, and careful analysis of similarity scores, statistical significance, the actual aligned residues, and the biological context are required. In the course of analyzing the G-protein-coupled receptor family, several proteins were found that, because of a high initn score and a low init1 score that increased almost 2-fold with optimization, appeared to be members of this family which were not previously recognized. RDF2 analysis showed borderline z values, and only a careful examination of the sequence alignments that focused on the conserved residues provided convincing evidence that the high scores were fortuitous. As sequence comparison methods become more powerful by becoming more sensitive, they become more likely to mislead, and even greater care is required.

COMPASS: A Tool for Comparison of Multiple Protein Alignments with Assessment of Statistical Significance

COMPASS server for remote homology inference

Rapid and Sensitive Protein Complex Alignment with Foldseek-Multimer

Multi-species comparative mapping in silico using the COMPASS strategy

Fold recognition by scoring protein maps using the congruence coefficient

Revealing divergent evolution, identifying circular permutations and detecting active-sites by protein structure comparison

Utilizing the Scale-Invariant Feature Transform Algorithm to Align Distance Matrices Facilitates Systematic Protein Structure Comparison

A Similarity Computing Algorithm for Proteins

COBALT: constraint-based alignment tool for multiple protein sequences

Finding Optimal Interaction Interface Alignments Between Biological Complexes

Multiple Alignment-Free Sequence Comparison

Rapid and sensitive sequence comparison with FASTP and FASTA

Rapid similarity search of proteins using alignments of domain arrangements

MRFalign: Protein Homology Detection through Alignment of Markov Random Fields

Protalign: a 3-dimensional protein alignment assessment tool

Combined alignments of sequences and domains characterize unknown proteins with remotely related protein search PSISearch2D

A comprehensive software suite for protein family construction and functional site prediction

Scoring profile-to-profile sequence alignments

CSA: comprehensive comparison of pairwise protein structure alignments

Sequence alignment using large protein structure alphabets improves sensitivity to remote homologs

Advancing the Accuracy of Protein Fold Recognition by Utilizing Profiles from Hidden Markov Models.