multiple sequence alignments of protein families are commonly represented as hidden markov models ( hmms ) ( 14 ) or position - specific scoring matrices ( pssms ) ( 5 ) . tremendous power of such protein profiles in enabling the detection of remote homologues is well known . for example , psi_blast ( 5 ) uses pssm generated , at the end of every round of search , as an input to the next round . such a use of dynamic pssms is known to be extremely effective in the detection of distant relatives of query in the sequence database . rps_blast and impala use a database of static pssms corresponding to homologous proteins and enables rapid match of the query sequence with various pssms in the database ( 6,7 ) . programs employing pssm matching algorithms are generally faster than the hmm matching programs and are commonly used in large - scale analyses ( 8,9 ) . recently we have focused on an important and sensitive feature of pssms ( 10 ) . it is the use of a reference sequence in the construction of a pssm starting from a multiple sequence alignment . reference sequence should be any one of the sequences involved in the multiple sequence alignment . basically , pssm integrates two kinds of information at every alignment position of a multiple sequence alignment . ( ii ) potential of replacement of the residue in the reference sequence with each one of the 20 residue types . the reference sequence corresponds to the query in the case of psi_blast . for rps_blast searches , reference sequence for generating a pssm hence , the pssms generated for a multiple sequence alignment with different homologues as reference sequences will be different . it has been shown convincingly that the sensitivity of the pssm is strongly dependent on the choice of the reference sequence . during rps_blast searches use of multiple pssms corresponding to different reference sequences for a given alignment results in remarkably improved specificity , sensitivity and error rate compared with the use of single pssm corresponding to an alignment and even hmm ( 10 ) . in this paper , we report the ready availability of multiple pssms for every multiple sequence alignment corresponding to large datasets of sequence and structural families . a web interface allows convenient use of rps_blast to search these datasets of multiple pssms and analysis of results . the set of protein domain sequence families corresponds to protein families database of alignments and hmms ( pfam version 17.0 ) , which consists of 7868 families ( 11,12 ) . the dataset of structural domain families has been obtained from the latest update of the database on phylogeny and alignment of homologous protein structures ( pali release 2.4 ) corresponding to 2625 families ( 13,14 ) . the list of protein domain families of known structures in pali release 2.4 has been obtained from 1.67 release of structural classification of proteins ( scop ) database ( 15,16 ) . the multiple sequence alignments in pali families have been derived on the basis of three - dimensional structural superposition of homologues using structural alignment of multiple proteins ( stamp , version 4.2 ) program ( 17 ) . homologues of yet unknown structure corresponding to the pali families have been identified by mining the universal protein resource ( uniprot ) database ( 18,19 ) using psi_blast ( e and h value cut - offs are 10 and 10 respectively ; cut - off for both query length coverage in the alignment and sequence identity are fixed at 30% ) . generation of multiple pssms with two very closely related sequences as references is unlikely to be worthwhile and it will unnecessarily add to the computer time while searching the pssm database . hence choice of reference sequences from a multiple sequence alignment should be optimal so that it maximizes the sensitivity during search in the pssm database and minimizes the computer time . hence we clustered sequences in every family of protein domains and identified a set of most disperse sequences within each family . blastclust was used for the generation of clusters of non - redundant sequences from the set of homologues of the family sequences ( 5 ) . in order to identify optimal level of clustering following exercise was performed : the clustering was performed at cut - offs of 30100% sequence identity for 250 randomly picked scop families . the sequences in pali database were then queried against the multi pssm database at each cut - off value . the number of true positives ( superfamily- or fold - related sequences ) and false positives ( class - related or across class connections ) were found for each cut - off value to assess the sensitivity and specificity of the profiles generated at each cut - off . it was found that at a sequence identity cut - off of 70% , the sensitivity of the method was equal to the sensitivity at 100% sequence identity cut - off whereas the number of profiles had decreased from 25  000 at 100% cut - off to 14  000 at 70% cut - off . at a sequence identity cut - off of 50% , the sensitivity is 75% of the sensitivity at 100% cut - off . apart from this point , considering the time taken for search that is dependent on the number of pssms in the database a cut - off of 50% was chosen as optimal . in our previous work , we have shown that in the multiple pssm approach , the false positive rate is 2% at an e - value cut - off of 10 ( 10 ) . the false positive rate remains the same after clustering the sequences at a cut - off of 50% sequence identity . for the structural families , the multiple sequence alignments were biased towards the structure - based alignment of the family . in principle , generation of profiles can be carried out using sequence - based alignments alone and it would make no difference to the quality of profiles if the sequences are not very divergent . in case of very divergent sequences , the structural alignments are more accurate than sequence alignments and thus in such cases , using the structural information will lead to more robust profiles . after identifying the reference sequences present in the multiple sequence alignment using the clustering algorithm , the following procedure is used for the generation of pssms . the multiple sequence alignment and the reference sequence are given as inputs to a psi_blast run to iterate against a database of sequences present in the input multiple sequence alignment . in such iterative searches , any hit identified is already present in the multiple sequence alignment fed to the search program . hence , the profile generated corresponds to the multiple sequence alignment fed to the psi_blast program . the current version of mulpssm database consists of 33  000 and 38  000 pssms corresponding to 7868 sequence and 2625 structural families , respectively . we believe , based on our earlier analysis ( 10 ) , that the use of multiple profiles for each family can lead to better annotation of genome sequences than is provided by hmm based or single profile based searches . for example , the protein gi|23508377 from plasmodium falciparum is annotated as a hypothetical protein in the ncbi genome database . using psi_blast search or pfam hmm based searches no known domain could be detected in the sequence with significant e - values . using the multiple pssm approach the e - value is 10 with a sequence identity of 22% over the alignment length of about 80 residues . in the mulpssm website ( ) , users can select datasets of pssms corresponding to either sequence families or structural families or both . although we recommend an e - value cut - off of 10 , the choice of e - value is left to the users . the results are represented in pictorial form , with ascii characters , to show the presence of various domains that may be present in the query sequence . links are provided to the alignment corresponding to the pssm with the best e - value , list of all the profiles within and across families that are hit to the query , links to the details of the families ( figure 1 ) . lists of all the sequence and structural families are provided in the mulpssm site with alignments and indication of reference sequences for each family . reference sequences used in every family are shown highlighted and the user can download all the pssms corresponding to a given family for use in a local machine in a rps_blast search . the complete set of all the pssms may be obtained from the authors upon request . this set is expected to be invaluable for rapid and sensitive genome - wide domain assignments . the main window gives the alignment of the families found as hits in a semi - graphical layout . the main window has a link to the details of the multiple hits for a family ( indicated by blue encircle ) which opens in a new window ( indicated by a blue arrow ) . similarly , the details of the family can be obtained by clicking on the name of the family ( indicated by a red encircle ) . the detailed view of the alignments and other features of a family hit ( e.g. the fraction of pssms of a family found as hits ) can help in assessing the accuracy of a hit .
representation of multiple sequence alignments of protein families in terms of position - specific scoring matrices ( pssms ) is commonly used in the detection of remote homologues . 
 a pssm is generated with respect to one of the sequences involved in the multiple sequence alignment as a reference . 
 we have shown recently that the use of multiple pssms corresponding to an alignment , with several sequences in the family used as reference , improves the sensitivity of the remote homology detection dramatically . 
 mulpssm contains pssms for a large number of sequence and structural families of protein domains with multiple pssms for every family . 
 the approach involves use of a clustering algorithm to identify most distinct sequences corresponding to a family . with each one of the distinct sequences as reference , 
 multiple pssms have been generated . 
 the current release of mulpssm contains 33  000 and 38  000 pssms corresponding to 7868 sequence and 2625 structural families . 
 a rps_blast interface allows sequence search against pssms of sequence or structural families or both . 
 an analysis interface allows display and convenient navigation of alignments and domain hits . 
 mulpssm can be accessed at .