BACKGROUND
sequencing of complete genomes has become commonplace in the last decade but traditional methods of protein structure determination cannot keep pace.
structural genomics  has emerged with the goal of providing a structure for every protein encoded by a genome using a high-throughput combination of experimental structure determination and homology modelling approaches.
the protein structure initiative  is a usa federal government, university, and industry effort that has now completed two phases  <cit> .
the first phase of the psi , from  <dig> to  <dig>  was dedicated to demonstrating the feasibility of high-throughput structure determination, solving unique protein structures, and developing the methodology and technology for a subsequent production phase.
the second phase, psi- <dig>  focused on implementing the high-throughput structure determination methods developed in psi- <dig>  as well as homology modelling and addressing bottlenecks like modelling membrane proteins.
there are four large-scale psi- <dig> centres and a number of smaller specialist centres.
the four large-scale centres are the joint center for structural genomics , the midwest center for structural genomics , the northeast structural genomics  consortium, and the new york sgx research center  for structural genomics.
when psi- <dig> began at the end of the 20th century there was much optimism that sg could provide structures to cover the whole of protein sequence space and together with improvements in homology modelling technology, it was only a matter of time before a reasonably accurate structure could be predicted for any protein in nature as soon as its sequence was known.
it has since become apparent that protein sequence space is much larger than was imagined at the time.
advancements in homology modelling have also not been as great as was hoped.
as a consequence sg has failed to deliver the hoped for level of coverage of sequence space and has been left with a collection of structures that much of the time were not targeted on the basis of their biological function.
this has led to a criticism of the psi that many of its structures are of proteins of unknown function and many do not have corresponding publications and therefore give little biological insight.
indeed, the third phase of the psi is called psi:biology, which is intended to reflect a new emphasis on the biological relevance of the work.
as psi- <dig> draws to a close, we attempt to partially address this problem by exploiting a wide range of bioinformatics tools to provide functional annotations for as many as possible of the protein structures experimentally determined so far by the mcsg.
an early analysis of structures solved worldwide by sg consortia assessed the new coverage of sequence and fold space using the cath  <cit>  and scop  <cit>  domain structure classifications  <cit> .
while sg was judged to be succeeding in structurally characterising new superfamilies, an early observation was that many of the sequences selected as likely to represent new folds were actually found to have existing folds.
in  <dig> xie and bourne adopted a different approach to measuring the impact of sg  <cit> .
they assessed functional coverage of the human genome by existing structures, structural genomics targets, and homology models.
using functional categories in the enzyme commission   <cit>  and gene ontology   <cit>  classifications they showed that, at the time, sg structures provided at least one domain that covered about a third of all the functional classes in the genome, and whole structure coverage for about a quarter of the genome.
even if all the registered sg targets were solved, which was a substantial number of targets even in  <dig>  then only about two thirds of the functional classes would be covered and there would be whole structure coverage for less than half of the genome.
they therefore argued in favour of target selection aimed at functional coverage, especially with a view to understanding human genetic disease.
also in  <dig> chandonia and brenner suggested the "pfam5000" strategy for target selection  <cit>  which involved selecting the  <dig> most important pfam families  <cit> .
this comprehensive collection of protein families is used extensively throughout the biological sciences, often for functional annotation of complete genomes, and it was felt that this strategy would be medically and biologically relevant, of good value, and tractable.
pfam <dig> was felt to be better value than the complete solution of several small to moderately sized bacterial proteomes or partial coverage of the human proteome since these would have limited impact on structural knowledge of other proteomes.
the jcsg, on the other hand, focused on selecting targets from thermotoga maritima leading to the generation of a three-dimensional reconstruction of the central metabolic network of this bacterium  <cit> .
the jcsg also championed the selection of sequences that were over-represented in the human gut metagenome  <cit>  and this strategy became part of the overall psi effort.
chandonia and brenner went on to analyze the novelty, cost, and impact of sg structures compared to non-sg structural biology  structures  <cit> .
determination of the first structure in a pfam family was established as a measure of success.
it should be noted, however, that a pfam family often contains proteins with a range of related functions and finer grained target selection would be necessary to truly achieve complete functional coverage.
in  <dig> about half of all first structures for pfam families were from sg rather than sb.
marsden et al.
<cit>  analyzed  <dig> complete genomes in the gene3d resource  <cit>  to provide new insights into protein family space.
the number of protein families was found to be continually expanding with time but a significant proportion of the proteomes could be assigned to relatively few large, well-characterized domain families while the vast majority of domain families were relatively rare and often species specific.
it was suggested that sg could provide structures for fewer than a thousand pfam families to achieve reasonable structural coverage of genomes.
within these large families it was suggested that multiple structures could be determined to reveal more about the evolution of the family and enable greater understanding of how function evolves.
similarly, the selection of targets from within large and diverse cath superfamilies was proposed to maximize modelling leverage  <cit> .
since many of the targets previously selected as likely to have a novel fold proved to have a known fold this could be seen as a more rational use of resources.
there was a coordinated strategy of target selection by the four large-scale centres during psi- <dig>  described by dessailly et al.
<cit> , involving four main approaches: targeting representatives from large, structurally uncharacterized protein domain superfamilies; targeting structurally uncharacterized subfamilies in very large and diverse superfamilies with incomplete structural coverage; community nominated targets; and biomedically important targets.
the first two approaches were primarily aimed at increasing the structural coverage of protein sequence space and many targets were selected regardless of their function.
metrics have been calculated in this and several other publications  <cit>  to show that this strategy has had the result that psi's per structure contribution to novel structural leverage has been much higher than that for sb, as would be expected since sb is not aiming to increase sequence coverage but rather provide biological insights.
the structural biology knowledgebase  <cit>  offers an easy way of keeping abreast of developments by the psi and by sg and sb in general.
it is a continually updated portal to research data and other resources from the psi.
metrics are regularly updated giving a summary of the total number of structures, the number of distinct and novel structures determined by the psi, and the modelling leverage that psi structures provide.
models are made available through the protein model portal  <cit> .
the jcsg has made a notable effort towards the annotation of psi structures with the open protein structure annotation network , a wiki-based portal for the dissemination of information for the broad biomedical community  <cit> .
the increased appearance of structures of proteins of unknown function in the pdb due to sg has stimulated the development of computational methods of function prediction from structure.
profunc, developed for the mcsg, combines a number of sequence-based and structure-based methods to gain clues about the function of a protein  <cit> .
the mcsg psi- <dig> structures were used to test and develop profunc  <cit> .
when a sequence match is weak and/or multiple functions are suggested, the structure-based methods within profunc can help select from the options that are presented and increase confidence in a particular putative function.
another recent and notable method that has been used for function prediction from psi structures is flora  <cit> .
this uses structural motifs associated with different functional sub-groups within functionally diverse cath domain superfamilies.
the mcsg selected targets from a broad range of pathogenic bacteria according to selection criteria that have evolved throughout the first two stages of the psi.
here we attempt to extend the functional annotations associated with mcsg structures by employing a diverse set of bioinformatics tools and resources.
we also analyse the structural novelty of protein structures solved by the mcsg.
we use sequence-based methods to assess the proportion of mcsg proteins that have a known function, a putative function, a possible function, and no known function.
this is followed up with profunc analysis to support some of the possible functional annotations and in some cases to suggest possible functions for the unknowns.
methods
a list of all protein structures released by the pdb  <cit> , their release date, and the source sg centre where applicable was downloaded from the pdb on the 3rd march  <dig>  in all subsequent analyses only data that were available on this date are used.
at the time  <dig>  mcsg protein structures for  <dig>  targets had been released by the pdb.
annotations
the source of most of the annotations used in this work is gene3d  <dig> .
<dig>  <cit> .
this database contains nearly  <dig> million protein sequences including all uniprot sequences  <cit>  and most complete genomes.
annotations are imported from multiple sources including ec numbers  <cit>  from uniprot, go terms  <cit> , kegg genes  <cit>  and protein family assignments from pfam  <cit> , tigrfams  <cit> , smart  <cit> , and panther  <cit> .
kegg currently has  <dig> reference pathway maps.
these are represented in a general way to be applicable to all organisms and thus are useful in the analysis of sg structural coverage.
each node in a reference pathway is represented by one or more manually curated kegg orthologs.
where there are multiple orthologs these may represent different subunits or different versions of the same enzyme.
also a node may belong to more than one pathway when different pathways interact.
kegg genes and orthologs are mapped to genomes in gene3d and pdb chains are mapped to uniprot entries using http://www.ebi.ac.uk/thornton-srv/databases/pdbsum/data/pdb_chain_sp_ec the data for which comes from the pdbe database  <cit> .
in addition to collecting these annotations from gene3d for mcsg structures, annotations are considered for inheritance by sequence-sequence comparison, sequence-profile comparison, and profile-profile comparison, such that increasingly remote sequence similarities are considered.
not all assignments to a protein family lead to a functional annotation.
some assignments are filtered out as being uninformative e.g.
assignments to a pfam family described as a duf .
the sequence of each mcsg structure is compared to all sequences in gene3d using blast  and five iterations of psi-blast , both with an e-value cut-off of  <dig>   <cit> .
ec numbers and molecular function go terms associated with matching sequences are considered for inheritance.
both of these types of annotation are hierarchic.
ec numbers all have four levels, the first level refers to the enzyme class , the second level refers to the type of bond or group that is acted on , and the next two levels give progressively more specific details of the catalyzed reaction and its substrates.
go molecular function terms are organized as a directed acyclic graph where there may be many nested levels with each child having an 'is a' relationship to its parent.
annotations associated with a sequence or sequence match are divided into deep annotations comprising level  <dig> or  <dig> ec numbers and go terms at level  <dig> or above, and the remainder which are generally less specific annotations.
hidden markov models  are built using the sequence of each mcsg structure as a seed for sam-t  <cit>  and then each of these hmms is compared to all pfam  <dig>  hmms  using prc  <cit> .
database annotations, inherited annotations, and annotations from the headers of the pdb files are presented together in a table at http://www.biochem.ucl.ac.uk/cgi-bin/dlee/mcsg_annotations used by the authors to manually assign an annotation status to each mcsg structure.
annotations are organized from left to right to reflect decreasing levels of both depth and reliability.
four levels of status are assigned; known, putative, possible, and unknown.
go provides evidence codes for functional annotations where the most reliable codes are the experimental evidence and author statement codes.
62% of mcsg solved targets have a molecular function go term associated with them but only 2% have such a high quality functional annotation, the remainder have annotations that are electronically inferred.
there are currently very few other annotation resources where evidence codes are available and sequence coverage is low.
we consider possession of one of the most reliable go evidence codes as too strict a criterion to define the 'known' status.
rather, possession of a deep database annotation and preferably agreement with at least one set of protein families is generally taken as evidence that the function is regarded as being known.
note that a go term level of  <dig> or higher is only a rough guide to the depth of an annotation and the association of a term such as go: <dig> for iron ion binding  is not considered to be sufficient in itself for a 'known' status.
this is one reason why manual assessment of status is so important.
the 'putative' status is generally assigned where there is no deep database annotation but there is a protein family assignment and the family description includes specific molecular functions.
the 'possible' status is generally assigned where only inherited annotations are available but without the added confidence provided by assignment to a family in a curated resource.
there may be a choice between a variety of deep annotations or perhaps only a very general, non-specific annotation is available presenting a broad range of possible deep annotations.
the 'unknown' status is applied where sequence methods are unable to provide any clue to the function of the protein.
profunc analysis  is linked to the online table for structures assigned a 'possible' or 'unknown' status.
go slims are cut-down versions of the go ontologies containing a subset of the terms in the whole go.
the go provides a generic go slim that gives a broad overview of all function categories without the detail of the specific fine grained terms.
this is useful for visualization of the broad functional coverage of mcsg structures and for comparison to the functional coverage of the pdb as a whole.
a sequence unique subset that approximately represents all pdb protein entries of known function is generated by selecting representatives annotated in gene3d with a molecular function go term at level six or higher.
go terms associated with these representatives are then mapped to the generic go slim using the go slimmer tool at http://amigo.geneontology.org/cgi-bin/amigo/slimmer.
similarly the uniprot entries for all mcsg solved targets with known function manually assigned as above are mapped to the generic go slim.
profunc analysis
the profunc web server at http://www.ebi.ac.uk/profunc <cit>  employs a number of complementary function-prediction methods, focusing, in particular, on methods based on 3d structural information.
this accepts a pdb file and runs a number of sequence- and structure-based analyses on it, listing any hits it finds to existing sequences and structures in descending order of significance.
these can help researchers identify any strong similarities which may be indicative of the protein's function.
of profunc's structure-based methods, the two that have been shown most accurately to suggest function  <cit>  are the ssm  fold-matching algorithm  <cit>  and the "reverse template" search  <cit> .
both these methods have a similar success rate and are able to find distant homologues when simple sequence-based searches fail.
the majority of their hits overlap, but occasionally one method picks up a match that the other fails to identify.
the advantage of the reverse template method over the fold-match is that its matches are local and the top-scoring hits tend to pick up sites of functional relevance.
to demonstrate the added value of 3d structure for proteins of uncertain or unknown function, we run profunc on the full data set of  <dig>  mcsg structures.
profunc results for targets assigned a 'possible' or 'unknown' status following sequence analysis are linked to the online table.
profunc results for targets of 'unknown' status are manually assessed to determine the improvement in functional coverage afforded by profunc for mcsg solved targets.
to further illustrate the value of profunc analysis we then focus on mcsg structures where no significant match to a known 3d structure in the pdb can be found from a search based on sequence alone, but is detectable using the structure-based reverse template method.
we investigate which parts of the matched structures the method identifies as the most significantly similar.
reverse templates
the reverse template method takes the query protein structure and breaks it up into a large set  of 3-residue templates.
each template consists of three neighbouring residues chosen such that the closest atom-atom distance between any two of the residues is no larger than  <dig> Å. templates containing more than one hydrophobic residue are rejected in order to bias the templates towards surface residues.
a scan of each template against a representative set of the structures in the pdb is then made using the jess algorithm  <cit>  to identify structures containing similar constellations of the three residues.
the representative set is a non-redundant list of protein chains in the pdb downloaded weekly from the pisces ftp server  <cit> .
hits are scored and assigned e-values as described in laskowski et al.
<cit> , the scoring being based on the local similarity of the equivalent regions when the structures are superposed on the matched residues.
the aim of the method is to pick up local structural similarity between two proteins and, in particular, to identify regions that have been conserved by evolution and hence are more likely to be functionally important.
thus, even though the sequences of two homologous proteins might have diverged considerably over evolutionary time, it is possible that, to maintain their function, the proteins' active sites have undergone less change.
the detection of such local conservation can provide strong support for the proteins having retained the same, or a similar, function.
<dig> residues and a  <dig> Å cut-off were found through empirical testing to give a good balance between the duration and success of a search .
identification of functionally important regions
to demonstrate this tendency to pick up functionally important regions, we select a subset of the mcsg data where a simple fasta  <cit>  search against the sequences in the pdb fails to return a significant result.
we use an e-value cut-off of  <dig> , which gives us  <dig> structures whose best match is at or below this threshold.
we then consider the top four hits returned by the reverse template search for each of these  <dig> structures to see if the template residues have any association with the protein's function.
we use the functional annotations given in pdbsum  <cit> , namely: residues belonging to a prosite  <cit>  pattern; catalytic residues, as determined from the catalytic site atlas   <cit> ; residues defined in the site records of the original pdb file; or residues in contact with a ligand, metal or dna/rna in the structure.
however, as the search data set used by the reverse template method is a representative one, it is possible that the pdb entry matched by the search has no functional annotation, yet a closely related pdb entry, excluded from the data set, does have some annotation.
there were nearly  <dig>  protein chains in the representative data set, compared with over  <dig>  structures in the pdb.
thus, to pick up such 'lost' annotations, we use the sequence annotated by structure  server  <cit>  to find additional functional annotations for each of the matched structures.
sas identifies closely related sequences from the pdb using fasta.
functional annotations are then imported from the resultant alignment where the sequence identity is at least 30%, the alignment overlap at least  <dig> residues, and the fasta e-value <  <dig> .
functional novelty
the release dates of all pdb entries that are associated with a go term  <cit> , are assigned to a pfam family  <cit> , or map to kegg orthologs  <cit>  are compared to find the earliest example released by the pdb for each annotation.
each structure that represents a first example is then categorized as being solved by a psi centre or as 'non-psi' in order to rate the performance of psi compared to the combined effort of all other laboratories, and to compare the performance of individual psi centres.
the dates that the first structures were solved are also divided into the two psi periods psi- <dig>  and psi- <dig>  a measure of novelty per pdb release is calculated and for comparison the novelty per non-redundant structure at 95% sequence identity following clustering using cd-hit is also given <cit> .
structural novelty
the release dates of all pdb entries containing a domain assigned to a superfamily in cath v <dig>  are compared to find the earliest entry in that superfamily  and fold .
structural novelty within a superfamily is determined by calculating a normalized rmsd .
a domain with a normrmsd of 5Å or more from all cath domains in structures previously released by the pdb is defined as structurally novel and belongs to a new structural sub-group .
both mcsg domains assigned to cath and unassigned mcsg pdb chains are compared to a representative library of cath domains  using cathedral  <cit>  as a fast filter for significant structural similarity.
the best matches by cathedral are then structurally aligned and superposed using the more computationally expensive and accurate method ssap  <cit>  and the rmsd is calculated.
the normrmsd is then calculated as follows:
where rmsd is the root mean square deviation of the superposition, max  is the length in residues of the longest domain in the superposition, and nmat is the number of aligned residue pairs   <cit> .
novel modelling leverage
novel modelling leverage is calculated according to the method of nair et al.
<cit> .
the sequences of all pdb protein chains, non-redundant at 100% sequence identity, are compared to uniref <dig> sequences downloaded on the 3rd of march  <dig> using psi-blast with three iterations and an e-value cut-off of 1e- <dig>  novel modelling leverage in residues is determined for each pdb entry  on the date of its release by the pdb.
the novel modelling leverage of the four large-scale psi centres is compared to each other and to the leverage of the combined non-psi laboratories.
human non-synonymous single nucleotide polymorphisms
the mcsg focused on solving structures of proteins from pathogenic bacteria and has only solved structures for two human proteins.
however, a number of mcsg structures may be used as templates for comparative modelling of human proteins.
all such structures identified in the modelling leverage analysis described above are used to analyze human non-synonymous nucleotide polymorphisms .
mcsg structures are included regardless of whether a closer template is available from another laboratory.
for each human modelling target identified above, variants are retrieved from uniprot, ensembl and omim and filtered for uniqueness.
uniprot variants and ensembl gene ids are listed for each protein in the uniprot text file available from uniprot.
the ensembl api is used to query their mysql database and retrieve all ensembl variants.
omim variants are linked to uniprot ids through the omim missense server at http://www.bioinf.org.uk/omim/.
each variant is inspected to see if it falls within a region having a template match and is thus modellable.
RESULTS
annotations
the functional annotation status of the mcsg solved targets is as follows: 31% known; 48% putative; 14% possible; and 7% unknown .
profunc analysis of the unknowns is summarised in figure  <dig> b. the majority of structures are therefore of direct relevance to biological investigations of protein function as well as fulfilling some of the goals of sg in expanding structural coverage of fold and function space as shown below.
the apparent discrepancy between the results obtained by looking at ec numbers in figure  <dig> and go terms in figure  <dig> needs explaining.
ec numbers are very high quality annotations but coverage of the pdb is less than that of go and some catalytic activity is not represented by ec numbers.
within the go catalytic activity category an over-representation of transferase activity and an under-representation of hydrolase activity is seen which does broadly agree with the ec analysis.
within the go binding category mcsg structures also show particularly significant under-representation of protein binding and calcium binding while nucleotide binding is actually a little over-represented.
profunc analysis
profunc's template searches are rated as certain, probable, possible, or long shots to indicate the likelihood that a prediction is correct.
profunc analysis is linked to all structures assigned a "possible" or "unknown" status in the online table at http://www.biochem.ucl.ac.uk/cgi-bin/dlee/mcsg_annotations profunc suggests functions for all of the  <dig> mcsg targets that are of unknown function following sequence analysis.
the overall likelihood of the best scoring profunc prediction being correct is shown in figure  <dig> b. one target has a prediction rated as being certain while another 17/ <dig>  targets have predictions rated as being probable.
these predictions are all made by the reverse templates method while two of the probable predictions are also made by ligand-binding templates.
profunc makes a valuable contribution to functional annotation coverage of mcsg structures.
profunc example
an example of a reverse template match is pdb entry 2aua which is the structure of bc <dig>  a protein from bacillus cereus that is classified as 'unknown' by the sequence analysis described above.
the reverse template search finds one 'certain' match, one 'probable' match and seven 'possible' matches.
figure  <dig> shows the template residues and surroundings for three of these matches, all being to the catalytic domains of bacterial toxins.
the first is diphtheria toxin , pdb code 1f0l  <cit> , and the other two are 'possible' matches: one to exotoxin a from pseudomonas aeruginosa, pdb code  <dig> × k <dig>  <cit>  and the second to cholix toxin from vibrio cholera, pdb entry 3ess  <cit> .
each of the templates includes one or more of the known catalytic residues, and the local region returned by the template match corresponds to the matched protein's known substrate binding site.
the catalytic domains correspond to ec  <dig> .
<dig> , which contains the nad+-diphthamide adp-ribosyltransferases.
all have the same alpha-beta fold, cath  <dig> .
<dig> , corresponding to domain  <dig> of diphtheria toxin.
these matches strongly suggest that the bc <dig> protein has a very similar function to these toxins.
bacillus cereus produces a number of known enterotoxins that are responsible for food poisoning  <cit> .
possibly, this represents another.
a second profunc example is presented  to show how profunc analysis can add weight and specific detail to a possible annotation suggested by sequence analysis.
functional novelty
so far  <dig>  out of  <dig>  pfams,  <dig>  out of  <dig>  go molecular function terms and  <dig>  out of  <dig>  terms at level  <dig> or higher have at least one representative with a solved structure.
note, however, that because it is generally only small pfam families that lack a structural representative, 88% of unique pfam assignments in complete genomes in gene3d belong to a pfam family with a structural representative.
table  <dig> gives an indication of the value of pdb releases in terms of functional novelty.
from an sg point of view the number of novel structures is given as a percentage of the total number of structures.
sb laboratories often solve multiple structures of the same protein with different ligands bound, at a different ph, or with single amino acid substitutions etc.
these investigations have different goals from those of sg but here the intention is to assess the performance of psi centres in terms of sg goals.
for an sb point of view the number of novel structures is given for comparison as a percentage of the number of non-redundant structures at 95% sequence identity.
although psi centres did not specifically target structurally uncharacterized go terms their overall performance is at least as good as that of all non-psi laboratories combined in terms of sg goals.
numbers of first structures and their percentage of total structures and non-redundant  structures at 95% sequence identity for molecular function go terms and pfam families during psi, psi- <dig>  and psi- <dig>
pfam families, on the other hand, were targeted by psi centres and the proportion of first structures solved out of all structures solved is significantly higher for psi than for non-psi laboratories as would be expected.
as would also be expected, the proportion of first structures for pfams are considerably higher than those for go terms since go terms are a much more general functional category compared to the evolutionarily related members of pfams.
the further back in time that we go the fewer is the number of go terms and pfams that were structurally characterized and thus there was a higher probability of solving the first structure for a term or family.
as structural coverage becomes more complete the probability of solving the first structure for a term or family diminishes.
this is one likely reason that the rate of solving first structures during psi- <dig> is higher than that during psi- <dig> for both psi and non-psi laboratories.
another possible explanation is that the easier functions have been solved and the remaining functions are more difficult to obtain a structure for.
perhaps we are seeing a combination of both factors.
by virtue of solving more structures than any other individual psi centre the mcsg also solved the highest number of first structures for both go terms and pfams.
in both cases the proportion of first structures is close to the average for all psi centres although they excelled at solving the first structures for pfams during psi- <dig>
39/ <dig> of the mcsg first structures for go molecular function terms correspond to a level  <dig> ec number and  <dig> of the remaining  <dig> also have catalytic activity e.g.
go: <dig> amino acid-transporting atpase activity, again showing a bias towards enzymes.
the ec and go annotation systems have different coverage from each other.
analysis of ec numbers in gene3d shows that the mcsg solved the first structure for  <dig> rather than  <dig> level  <dig> ec numbers.
<dig> of these ec numbers were not revealed by the go analysis so in total the mcsg solved the first structure for  <dig> level  <dig> ec numbers.
the novelty of these structures may also be considered in terms of providing structural coverage for metabolic pathways.
pathways
there are currently  <dig> kegg reference pathways.
these pathways have nodes, most but not all of which correspond to ec numbers.
kegg orthologs representing evolutionarily related groups of proteins map to many but not all of these nodes and many nodes have multiple orthologs.
currently,  <dig>  orthologs are assigned to the reference pathways and more are likely to be added in the future.
<dig>  out of a total of  <dig>  currently possible level  <dig> ec numbers map through kegg orthologs to kegg reference pathway nodes.
<dig> / <dig>   have at least one structural representative compared to  <dig> / <dig>   for the ec classification as a whole.
9/ <dig>  pathways have complete structural coverage of their orthologs but these are all small pathways, the largest having only  <dig> orthologs belonging to  <dig> nodes, ko <dig> synthesis and degradation of ketone bodies.
this pathway is also the largest pathway composed entirely of nodes corresponding to level  <dig> ec numbers with complete structural coverage.
the  <dig> nodes with ec numbers that the mcsg solved the first structure for belong to  <dig> kegg reference pathways but this number is a coincidence, with some of the nodes appearing in multiple pathways while multiple nodes also appear in single pathways.
the largest number of nodes within a single kegg reference pathway to gain a first structure from the mcsg is  <dig> for pathway ko <dig>  arginine and proline metabolism, but this is a large pathway with more than  <dig> nodes.
psi centres did not select targets with the aim of achieving complete coverage of metabolic pathways but this is a potential goal for the future.
overall the mcsg solved the first structure for  <dig> orthologs in  <dig> pathways.
<dig> of these orthologs belong to multiple pathways and thus join pathways together.
membership of multiple pathways is, however, not rare with about a half of all kegg orthologs  belonging to more than one reference pathway.
structural novelty
of the  <dig>  mcsg pdb entries analyzed in this work only  <dig> are so far classified in cath v <dig>  and so the analysis of structural novelty is currently incomplete.
also, since a special effort was made by the curators of cath to update the classification of mcsg structures as a priority over other structures, mcsg structural novelty cannot be compared to that of other psi centres or sb in general.
mcsg structural novelty is shown in table  <dig>
number of novel cath structural domains solved by the mcsg.
three levels of novelty are defined: a domain with a novel fold is the first example released by the pdb of a fold in cath ; a domain belonging to a novel superfamily is the first example released by the pdb of a superfamily in cath ; a domain belonging to a novel structural sub-group  has a normalized rmsd of at least 5Å from all previously existing domains in cath.
by definition all new folds are also new superfamilies and all new superfamilies are also new structural subgroups .
however, in this analysis not all new ssgs are classified in cath or necessarily belong to existing superfamilies or folds.
of the  <dig> domains that belong to new ssgs but are not classified in cath as new folds or superfamilies,  <dig> are found to be classified into pre-existing cath superfamilies while  <dig> are not classified in cath v <dig> .
<dig> of the  <dig> not classified achieve a score > <dig> from the structural comparison program ssap following the cathedral scan suggesting that they probably do belong to pre-existing superfamilies.
the remaining  <dig> domains could potentially represent new superfamilies or folds.
novel modelling leverage
the average novel modelling leverage in residues per structure is considerably higher during psi- <dig> compared to psi- <dig> for both psi and non-psi .
this is despite total modelling leverage per structure actually being higher for both psi and non-psi during psi- <dig> compared to during psi- <dig> .
the reason is that as the cumulative total and diversity of solved structures increases and targets with good modelling leverage are preferentially selected in the early stages of psi, any new structures that are solved have less chance of providing novel modelling leverage for a fixed set of protein sequences.
during both psi- <dig> and psi- <dig>  psi substantially outperformed non-psi in terms of novel modelling leverage per structure as would be expected since non-psi often target multiple forms of a given target protein to better characterise the biology.
with the increased output of psi during psi- <dig>  the contribution of psi to total worldwide novel modelling leverage more than doubled from  <dig> % during psi- <dig> to  <dig> % during psi- <dig>  during both psi- <dig> and psi- <dig>  the mcsg substantially outperformed the other three large-scale centres in terms of both total novel modelling leverage and average novel modelling leverage per structure .
novel modelling leverage for residues in uniref <dig> during psi, psi- <dig>  and psi- <dig>
human non-synonymous single nucleotide polymorphisms
in the modelling leverage analysis  <dig> mcsg structures are identified as templates for  <dig> human proteins in uniprot.
in this analysis of modellable nssnps all modelling leverage of mcsg structures is considered, not just novel leverage, and the mcsg structure is not necessarily the closest available template.
within these human protein sequences  <dig>  unique ensembl nssnps,  <dig>  unique uniprot nssnps, and  <dig> unique omim nssnps are identified.
<dig>  of these nssnps are within a modellable region when using an mcsg structure as a template.
a good example of an mcsg structure helping to explain human disease is illustrated by pdb entry 2hma.
human mitochondrial trna-specific 2-thiouridylase  <dig>  matches 2hma with 42% sequence identity in a blast alignment.
the enzyme has been implicated in aggravating mitochondrial 12s ribosomal rna aminoglycoside-induced and non-syndromic deafness.
it catalyzes the 2-thiolation of uridine at the wobble position  of mitochondrial trna, trna and trna.
a few natural variants have been identified and one  has been linked to a decrease in enzymatic activity.
to attempt to explain this, a homology model is built using 2hma as a template which includes the substrate, s-adenosyl methionine .
when inspecting the model, it is clear that ala  <dig> is part of the sam binding pocket of the enzyme .
in 2hma the adenosine group fits neatly into this pocket where ala  <dig> is located.
in the ala10ser mutation, an additional hydrogen bond between the protein and sam is probably introduced.
this increases the binding affinity between sam and the protein and thus slows down the release of sam once the sulphur has been transferred.
this agrees with observations made by guan et al.
<cit>  that the ala10ser mutation reduces activity.
CONCLUSIONS
the mcsg has performed well during the first two phases of the psi in terms of the goals that were established during this period.
whilst there has been concern that the psi produced too many structures of unknown function our analysis reveals that by using a range of bioinformatics tools and resources we are able to provide functional annotations for more than 90% of the structures solved by the mcsg.
structures of unknown function have helped stimulate the development of methods such as profunc that can predict function from structure.
profunc analysis complements sequence analysis by both adding weight and specific detail to predicted function and by suggesting function where sequence methods have failed to do so.
authors' contributions
rl carried out the profunc analyses, tdeb carried out the nssnp analyses, and dl carried out all other analyses and drafted the manuscript.
all authors participated in the design of the study, helped draft the manuscript, and read and approved the final manuscript.
supplementary material
additional file 1
additional profunc example.
results and discussion, figure, and reference showing how profunc can be used to refine, add detail to, and support a protein function prediction.
click here for file
acknowledgements
this work was funded by the nih as part of the protein structure initiative.

 sequencing of complete genomes has become commonplace in the last decade but traditional methods of protein structure determination cannot keep pace. structural genomics has emerged with the goal of providing a structure for every protein encoded by a genome using a high-throughput combination of experimental structure determination and homology modelling approaches. the protein structure initiative is a usa federal government, university, and industry effort that has now completed two phases. the first phase of the psi, from to was dedicated to demonstrating the feasibility of high-throughput structure determination, solving unique protein structures, and developing the methodology and technology for a subsequent production phase. the second phase, psi- focused on implementing the high-throughput structure determination methods developed in psi- as well as homology modelling and addressing bottlenecks like modelling membrane proteins. there are four large-scale psi- centres and a number of smaller specialist centres. the four large-scale centres are the joint center for structural genomics, the midwest center for structural genomics, the northeast structural genomics consortium, and the new york sgx research center for structural genomics. when psi- began at the end of the 20th century there was much optimism that sg could provide structures to cover the whole of protein sequence space and together with improvements in homology modelling technology, it was only a matter of time before a reasonably accurate structure could be predicted for any protein in nature as soon as its sequence was known. it has since become apparent that protein sequence space is much larger than was imagined at the time. advancements in homology modelling have also not been as great as was hoped. as a consequence sg has failed to deliver the hoped for level of coverage of sequence space and has been left with a collection of structures that much of the time were not targeted on the basis of their biological function. this has led to a criticism of the psi that many of its structures are of proteins of unknown function and many do not have corresponding publications and therefore give little biological insight. indeed, the third phase of the psi is called psi:biology, which is intended to reflect a new emphasis on the biological relevance of the work. as psi- draws to a close, we attempt to partially address this problem by exploiting a wide range of bioinformatics tools to provide functional annotations for as many as possible of the protein structures experimentally determined so far by the mcsg. an early analysis of structures solved worldwide by sg consortia assessed the new coverage of sequence and fold space using the cath and scop domain structure classifications. while sg was judged to be succeeding in structurally characterising new superfamilies, an early observation was that many of the sequences selected as likely to represent new folds were actually found to have existing folds. in xie and bourne adopted a different approach to measuring the impact of sg. they assessed functional coverage of the human genome by existing structures, structural genomics targets, and homology models. using functional categories in the enzyme commission and gene ontology classifications they showed that, at the time, sg structures provided at least one domain that covered about a third of all the functional classes in the genome, and whole structure coverage for about a quarter of the genome. even if all the registered sg targets were solved, which was a substantial number of targets even in then only about two thirds of the functional classes would be covered and there would be whole structure coverage for less than half of the genome. they therefore argued in favour of target selection aimed at functional coverage, especially with a view to understanding human genetic disease. also in chandonia and brenner suggested the "pfam5000" strategy for target selection which involved selecting the most important pfam families. this comprehensive collection of protein families is used extensively throughout the biological sciences, often for functional annotation of complete genomes, and it was felt that this strategy would be medically and biologically relevant, of good value, and tractable. pfam was felt to be better value than the complete solution of several small to moderately sized bacterial proteomes or partial coverage of the human proteome since these would have limited impact on structural knowledge of other proteomes. the jcsg, on the other hand, focused on selecting targets from thermotoga maritima leading to the generation of a three-dimensional reconstruction of the central metabolic network of this bacterium. the jcsg also championed the selection of sequences that were over-represented in the human gut metagenome and this strategy became part of the overall psi effort. chandonia and brenner went on to analyze the novelty, cost, and impact of sg structures compared to non-sg structural biology structures. determination of the first structure in a pfam family was established as a measure of success. it should be noted, however, that a pfam family often contains proteins with a range of related functions and finer grained target selection would be necessary to truly achieve complete functional coverage. in about half of all first structures for pfam families were from sg rather than sb. marsden et al. analyzed complete genomes in the gene3d resource to provide new insights into protein family space. the number of protein families was found to be continually expanding with time but a significant proportion of the proteomes could be assigned to relatively few large, well-characterized domain families while the vast majority of domain families were relatively rare and often species specific. it was suggested that sg could provide structures for fewer than a thousand pfam families to achieve reasonable structural coverage of genomes. within these large families it was suggested that multiple structures could be determined to reveal more about the evolution of the family and enable greater understanding of how function evolves. similarly, the selection of targets from within large and diverse cath superfamilies was proposed to maximize modelling leverage. since many of the targets previously selected as likely to have a novel fold proved to have a known fold this could be seen as a more rational use of resources. there was a coordinated strategy of target selection by the four large-scale centres during psi- described by dessailly et al., involving four main approaches: targeting representatives from large, structurally uncharacterized protein domain superfamilies; targeting structurally uncharacterized subfamilies in very large and diverse superfamilies with incomplete structural coverage; community nominated targets; and biomedically important targets. the first two approaches were primarily aimed at increasing the structural coverage of protein sequence space and many targets were selected regardless of their function. metrics have been calculated in this and several other publications to show that this strategy has had the result that psi's per structure contribution to novel structural leverage has been much higher than that for sb, as would be expected since sb is not aiming to increase sequence coverage but rather provide biological insights. the structural biology knowledgebase offers an easy way of keeping abreast of developments by the psi and by sg and sb in general. it is a continually updated portal to research data and other resources from the psi. metrics are regularly updated giving a summary of the total number of structures, the number of distinct and novel structures determined by the psi, and the modelling leverage that psi structures provide. models are made available through the protein model portal. the jcsg has made a notable effort towards the annotation of psi structures with the open protein structure annotation network, a wiki-based portal for the dissemination of information for the broad biomedical community. the increased appearance of structures of proteins of unknown function in the pdb due to sg has stimulated the development of computational methods of function prediction from structure. profunc, developed for the mcsg, combines a number of sequence-based and structure-based methods to gain clues about the function of a protein. the mcsg psi- structures were used to test and develop profunc. when a sequence match is weak and/or multiple functions are suggested, the structure-based methods within profunc can help select from the options that are presented and increase confidence in a particular putative function. another recent and notable method that has been used for function prediction from psi structures is flora. this uses structural motifs associated with different functional sub-groups within functionally diverse cath domain superfamilies. the mcsg selected targets from a broad range of pathogenic bacteria according to selection criteria that have evolved throughout the first two stages of the psi. here we attempt to extend the functional annotations associated with mcsg structures by employing a diverse set of bioinformatics tools and resources. we also analyse the structural novelty of protein structures solved by the mcsg. we use sequence-based methods to assess the proportion of mcsg proteins that have a known function, a putative function, a possible function, and no known function. this is followed up with profunc analysis to support some of the possible functional annotations and in some cases to suggest possible functions for the unknowns. methods a list of all protein structures released by the pdb, their release date, and the source sg centre where applicable was downloaded from the pdb on the 3rd march in all subsequent analyses only data that were available on this date are used. at the time mcsg protein structures for targets had been released by the pdb. annotations the source of most of the annotations used in this work is gene3d.. this database contains nearly million protein sequences including all uniprot sequences and most complete genomes. annotations are imported from multiple sources including ec numbers from uniprot, go terms, kegg genes and protein family assignments from pfam, tigrfams, smart, and panther. kegg currently has reference pathway maps. these are represented in a general way to be applicable to all organisms and thus are useful in the analysis of sg structural coverage. each node in a reference pathway is represented by one or more manually curated kegg orthologs. where there are multiple orthologs these may represent different subunits or different versions of the same enzyme. also a node may belong to more than one pathway when different pathways interact. kegg genes and orthologs are mapped to genomes in gene3d and pdb chains are mapped to uniprot entries using http://www.ebi.ac.uk/thornton-srv/databases/pdbsum/data/pdb_chain_sp_ec the data for which comes from the pdbe database. in addition to collecting these annotations from gene3d for mcsg structures, annotations are considered for inheritance by sequence-sequence comparison, sequence-profile comparison, and profile-profile comparison, such that increasingly remote sequence similarities are considered. not all assignments to a protein family lead to a functional annotation. some assignments are filtered out as being uninformative e.g. assignments to a pfam family described as a duf. the sequence of each mcsg structure is compared to all sequences in gene3d using blast and five iterations of psi-blast, both with an e-value cut-off of. ec numbers and molecular function go terms associated with matching sequences are considered for inheritance. both of these types of annotation are hierarchic. ec numbers all have four levels, the first level refers to the enzyme class, the second level refers to the type of bond or group that is acted on, and the next two levels give progressively more specific details of the catalyzed reaction and its substrates. go molecular function terms are organized as a directed acyclic graph where there may be many nested levels with each child having an'is a' relationship to its parent. annotations associated with a sequence or sequence match are divided into deep annotations comprising level or ec numbers and go terms at level or above, and the remainder which are generally less specific annotations. hidden markov models are built using the sequence of each mcsg structure as a seed for sam-t and then each of these hmms is compared to all pfam hmms using prc. database annotations, inherited annotations, and annotations from the headers of the pdb files are presented together in a table at http://www.biochem.ucl.ac.uk/cgi-bin/dlee/mcsg_annotations used by the authors to manually assign an annotation status to each mcsg structure. annotations are organized from left to right to reflect decreasing levels of both depth and reliability. four levels of status are assigned; known, putative, possible, and unknown. go provides evidence codes for functional annotations where the most reliable codes are the experimental evidence and author statement codes. 62% of mcsg solved targets have a molecular function go term associated with them but only 2% have such a high quality functional annotation, the remainder have annotations that are electronically inferred. there are currently very few other annotation resources where evidence codes are available and sequence coverage is low. we consider possession of one of the most reliable go evidence codes as too strict a criterion to define the'known' status. rather, possession of a deep database annotation and preferably agreement with at least one set of protein families is generally taken as evidence that the function is regarded as being known. note that a go term level of or higher is only a rough guide to the depth of an annotation and the association of a term such as go: for iron ion binding is not considered to be sufficient in itself for a'known' status. this is one reason why manual assessment of status is so important. the'putative' status is generally assigned where there is no deep database annotation but there is a protein family assignment and the family description includes specific molecular functions. the'possible' status is generally assigned where only inherited annotations are available but without the added confidence provided by assignment to a family in a curated resource. there may be a choice between a variety of deep annotations or perhaps only a very general, non-specific annotation is available presenting a broad range of possible deep annotations. the'unknown' status is applied where sequence methods are unable to provide any clue to the function of the protein. profunc analysis is linked to the online table for structures assigned a'possible' or'unknown' status. go slims are cut-down versions of the go ontologies containing a subset of the terms in the whole go. the go provides a generic go slim that gives a broad overview of all function categories without the detail of the specific fine grained terms. this is useful for visualization of the broad functional coverage of mcsg structures and for comparison to the functional coverage of the pdb as a whole. a sequence unique subset that approximately represents all pdb protein entries of known function is generated by selecting representatives annotated in gene3d with a molecular function go term at level six or higher. go terms associated with these representatives are then mapped to the generic go slim using the go slimmer tool at http://amigo.geneontology.org/cgi-bin/amigo/slimmer. similarly the uniprot entries for all mcsg solved targets with known function manually assigned as above are mapped to the generic go slim. profunc analysis the profunc web server at http://www.ebi.ac.uk/profunc employs a number of complementary function-prediction methods, focusing, in particular, on methods based on 3d structural information. this accepts a pdb file and runs a number of sequence- and structure-based analyses on it, listing any hits it finds to existing sequences and structures in descending order of significance. these can help researchers identify any strong similarities which may be indicative of the protein's function. of profunc's structure-based methods, the two that have been shown most accurately to suggest function are the ssm fold-matching algorithm and the "reverse template" search. both these methods have a similar success rate and are able to find distant homologues when simple sequence-based searches fail. the majority of their hits overlap, but occasionally one method picks up a match that the other fails to identify. the advantage of the reverse template method over the fold-match is that its matches are local and the top-scoring hits tend to pick up sites of functional relevance. to demonstrate the added value of 3d structure for proteins of uncertain or unknown function, we run profunc on the full data set of mcsg structures. profunc results for targets assigned a'possible' or'unknown' status following sequence analysis are linked to the online table. profunc results for targets of'unknown' status are manually assessed to determine the improvement in functional coverage afforded by profunc for mcsg solved targets. to further illustrate the value of profunc analysis we then focus on mcsg structures where no significant match to a known 3d structure in the pdb can be found from a search based on sequence alone, but is detectable using the structure-based reverse template method. we investigate which parts of the matched structures the method identifies as the most significantly similar. reverse templates the reverse template method takes the query protein structure and breaks it up into a large set of 3-residue templates. each template consists of three neighbouring residues chosen such that the closest atom-atom distance between any two of the residues is no larger than Å. templates containing more than one hydrophobic residue are rejected in order to bias the templates towards surface residues. a scan of each template against a representative set of the structures in the pdb is then made using the jess algorithm to identify structures containing similar constellations of the three residues. the representative set is a non-redundant list of protein chains in the pdb downloaded weekly from the pisces ftp server. hits are scored and assigned e-values as described in laskowski et al., the scoring being based on the local similarity of the equivalent regions when the structures are superposed on the matched residues. the aim of the method is to pick up local structural similarity between two proteins and, in particular, to identify regions that have been conserved by evolution and hence are more likely to be functionally important. thus, even though the sequences of two homologous proteins might have diverged considerably over evolutionary time, it is possible that, to maintain their function, the proteins' active sites have undergone less change. the detection of such local conservation can provide strong support for the proteins having retained the same, or a similar, function. residues and a Å cut-off were found through empirical testing to give a good balance between the duration and success of a search. identification of functionally important regions to demonstrate this tendency to pick up functionally important regions, we select a subset of the mcsg data where a simple fasta search against the sequences in the pdb fails to return a significant result. we use an e-value cut-off of, which gives us structures whose best match is at or below this threshold. we then consider the top four hits returned by the reverse template search for each of these structures to see if the template residues have any association with the protein's function. we use the functional annotations given in pdbsum, namely: residues belonging to a prosite pattern; catalytic residues, as determined from the catalytic site atlas; residues defined in the site records of the original pdb file; or residues in contact with a ligand, metal or dna/rna in the structure. however, as the search data set used by the reverse template method is a representative one, it is possible that the pdb entry matched by the search has no functional annotation, yet a closely related pdb entry, excluded from the data set, does have some annotation. there were nearly protein chains in the representative data set, compared with over structures in the pdb. thus, to pick up such'lost' annotations, we use the sequence annotated by structure server to find additional functional annotations for each of the matched structures. sas identifies closely related sequences from the pdb using fasta. functional annotations are then imported from the resultant alignment where the sequence identity is at least 30%, the alignment overlap at least residues, and the fasta e-value <. functional novelty the release dates of all pdb entries that are associated with a go term, are assigned to a pfam family, or map to kegg orthologs are compared to find the earliest example released by the pdb for each annotation. each structure that represents a first example is then categorized as being solved by a psi centre or as'non-psi' in order to rate the performance of psi compared to the combined effort of all other laboratories, and to compare the performance of individual psi centres. the dates that the first structures were solved are also divided into the two psi periods psi- and psi- a measure of novelty per pdb release is calculated and for comparison the novelty per non-redundant structure at 95% sequence identity following clustering using cd-hit is also given. structural novelty the release dates of all pdb entries containing a domain assigned to a superfamily in cath v are compared to find the earliest entry in that superfamily and fold. structural novelty within a superfamily is determined by calculating a normalized rmsd. a domain with a normrmsd of 5Å or more from all cath domains in structures previously released by the pdb is defined as structurally novel and belongs to a new structural sub-group. both mcsg domains assigned to cath and unassigned mcsg pdb chains are compared to a representative library of cath domains using cathedral as a fast filter for significant structural similarity. the best matches by cathedral are then structurally aligned and superposed using the more computationally expensive and accurate method ssap and the rmsd is calculated. the normrmsd is then calculated as follows: where rmsd is the root mean square deviation of the superposition, max is the length in residues of the longest domain in the superposition, and nmat is the number of aligned residue pairs. novel modelling leverage novel modelling leverage is calculated according to the method of nair et al.. the sequences of all pdb protein chains, non-redundant at 100% sequence identity, are compared to uniref sequences downloaded on the 3rd of march using psi-blast with three iterations and an e-value cut-off of 1e- novel modelling leverage in residues is determined for each pdb entry on the date of its release by the pdb. the novel modelling leverage of the four large-scale psi centres is compared to each other and to the leverage of the combined non-psi laboratories. human non-synonymous single nucleotide polymorphisms the mcsg focused on solving structures of proteins from pathogenic bacteria and has only solved structures for two human proteins. however, a number of mcsg structures may be used as templates for comparative modelling of human proteins. all such structures identified in the modelling leverage analysis described above are used to analyze human non-synonymous nucleotide polymorphisms. mcsg structures are included regardless of whether a closer template is available from another laboratory. for each human modelling target identified above, variants are retrieved from uniprot, ensembl and omim and filtered for uniqueness. uniprot variants and ensembl gene ids are listed for each protein in the uniprot text file available from uniprot. the ensembl api is used to query their mysql database and retrieve all ensembl variants. omim variants are linked to uniprot ids through the omim missense server at http://www.bioinf.org.uk/omim/. each variant is inspected to see if it falls within a region having a template match and is thus modellable.  annotations the functional annotation status of the mcsg solved targets is as follows: 31% known; 48% putative; 14% possible; and 7% unknown. profunc analysis of the unknowns is summarised in figure b. the majority of structures are therefore of direct relevance to biological investigations of protein function as well as fulfilling some of the goals of sg in expanding structural coverage of fold and function space as shown below. the apparent discrepancy between the results obtained by looking at ec numbers in figure and go terms in figure needs explaining. ec numbers are very high quality annotations but coverage of the pdb is less than that of go and some catalytic activity is not represented by ec numbers. within the go catalytic activity category an over-representation of transferase activity and an under-representation of hydrolase activity is seen which does broadly agree with the ec analysis. within the go binding category mcsg structures also show particularly significant under-representation of protein binding and calcium binding while nucleotide binding is actually a little over-represented. profunc analysis profunc's template searches are rated as certain, probable, possible, or long shots to indicate the likelihood that a prediction is correct. profunc analysis is linked to all structures assigned a "possible" or "unknown" status in the online table at http://www.biochem.ucl.ac.uk/cgi-bin/dlee/mcsg_annotations profunc suggests functions for all of the mcsg targets that are of unknown function following sequence analysis. the overall likelihood of the best scoring profunc prediction being correct is shown in figure b. one target has a prediction rated as being certain while another 17/ targets have predictions rated as being probable. these predictions are all made by the reverse templates method while two of the probable predictions are also made by ligand-binding templates. profunc makes a valuable contribution to functional annotation coverage of mcsg structures. profunc example an example of a reverse template match is pdb entry 2aua which is the structure of bc a protein from bacillus cereus that is classified as'unknown' by the sequence analysis described above. the reverse template search finds one'certain' match, one'probable' match and seven'possible' matches. figure shows the template residues and surroundings for three of these matches, all being to the catalytic domains of bacterial toxins. the first is diphtheria toxin, pdb code 1f0l, and the other two are'possible' matches: one to exotoxin a from pseudomonas aeruginosa, pdb code × k and the second to cholix toxin from vibrio cholera, pdb entry 3ess. each of the templates includes one or more of the known catalytic residues, and the local region returned by the template match corresponds to the matched protein's known substrate binding site. the catalytic domains correspond to ec., which contains the nad+-diphthamide adp-ribosyltransferases. all have the same alpha-beta fold, cath., corresponding to domain of diphtheria toxin. these matches strongly suggest that the bc protein has a very similar function to these toxins. bacillus cereus produces a number of known enterotoxins that are responsible for food poisoning. possibly, this represents another. a second profunc example is presented to show how profunc analysis can add weight and specific detail to a possible annotation suggested by sequence analysis. functional novelty so far out of pfams, out of go molecular function terms and out of terms at level or higher have at least one representative with a solved structure. note, however, that because it is generally only small pfam families that lack a structural representative, 88% of unique pfam assignments in complete genomes in gene3d belong to a pfam family with a structural representative. table gives an indication of the value of pdb releases in terms of functional novelty. from an sg point of view the number of novel structures is given as a percentage of the total number of structures. sb laboratories often solve multiple structures of the same protein with different ligands bound, at a different ph, or with single amino acid substitutions etc. these investigations have different goals from those of sg but here the intention is to assess the performance of psi centres in terms of sg goals. for an sb point of view the number of novel structures is given for comparison as a percentage of the number of non-redundant structures at 95% sequence identity. although psi centres did not specifically target structurally uncharacterized go terms their overall performance is at least as good as that of all non-psi laboratories combined in terms of sg goals. numbers of first structures and their percentage of total structures and non-redundant structures at 95% sequence identity for molecular function go terms and pfam families during psi, psi- and psi- pfam families, on the other hand, were targeted by psi centres and the proportion of first structures solved out of all structures solved is significantly higher for psi than for non-psi laboratories as would be expected. as would also be expected, the proportion of first structures for pfams are considerably higher than those for go terms since go terms are a much more general functional category compared to the evolutionarily related members of pfams. the further back in time that we go the fewer is the number of go terms and pfams that were structurally characterized and thus there was a higher probability of solving the first structure for a term or family. as structural coverage becomes more complete the probability of solving the first structure for a term or family diminishes. this is one likely reason that the rate of solving first structures during psi- is higher than that during psi- for both psi and non-psi laboratories. another possible explanation is that the easier functions have been solved and the remaining functions are more difficult to obtain a structure for. perhaps we are seeing a combination of both factors. by virtue of solving more structures than any other individual psi centre the mcsg also solved the highest number of first structures for both go terms and pfams. in both cases the proportion of first structures is close to the average for all psi centres although they excelled at solving the first structures for pfams during psi- 39/ of the mcsg first structures for go molecular function terms correspond to a level ec number and of the remaining also have catalytic activity e.g. go: amino acid-transporting atpase activity, again showing a bias towards enzymes. the ec and go annotation systems have different coverage from each other. analysis of ec numbers in gene3d shows that the mcsg solved the first structure for rather than level ec numbers. of these ec numbers were not revealed by the go analysis so in total the mcsg solved the first structure for level ec numbers. the novelty of these structures may also be considered in terms of providing structural coverage for metabolic pathways. pathways there are currently kegg reference pathways. these pathways have nodes, most but not all of which correspond to ec numbers. kegg orthologs representing evolutionarily related groups of proteins map to many but not all of these nodes and many nodes have multiple orthologs. currently, orthologs are assigned to the reference pathways and more are likely to be added in the future. out of a total of currently possible level ec numbers map through kegg orthologs to kegg reference pathway nodes. / have at least one structural representative compared to / for the ec classification as a whole. 9/ pathways have complete structural coverage of their orthologs but these are all small pathways, the largest having only orthologs belonging to nodes, ko synthesis and degradation of ketone bodies. this pathway is also the largest pathway composed entirely of nodes corresponding to level ec numbers with complete structural coverage. the nodes with ec numbers that the mcsg solved the first structure for belong to kegg reference pathways but this number is a coincidence, with some of the nodes appearing in multiple pathways while multiple nodes also appear in single pathways. the largest number of nodes within a single kegg reference pathway to gain a first structure from the mcsg is for pathway ko arginine and proline metabolism, but this is a large pathway with more than nodes. psi centres did not select targets with the aim of achieving complete coverage of metabolic pathways but this is a potential goal for the future. overall the mcsg solved the first structure for orthologs in pathways. of these orthologs belong to multiple pathways and thus join pathways together. membership of multiple pathways is, however, not rare with about a half of all kegg orthologs belonging to more than one reference pathway. structural novelty of the mcsg pdb entries analyzed in this work only are so far classified in cath v and so the analysis of structural novelty is currently incomplete. also, since a special effort was made by the curators of cath to update the classification of mcsg structures as a priority over other structures, mcsg structural novelty cannot be compared to that of other psi centres or sb in general. mcsg structural novelty is shown in table number of novel cath structural domains solved by the mcsg. three levels of novelty are defined: a domain with a novel fold is the first example released by the pdb of a fold in cath; a domain belonging to a novel superfamily is the first example released by the pdb of a superfamily in cath; a domain belonging to a novel structural sub-group has a normalized rmsd of at least 5Å from all previously existing domains in cath. by definition all new folds are also new superfamilies and all new superfamilies are also new structural subgroups. however, in this analysis not all new ssgs are classified in cath or necessarily belong to existing superfamilies or folds. of the domains that belong to new ssgs but are not classified in cath as new folds or superfamilies, are found to be classified into pre-existing cath superfamilies while are not classified in cath v. of the not classified achieve a score > from the structural comparison program ssap following the cathedral scan suggesting that they probably do belong to pre-existing superfamilies. the remaining domains could potentially represent new superfamilies or folds. novel modelling leverage the average novel modelling leverage in residues per structure is considerably higher during psi- compared to psi- for both psi and non-psi. this is despite total modelling leverage per structure actually being higher for both psi and non-psi during psi- compared to during psi-. the reason is that as the cumulative total and diversity of solved structures increases and targets with good modelling leverage are preferentially selected in the early stages of psi, any new structures that are solved have less chance of providing novel modelling leverage for a fixed set of protein sequences. during both psi- and psi- psi substantially outperformed non-psi in terms of novel modelling leverage per structure as would be expected since non-psi often target multiple forms of a given target protein to better characterise the biology. with the increased output of psi during psi- the contribution of psi to total worldwide novel modelling leverage more than doubled from % during psi- to % during psi- during both psi- and psi- the mcsg substantially outperformed the other three large-scale centres in terms of both total novel modelling leverage and average novel modelling leverage per structure. novel modelling leverage for residues in uniref during psi, psi- and psi- human non-synonymous single nucleotide polymorphisms in the modelling leverage analysis mcsg structures are identified as templates for human proteins in uniprot. in this analysis of modellable nssnps all modelling leverage of mcsg structures is considered, not just novel leverage, and the mcsg structure is not necessarily the closest available template. within these human protein sequences unique ensembl nssnps, unique uniprot nssnps, and unique omim nssnps are identified. of these nssnps are within a modellable region when using an mcsg structure as a template. a good example of an mcsg structure helping to explain human disease is illustrated by pdb entry 2hma. human mitochondrial trna-specific 2-thiouridylase matches 2hma with 42% sequence identity in a blast alignment. the enzyme has been implicated in aggravating mitochondrial 12s ribosomal rna aminoglycoside-induced and non-syndromic deafness. it catalyzes the 2-thiolation of uridine at the wobble position of mitochondrial trna, trna and trna. a few natural variants have been identified and one has been linked to a decrease in enzymatic activity. to attempt to explain this, a homology model is built using 2hma as a template which includes the substrate, s-adenosyl methionine. when inspecting the model, it is clear that ala is part of the sam binding pocket of the enzyme. in 2hma the adenosine group fits neatly into this pocket where ala is located. in the ala10ser mutation, an additional hydrogen bond between the protein and sam is probably introduced. this increases the binding affinity between sam and the protein and thus slows down the release of sam once the sulphur has been transferred. this agrees with observations made by guan et al. that the ala10ser mutation reduces activity.  the mcsg has performed well during the first two phases of the psi in terms of the goals that were established during this period. whilst there has been concern that the psi produced too many structures of unknown function our analysis reveals that by using a range of bioinformatics tools and resources we are able to provide functional annotations for more than 90% of the structures solved by the mcsg. structures of unknown function have helped stimulate the development of methods such as profunc that can predict function from structure. profunc analysis complements sequence analysis by both adding weight and specific detail to predicted function and by suggesting function where sequence methods have failed to do so. authors' contributions rl carried out the profunc analyses, tdeb carried out the nssnp analyses, and dl carried out all other analyses and drafted the manuscript. all authors participated in the design of the study, helped draft the manuscript, and read and approved the final manuscript. supplementary material additional file 1 additional profunc example. results and discussion, figure, and reference showing how profunc can be used to refine, add detail to, and support a protein function prediction. click here for file acknowledgements this work was funded by the nih as part of the protein structure initiative.
structures of unknown function have stimulated the development of methods for function prediction from structure. we show that the mcsg has met the stated goals of the psi and use online resources and readily available function prediction methods to provide functional annotations for more than 90% of the mcsg structures. the structure-to-function prediction method profunc provides likely functions for many of the mcsg structures that cannot be annotated by sequence-based methods. the midwest center for structural genomics is one of the large-scale centres of the protein structure initiative. although the focus of the psi was structural coverage, many of the structures solved by the mcsg can also be associated with functional classes and biological roles of possible biomedical value. a criticism of structural genomics is that target selection strategies mean that some structures are solved without having a known function and thus are of little biomedical significance. during the first two phases of the psi the mcsg has solved over a thousand protein structures.
 the midwest center for structural genomics is one of the large-scale centres of the protein structure initiative. during the first two phases of the psi the mcsg has solved over a thousand protein structures. a criticism of structural genomics is that target selection strategies mean that some structures are solved without having a known function and thus are of little biomedical significance. structures of unknown function have stimulated the development of methods for function prediction from structure.  we show that the mcsg has met the stated goals of the psi and use online resources and readily available function prediction methods to provide functional annotations for more than 90% of the mcsg structures. the structure-to-function prediction method profunc provides likely functions for many of the mcsg structures that cannot be annotated by sequence-based methods.  although the focus of the psi was structural coverage, many of the structures solved by the mcsg can also be associated with functional classes and biological roles of possible biomedical value. 