the central paradigm of protein science suggests that protein functions are directly controlled by protein structures . with the increasing number of solved protein structures a number of semi or fully automated structure comparison methods have been developed based on methodologies like alignment of secondary structure elements ( 13 ) , environmental profiles ( 4 ) and distance measure matrices ( 5 ) . most of these methods use regular secondary structure information in their algorithms . by analyzing local protein structures , many groups have found recurring short structural motifs also called structural alphabet ( sa ) spanning structural space ( 68 ) . these short motifs represent local structure variations in protein space upon which backbone model of most proteins can be built . they have been shown to be informative to analyze protein structures ( 9 ) and have been used in structure prediction ( 10 ) , backbone reconstruction ( 11,12 ) and loop modeling ( 13 ) . we present a web - based service called protein block expert ( pbe ) for protein structure comparison and analysis using a sa of 16 pentapeptide structural motifs known as  protein blocks ' ( pbs ) ( 14,15 ) . a protein structure can be encoded into sequence of pbs by sliding an overlapping window of five residues . hence , simplified 1d representation of protein structure can be used just like amino acid sequence analysis to find similarity , dissimilarity and relationship among proteins in terms of structure . its concept is similar to sa - search ( 18 ) web server , but differs greatly as it uses a specialized sa substitution matrix derived on the basis of aligned homologous proteins present in a large database of phylogeny and alignment of homologous protein structures ( pali ) ( 19,20 ) . applications and validation of such a matrix have been demonstrated ( m. tyagi , v. s. gowri , n. srinivasan , a. g. de brevern , b. offmann , manuscript submitted ) . pbe is not only a service to find structural similarities between proteins or a mining tool for recognizing the fold of a protein structure , it also provides an interface to a database to study proteins in terms of pbs at the levels of superfamily and family . pbe provides the following features to the user : 
 a tool to encode protein structure into pbs sequence.structure comparison between a pair of proteins using pb description using both local and global alignment algorithms.mining a databank , which is derived from scop , for proteins with similar fold.access to a database of preprocessed pb sequences and pairwise alignments at family and superfamily levels . a tool to encode protein structure into pbs sequence . structure comparison between a pair of proteins using pb description using both local and global alignment algorithms . mining a databank , which is derived from scop , for proteins with similar fold . access to a database of preprocessed pb sequences and pairwise alignments at family and superfamily levels . the set of pbs consists of 16 structural motifs of each five residue long ( 14,15 ) . each of the pbs is represented by a vector of eight , dihedral angles associated with five consecutive c atoms and the 16 pbs are denoted by the letters a , b ,  , p. encoding of protein 3d structure into sequence of pbs as implemented in our server is a two - step process . first , protein backbone is encoded into sequence of ( , ) angles calculated from backbone atomic positions . second , an overlapping window of five c atoms , i.e. vector of eight ( , ) angles is moved along the backbone . pbs for each window is assigned on the basis of smallest dissimilarity measure called root mean square deviation on angular values or r.m.s.d.a . ( 21 ) between observed ( , ) values in the window and the standard dihedral angles for various pbs . it accepts a structure as an input and lists the sequence of pbs in the structure . analysis of sequence of pbs using classical amino acid sequence alignment algorithms allows us to explore possibility of finding structural similarities between two proteins using reduced complexity of protein structure . it allows user to compare two proteins using simple dynamic programming ( dp ) algorithm by aligning two pb sequences using our pb substitution matrix . the substitution table used in our study was derived by re - encoding in terms of pbs the structurally aligned homologous proteins present in the pali database ( 19,20 ) . the detailed description of calculation , discussion on pb substitution matrix and proposed applications are reported elsewhere ( m. tyagi , v. s. gowri , n. srinivasan , a. g. de brevern , b. offmann , manuscript submitted ) . indeed , local structural similarities between two uploaded protein structures are found using pb sequence alignment . this approach has already been successfully benchmarked and compared to standard flexible alignment methods like dali ( 5 ) and rigid body superposition methods like stamp ( 25 ) where > 75% of structurally equivalent residues in our pb alignment method overlapped with those identified with these standard methods . moreover , careful inspection of aligned coordinates from pb - alignc after aligning pbs indeed shows identical and in some cases , lower r.m.s.d . values than dali ( m. tyagi , v. s. gowri , n. srinivasan , a. g. de brevern , b. offmann , manuscript submitted ) . these results are expected to improve by using more robust dp algorithm combined with optimized gap penalty . interestingly , in the same study we have shown how pb alignment method is able to pick up subtle similarities at local level between two proteins . hence pb alignment is providing both local and global flavors of dp algorithms . in pbe - alignc , the user is required to upload two protein structures in pdb format . after transforming the 3d structures into 1d pb sequences and the latter if the uploaded protein structures have more than one chain , option to select any one of the possible pair for alignment is presented . the output displays the aligned pb sequences along with the information like length of proteins , alignment length , best fit superposition r.m.s.d . pb alignment is transformed into amino acid alignment to define equivalent regions required by profit and further iterations are done to obtain best fit r.m.s.d . the server provides the possibility to download the initial pb and corresponding amino acid alignments in fasta format as well as the superimposed coordinates between the two structures . as pbe requires only backbone atoms to generate pb sequence and is independent of residues , the user can upload anonymous protein structures by changing all residues to any one kind and giving only coordinates of backbone atoms in the pdb format . hence newly solved structures can be easily analyzed without making them public . pbe - alignm uses a database of protein domains derived from scop database ( 27 ) . protein structures were extracted from scop 1.65 via the astral ( 28 ) server using a sequence identity cutoff of 95% ( scop95 ) with 9392 domains . these domains were encoded into pb sequences and are made available for user to query at family and superfamily levels in pbe - sadb database . further , an extensive all - against - all pairwise pb sequence alignments between all 7195 domains were generated using dp and our pb substitution matrix . protein domains in scop95 having any chain breaks were not considered for pb sequence alignment process . pairwise alignments within each of seven major structural classes from scop95 , which amounts 5 405 433 alignments , are featured in pbe - sadb database where option is provided to the user to view / download pairwise pb alignments at the level of family or superfamily . each pb alignment in the generated databank had raw score given by dp algorithm . to remove the dependence of this value on the length of the two proteins , the score is also provided in the normalized form by dividing the raw score by the length of the alignment including gaps . this normalized score from global alignment algorithm is used to rank alignments during the following analysis . in the first study we analyzed the efficiency of our method to discriminate between various scop classes or in other words what is the overlap of scores between classes based on pb sequence alignment . this question is important since 1d representation of protein structure using pb sequence lacks topological information , which can create confusion due to identical linear sequence of regular structures in two proteins having different topologies . a dataset of 1500 protein domains was selected randomly from scop at 95% keeping the relative proportion of seven major classes same as in original databank . each time one domain was selected and was queried against the databank to find top 10 ranking pb alignments against the given query and statistics was calculated for true hits at each rank position . analysis of the distribution of true hits shows that 85.9% of them are at first rank and a hit rate of 98.2% is achieved when first 10 ranking alignments are considered ( data not shown ) . it should be noted that the value increases from 85.9 to 93% when the same analysis is performed on the 7267  7267 pairwise alignments . a confusion matrix between seven scop classes matrix is populated simply based on criteria if query protein and first rank protein have same class or not . alpha plus beta class was most confused class with 76.2% cases finding itself at first rank . beta class was most well - behaved class with accuracy of 94.4% followed by alpha beta and alpha class . low accuracy rate of multi - domain and membrane class can be attributed to very low number of proteins of this present in our dataset . in a second study , we assessed how well a pb alignment can extract protein of similar fold from a databank within given a class . a jackknife approach ( as done in previous analysis , cf . infra ) was applied to calculate statistics for identifying true fold of a protein as defined by scop at various levels . it is noteworthy that 81.3% of true hits are from first rank while 89.3% true hits are within top 10 ranking alignments . further efficiency of mining similar folds within each seven major classes was studied and results are reported in table 2 . for each class , hit rate was calculated at three different levels , top10 , top5 and top1 where first 10 , 5 and first ranking alignments were considered , respectively . the ability of our method to extract same scop fold within top10 level vary from a hit rate of 86.1% for alpha class to 93.6% for alpha / beta class . similarly , at top1 level , the hit rate varies from 70% ( small protein class ) to 88.4% ( alpha / beta class ) . consistent good level of hit rates across various classes to mine similar fold using pb alignment method gives support to basic ability of the method and quality of the substitution matrix . these results hence illustrates that the use of pb substitution matrix with simple dp algorithm along with nave scoring function is efficient to extract proteins sharing structural similarities from large dataset . pbe - alignm provides this facility for mining structural similarities from a databank using a reduced representation of protein structures . user can upload a protein structure in pdb format and can decide against which databank the structure is to be queried . pbe gives option to select local or global alignment algorithm , setting up parameters like minimum length of proteins against which query should be aligned . option is also given to decide if the user would like to align against whole databank or with some specific scop class proteins . we have created two databases of protein structures and are grouped under the pbe - sadb facility . first is a database of 9392 protein domains extracted from scop95 that were translated into pb sequences . second is a database of all possible pairwise pb sequence alignments within each scop class . an interface gives option for querying any one of these two databases at superfamily or family level by entering appropriate scop code . list of all family and superfamily codes and their description present in our database is available in the help section . in addition pb sequences or alignments can also be accessed by specifying a pdb i d of a protein . because pdb was filtered for 95% sequence identity cutoff , the list of the available pdb structures can also be checked in help section . facility to query and download pb sequences or pb alignments at family or superfamily level is expected to be of great help in studying protein structure conservation . this can also aid studies on variations in homologous proteins in terms of structural alphabets , which might provide better insight into sequence to structure relationship . analysis of pb alignments to study conservation or variability of local structures is expected to provide better understanding of relationship between structure and function of homologous proteins . pairwise pb alignments for each given class were calculated using 32 processor ibm aix52 machine . web server front end and back end processing are handled using html , cgi and perl scripting along with java ( pbs encoding of protein ) and c ( dp ) programs . job requests in pbe - alignm are queued and provided a randomly generated job - id that guarantees the inaccessibility of jobs to other users of the servers . pbe server is maintained on a linux - based single processor machine and is accessible at . decrease in the complexity of protein space from three dimensions to one dimension with the combination of sequence analysis methods to study protein structure has opened a simple and exciting way of looking at protein structure space . initial results of mining similar fold structures in databases and finding local structural similarities between proteins has been a promising start though far from exploiting full potential of such methodology . low confusion rate across various scop classes and high efficiency rate to mine similar fold proteins from large database based on nave scoring scheme indicates that pb alignment method is efficient enough to discriminate between different topologies despite lack of topological information . this success can be attributed to both , efficiency of pbs to represent local structural properties in more refined form compared to simple sse representation and quality of substitution matrix . sequence of pbs between regular sses and their alignment or misalignment might be playing important role in discriminating between correct and incorrect hits . pairwise comparison of proteins using pbe - alignc performed reasonably well when compared to standard methods like dali and this was further validated on a large - scale basis here ( 7195  7195 pairwise alignments ) . structure alignment using pbs has also shown its efficiency to locate subtle similarities at local level and to very efficiently mine for local structural similarities from large structural databases . though , pb alignment is expected to be very advantageous in cases of distantly related proteins where residue finally , because prediction of protein backbone in terms of pb sequences is possible from amino acid sequence ( 14 ) , this work opens up interesting perspectives for large - scale structural annotation of genomic data . if a given query and extracted alignment have same fold , a hit is counted at that position . mining scop for similar structures using pb alignment confusion matrix between true ( vertical ) and predicted ( horizontal ) scop classes . mining scop for similar structures using pb alignment hit rates ( in percentage ) for identifying similar fold within each scop classes . are given rates that take into account top10 , top5 and top1 ranking alignments .