Subset Of Residues Research Articles

Crystal Structure of Helicobacter pylori Formamidase AmiF Reveals a Cysteine-Glutamate-Lysine Catalytic Triad

Helicobacter pylori AmiF formamidase that hydrolyzes formamide to produce formic acid and ammonia belongs to a member of the nitrilase superfamily. The crystal structure of AmiF was solved to 1.75A resolution using single-wavelength anomalous dispersion methods. The structure consists of a homohexamer related by 3-fold symmetry in which each subunit has an alpha-beta-beta-alpha four-layer architecture characteristic of the nitrilase superfamily. One exterior alpha layer faces the solvent, whereas the other one associates with that of the neighbor subunit, forming a tight alpha-beta-beta-alpha-alpha-beta-beta-alpha dimer. The apo and liganded crystal structures of an inactive mutant C166S were also determined to 2.50 and 2.30 A, respectively. These structures reveal a small formamide-binding pocket that includes Cys(166), Glu(60), and Lys(133) catalytic residues, in which Cys(166) acts as a nucleophile. Analysis of the liganded AmiF and N-carbamoyl d-amino acid amidohydrolase binding pockets reveals a common Cys-Glu-Lys triad, another conserved glutamate, and different subsets of ligand-binding residues. Molecular dynamic simulations show that the conserved triad has minimal fluctuations, catalyzing the hydrolysis of a specific nitrile or amide in the nitrilase superfamily efficiently.

Electron Paramagnetic Resonance Spectroscopy of Site-directed Spin Labels Reveals the Structural Heterogeneity in the N-terminal Domain of ApoA-I in Solution

Apolipoprotein A-I (apoA-I) is the major protein constituent of high density lipoprotein (HDL) and plays a central role in phospholipid and cholesterol metabolism. This 243-residue long protein is remarkably flexible and assumes numerous lipid-dependent conformations. Consequently, definitive structural determination of lipid-free apoA-I in solution has been difficult. Using electron paramagnetic spectroscopy of site-directed spin labels in the N-terminal domain of apoA-I (residues 1-98) we have mapped a mixture of secondary structural elements, the composition of which is consistent with findings from other in-solution methods. Based on side chain mobility and their accessibility to polar and non-polar spin relaxers, the precise location of secondary elements for amino acids 14-98 was determined for both lipid-free and lipid-bound apoA-I. Based on intermolecular dipolar coupling at positions 26, 44, and 64, these secondary structural elements were arranged into a tertiary fold to generate a structural model for lipid-free apoA-I in solution.

Structure and Dynamics of Pin1 During Catalysis by NMR

The link between internal enzyme motions and catalysis is poorly understood. Correlated motions in the microsecond-to-millisecond timescale may be critical for enzyme function. We have characterized the backbone dynamics of the peptidylprolyl isomerase (Pin1) catalytic domain in the free state and during catalysis. Pin1 is a prolyl isomerase of the parvulin family and specifically catalyzes the isomerization of phosphorylated Ser/Thr-Pro peptide bonds. Pin1 has been shown to be essential for cell-cycle progression and to interact with the neuronal tau protein inhibiting its aggregation into fibrillar tangles as found in Alzheimer's disease. 15N relaxation dispersion measurements performed on Pin1 during catalysis reveal conformational exchange processes in the microsecond timescale. A subset of active site residues undergo kinetically similar exchange processes even in the absence of a substrate, suggesting that this area is already “primed” for catalysis. Furthermore, structural data of the turning-over enzyme were obtained through inter- and intramolecular nuclear Overhauser enhancements. This analysis together with a characterization of the substrate concentration dependence of the conformational exchange allowed the distinguishing of regions of the enzyme active site that are affected primarily by substrate binding versus substrate isomerization. Together these data suggest a model for the reaction trajectory of Pin1 catalysis.

Characterization of Folding Intermediates of a Domain-Swapped Protein by Solid-State NMR Spectroscopy

We have employed two-dimensional solid-state NMR to study structure and dynamics of insoluble folding states of the domain-swapped protein Crh. Starting from the protein precipitated at its pI, conformational changes due to a modest temperature increase were investigated at the level of individual residues and in real-time. As compared to the crystalline state, Crh pI-precipitates exhibited a higher degree of molecular mobility for several regions of the protein. A rigidly intact center was observed including a subset of residues of the hydrophobic core. Raising the temperature by 13 K to 282 K created a partially unfolded intermediate state that was converted into beta-sheet-rich aggregates that are mostly of spherical character according to electron microscopy. Residue-by-residue analysis indicated that two out of three alpha-helices in aggregated Crh underwent major structural rearrangements while the third helix was preserved. Residues in the hinge region exhibited major chemical-shift changes, indicating that the domain swap was not conserved in the aggregated form. Our study provides direct evidence that protein aggregates of a domain-swapped protein retain a significant fraction of native secondary structure and demonstrates that solid-state NMR can be used to directly monitor slow molecular folding events.

Exhaustive assignment of compositional bias reveals universally prevalent biased regions: analysis of functional associations in human and Drosophila

BackgroundCompositionally biased (CB) regions are stretches in protein sequences made from mainly a distinct subset of amino acid residues; such regions are frequently associated with a structural role in the cell, or with protein disorder.ResultsWe derived a procedure for the exhaustive assignment and classification of CB regions, and have applied it to thirteen metazoan proteomes. Sequences are initially scanned for the lowest-probability subsequences (LPSs) for single amino-acid types; subsequently, an exhaustive search for lowest probability subsequences (LPSs) for multiple residue types is performed iteratively until convergence, to define CB region boundaries. We analysed > 40,000 CB regions with > 20 million residues; strikingly, nine single-/double- residue biases are universally abundant, and are consistently highly ranked across both vertebrates and invertebrates. To home in subpopulations of CB regions of interest in human and D. melanogaster, we analysed CB region lengths, conservation, inferred functional categories and predicted protein disorder, and filtered for coiled coils and protein structures. In particular, we found that some of the universally abundant CB regions have significant associations to transcription and nuclear localization in Human and Drosophila, and are also predicted to be moderately or highly disordered. Focussing on Q-based biased regions, we found that these regions are typically only well conserved within mammals (appearing in 60–80% of orthologs), with shorter human transcription-related CB regions being unconserved outside of mammals; they are also preferentially linked to protein domains such as the homeodomain and glucocorticoid-receptor DNA-binding domain. In general, only ~40–50% of residues in these human and Drosophila CB regions have predicted protein disorder.ConclusionThis data is of use for the further functional characterization of genes, and for structural genomics initiatives.

Turnover of the Active Fraction of IRS1 Involves Raptor-mTOR- and S6K1-Dependent Serine Phosphorylation in Cell Culture Models of Tuberous Sclerosis

The TSC1-TSC2/Rheb/Raptor-mTOR/S6K1 cell growth cassette has recently been shown to regulate cell autonomous insulin and insulin-like growth factor I (IGF-I) sensitivity by transducing a negative feedback signal that targets insulin receptor substrates 1 and 2 (IRS1 and -2). Using two cell culture models of the familial hamartoma syndrome, tuberous sclerosis, we show here that Raptor-mTOR and S6K1 are required for phosphorylation of IRS1 at a subset of serine residues frequently associated with insulin resistance, including S307, S312, S527, S616, and S636 (of human IRS1). Using loss- and gain-of-function S6K1 constructs, we demonstrate a requirement for the catalytic activity of S6K1 in both direct and indirect regulation of IRS1 serine phosphorylation. S6K1 phosphorylates IRS1 in vitro on multiple residues showing strong preference for RXRXXS/T over S/T,P sites. IRS1 is preferentially depleted from the high-speed pellet fraction in TSC1/2-deficient mouse embryo fibroblasts or in HEK293/293T cells overexpressing Rheb. These studies suggest that, through serine phosphorylation, Raptor-mTOR and S6K1 cell autonomously promote the depletion of IRS1 from specific intracellular pools in pathological states of insulin and IGF-I resistance and thus potentially in lesions associated with tuberous sclerosis.

A simple and reliable approach to docking protein–protein complexes from very sparse NOE-derived intermolecular distance restraints

A simple and reliable approach for docking protein-protein complexes from very sparse NOE-derived intermolecular distance restraints (as few as three from a single point) in combination with a novel representation for an attractive potential between mapped interaction surfaces is described. Unambiguous assignments of very sparse intermolecular NOEs are obtained using a reverse labeling strategy in which one the components is fully deuterated with the exception of selective protonation of the delta-methyl groups of isoleucine, while the other component is uniformly (13)C-labeled. This labeling strategy can be readily extended to selective protonation of Ala, Leu, Val or Met. The attractive potential is described by a 'reduced' radius of gyration potential applied specifically to a subset of interfacial residues (those with an accessible surface area > or = 50% in the free proteins) that have been delineated by chemical shift perturbation. Docking is achieved by rigid body minimization on the basis of a target function comprising the sparse NOE distance restraints, a van der Waals repulsion potential and the 'reduced' radius of gyration potential. The method is demonstrated for two protein-protein complexes (EIN-HPr and IIA(Glc)-HPr) from the bacterial phosphotransferase system. In both cases, starting from 100 different random orientations of the X-ray structures of the free proteins, 100% convergence is achieved to a single cluster (with near identical atomic positions) with an overall backbone accuracy of approximately 2 A. The approach described is not limited to NMR, since interfaces can also be mapped by alanine scanning mutagenesis, and sparse intermolecular distance restraints can be derived from double cycle mutagenesis, cross-linking combined with mass spectrometry, or fluorescence energy transfer.

Polypeptide binding proteins: what remains to be discovered?

Many proteins have developed the potential to sequester a client polypeptide chain in its various folding states as a specific intermolecular ligand and, thus, exhibit the properties of a holding chaperone. The resulting complexes can be of a diverse nature in terms of structure and reaction dynamics and are characterized on the basis of various microscopic properties including formation and decay of encounter and Michaelis complexes as well as reactant and product stability. Interpretation of the functional consequences of complex formation in the cell generally tends to be rather complicated, with notable exceptions including complexes formed during the reaction pathways of proteases, protein kinases and protein phosphatases. Peptide bond cis/trans isomerases take up an intermediate position among the poly(oligo)peptide binding proteins because, although the relationship between chain sequestration and catalysis of isomerization can easily be delineated in vitro, it is sometimes difficult to resolve in the cell. Time-resolved studies on interactions involving peptide bond cis/trans isomerases have led to the establishment of generally applicable methods for studying protein-poly(oligo)peptide interactions that are capable of identifying new types of biocatalysis.

Long-range cooperative binding effects in a T cell receptor variable domain

Although cellular processes depend on protein-protein interactions, our understanding of molecular recognition between proteins remains far from comprehensive. Protein-protein interfaces are structural and energetic mosaics in which a subset of interfacial residues, called hot spots, contributes disproportionately to the affinity of the complex. These hot-spot residues can be further clustered into hot regions. It has been proposed that binding energetics between residues within a hot region are cooperative, whereas those between hot regions are strictly additive. If this idea held true for all protein-protein interactions, then energetically significant long-range conformational effects would be unlikely to occur. In the present study, we show cooperative binding energetics between distinct hot regions that are separated by >20 A. Using combinatorial mutagenesis and surface plasmon resonance binding analysis to dissect additivity and cooperativity in a complex formed between a variable domain of a T cell receptor and a bacterial superantigen, we find that combinations of mutations from each of two hot regions exhibited significant cooperative energetics. Their connecting sequence is composed primarily of a single beta-strand of the T cell receptor variable Ig domain, which has been observed to undergo a strand-switching event and does not form an integral part of the stabilizing core of this Ig domain. We propose that these cooperative effects are propagated through a dynamic structural network. Cooperativity between hot regions has significant implications for the prediction and inhibition of protein-protein interactions.

Gaussian-Weighted RMSD Superposition of Proteins: A Structural Comparison for Flexible Proteins and Predicted Protein Structures

Many proteins contain flexible structures such as loops and hinged domains. A simple root mean square deviation (RMSD) alignment of two different conformations of the same protein can be skewed by the difference between the mobile regions. To overcome this problem, we have developed a novel method to overlay two protein conformations by their atomic coordinates using a Gaussian-weighted RMSD (wRMSD) fit. The algorithm is based on the Kabsch least-squares method and determines an optimal transformation between two molecules by calculating the minimal weighted deviation between the two coordinate sets. Unlike other techniques that choose subsets of residues to overlay, all atoms are included in the wRMSD overlay. Atoms that barely move between the two conformations will have a greater weighting than those that have a large displacement. Our superposition tool has produced successful alignments when applied to proteins for which two conformations are known. The transformation calculation is heavily weighted by the coordinates of the static region of the two conformations, highlighting the range of flexibility in the overlaid structures. Lastly, we show how wRMSD fits can be used to evaluate predicted protein structures. Comparing a predicted fold to its experimentally determined target structure is another case of comparing two protein conformations of the same sequence, and the degree of alignment directly reflects the quality of the prediction.

Common Motifs and Topological Effects in the Protein Folding Transition State

Through extensive experiment, simulation, and analysis of protein S6 (1RIS), we find that variations in nucleation and folding pathway between circular permutations are determined principally by the restraints of topology and specific nucleation, and affected by changes in chain entropy. Simulations also relate topological features to experimentally measured stabilities. Despite many sizable changes in phi values and the structure of the transition state ensemble that result from permutation, we observe a common theme: the critical nucleus in each of the mutants share a subset of residues that can be mapped to the critical nucleus residues of the wild-type. Circular permutations create new N and C termini, which are the location of the largest disruption of the folding nucleus, leading to a decrease in both phi values and the role in nucleation. Mutant nuclei are built around the wild-type nucleus but are biased towards different parts of the S6 structure depending on the topological and entropic changes induced by the location of the new N and C termini.

The Differential Engagement of Arrestin Surface Charges by the Various Functional Forms of the Receptor

G-protein-coupled receptor signaling is terminated by arrestin proteins that preferentially bind to the activated phosphorylated form of the receptor. Arrestins also bind active unphosphorylated and inactive phosphorylated receptors. Binding to the non-preferred forms of the receptor is important for visual arrestin translocation in rod photoreceptors and the regulation of receptor signaling and trafficking by non-visual arrestins. Given the importance of arrestin interactions with the various functional forms of the receptor, we performed an extensive analysis of the receptor-binding surface of arrestin using site-directed mutagenesis. The data indicated that a large number of surface charges are important for arrestin interaction with all forms of the receptor. Arrestin elements involved in receptor binding are differentially engaged by the various functional forms of the receptor, each requiring a unique subset of arrestin residues in a specific spatial configuration. We identified several additional phosphate-binding elements in the N-domain and demonstrated for the first time that the active receptor preferentially engages the arrestin C-domain. We also found that the interdomain contact surface is important for arrestin interaction with the non-preferred forms of the receptor and that residues in this region play a role in arrestin transition into its high affinity receptor binding state.

Overlapping, Nonidentical Binding Sites of Different Classes of Nonpeptide Antagonists for the Human Gonadotropin-Releasing Hormone Receptor

Peptide agonists and antagonists of the human gonadotropin-releasing hormone receptor (GnRH-R) are widely used to treat a range of reproductive hormone related diseases. Recently, nonpeptide, orally available GnRH-R antagonists have emerged from several chemical classes. To understand how a relatively large peptide-binding pocket can recognize numerous nonpeptide ligands, we undertook a systematic mapping of GnRH-R residues involved in the binding of three nonpeptide antagonists. A region composed of the extracellular portions of transmembrane helices 6 and 7, extracellular loop 3, and the N-terminal domain significantly contributed to nonpeptide antagonist binding. However, each molecule was affected by a different subset of residues in these regions, indicating that each appears to occupy distinct, partially overlapping subregions within the more extensive peptide-binding pocket. Moreover, the resulting receptor interaction maps provide a basis to begin to reconcile structure-activity relationships between various nonpeptide and peptide series and facilitate the design of improved therapeutic agents.

Query3d: a new method for high-throughput analysis of functional residues in protein structures

BackgroundThe identification of local similarities between two protein structures can provide clues of a common function. Many different methods exist for searching for similar subsets of residues in proteins of known structure. However, the lack of functional and structural information on single residues, together with the low level of integration of this information in comparison methods, is a limitation that prevents these methods from being fully exploited in high-throughput analyses.ResultsHere we describe Query3d, a program that is both a structural DBMS (Database Management System) and a local comparison method. The method conserves a copy of all the residues of the Protein Data Bank annotated with a variety of functional and structural information. New annotations can be easily added from a variety of methods and known databases. The algorithm makes it possible to create complex queries based on the residues' function and then to compare only subsets of the selected residues. Functional information is also essential to speed up the comparison and the analysis of the results.ConclusionWith Query3d, users can easily obtain statistics on how many and which residues share certain properties in all proteins of known structure. At the same time, the method also finds their structural neighbours in the whole PDB. Programs and data can be accessed through the PdbFun web interface.

Similarity networks of protein binding sites

An increasing attention has been dedicated to the characterization of complex networks within the protein world. This work is reporting how we uncovered networked structures that reflected the structural similarities among protein binding sites. First, a 211 binding sites dataset has been compiled by removing the redundant proteins in the Protein Ligand Database (PLD) (http://www-mitchell.ch.cam.ac.uk/pld/). Using a clique detection algorithm we have performed all-against-all binding site comparisons among the 211 available ones. Within the set of nodes representing each binding site an edge was added whenever a pair of binding sites had a similarity higher than a threshold value. The generated similarity networks revealed that many nodes had few links and only few were highly connected, but due to the limited data available it was not possible to definitively prove a scale-free architecture. Within the same dataset, the binding site similarity networks were compared with the networks of sequence and fold similarity networks. In the protein world, indications were found that structure is better conserved than sequence, but on its own, sequence was better conserved than the subset of functional residues forming the binding site. Because a binding site is strongly linked with protein function, the identification of protein binding site similarity networks could accelerate the functional annotation of newly identified genes. In view of this we have discussed several potential applications of binding site similarity networks, such as the construction of novel binding site classification databases, as well as the implications for protein molecular design in general and computational chemogenomics in particular.

The Binding of W7, an Inhibitor of Striated Muscle Contraction, to Cardiac Troponin C

W7 is a well-characterized calmodulin antagonist. It decreases the maximal tension and rate of ATP hydrolysis in cardiac muscle fibers. Cardiac troponin C (cTnC) has been previously implicated as the mechanistically significant target for W7 in the myofilament. Two-dimensional NMR spectra ({1H,15N}- and {1H,13C}-HSQCs) were used to monitor the Ca2+-dependent binding of W7 to cTnC. Titration of cTnC x 3Ca2+ with W7 indicated binding to both domains of the protein. We examined the binding of W7 to the separated domains of cTnC to simplify the spectral analysis. In the titration of the C-terminal domain (cCTnC x 2Ca2+), the spectral peaks originating from a subset of residues changed nonuniformly, and could not be well-described as single-site binding. A global fit of the cCTnC x 2Ca2+ titration data to a two-site, sequential binding model (47 residues simultaneously fit) yielded a dissociation constant (Kd1) of 0.85-0.91 mM for the singly bound state, with the second dissociation constant fit to 3.40-3.65 mM (> or = 4 x Kd1). The titration data for the N-terminal domain (cNTnC x Ca2+) was globally fit to single-site binding model with a Kd of 0.15-0.30 mM (41 residues fit). The data are consistent with W7 binding to each domain's major hydrophobic pocket, coordinating side chains responsible for liganding cTnI. When in muscle fibers, W7 may compete with cTnI for target sites on cTnC.

Dpp-responsive Silencers Are Bound by a Trimeric Mad-Medea Complex

Transcriptional regulation by transforming growth factor-beta signaling is mediated by the Smad family of transcription factors. It is generally accepted that Smads must interact with other transcription factors to bind to their targets. However, recently it has been shown that a complex of the Drosophila Smad proteins, Mad and Medea, binds with high affinity to silencer elements that repress brinker and bag of marbles in response to Dpp signaling. Here we report that these silencers are bound by a heterotrimer containing two Mad subunits and one Medea subunit. We found that the MH1 domains of all three subunits contributed directly to sequence-specific DNA contact, thus accounting for the exceptionally high stability of the Smad-silencer complex. The Medea MH1 domain binds to a canonical Smad box (GTCT), whereas the Mad MH1 domains bind to a GC-rich sequence resembling Mad binding sites previously identified in Dpp-responsive enhancer elements. The consensus for this sequence, GRCGNC, differs from that of the canonical Smad box, but we found that Mad binding nonetheless required the same beta-hairpin amino acids that mediate base-specific contact with GTCT. Binding was also affected by alanine substitutions in Mad and Med at a subset of basic residues within and flanking helix 2, indicating a contribution to binding of the GRCGNC and GTCT sites. The slight alteration of the Dpp silencers caused them to activate transcription in response to Dpp signaling, indicating that the potential for Smad complexes to recognize specific targets need not be limited to repression.

A Conserved Glutamate Residue in Transmembrane Helix 10 Influences Substrate Specificity of Rabbit OCT2 (SLC22A2)

OCT1 and OCT2 are involved in renal secretion of cationic drugs. Although they have similar selectivity for some substrates (e.g. tetraethylammonium (TEA)), they have distinct selectivities for others (e.g. cimetidine). We postulated that "homolog-specific residues," i.e. the 24 residues that are conserved in OCT1 orthologs as one amino acid and in OCT2 as a different one, influence homolog-specific selectivity and examined the influence on substrate binding of three of these conserved residues that are found in the C-terminal half of the rabbit orthologs of OCT1/2. The N353L and R403I substitutions (OCT2 to OCT1) did not significantly change the properties of OCT2. However, the E447Q replacement shifted substrate selectivity toward an OCT1-like phenotype. Substitution of glutamate with cationic amino acids (E447K and E447R) abolished transport activity, and the E447L mutant displayed markedly reduced transport of TEA and cimetidine while retaining transport of 1-methyl-4-phenylpyridinium. In a novel homology model of the three-dimensional structure of OCT2, Glu(447) was found in a putative docking region within a hydrophilic cleft of the protein. In addition, six residues identified in separate studies as exerting significant effects on OCT binding were also found within the putative cleft region. There was a significant correlation (r(2) = 0.82) between the IC(50) values for inhibition of TEA transport by 14 different compounds and their calculated K(D) values for binding to the model of rabbit OCT2. The results suggest that homology modeling offers an opportunity to direct future site-directed studies of OCT/substrate interaction.

NMR Investigation of Main-Chain Dynamics of the H80E Mutant of Bovine Neurophysin-I: Demonstration of Dimerization-Induced Changes at the Hormone-Binding Site

Neurophysins are hormone-binding proteins composed of two partially homologous domains. Ligand-binding (localized to the amino domain) and dimerization (involves both domains) are cooperatively linked by an as yet undefined allosteric mechanism. To help define this mechanism, we investigated the backbone dynamics of the unliganded monomeric state of the H80E mutant of bovine neurophysin-I by (15)N NMR. Model-free analysis of the NMR relaxation parameters indicated significantly greater flexibility in the carboxyl domain than in the amino domain, particularly at their dimerization interface segments. Amino domain residues critical to hormone binding were highly structured, constraining potential allosteric mechanisms. Model-free analysis additionally demonstrated chemical exchange effects, manifest as R(ex) terms, in 16 residues, 14 of which are located in the amino domain at, or immediately adjacent to, either the dimerization interface or the hormone-binding site. The chemical exchange process was further characterized using relaxation-compensated CPMG measurements, the results allowing assignment of the process to monomer-dimer exchange and calculation of the exchange kinetics, which were slow on the NMR time scale. An apparently different concentration-dependent process, distinguished from normal dimerization by its fast exchange behavior and pH-independence, also principally involved a subset of residues at and immediately adjacent to either the hormone-binding site or the amino domain dimerization interface. The data represent the first direct demonstration of an effect of dimerization in the unliganded state on neurophysin's hormone-binding site, the effect particularly involving residues that interact with hormone residue 2, and specifically identify Ser25 and Ile26 as likely intermediaries between the sites of dimerization and of hormone binding. Consistent with recent views of the role of anchor residues in protein interactions, we propose that dimerization proceeds by a fast pH-independent association of the well-structured amino domain interface that is rapidly communicated to the binding site for hormone residue 2, followed by a rate-determining pH-dependent interaction of the less structured carboxyl domain interface.

Functional annotation by identification of local surface similarities: a novel tool for structural genomics

BackgroundProtein function is often dependent on subsets of solvent-exposed residues that may exist in a similar three-dimensional configuration in non homologous proteins thus having different order and/or spacing in the sequence. Hence, functional annotation by means of sequence or fold similarity is not adequate for such cases.ResultsWe describe a method for the function-related annotation of protein structures by means of the detection of local structural similarity with a library of annotated functional sites. An automatic procedure was used to annotate the function of local surface regions. Next, we employed a sequence-independent algorithm to compare exhaustively these functional patches with a larger collection of protein surface cavities. After tuning and validating the algorithm on a dataset of well annotated structures, we applied it to a list of protein structures that are classified as being of unknown function in the Protein Data Bank. By this strategy, we were able to provide functional clues to proteins that do not show any significant sequence or global structural similarity with proteins in the current databases.ConclusionThis method is able to spot structural similarities associated to function-related similarities, independently on sequence or fold resemblance, therefore is a valuable tool for the functional analysis of uncharacterized proteins. Results are available at