IsarPipeline: Combining MMseqs2 and PSI-BLAST to Quickly Generate Extensive Protein Sequence Alignment Profiles

Many of the machine learning (ML) models used in the field of bioinformatics and computational biology to predict either function or structure of proteins rely on the evolutionary information as summarized in multiple-sequence alignments (MSAs) or the resulting position-specific scoring matrices (PSSMs), as generated by PSI-BLAST. The current procedure used in protein structure and function prediction is computationally exhaustive and time-consuming. The main issue relies on the PSI-BLAST software being forced to load the current database of sequences (about 220 GB) in batches and search for similar sequence alignments to a query sequence. This leads to an average runtime of about 40-60 min for a medium-sized (450 Amino Acids) query protein. This average runtime is strictly dependent on the hardware used to run the software. The issue is becoming more problematic since the bio-sequence data pools are increasing in size exponentially over time, hence raising PSI-BLAST runtime as well. A prominent solution claims to speed up the current process by 100 folds. The MMseqs2 method, given enough memory, will load the whole database in memory and apply certain heuristics to retrieve the relevant set of aligned sequences. However, this solution cannot be used directly to generate the final output in the desired PSI-BLAST alignment and PSSM profile data format. In this research project, we analyzed the runtime performance of each tool separately. Furthermore, we built a pipeline that combines both MMseqs2 and PSI-BLAST to obtain a robust, optimized and very fast hybrid alignment tool, faster than PSI-BLAST by two orders of magnitude. It is implemented in C++ and is freely available under the MIT license at https://github.com/issararab/IsarPipeline. The output of our pipeline was evaluated on two previously built predictive models.