Method for prediction of protein function from sequence using the sequence-to-structure-to-function paradigm with application to Glutaredoxins/Thioredoxins and T 1 Ribonucleases 1 1Edited by F. Cohen
Journal of Molecular Biology
The practical exploitation of the vast numbers of sequences in the genome sequence databases is crucially dependent on the ability to identify the function of each sequence. Unfortunately, current methods, including global sequence alignment and local sequence motif identi®cation, are limited by the extent of sequence similarity between sequences of unknown and known function; these methods increasingly fail as the sequence identity diverges into and beyond the twilight zone of sequence
... . To address this problem, a novel method for identi®cation of protein function based directly on the sequence-to-structure-to-function paradigm is described. Descriptors of protein active sites, termed"fuzzy functional forms" or FFFs, are created based on the geometry and conformation of the active site. By way of illustration, the active sites responsible for the disul®de oxidoreductase activity of the glutaredoxin/ thioredoxin family and the RNA hydrolytic activity of the T 1 ribonuclease family are presented. First, the FFFs are shown to correctly identify their corresponding active sites in a library of exact protein models produced by crystallography or NMR spectroscopy, most of which lack the speci-®ed activity. Next, these FFFs are used to screen for active sites in lowto-moderate resolution models produced by ab initio folding or threading prediction algorithms. Again, the FFFs can speci®cally identify the functional sites of these proteins from their predicted structures. The results demonstrate that low-to-moderate resolution models as produced by state-of-the-art tertiary structure prediction algorithms are suf®cient to identify protein active sites. Prediction of a novel function for the gamma subunit of a yeast glycosyl transferase and prediction of the function of two hypothetical yeast proteins whose models were produced via threading are presented. This work suggests a means for the large-scale functional screening of genomic sequence databases based on the prediction of structure from sequence, then on the identi®cation of functional active sites in the predicted structure.