PSI-Search: iterative HOE-reduced profile SSEARCH searching

Summary: Iterative similarity searches with PSI-BLAST position-specific score matrices (PSSMs) find many more homologs than single searches, but PSSMs can be contaminated when homologous alignments are extended into unrelated protein domains—homologous over-extension (HOE). PSI-Search combines an optimal Smith–Waterman local alignment sequence search, using SSEARCH, with the PSI-BLAST profile construction strategy. An optional sequence boundary-masking procedure, which prevents alignments from being extended after they are initially included, can reduce HOE errors in the PSSM profile. Preventing HOE improves selectivity for both PSI-BLAST and PSI-Search, but PSI-Search has ~4-fold better selectivity than PSI-BLAST and similar sensitivity at 50% and 60% family coverage. PSI-Search is also produces 2- for 4-fold fewer false-positives than JackHMMER, but is ~5% less sensitive.

Availability and implementation: PSI-Search is available from the authors as a standalone implementation written in Perl for Linux-compatible platforms. It is also available through a web interface (www.ebi.ac.uk/Tools/sss/psisearch) and SOAP and REST Web Services (www.ebi.ac.uk/Tools/webservices).

Contact:pearson@virginia.edu; rodrigo.lopez@ebi.ac.uk

Accelerating pairwise statistical significance estimation for local alignment by harvesting GPU's power

Background

Pairwise statistical significance has been recognized to be able to accurately identify related sequences, which is a very important cornerstone procedure in numerous bioinformatics applications. However, it is both computationally and data intensive, which poses a big challenge in terms of performance and scalability.

Results

We present a GPU implementation to accelerate pairwise statistical significance estimation of local sequence alignment using standard substitution matrices. By carefully studying the algorithm's data access characteristics, we developed a tile-based scheme that can produce a contiguous data access in the GPU global memory and sustain a large number of threads to achieve a high GPU occupancy. We further extend the parallelization technique to estimate pairwise statistical significance using position-specific substitution matrices, which has earlier demonstrated significantly better sequence comparison accuracy than using standard substitution matrices. The implementation is also extended to take advantage of dual-GPUs. We observe end-to-end speedups of nearly 250 (370) × using single-GPU Tesla C2050 GPU (dual-Tesla C2050) over the CPU implementation using Intel^© Core™i7 CPU 920 processor.

Conclusions

Harvesting the high performance of modern GPUs is a promising approach to accelerate pairwise statistical significance estimation for local sequence alignment.

GTP-dependent RNA 3'-terminal phosphate cyclase from the hyperthermophilic archaeon Pyrococcus furiosus

We discovered that the PF1549 gene in Pyrococcus furiosus encodes a very heat-stable RNA 3'-terminal phosphate cyclase (Pf-Rtc). Although all previously reported Rtc proteins are ATP-dependent enzymes, we found that Pf-Rtc requires GTP for its cyclase activity at 95 °C. Low-level activation of the enzyme was also observed in the presence of dGTP but not other dNTPs, indicating that the guanine base is very important for Pf-Rtc activity. We analyzed a series of GTP analogues and found that the conversion from GTP to GMP is important for Pf-Rtc activity and that an excess of GMP inhibits this activity. Gel-shift analysis clearly showed that the RNA-binding activity of Pf-Rtc is totally dependent on the linear form of the 3'-terminal phosphate, with an apparent K(d) value of 20 nm at 95°C. Furthermore, we found that Pf-Rtc may contribute to GTP-dependent RNA ligation activity through the PF0027 protein (a 2'-5' RNA ligase-like protein in P. furiosus). The possible roles of Pf-Rtc and the importance of terminal phosphate structures in RNA are discussed.