The protein structure prediction problem could be solved using the current PDB library Yang Zhang and Jeffrey Skolnick * Center of Excellence in Bioinformatics , University at Buffalo , 901 Washington Street , Buffalo , NY 14203 Communicated by R. Stephen Berry , University of Chicago , Chicago , IL , September 27 , 2004 -LRB- received for review November 10 , 2003 -RRB- For single-domain proteins , we examine the completeness of the structures in the current Protein Data Bank -LRB- PDB -RRB- library for use in full-length model construction of unknown sequences . To address this issue , we employ a comprehensive benchmark set of 1,489 medium-size proteins that cover the PDB at the level of 35 % sequence identity and identify templates by structure alignment . With homologous proteins excluded , we can always find similar folds to native with an average rms deviation -LRB- RMSD -RRB- from native of 2.5 with 82 % alignment coverage . These template structures often contain a significant number of insertions deletions . The TASSER algorithm was applied to build full-length models , where continuous fragments are excised from the top-scoring templates and reassembled under the guide of an optimized force field , which includes consensus restraints taken from the templates and knowledge-based statistical potentials . For almost all targets -LRB- except for 2 1,489 -RRB- , the resultant full-length models have an RMSD to native below 6 -LRB- 97 % of them below 4 -RRB- . On average , the RMSD of full-length models is 2.25 , with aligned regions improved from 2.5 to 1.88 , comparable with the accuracy of low-resolution experimental structures . Furthermore , starting from state-of-theart structural alignments , we demonstrate a methodology that can consistently bring template-based alignments closer to native . These results are highly suggestive that the protein-folding problem can in principle be solved based on the current PDB library by developing efficient fold recognition algorithms that can recover such initial alignments . s of December 30 , 2003 , 23,000 solved protein structures have been deposited in the Brookhaven Protein Data Bank -LRB- PDB -RRB- -LRB- 1 -RRB- . This number keeps increasing , with 300 new entries added each month . The size and completeness of the PDB is essential to the success of template-based approaches to protein structure prediction . These methods include comparative modeling -LRB- 2 , 3 -RRB- and threading -LRB- 4 7 -RRB- , which are designed to infer an unknown sequence 's structure based on solved , similarly folded protein structures in the PDB . Because an accurate theory for the prediction of protein structure on the basis of physical principles does not yet exist , comparative modeling threading approaches are the only reliable strategy for high-resolution tertiary structure prediction -LRB- 8 10 -RRB- . On the other hand , the percentage of new folds in these new entries , the topology of which has never been seen in the current PDB library , keeps decreasing -LRB- e.g. , the percentage of new folds was 27 % in 1995 but 5 % in 2001 -RRB- . The apparent saturation of new folds immediately raises an important question : At least for singledomain proteins , is the current structure library already complete enough to in principle solve the protein tertiary structure prediction problem at low-to-moderate resolutions ? By means of a variety of structure comparison tools -LRB- 11 14 -RRB- , this issue has been partially addressed by many authors -LRB- 15 20 -RRB- . It was demonstrated through systematic protein structure classification -LRB- 15 17 -RRB- that protein fold space is highly dense , and all solved PDB structures can be grouped into a very limited number of hierarchical families . Several authors -LRB- 17 20 -RRB- found that many proteins from different fold families share common substructures motifs and the protein fold space tends to be continuous . Especially , Kihara and Skolnick -LRB- 20 -RRB- showed that -LRB- after excluding homologues -RRB- for 90 % of single-domain proteins below 200 residues , there exists at least www.pnas.org cgi doi 10.1073 pnas .0407152101 A one structure -LRB- actually many -RRB- in the PDB having an rms deviation -LRB- RMSD -RRB- root to native below 4 with 79 % alignment coverage . This finding suggests that , at least at the level of structure alignments , the current PDB is almost a complete set of single-domain protein structures . However , the alignments identified from structure superposition usually contain numerous gaps . Starting from such alignments , it is still unknown in how many cases one could successfully build full-length models useful for biological annotation -LRB- 21 23 -RRB- . It could be that such models , while bearing a structural resemblance to the native state , might be sufficiently distorted that they could not be used as starting templates to build physically reasonable structures . Then our prior conclusion about the completeness of the PDB would be of purely academic interest , without practical ramifications . If , however , biologically useful models could be built , then the observation of the completeness of the PDB would have immediate practical value , not the least being that the protein structure prediction problem could in principle be solved on the basis of the current PDB library , if a sufficiently powerful fold recognition algorithm could be developed to recover such alignments -LRB- 21 , 23 , 25 -RRB- . It is the desire to explore this issue that provided the impetus of the work described here . Certainly , the requisite resolution required for different aspects of functional analysis -LRB- ligand docking , active site identification , etc. -RRB- may vary -LRB- 21 -RRB- . Recent investigations show that the active sites of enzymes can be successfully identified in about one-third of decoy structures of 3 - to 4 - RMSD -LRB- 24 -RRB- . Another important issue faced by protein structure modelers is related to the template refinement process . Until recently , proteinmodeling procedures usually drive the models farther from native , compared with initial template alignments -LRB- 8 10 , 26 -RRB- . It is a hard challenge to start from the structural -LRB- as opposed to threadingbased -RRB- alignments and improve upon them . To date , there has been no systematic demonstration that this was possible . The exploration of this issue provides the second motivation for this work . In this work , using a recently developed modeling algorithm , TASSER -LRB- 27 -RRB- , we examine both issues , by building full-length models from the templates identified by a state-of-the-art structure alignment algorithm -LRB- 20 -RRB- , and by demonstrating that in many cases the initial alignments are improved . To assess the generality of the conclusion , the strategy will be applied to a comprehensive , largescale PDB test set , with homologous proteins removed from the template library . BIOPHYSICS Methods The protein structure modeling procedure used in this work consists of template identification , force field construction , fragment assembly , and model selection . A flowchart is presented in Fig. 1 . Template Identification . Templates are identified from the solved structures in the PDB by structurally aligning the native of query proteins to templates by using our Structure Alignment -LRB- SAL algorithm -LRB- 20 -RRB- . The alignment is performed by a Needleman Abbreviations : RMSD , rms deviation ; SG , side-chain group ; FJC , freely jointed chain . * To whom correspondence should be addressed . E-mail : skolnick@buffalo.edu . 2005 by The National Academy of Sciences of the USA PNAS January 25 , 2005 vol . 102 no. 4 1029 1034 Fig. 2 . RMSD to native of the templates identified by the structure alignment program SAL -LRB- 20 -RRB- versus the alignment coverage . Fig. 1 . Overview of the TASSER structure prediction methodology that consists of template identification -LRB- here by structure alignment -RRB- , force field construction , structure assembly , and model selection . Results and Discussion Benchmark of Targets and Templates . For test proteins , we develop Wunsch dynamic program -LRB- 28 -RRB- with the score matrix defined as -LRB- 29 -RRB- 3 score -LRB- i , j -RRB- 20 -LRB- 1 dij 5 -RRB- , where dij is the spatial distance between the ith and jth residues after an initial guessed superposition . A number of iterations are performed until the structure alignment converges . Various gap penalties are implemented , and the best alignment is selected on the basis of the Z-score of the relative RMSD of two aligned proteins -LRB- 30 -RRB- . Force Field Construction . The force field used in TASSER includes four classes of terms : -LRB- i -RRB- C and side-chain group -LRB- SG -RRB- regularities correlations from the statistics of the PDB , -LRB- ii -RRB- propensities for predicted secondary structure from PSIPRED -LRB- 31 -RRB- , -LRB- iii -RRB- tertiary consensus contact distance restraints , and -LRB- iv -RRB- a protein-specific SG pair potential , both extracted from the identified multiple templates . The construction and implementation of the potentials in i and ii have been described -LRB- 32 , 33 -RRB- . The tertiary restraints in iii are constructed as done by our threading program PROSPECTOR 3 -LRB- 6 -RRB- ; and the details of the new pair potential in iv are in the Appendix . Having all of the energy terms , optimization for the combination weight factors was performed based on 100 training proteins -LRB- outside the benchmark test set used here -RRB- , each with 60,000 structure decoys , where we maximize the correlation between the energy and RMSD from native to the decoys -LRB- 32 -RRB- . Structure Assembly . Full-length models are constructed by reassem - a representative benchmark set of all single-domain structures in PDB with 41 200 aa . This target set contains 1,489 nonhomologous proteins having 448 , 434 , and 550 - , - , and - proteins , respectively -LRB- the other 57 are C - only targets or have irregular secondary structures -RRB- . The template library consists of 3,575 representative proteins from the PDB with a maximum 35 % pairwise sequence identity to each other ; all templates are taken from this library . Fig. 2 shows a summary of the resulting templates that have the highest RMSD Z-score obtained from SAL -LRB- 20 -RRB- for all 1,489 test proteins , where all templates with 25 % sequence identity to the target protein are excluded . The majority of targets have 70 % coverage and 4 - RMSD to native , which shows the significant denseness of protein structure space . The average sequence identity between template and target is 13 % in the aligned regions . Summary of Folding Results . Table 1 presents a summary of the folding results , where , for each protein , the two templates with the highest RMSD Z-score are used in TASSER . A detailed list of templates and final models , as well as the simulation trajectories , can be obtained by contacting J.S. . In Table 1 , the second column shows the best templates in the top two with the lowest RMSD to native in the aligned regions . On average , the RMSD to native is 2.51 with 82 % coverage . The final models show improvement over the initial template alignments . Over the same aligned regions , the average RMSD to native is reduced to 1.88 . Many low-resolution templates have been shifted by TASSER refinement to structures with an acceptable resolution for biochemical function annotation -LRB- 24 -RRB- . For example , there are 358 targets shifted from 3 - to 3 - RMSD to native , and 424 targets shifted from above to below 2 . For full-length models , almost all targets -LRB- except for 1a2kC and 1k5dB -RRB- have an RMSD to native below 6 for the best of the top five models , with an average rank of 1.7 ; 97 % have an RMSD 4 . The average RMSD to native is 2.25 , comparable with the accuracy of a low-resolution NMR or x-ray structure -LRB- 25 , 36 , 37 -RRB- . For the rank one cluster that has the highest structure density , the average RMSD to native is 2.35 . As a reference , we also list in the right hand columns of Table 1 the results of refined models from the publicly available comparative modeling program MODELLER -LRB- 2 , 3 -RRB- , using the same templates from SAL . The average RMSD of the best of top five models by MODELLER is 3.74 -LRB- average rank of 2.9 ; here , the rank of MODELLER models is decided by their objective function -RRB- , with 2.71 in the aligned regions . In general , successful Zhang and Skolnick bling the continuous fragments excised from the templates under the guide of the optimized force field . These fragments are elemental building blocks with internal conformations kept invariant during modeling . Residues in gapped regions are generated from an ab initio lattice modeling approach -LRB- 32 -RRB- . These regions also serve as linkage points for the rigid block movements . Conformational space is searched by using parallel hyperbolic Monte Carlo sampling -LRB- 34 -RRB- , where the tertiary topology varies by continuous translations and rotations of the rigid blocks . The magnitude of the move is scaled by the size of the blocks . Forty to fifty replicas are used in the simulations depending on protein size , and the trajectories in the 14 lowest temperature replicas are submitted to SPICKER -LRB- 35 -RRB- for clustering . The final models are combined from members of the structure clusters , ranked on the basis of cluster structure density . 1030 www.pnas.org cgi doi 10.1073 pnas .0407152101 Table 1 . Summary of templates by SAL SAL -LRB- 20 -RRB- and models built by TASSER TASSER or MODELLER -LRB- 2,3 -RRB- MODELLER Best * RMSD , Coverage , % NRMSD 6 NRMSD 5 NRMSD 4 NRMSD 3 NRMSD 2 NRMSD 1 2.510 82 1,489 1,472 1,369 1,064 498 46 Align 1.877 82 1,489 1,489 1,488 1,422 922 83 Top five 2.246 100 1,487 1,481 1,447 1,259 623 52 Top one 2.352 100 1,481 1,464 1,423 1,206 582 49 Align 2.708 82 1,462 1,395 1,255 1,008 520 37 Top five 3.740 100 1,326 1,195 984 647 300 20 Top one 4.318 100 1,202 1,060 841 551 244 15 * The template of the lowest RMSD to native . The best model in top five by TASSER and MODELLER with the RMSD calculated in the aligned region . The best model in top five where the RMSD is calculated for the entire chain . The first model where the RMSD is calculated for the entire chain . No . of targets with RMSD below the specified threshold -LRB- -RRB- . modeling in MODELLER has a stronger dependence on the template coverage than in TASSER . For example , if we look at those targets with 90 % coverage -LRB- 437 in total -RRB- , the average RMSDs of the full-length models by TASSER and MODELLER are fairly close , i.e. , 1.56 and 2.19 , respectively . However , for targets with initial alignment coverage 75 % -LRB- 386 in total -RRB- , the average RMSDs from native to models using TASSER and MODELLER are 2.92 and 6.05 , respectively , a significant difference . Overall , in 1,120 targets , TASSER models have a lower RMSD to native , where the alignment coverage in those targets is on average 81 % . For 102 targets , MODELLER does better , where the coverage is on average 92 % . Nevertheless , both procedures lead to the conclusion that the structure alignments produced by SAL can produce buildable models . Improvements of Initial Alignment . In Fig. 3 , we plot a detailed comparison of the final models with respect to the template in the aligned regions . In the majority of cases , TASSER models show obvious improvement -LRB- see Fig. 3a -RRB- . As shown in Fig. 3c , for targets having initial templates of aligned regions with an RMSD ranging from 2 to 3 , for 61 % of these cases , the models have at least a 0.5 - improvement ; and for initial alignments with an initial RMSD from 3 to 4 , in 49 % of cases , the final models improve by at least 1.0 . This result is partly because the force field takes consensus information from multiple templates , which can have higher accuracy confidence than that from individual templates . In TASSER , the local fragments from individual templates are rearranged under the guide of the force field , and the global topology can therefore move closer to native . This consensus information is further reinforced during the final model combination procedure of SPICKER clustering -LRB- 35 -RRB- . Another factor that contributes to the improvement is the protein-like energy terms , representing an optimal combination of statistical potentials , hydrogen bonds , and secondary structure predictions that lead to better side-chain and backbone structure packing than in the initial template-based alignments -LRB- 27 , 32 -RRB- . In Fig. 3 b and d , we also show the comparison between the models generated by MODELLER and the initial template alignments . In the majority of cases , MODELLER keeps the topology of models near the template , which is understandable because it was designed to optimally satisfy the spatial restraints from templates for homologous proteins -LRB- 2 , 3 -RRB- . However , in a few cases -LRB- 10 % -RRB- , the MODELLER models can be 1 worse than the initial templates . Modeling Unaligned Loop Regions . Because there is no spatial information provided from the template alignments , modeling the unaligned or loop regions is a hard , unsolved problem -LRB- 3 , 38 , 39 -RRB- . Here , we define an unaligned or loop region -LRB- including tails -RRB- as a piece of continuous sequence that has no coordinate assignment in the SAL template alignments . For each piece of those sequences , Zhang and Skolnick two types of accuracy are calculated -LRB- 3 -RRB- : RMSDlocal denotes the RMSD between native and the modeled loop with direct superposition of the unaligned region ; RMSDglobal is the RMSD between native and modeled loop after superposition of up to five neighboring stem residues on each side of the loop -LRB- for tails , the supposition is done on the side including five stem residues -RRB- . RMSDlocal measures the modeling accuracy of the local conformation , whereas RMSDglobal measures both the accuracy of the local conformation and the global orientation . There are in total 11,380 unaligned loop regions with size from 1 84 residues in the 1,489 targets . In Fig. 4 , we show the average values of RMSDlocal and RMSDglobal of TASSER and MODELLER models versus loop length L . In both cases , the accuracy of loop Fig. 3 . Comparison of initial and final alignments to the target structure . -LRB- a -RRB- Scatter plot of RMSD from native of the final models built by TASSER refinements versus RMSD from native of the initial template alignments identified by SAL . The same aligned regions are used in both RMSD calculations . -LRB- b -RRB- Similar data as in a , but the models are from MODELLER refinements . -LRB- c -RRB- Fraction of targets with an RMSD improvement `` d' ' by TASSER greater than some threshold value . Here d -LRB- RMSD of template -RRB- -LRB- RMSD of final model -RRB- , where both RMSDs are calculated over aligned regions . Each point in c is calculated with a bin width of 1 . -LRB- d -RRB- Similar data as in c , but the models are from MODELLER . PNAS January 25 , 2005 vol . 102 no. 4 1031 BIOPHYSICS Table 2 . Result of modeling of unaligned loop regions -LRB- > 4 residues -RRB- by TASSER and MODELLER TASSER MODELLER RMSDlocal * NRMSD 6 NRMSD 5 NRMSD 4 NRMSD 3 NRMSD 2 NRMSD 1 RMSD , 1,670 1,664 1,631 1,527 1,173 519 1.62 RMSDglobal * 1,386 1,199 924 549 193 18 4.34 RMSDlocal * 1,633 1,603 1,528 1,342 1,009 498 2.02 RMSDglobal * 1,011 800 507 244 64 10 6.59 * See text for definitions of RMSDlocal and RMSDglobal . No . of targets with a RMSD below the specified threshold -LRB- -RRB- . New Fold Targets in CASP5 . We revisit the new fold targets in the Fig. 4 . RMSDlocal -LRB- a -RRB- and RMSDglobal -LRB- b -RRB- of unaligned loop regions as a function of loop length -LRB- L -RRB- . TASSER and MODELLER models are denoted by triangles and circles , respectively . The lines connecting the points serve to guide the eye . The dashed line in b denotes an RMSDglobal cutoff of 7 . modeling decreases with increasing loop size . However , for all size ranges , the loops in TASSER models have lower average RMSDlocal and RMSDglobal . If we make a cutoff of RMSDglobal 7 in Fig. 4b , MODELLER generates reasonable models for unaligned loop regions of length up to 10 residues ; TASSER can have the same accuracy cutoff for the loops up to 28 residues . If using a lower RMSDglobal cutoff , the acceptable loop size in both approaches will decrease , and the difference between MODELLER and TASSER becomes smaller . Most unaligned regions loops in SAL alignments are of small size , which are relatively easier to model because of the limited configuration entropy . If we focus only on the unaligned loops of length greater than or equal to four residues , there are 1,675 cases with an average length of 8.8 residues . The distribution of modeling accuracy is summarized in Table 2 . Consistent with Fig. 4a , the distribution of RMSDlocal is quite close using TASSER and MODELLER . However , TASSER shows an obviously better control of loop orientations . For example , in one-third -LRB- 549 1,675 -RRB- of the cases , including loops and tails , TASSER generates models of RMSDglobal 3 , whereas the fraction of MODELLER models having RMSDglobal 3 is around one-seventh -LRB- 244 1,675 -RRB- . Representative Examples . In Fig. 5 , three representative examples of TASSER modeling results are provided : 1jm7A -LRB- an - protein -RRB- , 1b2iA -LRB- a - protein -RRB- , and 1xer -LRB- an - protein -RRB- . In all three , the template topologies in the core identified by SAL are quite similar to native -LRB- 5 -RRB- ; however , the local packing of the fragments and sometimes the termini are misoriented . Rearrangement using the TASSER force field results in a 2 improvement in the aligned region . In Fig. 6 , we show the predicted structure of 1k5dB , one of the two cases where TASSER failed to generate models with an RMSD to native below 6 . This is a Ran-binding protein -LRB- Ran-BP1 -RRB- , i.e. , chain B of the Ran-RanBP1 RanGAP protein complex -LRB- 40 -RRB- , which has a long tail interacting with chain A -LRB- Ran -RRB- . Because interactions with partner chains are not included , we failed to model the configuration of the tail ; this case results in a full-length RMSD of 7.8 to native . If we cut the first 22 residues in the N terminus associated with intermolecular interactions , the core region of the first model has a 1.4 - RMSD -LRB- Fig. 6c -RRB- . CASP5 folding experiment -LRB- 10 -RRB- , because , by definition , those targets putatively adopt a novel tertiary topology -LRB- 1 -RRB- . However , as shown in Table 3 , there are still more or less similar folds found in PDB , although the sequence identity is very low -LRB- 7 % on average -RRB- . The new entries in PDB released later than the CASP5 prediction season are not included in our template library . On average , these templates have shorter alignments -LRB- 73 % coverage -RRB- with higher RMSD -LRB- 3.8 -RRB- than those identified for the benchmark proteins . Still , acceptable models can be built from the initial template alignments using TASSER , with an average RMSD from native of 2.87 for the first predicted model . Fig. 5 . Representative examples showing the improvement of final models with respect to their initial template alignments . The left column is the superposition of the template alignments and native , and the right column is that for refined models and native . The thin lines are native structures ; the thick lines signify templates or final models . Blue to red runs from N to C terminus . The numbers in parentheses are the coverage of the templates or models on which the denoted RMSD has been calculated . 1032 www.pnas.org cgi doi 10.1073 pnas .0407152101 Zhang and Skolnick Fig. 6 . A representative example of targets where the final model has an RMSD to native 6 . The native structure of 1k5dB is shown in white with thin backbones ; the predicted model of the highest cluster density in blue with thick backbones . The red wire-frame in a denotes the partner chains in the Ran-RanBP1 RanGAP complex . -LRB- a -RRB- 1k5dB in the entire complex . -LRB- b -RRB- Native model superposition . -LRB- c -RRB- As in b but with the tail cut off . Concluding Remarks In this article , we examined the issues of whether all single-domain proteins are foldable based on the set of solved structures currently deposited in PDB -LRB- 1 -RRB- and whether the templates can be further improved by rearranging the fragments . We used our structure alignment algorithm , SAL -LRB- 20 -RRB- , to identify the best possible target template pairs , and then attempted to build and refine the fulllength model using the template assembly refinement algorithm TASSER -LRB- 27 -RRB- . This strategy was applied to a comprehensive PDB benchmark set of 1,489 medium-size , single-domain proteins . With homologous proteins excluded , similar folds can be found for all benchmark proteins , and the majority have a RMSD to native 4 over 70 % of their sequence . On average , the RMSD between template and native is 2.51 with 82 % alignment coverage . After TASSER , the average RMSD in the aligned region improves to 1.88 . The average global RMSD of unaligned loop tail regions -LRB- 4 residues -RRB- generated by TASSER is 4.3 . Almost all targets have at least one full-length model in the top five with an RMSD to native below 6 -LRB- 97 % are below 4 -RRB- . The average RMSD to native is 2.25 , comparable with the accuracy of low-to-moderateresolution experimental structures . In this sense , the answer to the question of completeness of the current PDB library for model construction of single-domain proteins is quite positive . Not only can physically reasonable models be built , but , starting from structural alignments , there is a significant improvement in many models ; 349 1,489 targets have a RMSD improvement of 1.0 in the aligned region . Thus , it is suggestive that the barrier to structure refinement noted in CASP5 -LRB- 10 -RRB- has been broken . In contrast to previous approaches , there are several reasons that contribute to the improvement of model quality compared with initial templates . First , the force field includes multiple sources of knowledge-based potentials and consensus tertiary restraints from multiple templates . The consensus spatial information usually has higher accuracy confidence than that of individual templates . Second , the combination of the different types of energy terms was optimized on the basis of a large-scale set of structure decoys -LRB- including 100 60,000 extrinsic targets structures -RRB- to yield an optimized potential that can provide better packing of the sidechains and peptides . This improvement occurs because of a better correlation between model quality and energy -LRB- the correlation coefficient between RMSD and the combined potential for test cases is 0.7 -RRB- . Finally , templates usually contain unphysical alignments because chain connectivity was not considered in the initial alignments . The reassembly procedure of TASSER that converts these unphysical alignments into physical models also contributes to the improvement in model quality relative to that of the initial template alignment . Unlike many other comparative modeling approaches , e.g. , MODELLER -LRB- 2 -RRB- , whose goal is to optimally satisfy the spatial restraints of an initial template , the relative orientation of template fragments , and therefore global topology in TASSER models , can change . On the other hand , the local conformation of the continuous fragments is kept rigid during the modeling procedure , which helps the models retain the accuracy of well-aligned regions from native and reduces the conformational entropy . For the more realistic situation where templates alignments are identified by using our threading program PROSPECTOR 3 -LRB- 6 -RRB- , the success rate for the same benchmark set of targets proteins is about two-thirds -LRB- where a foldable case is defined if one of top five full-length models has an RMSD below 6.5 to native -RRB- -LRB- 27 -RRB- . The results reported here highlight the urgent need to develop more efficient fold recognition algorithms that can provide acceptable templates for the remaining one-third of proteins , as well as better alignments to improve the overall quality of the predictions . In previous work -LRB- 27 -RRB- , we also observed an improvement in the models relative to their initial template alignments , but because threading models tend to be of poorer quality than those obtained from structural alignments , it could be argued that the results are not that significant -LRB- i.e. , the predicted models might be poorer than the best structural alignments -RRB- . Here we have demonstrated that even when Table 3 . Folding results for the new fold targets in CASP5 PDB ID 1h40 1iznC 1m6yB 1nyn 1o0uB 1o12C CASP5 ID T0170 T0162 T0172 T0181 T0187 T0186 Average 3 2 1 3 Length 69 168 101 111 165 36 108.3 RTali , * -LRB- % -RRB- 2.81 -LRB- 83 -RRB- 5.82 -LRB- 61 -RRB- 3.31 -LRB- 71 -RRB- 5.10 -LRB- 74 -RRB- 3.61 -LRB- 56 -RRB- 2.16 -LRB- 94 -RRB- 3.80 -LRB- 73 -RRB- RMali , 1.57 3.28 2.62 3.88 3.10 1.71 2.69 R5 , 1.70 -LRB- 2 -RRB- 3.05 -LRB- 2 -RRB- 2.79 -LRB- 1 -RRB- 3.94 -LRB- 1 -RRB- 3.56 -LRB- 1 -RRB- 1.82 -LRB- 1 -RRB- 2.81 -LRB- 1.3 -RRB- R1 , ŧ 1.79 3.31 2.79 3.94 3.56 1.82 2.87 ali , aligned . * RMSD to native of the best template -LRB- RT -RRB- and the alignment coverage . RMSD from the final model -LRB- RM -RRB- to native over initial aligned regions . RMSD to native for the best model in the top five . The number in parentheses is the rank of the best model . RMSD to native for the rank-one model . Zhang and Skolnick PNAS January 25 , 2005 vol . 102 no. 4 1033 BIOPHYSICS recent work -LRB- unpublished results -RRB- , we found that using the alignment from other software -LSB- in particular , CE -LRB- 14 -RRB- -RSB- as the initial alignment in SAL results in structure alignments of longer coverage and lower RMSD to native . Better structure alignment algorithms only serve to identify better templates from the PDB , which should result in better final models than reported here . Nevertheless , even now , it seems that the library of the solved protein structures is complete at the level of single-domain proteins . This structure completeness should have significant implications for both protein structure prediction and structural genomics -LRB- 21 , 23 -RRB- . Appendix : Multiple-Template-Based SG Pair Potential The protein-specific pair potential in our force field -LRB- term iv in the potential -RRB- is calculated from the identified multiple templates by using V i , j ln 0 if i , j are in aligned regions if i , j are in gapped regions , -LSB- 1 -RSB- q i , j Q i , j q 0 i , j , N L ln q i , j Q i , j q 0 i , j , N L Fig. 7 . Side-chain contact probability for random chains as a function of chain length . -LRB- a -RRB- Contact probability q0 -LRB- i , j , N L -RRB- for the ith and jth residues at given scaled contact-numbers , N L , versus the distance i j along the protein chain , from simulations of an FJC of 50 linkers with excluded volume . The Kuhn-length of the FJC is taken to be 11.4 , according to experimental single protein molecule stretching data -LRB- 43 -RRB- . The definition of excluded volume is introduced on the basis of the minimal observed C distance in the PDB , i.e. , no pair of linkers of the FJC could go closer than 4.5 . Contacts are defined based on a weighted average of C distances for contacting residues in the PDB , i.e. , a contact occurs if two linkers are closer than 7.27 . The curves are the least square fits to the power-law at each given N L. -LRB- b -RRB- The power versus the scaled contact-number N L. Data are from the FJC of different lengths , i.e. , 50 - , 100 - , 200 - , and 300-residue linkers , which correspond to the protein lengths of 150 , 300 , 600 , and 900 residues , respectively , because one Kuhnlength here corresponds approximately to three C C backbone lengths . The error bars denote the power-law fitting errors of a . The curve denotes the least 2.8 2.5 N L. square fit equation : the best structural alignments are used , we can often improve the models . This demonstration represents significant progress in the field . Because the average sequence identity between the target proteins and the best templates identified here is only 13 % , much lower than the `` twilight zone '' of sequence identity , correctly aligning the sequence to these templates will be a major challenge . This result is certainly true for our threading program , PROSPECTOR 3 , where for 90 % of targets , at least one correct fold can be identified in a large scale test ; however , only around 62 % were aligned correctly -LRB- 6 -RRB- . The results reported here provide a lower bound to the completeness and utility of the current PDB library . Certainly , the structure alignment program SAL is not perfect . It is not guaranteed to find the best structural alignment because the final alignment in this algorithm is sensitive to the initial guessed superposition . In 1 . Berman , H. M. , Westbrook , J. , Feng , Z. , Gilliland , G. , Bhat , T. N. , Weissig , H. , Shindyalov , I. N. & Bourne , P. E. -LRB- 2000 -RRB- Nucleic Acids Res . 28 , 235 242 . 2 . Sali , A. & Blundell , T. L. -LRB- 1993 -RRB- J. Mol . Biol . 234 , 779 815 . 3 . Fiser , A. , Do , R. K. & Sali , A. -LRB- 2000 -RRB- Protein Sci . 9 , 1753 1773 . 4 . Bowie , J. U. , Luthy , R. & Eisenberg , D. -LRB- 1991 -RRB- Science 253 , 164 170 . 5 . Jones , D. T. -LRB- 1999 -RRB- J. Mol . Biol . 287 , 797 815 . 6 . Skolnick , J. , Kihara , D. & Zhang , Y. -LRB- 2004 -RRB- Proteins 56 , 502 518 . 7 . Bryant , S. H. & Altschul , S. F. -LRB- 1995 -RRB- Curr . Opin . Struct . Biol . 5 , 236 244 . 8 . Moult , J. , Hubbard , T. , Fidelis , K. & Pedersen , J. T. -LRB- 1999 -RRB- Proteins 37 , Suppl . 3 , 2 6 . 9 . Moult , J. , Fidelis , K. , Zemla , A. & Hubbard , T. -LRB- 2001 -RRB- Proteins 45 , Suppl . 5 , 2 7 . 10 . Moult , J. , Fidelis , K. , Zemla , A. & Hubbard , T. -LRB- 2003 -RRB- Proteins 53 , Suppl . 6 , 334 339 . 11 . Taylor , W. R. , Flores , T. P. & Orengo , C. A. -LRB- 1994 -RRB- Protein Sci . 3 , 1858 1870 . 12 . Holm , L. & Sander , C. -LRB- 1995 -RRB- Trends Biochem . Sci . 20 , 478 480 . 13 . Gibrat , J. F. , Madej , T. & Bryant , S. H. -LRB- 1996 -RRB- Curr . Opin . Struct . Biol . 6 , 377 385 . 14 . Shindyalov , I. N. & Bourne , P. E. -LRB- 1998 -RRB- Protein Eng . 11 , 739 747 . 15 . Holm , L. & Sander , C. -LRB- 1996 -RRB- Science 273 , 595 603 . 16 . Murzin , A. G. , Brenner , S. E. , Hubbard , T. & Chothia , C. -LRB- 1995 -RRB- J. Mol . Biol . 247 , 536 540 . 17 . Orengo , C. A. , Michie , A. D. , Jones , S. , Jones , D. T. , Swindells , M. B. & Thornton , J. M. -LRB- 1997 -RRB- Structure 5 , 1093 1108 . 18 . Yang , A. S. & Honig , B. -LRB- 2000 -RRB- J. Mol . Biol . 301 , 665 678 . 19 . Harrison , A. , Pearl , F. , Mott , R. , Thornton , J. & Orengo , C. -LRB- 2002 -RRB- J. Mol . Biol . 323 , 909 926 . 20 . Kihara , D. & Skolnick , J. -LRB- 2003 -RRB- J. Mol . Biol . 334 , 793 802 . 21 . Skolnick , J. , Fetrow , J. S. & Kolinski , A. -LRB- 2000 -RRB- Nat . Biotechnol . 18 , 283 287 . 22 . Baxter , S. M. & Fetrow , J. S. -LRB- 2001 -RRB- Curr . Opin . Drug Discov . Devel . 4 , 291 295 . where q -LRB- i , j -RRB- is the number of SG contacts between residues i and j in all of the templates ; Q -LRB- i , j -RRB- is the number of templates that have both residue i and j aligned ; and q0 -LRB- i , j , N L -RRB- is the expected probability of contacts between residues i and j for a random chain of size L having a given total contact number N . The average ... is over all aligned pairs of -LRB- i , j -RRB- ; the shift sets the potential in gapped regions equal to the average magnitude of that in the aligned regions . To calculate q0 -LRB- i , j , N L -RRB- , we performed a Monte Carlo calculation of the freely jointed chain -LRB- FJC -RRB- model -LRB- 41 , 42 -RRB- with excluded volume . Fig. 7a shows the result for q0 -LRB- i , j , N L -RRB- from a random chain of 50 linkers . The contact probability of the FJC follows a power-law over more than three orders of magnitude : q0 -LRB- i , j , N L -RRB- i j -LRB- N/L -RRB- . Similar results are also obtained from the chains with 100 , 200 , and 300 linkers -LRB- data not shown -RRB- . In Fig. 7b , as a function of different scaled contact-numbers , -LRB- N L -RRB- , is presented . The data are well fit by : 2.8 2.5 N L. Integration of the contact probability results in a power-law : P -LRB- i , j -RRB- i j 1.84 , which coincides with the N/L q0 -LRB- i , j , N L -RRB- estimate from a Gaussian random chain with excluded volume that has a power of 9 5 -LRB- 41 -RRB- . This research was supported in part by National Institutes of Health Grants GM-37408 and GM-48835 . 23 . Baker , D. & Sali , A. -LRB- 2001 -RRB- Science 294 , 93 96 . 24 . Arakaki , A. K. , Zhang , Y. & Skolnick , J. -LRB- 2004 -RRB- Bioinformatics 20 , 1087 1096 . 25 . Marti-Renom , M. A. , Stuart , A. C. , Fiser , A. , Sanchez , R. , Melo , F. & Sali , A. -LRB- 2000 -RRB- Annu . Rev. Biophys . Biomol . Struct . 29 , 291 325 . 26 . Tramontano , A. & Morea , V. -LRB- 2003 -RRB- Proteins 53 , Suppl . 6 , 352 368 . 27 . Zhang , Y. & Skolnick , J. -LRB- 2004 -RRB- Proc . Natl. Acad . Sci . USA 101 , 7594 7599 . 28 . Needleman , S. B. & Wunsch , C. D. -LRB- 1970 -RRB- J. Mol . Biol . 48 , 443 453 . 29 . Gerstein , M. & Levitt , M. -LRB- 1996 -RRB- Proc . Int . Conf . Intell . Syst . Mol . Biol . 4 , 59 67 . 30 . Betancourt , M. R. & Skolnick , J. -LRB- 2001 -RRB- J. Comp . Chem . 22 , 339 353 . 31 . Jones , D. T. -LRB- 1999 -RRB- J. Mol . Biol . 292 , 195 202 . 32 . Zhang , Y. , Kolinski , A. & Skolnick , J. -LRB- 2003 -RRB- Biophys . J. 85 , 1145 1164 . 33 . Kolinski , A. & Skolnick , J. -LRB- 1998 -RRB- Proteins 32 , 475 494 . 34 . Zhang , Y. , Kihara , D. & Skolnick , J. -LRB- 2002 -RRB- Proteins 48 , 192 201 . 35 . Zhang , Y. & Skolnick , J. -LRB- 2004 -RRB- J. Comput . Chem . 25 , 865 871 . 36 . Glusker , J. P. -LRB- 1994 -RRB- Methods Biochem . Anal . 37 , 1 72 . 37 . Wagner , G. , Hyberts , S. G. & Havel , T. F. -LRB- 1992 -RRB- Annu . Rev. Biophys . Biomol . Struct . 21 , 167 198 . 38 . Mosimann , S. , Meleshko , R. & James , M. N. -LRB- 1995 -RRB- Proteins 23 , 301 317 . 39 . Martin , A. C. , MacArthur , M. W. & Thornton , J. M. -LRB- 1997 -RRB- Proteins 29 , Suppl . 1 , 14 28 . 40 . Seewald , M. J. , Korner , C. , Wittinghofer , A. & Vetter , I. R. -LRB- 2002 -RRB- Nature 415 , 662 666 . 41 . Doi , M. & Edwards , S. F. -LRB- 1986 -RRB- The Theory of Polymer Dynamics -LRB- Clarendon Press , Oxford -RRB- . 42 . Zhang , Y. , Zhou , H. J. & Ouyang , Z. C. -LRB- 2001 -RRB- Biophys . J. 81 , 1133 1143 . 43 . Rief , M. , Fernandez , J. M. & Gaub , H. E. -LRB- 1998 -RRB- Phys . Rev. Lett . 81 , 4764 4767 . 1034 www.pnas.org cgi doi 10.1073 pnas .0407152101 Zhang and Skolnick

O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O B O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O B O B O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O B O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O B O O O O O O O O O O B O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O B O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O B O O O O O O O B B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O