Abstract:Motivation There are now 225 million sequences in the UniProtKB database, as of January 2022 and 451 million protein sequences in the NCBI non-redundant database. This huge sequence data is ripe for analysis and can be extremely informative about biological function when analyzed with the appropriate methods. Evolutionary information such as the relationship among protein sequences is key to performing sequence analyses. Since sequence matching is one of the primary ways that annotations are found, higher-quality sequence matches yield a larger number of identified homologs. Thus, there is an essential need for a faster and more accurate homolog detection method to process the huge amount of rapidly growing biological sequences. Method Recently, we have seen major improvements in various predictive computational tasks such as structure prediction from the ever-improving artificial intelligence methods. One such approach has been to use language models to represent proteins numerically in a representation matrix (embeddings) while retaining context-dependent biochemical, biophysical, and evolutionary information. Computational transformer architectures that utilize attention neural networks can generate these context-aware numerical representations in an unsupervised fashion. One such use for these protein embeddings is remote homolog detection. In this work, we utilize protein language models and then apply discrete cosine transforms to extract the essential part of these embeddings, resulting in a significantly smaller fixed-size matrix for each sequence. This allows us to numerically and efficiently calculate the distance between all pairs of proteins resulting in homolog detection. Results Our Protein LAnguage model Search Tool (PLAST) is significantly faster, with linear runtimes in the number of sequences within the query database. With only one CPU core, it can scan a million sequences in less than a second. It essentially removes the noise in the sequence data and leads to significant improvements. PLAST is more accurate in the benchmarks tested from the PFAM, SCOP, and CATH databases than other approaches. When benchmarked with the PFAM database, the increase in the area under the receiver operating characteristic curve (AUROC) 3.1% when compared with NCBI-BLAST. The number of remote homologs that are detectable now is significantly larger and pushes sequence matches deeply into the usual twilight zone. Compared with the state-of-the-art profile-based homology search tools like CSBLAST, the increase was still 2.0%. PLAST can find remote homologs for a significant number of proteins that had been thought to be unique due to homolog detection failure. These homologs that are found usually have less than 20% sequence identity making them indistinguishable from noise with most other sequence matching methods. Conclusion PLAST is an accurate and fast homolog detection tool essential for easy and rapid progress to utilize the vast amount of data generated by next-generation sequencing methods. Quantization of sequence embeddings into highly-compressed noise-free representations with the use of direct cosine transforms allows for the efficient and accurate detection of normal homologs and remote ones that are un-detectable by other sequence similarity methods. The PLAST web server is accessible from https://mesihk.github.io/plast .

Rapid multiple protein sequence search by parallel and heterogeneous computation

Gene Sequence Alignment on a Public Computing Platform

An efficient parallel algorithm for multiple sequence similarities calculation using a low complexity method.

A Multiple Criteria Framework for 3D Protein Structure Similarity Retrieval

fastMSA: Accelerating Multiple Sequence Alignment with Dense Retrieval on Protein Language

Protein language model powers accurate and fast sequence search for remote homology

Extreme-scale many-against-many protein similarity search

CUDAMPF++: A Proactive Resource Exhaustion Scheme for Accelerating Homologous Sequence Search on CUDA-enabled GPU

Rapid and Sensitive Protein Complex Alignment with Foldseek-Multimer

GPU-accelerated homology search with MMseqs2

A high-throughput gene sequence alignment strategy using parallel computing

HiCOPS: High Performance Computing Framework for Tera-Scale Database Search of Mass Spectrometry based Omics Data

PLMSearch: Protein language model powers accurate and fast sequence search for remote homology

Accurate Multiple Sequence Alignment of Ultramassive Genome Sets

A PLMs based protein retrieval framework

A Parallel Implementation for Large-Scale TSR-based 3D Structural Comparisons of Protein and Amino Acid

Scalable Protein Sequence Similarity Search using Locality-Sensitive Hashing and MapReduce

Exploring large protein sequence space through homology- and representation-based hierarchical clustering

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

Protein Language Model Performs Efficient Homology Detection

CRAFT: Compact genome Representation toward large-scale Alignment-Free daTabase