36 research outputs found
Accelerated Profile HMM Searches
Profile hidden Markov models (profile HMMs) and probabilistic inference methods have made important contributions to the theory of sequence database homology search. However, practical use of profile HMM methods has been hindered by the computational expense of existing software implementations. Here I describe an acceleration heuristic for profile HMMs, the “multiple segment Viterbi” (MSV) algorithm. The MSV algorithm computes an optimal sum of multiple ungapped local alignment segments using a striped vector-parallel approach previously described for fast Smith/Waterman alignment. MSV scores follow the same statistical distribution as gapped optimal local alignment scores, allowing rapid evaluation of significance of an MSV score and thus facilitating its use as a heuristic filter. I also describe a 20-fold acceleration of the standard profile HMM Forward/Backward algorithms using a method I call “sparse rescaling”. These methods are assembled in a pipeline in which high-scoring MSV hits are passed on for reanalysis with the full HMM Forward/Backward algorithm. This accelerated pipeline is implemented in the freely available HMMER3 software package. Performance benchmarks show that the use of the heuristic MSV filter sacrifices negligible sensitivity compared to unaccelerated profile HMM searches. HMMER3 is substantially more sensitive and 100- to 1000-fold faster than HMMER2. HMMER3 is now about as fast as BLAST for protein searches
Recommended from our members
Nuclear physics at PEP: First test and future plans
A test run of internal target nuclear physics at the PEP storage ring is described. The Time Projection Chamber (TPC-2..gamma.. detector) was used to detect the inelastically scattered electron and complete hadronic final state in the interaction of 14.5 GeV electrons with D/sub 2/, Ar and Xe gas targets. The data comprise mostly low-x low-Q/sup 2/ events, but some deep inelastic scattering as well. The future possibilities of a dedicated nuclear physics program at PEP are outlined. 15 refs., 25 figs
Recommended from our members
Exploratory gas-target experiment at PEP using the TPC/2. gamma. facility
A pilot nuclear-physics experiment at the PEP ring was performed by bleeding deuterium, argon, and xenon gases into the interaction region containing the TPC/2Y facility. The purpose of the experiment was to obtain information for the design of a dedicated nuclear-physics interaction region at PEP, as well as to obtain a first view of the physics that could be studied at such a facility. 1 ref., 4 figs