Introduction {#s1}
============

The cupin superfamily of proteins possesses remarkable functional diversity with representatives found in Archaea, Eubacteria, and Eukaryota [@pone.0074477-Dunwell1], [@pone.0074477-Dunwell2], [@pone.0074477-Dunwell3], [@pone.0074477-Dunwell4], [@pone.0074477-Dunwell5]. The identification of the cupin superfamily was originally based on the recognition that the wheat protein germin shared a nine amino acid sequence with another protein, spherulin, produced by the slime mold *Physarum polycephalun* during starvation [@pone.0074477-Dunwell3]. This sequence similarity was also observed in a number of seed storage proteins called germin-like proteins (GLPs). Knowledge of the three dimensional structures of these proteins led to the collective name "cupin" on the basis of their β-barrel shape ("cupa" means "small barrel" in Latin) [@pone.0074477-Dunwell4]. Characteristic features of proteins with this fold include high thermal stability and resistance to proteases. These features are consistent with their high degree of subunit contacts, hydrophobic interactions, and short loops. The cupin domain was originally described as two conserved motifs, each composed of two β-strands [@pone.0074477-Dunwell2], [@pone.0074477-Dunwell5]. Motif 1 was designated as G(X)~5~HXH(X)~3,4~E(X)~6~G and Motif 2 as G(X)~5~PXG(X)~2~H(X)~3~N. The two motifs are separated by an intermotif region which ranges from 15 to 50 amino acids long. With more sequences analyzed, it became clear that the primary sequence of the two motifs is not as highly conserved as previously thought [@pone.0074477-Dunwell1], [@pone.0074477-Woo1].

Genomic enzymology; large-scale sequence, functional, and structural databases; and inferences from molecular evolution have informed family, superfamily, and suprafamily designations [@pone.0074477-Anantharaman1], [@pone.0074477-Gerlt1], [@pone.0074477-Glasner1], [@pone.0074477-Galperin1]. Enzyme families are relatively recently-diverged groups of enzymes that share similar three-dimensional structures and functions. Enzyme superfamilies consist of enzymes that diverged earlier and possess fewer common elements but share a common mechanistic attribute, while suprafamilies share conserved residues but do not share a common mechanistic attribute [@pone.0074477-Gerlt1], [@pone.0074477-Allewell1]. Inconsistencies exist in the usage of the term 'cupin'. According to the Structural Classification of Proteins (SCOP) database [@pone.0074477-Murzin1], [@pone.0074477-LoConte1], [@pone.0074477-Andreeva1], cupin proteins are members of the 'RmlC-like Cupins' superfamily within the double-stranded β-helix (DSBH) multicatalytic fold [@pone.0074477-Agarwal1]. The term 'cupin superfamily' has often been used to refer to those proteins defined by the SCOP database as well as the 2-oxoglutarate-Fe^2+^-dependent dioxygenase superfamily that also possesses the DSBH fold [@pone.0074477-Dunwell1], [@pone.0074477-Dunwell5]. However defined, the cupin superfamily is extremely diverse and includes catalytically inactive seed storage and sugar-binding metal-independent proteins as well as metal-dependent enzymes possessing dioxygenase, decarboxylase, and other activities [@pone.0074477-Galperin1]. [Figure 1](#pone-0074477-g001){ref-type="fig"} (A--J) shows the structures of 10 representative proteins designated as cupins in the Pfam database. Although the majority of enzymatic cupins contain iron as an active site metal, other members may contain nickel, zinc, manganese, cobalt, or copper. Each cofactor allows a different type of chemistry to occur within the conserved tertiary structure. Proposed reaction mechanisms of the metal-dependent cupins generally involve sequential binding of the substrate and dioxygen to the catalytic metal cation [@pone.0074477-McCoy1]. [Figure 2](#pone-0074477-g002){ref-type="fig"} (A-J) shows the metal-binding sites of 8 representative members of the enzymatic cupins.

![Structures of representative of members of the cupin superfamily.\
A. oxalate oxidase (PDB code: 2et1) [@pone.0074477-Opaleye1], B. oxalate decarboxylase (PDB code: 1uw8) [@pone.0074477-Just1], C. seed storage protein Ara h (PDB code: 3s7i) [@pone.0074477-Chruszcz1], D. NovW, a 4-keto-6-deoxy sugar epimerase (PDB code: 2c0z) [@pone.0074477-Jakimowicz1], E. cysteine dioxygenase (PDB code: 2q4s) [@pone.0074477-Levin1], F. phosphomannose isomerase (PDB code: 1 pmi) [@pone.0074477-Cleasby1] G. acireductone dioxygenase (PDB code: 1 zrr) [@pone.0074477-Pochapsky2] H. taurine/alpha-ketoglutarate dioxygenase (PDB code: 1os7) [@pone.0074477-OBrien1], I. hypoxia-inducible factor 1-alpha inhibitor (PBD code: 2y0i) [@pone.0074477-Yang1], J. lysine-specific demethylase 6B (PDB code: 2 xue) [@pone.0074477-Kruidenier1]. β-sheets are shown in green, α-helices are shown in red, and random coils are shown in grey. Spheres represent bound metal ions. Figures were generated using Pymol (The PyMOL Molecular Graphics System, Schrödinger, LLC).](pone.0074477.g001){#pone-0074477-g001}

![Metal-binding sites of representative members of the cupin superfamily.\
A. Mn ion of oxalate oxidase coordinated by His88, His90, Glu95, and His137 (PDB code: 2et1) [@pone.0074477-Opaleye1], B. N-terminal Mn ion of oxalate decarboxylase coordinated by His95, His97, Glu101, and His140 (PDB code: 1uw8) [@pone.0074477-Just1], C. C-terminal Mn ion of oxalate decarboxylase coordinated by His273, His275, Glu280, and His319 (PDB code: 1uw8) [@pone.0074477-Just1], D. Ni ion of cysteine dioxygenase coordinated by His86, His88, and His140 (PDB code: 2q4s) [@pone.0074477-Levin1], E. Zn ion of phosphomannose isommerase coordinated by Gln111, His113, Glu138, and His285 (PDB code: 1 pmi) [@pone.0074477-Cleasby1] F. Ni ion of acireductone dioxygenase coordinated by His96, His98, Glu102, and His140 (PDB code: 1 zrr) [@pone.0074477-Pochapsky2], G. Fe ion of taurine/alphaketoglutarate dioxygenase coordinated by His99, Asp101, and His255 (PDB code: 1os7) [@pone.0074477-OBrien1], H. Fe ion of hypoxia-inducible factor 1-alpha inhibitor coordinated by His199, Asp201, and His279 (PBD code: 2y0i) [@pone.0074477-Yang1], I. Fe ion of lysine-specific demethylase 6B coordinated by His1390, Glu1392, and His1470 (PDB code: 2 xue) [@pone.0074477-Kruidenier1], J. Zn ion of lysine-specific demethylase 6B coordinated by Cys1575, Cys1578, Cys1602, and Cys1605 (not part of the cupin domain) (PDB code: 2 xue) [@pone.0074477-Kruidenier1]. β-sheets are shown in green, α-helices are shown in red, and random coils are shown in grey. Spheres represent bound metal ions. Figures were generated using Pymol (The PyMOL Molecular Graphics System, Schrödinger, LLC).](pone.0074477.g002){#pone-0074477-g002}

The DSBH (also referred to as jelly roll) topology is most often composed of eight β-strands that form a β-sandwich structure comprised of two four-stranded antiparallel β-sheets. Although four forms are possible (two hands of the helix and two directions to trace the structure), only one form (right-handed class I) is prevalent in nature [@pone.0074477-Stirk1]. Ancestral cupins can be evolutionarily reconstructed as simple, small molecule-binding domains that likely bound sugars and cyclic nucleotides [@pone.0074477-Dunwell5], [@pone.0074477-Anantharaman1], [@pone.0074477-Anantharaman2]. These sugar-binding domains later gave rise to sugar-modifying domains such as isomerases and epimerases [@pone.0074477-Woo2]. Analyses of the evolution of the fold suggest that a set of conserved histidine residues employed in sugar-binding in the ancestral non-enzymatic domain evolved into the metal-coordinating histidine residues observed in oxalate oxidase ([Figures 1A](#pone-0074477-g001){ref-type="fig"} and [2A](#pone-0074477-g002){ref-type="fig"}) [@pone.0074477-Woo2] and oxalate decarboxylase ([Figures 1B](#pone-0074477-g001){ref-type="fig"}, [2B, and 2C](#pone-0074477-g002){ref-type="fig"}) [@pone.0074477-Anand1]. Another lineage of DSBH domains acquired a new set of conserved residues with the ability to bind 2-oxoglutarate which gave rise to the 2-oxoglutarate-Fe^2+^-dependent dioxygenases [@pone.0074477-Anantharaman1], [@pone.0074477-Iyer1].

The exponential growth of structural information for proteins provides abundant material for the analysis of how protein structure informs biological function and chemical reactivity. Babbitt *et al.* recognized the need to associate structure and sequence information with biological function in ways that are accessible to both experimental and computational biologists. They presented protein similarity networks (PSNs) to fulfill this need [@pone.0074477-Atkinson1]. PSNs have contributed to our understanding of a number of large groups of proteins including the enolase superfamily [@pone.0074477-Gerlt2], the ePK-like superfamily [@pone.0074477-Brown1], glutathione transferases [@pone.0074477-Atkinson2], [@pone.0074477-Atkinson3], strictosidine synthase-like proteins [@pone.0074477-Hicks1], cysteine peptidases [@pone.0074477-Atkinson4], and proteins used in algal metal transport [@pone.0074477-BlabyHaas1]. These studies have yielded meaningful insights, validated PSN methodology, and provided an understanding of the caveats and limitations of PSNs. PSNs are complementary to phylogenetic studies and provide different and new information compared to other methods relating structures and sequences. It has been noted that protein similarity networks are most compelling when painted with functional or structural information that is orthogonal to the data used to generate the networks [@pone.0074477-Atkinson1].

The Pfam database [@pone.0074477-Punta1] lists 112,082 cupin sequences represented in 6529 species and 945 associated protein structures. This represents a greater than 10-fold increase in the number of sequences in only four years [@pone.0074477-Agarwal1]. Protein similarity networks are a way to visualize large-scale computational analyses of sequence and structure among a given set of proteins [@pone.0074477-Atkinson1], [@pone.0074477-Brown1] and have been used to guide experimental design and data interpretation [@pone.0074477-Hicks1], [@pone.0074477-Pieper1]. In this work we describe the first visualization of sequence and structure data of the cupin superfamily using PSNs. The PSNs were built using open source software programs Pythoscape [@pone.0074477-Barber1] and Cytoscape [@pone.0074477-Smoot1]. In these networks the protein structures (or sequences) are represented as nodes, and their similarities to each other are represented as edges. A value of this work lies in displaying nodes by particular attributes. Node attributes that can be gathered automatically through Pythoscape include sequence, organism, description, and identification codes. Other manually-retrieved attributes may be inputted, such as the presence of a catalytic motif or reaction intermediates. Overlaying attribute information onto clusters allows visualization of the relationship between structure and function. There exists a large amount of structure and sequence information for cupin superfamily members for which there has been no experimental work conducted. By facilitating the prediction of the functions of these proteins, PSNs may be used to guide experimental inquiry into the cupin superfamily.

Results and Discussion {#s2}
======================

The Pfam database [@pone.0074477-Punta1] contains 484 unique cupin structures. For visual simplicity, only a single representative structure of a protein with multiple structures was used in the networks. This reduced the networks to 183 structures, which are listed in Supporting Information ([Tables S1--S7 in File S1](#pone.0074477.s001){ref-type="supplementary-material"}). Of the 183 PDB structures representing unique members of the cupin superfamily, 76 bind no metal, 49 bind iron, 18 bind nickel, 16 bind zinc, 10 bind manganese, 1 binds cadmium (protein yhhW, PDB code: 1tq5) [@pone.0074477-Adams1], 1 binds copper (quercetin dioxygenase from *Aspergillus japonicus*, PDB code:1 juh) [@pone.0074477-Fusetti1], 1 binds mercury (protein YhcH, PDB code: 1 jop) [@pone.0074477-Teplyakov1], and 11 bind multiple metals.

[Figure 3](#pone-0074477-g003){ref-type="fig"} shows pairwise similarities for this non-redundant set of structures calculated using TM-align [@pone.0074477-Zhang1] (See Methods) at two different stringency thresholds. This permits alternative views of the same structural relationships. These networks are painted according to the identity of the bound metal in the structure. In the less stringent network ([Figure 3A](#pone-0074477-g003){ref-type="fig"}), all but four structures have connections to a single large constellation of structures. This constellation partitions off into smaller clusters of structures in the more stringent network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}). It can be observed that often proteins of the same function cluster together. For example, at the higher stringency the catalytically inactive seed storage proteins cluster together with two metal-binding cupins, oxalate decarboxylase and the MncA protein. Other clusters analyzed in this work include the cysteine dioxygenases (CDOs), the Jumonji C domain-containing proteins (JmjC), and the RmlC epimerases.

![Structure similarity networks of the cupin protein stuctures colored by metal ligand.\
Pairwise similarities for a non-redundant set of 183 structures from the Pfam cupin clan (CL0029) were calculated using TM-align [@pone.0074477-Zhang1]. Each node represents a structure. Nodes were arranged using the yfiles organic layout of of Cytoscape version 3.0. A. Edges between nodes were drawn only if the average TM-score \>0.53 for that edge. At this cutoff, the average r.m.s.d. is 2.91 Å with an average of 158.0 Cα atoms aligned. B. Edges between nodes were drawn only if the average TM-score \>0.65 for that edge. At this cutoff, the average r.m.s.d. is 2.44 Å with an average of 185.4 Cα atoms aligned.](pone.0074477.g003){#pone-0074477-g003}

Seed Storage Proteins Cluster Together {#s2a}
--------------------------------------

The seed storage proteins cluster together at the lower threshold ([Figure 3A](#pone-0074477-g003){ref-type="fig"}), but share edges with other protein structures. At the higher threshold ([Figure 3B](#pone-0074477-g003){ref-type="fig"}), however, these proteins partition exclusively. Members of this cluster contain two cupin domains and are, therefore, classified as bicupins. The majority of the members of this group contain no metal and include procruciferin (PDB code: 3 kgl) [@pone.0074477-TandangSilvas1], Ara h (PDB code: 3s7i) [@pone.0074477-Chruszcz1], phaseolin (PDB code: 2 phl) [@pone.0074477-Lawrence1], beta-conglycinin (PDB code: 1 uij) [@pone.0074477-Maruyama1], and proglycinin (PDB code: 1 fxz) [@pone.0074477-Adachi1]. The two metal-containing members of this cluster are oxalate decarboxylase (PDB code: 1uw8) [@pone.0074477-Just1] and MncA (PDB code: 2 vqa) [@pone.0074477-Tottey1] with both incorporating manganese. The architecture of this group is represented in [Figure 1B](#pone-0074477-g001){ref-type="fig"} by oxalate decarboxylase (OxDC) and in [Figure 1C](#pone-0074477-g001){ref-type="fig"} by Ara h, of interest as a major peanut allergen. The two nearly identical Mn-binding sites of OxDC are shown in [Figures 2A and 2B](#pone-0074477-g002){ref-type="fig"}. Previous work with protein similarity networks in general had indicated that the clustering of proteins does not change dramatically whether the domain in common is isolated or is a component of a larger multi-domain complex [@pone.0074477-Atkinson1], and this is our observation also. However, we were concerned that monocupins were generally segregated from bicupins in the network built using entire structures. In an effort to explore this possible limitation, we constructed networks of isolated cupin domains for the structures used in [Figure 3](#pone-0074477-g003){ref-type="fig"}. [Figures 4A and 4B](#pone-0074477-g004){ref-type="fig"} show the cupin domain networks at the same stringencies as those in [Figures 3A and 3B](#pone-0074477-g003){ref-type="fig"}, respectively, also colored by metal ligand. At both stringencies, monocupin oxalate oxidase (OxOx) ([Figures 1A](#pone-0074477-g001){ref-type="fig"} and [2A](#pone-0074477-g002){ref-type="fig"}) clusters with both the N-terminal and C-terminal cupin domains of bicupin OxDC. At the higher stringency ([Figure 4B](#pone-0074477-g004){ref-type="fig"}), OxDC (both domains), MncA (both domains), and OxOx partition away from the group labeled "seed storage cluster." Interestingly, many of the seed storage proteins such as Ara h split one domain to the "seed storage cluster" and one domain remains in the group containing OxDC, MncA, and OxOx. This observation is consistent with the conclusions from phylogenetic analyses that the N-and C-terminal domains of many bicupins arose through independent evolutionary events [@pone.0074477-Agarwal1] whereas others such as OxDC arose through gene duplication events [@pone.0074477-Dunwell2].

![Structure similarity networks of cupin domains colored by metal ligand.\
Pairwise similarities for a non-redundant set of 213 domains from the Pfam cupin clan (CL0029) were calculated using TM-align [@pone.0074477-Zhang1]. Each node represents a domain. Nodes were arranged using the yfiles organic layout of of Cytoscape version 3.0. A. Edges between nodes were drawn only if the average TM-score \>0.53 for that edge. At this cutoff, the average r.m.s.d. is 1.73 Å with 74.2 Cα atoms aligned. B. Edges between nodes were drawn only if the average TM-score \>0.65 for that edge. At this cutoff, the average r.m.s.d. is 1.42 Å with 80.1 Cαatoms aligned.](pone.0074477.g004){#pone-0074477-g004}

[Figure 5](#pone-0074477-g005){ref-type="fig"} reproduces the same network as [Figure 3B](#pone-0074477-g003){ref-type="fig"} but is painted by kingdom ([Figure 5A](#pone-0074477-g005){ref-type="fig"}) and function ([Figure 5B](#pone-0074477-g005){ref-type="fig"}). [Figure 5A](#pone-0074477-g005){ref-type="fig"} illustrates that OxDC and MncA are bacterial proteins and the only non-plant members of the seed storage cluster. OxDC is a Mn-dependent enzyme that catalyzes the carbon-carbon bond cleavage of oxalate to yield carbon dioxide and formate in a reaction with no net change in oxidation state (the only lyase represented in [Figure 5B](#pone-0074477-g005){ref-type="fig"}) [@pone.0074477-Tanner1]. MncA is the most abundant Mn-containing protein in cyanobateria *Synechocystis* PCC 6803 and played a key role in a study that elucidated a mechanism whereby the compartment in which a protein is folded overrides its binding preference to control its metal content [@pone.0074477-Tottey1]. MncA is first neighbors with (shares edges) in the structure networks ([Figures 3A, 3B](#pone-0074477-g003){ref-type="fig"}, [5A, and 5B](#pone-0074477-g005){ref-type="fig"}) with all members of the seed storage cluster as well as BacB (PDB code: 3h7j) [@pone.0074477-Rajavel1] in an adjacent cluster. BacB binds both Fe and Co, has been shown to play a role in the synthesis of bacilysin, and clusters with pirins and quercetin dioxygenases (see below). Networks based on sequence ([Figures 6A and 6B](#pone-0074477-g006){ref-type="fig"}) similarly cluster the seed storage proteins which share edges with OxDC and MncA.

![Structure similarity networks colored by species and function.\
This is the same network as in 3B but colored according to A. species and B. function.](pone.0074477.g005){#pone-0074477-g005}

![Sequence similarity networks colored by metal ligand.\
Networks were generated by all-by all BLAST comparisons of the 183 sequences corresponding to the unique cupin structures shown in [Figure 3](#pone-0074477-g003){ref-type="fig"}. Nodes were arranged using the yfiles organic layout of of Cytoscape version 3.0. A. Edges between nodes are drawn only if the E-value is better than of 1E-3.5. At this cutoff, edges at this threshold represent alignments with a median 32.1% identity over 93 residues. B. Edges between nodes are drawn only if the E-value is better than of 1E-6.0. At this cutoff, edges at this threshold represent alignments with a median 36.2% identity over 185 residues.](pone.0074477.g006){#pone-0074477-g006}

Enzymes Cluster Based on Function {#s2b}
---------------------------------

The structure network at a more permissive threshold ([Figure 3A](#pone-0074477-g003){ref-type="fig"}) clearly shows clustering correlated with function. For example, the four cysteine dioxygenases, which catalyzes the oxidation of the L-cysteine to its sulfinic acid by incorporating both of the oxygen atoms from O~2~ into the product, are part of a larger cluster also containing oxalate oxidase (OxOx) and auxin-binding protein at the lower stringency ([Figure 3A](#pone-0074477-g003){ref-type="fig"}). In the more stringent network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}), the CDO structures (*Homo sapiens*, PDB code: 2ic1 [@pone.0074477-Ye1]; *Ralstonia eutropha,* PDB code: 2gm6; *Mus musculus* PDB code: 2q4s [@pone.0074477-Levin1]; and *Pseudomonas aeruginosa* PDB code: 3 uss) share first neighbor status with each other only and are connected to the larger cluster through the putative CDO from *Bacillus subtilis* (PDB code: 3 eqe). These monocupins are represented in [Figure 1E](#pone-0074477-g001){ref-type="fig"} by the *Mus musculus* (mouse) enzyme. The metal ion in this structure is Ni(II) and is coordinated by three conserved residues (His 86, His88, and His140) as shown in [Figure 2D](#pone-0074477-g002){ref-type="fig"}. Biochemical characterization of this protein has shown that iron is required for catalytic activity [@pone.0074477-McCoy1]. CDOs cluster by function when the entire structures are compared ([Figure 3](#pone-0074477-g003){ref-type="fig"}), but remain separate due to the smaller region of alignment when compared at the domain level at the same stringency ([Figure 4](#pone-0074477-g004){ref-type="fig"}, data not indicated with arrows). The majority of cupins for which there is structural information are from bacterial organisms. The CDOs cluster by function in the sequence networks ([Figure 6](#pone-0074477-g006){ref-type="fig"}). These networks confirm and extend previous observations that members of the cupin superfamily cluster by function when the function is defined at a fine level, such as the oxidation of cysteine but do not do so when function is defined more broadly [@pone.0074477-Agarwal1]. For example, the Fe-containing 2-oxoglutarate-dependent dioxygenases such as phytanoyl-CoA dioxygenase, taurine dioxygenase (TauD) ([Figures 1H](#pone-0074477-g001){ref-type="fig"} and [2G](#pone-0074477-g002){ref-type="fig"}), and lysine-specific demethylases (JmjC cluster, see below) are separated in both the structure and sequence networks. Furthermore, the networks painted by metal ion ([Figures 3](#pone-0074477-g003){ref-type="fig"}, [4](#pone-0074477-g004){ref-type="fig"}, and [6](#pone-0074477-g006){ref-type="fig"}) show that neither structure-based nor sequence-based clustering correlates with the identity of the active site metal.

Structures are available for 3-hydroxyanthranilate-3,4-dioxygenase (HAD) from four species. HAD from *Cupriavidus metallidurans* (formerly *Ralstonia metallidurans*) (PDB code:1 yfu) [@pone.0074477-Zhang2], [@pone.0074477-Colabroy1] and *Bos taurus* (PDB code: 3fe5) [@pone.0074477-Dilovic1] contains iron, while HAD from *Saccharomyces cerevisiae* (PDB code:1 zvf) and *Homo sapiens* (PDB code: 2 qnk) contains nickel. These four structures group together by function in the network generated using the whole protein ([Figure 3](#pone-0074477-g003){ref-type="fig"}) but not when the cupin domain alone is used to generate the network ([Figure 4](#pone-0074477-g004){ref-type="fig"}). The HADs cluster by function in the sequence networks ([Figure 6](#pone-0074477-g006){ref-type="fig"}). Similarly, four structures are available for mannose-6-phosphate isomerases (MPIs) ([Figures 1F](#pone-0074477-g001){ref-type="fig"} and [2E](#pone-0074477-g002){ref-type="fig"}). These structures segregate exclusively at the more stringent thresholds of the whole protein structure network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}) and sequence network ([Figure 6B](#pone-0074477-g006){ref-type="fig"}), but not in the domain networks ([Figure 4](#pone-0074477-g004){ref-type="fig"}). Of these, three contain Zn (*Bacillus subtilis*, PDB code: 1 qwr; *Salmonella typhimurium*, PDB code: 3h1m [@pone.0074477-Sagurthi1]; and *Candida albicans*, PDB code: 1 pmi [@pone.0074477-Cleasby1]).One putative PMI was crystallized with no metal present (*Archaeoglobus fulgidus*, PDB code:1zx5).

Quercetin 2,3-dioxygenase (QueD) catalyzes the insertion of molecular oxygen into polyphenolic flavonols. Structural information is available for quercetin 2,3-dioxygenases from *Aspergillus japonicus* (PDB code: 1 juh, [@pone.0074477-Fusetti1]) and *Bacillus subtilis*, (PDB code: 1y3t [@pone.0074477-Gopal1]). QueD is able to employ a variety of metals in catalysis [@pone.0074477-Gopal1], [@pone.0074477-Schaab1], [@pone.0074477-Bowater1]. Incubation of apoprotein with transition metal ions has been done to examine the effects of different metal ions on enzymatic activity. Results yielded an activity profile with trends that were consistent with the Irving-Williams metal ion series based on the stability of metal ion complexes [@pone.0074477-Irving1]. Data suggest that Mn(II) is the preferred cofactor for the enzyme [@pone.0074477-Schaab1]. Pirins have been shown to possess QueD activity [@pone.0074477-Adams1], and structural information is available for pirins from *Escherichia coli* (PDB code: 1tq5) [@pone.0074477-Adams1], *Homo sapiens* (PDB code: 1j1l), and *Geobacillus kaustophilus* (PDB code: 2p17). The QueDs and pirins are members of the same larger cluster in the structure-based networks ([Figure 3](#pone-0074477-g003){ref-type="fig"}). Consistent with previous analyses and observations [@pone.0074477-Agarwal1], the pirin N- and C- terminal domains do not share edges with each other in the domain-based networks, suggesting an independent evolution of these domains in bicupins. This is in contrast to the N- and C- terminal domain clustering of bicupins such as OxDC, QueD, and BacB from *Bacillis subilis* where gene duplication events are proposed to play a large role in i) increasing the size of the genome and ii) producing bicupin architectures [@pone.0074477-Dunwell2].

Acireductone dioxygenase (ARD) is a particularly interesting example of how the cupin scaffold can be used to catalyze different reactions in that this enzyme catalyzes different reactions depending on the type of metal ion bound in the active site [@pone.0074477-Dai1], [@pone.0074477-Dai2]. ARD serves as a branch point in the methionine salvage pathway. This pathway returns the γ-thiomethyl group of methylthioadenosine to methionine [@pone.0074477-Pochapsky1]. The overexpression of the ARD gene in *E. coli* yields two enzymes with different activities that are able to be separated by conventional chromatographic techniques of ion exchange and hydrophobic interaction resins. One enzyme catalyzes the on-pathway oxidation of acireductone to ketoacid and formate; the other enzyme catalyzes an off-pathway oxidation with formation of CO. The only difference between the two enzymes (Fe-ARD and Ni-ARD) was the metal bound in the active sites [@pone.0074477-Dai1]. There is structural information available for ARD from *Mus musculus* (PDB code: 1vr3) [@pone.0074477-Xu1] and *Klebsiella pneumonia* (PDB code: 1 zrr) ([Figures 1G](#pone-0074477-g001){ref-type="fig"} and [2F](#pone-0074477-g002){ref-type="fig"}) [@pone.0074477-Pochapsky2]. In the more permissive structure network ([Figure 3A](#pone-0074477-g003){ref-type="fig"}) both ARD structures, like the RmlC epimerases, are connected to the larger cluster through dTDP-6-deoxy-3,4-keto-hexulose (PDB code: 2pa7) [@pone.0074477-Davis1] but lose all edges in the more stringent network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}). Similarly, in the more permissive sequence network ([Figure 6A](#pone-0074477-g006){ref-type="fig"}) the ARD structures share an edge with each other as well as other proteins but partition off into an independent doublet (sharing a single edge) under a more stringent threshold ([Figure 6B](#pone-0074477-g006){ref-type="fig"}).

That members of the cupin superfamily cluster by function when the function is defined at a fine level allows one to make functional inferences which may provide a starting point for biochemical investigations. For example, nearest neighbor analysis of the protein product of AT3G21360 of *Arabidopsis thaliana* (PBD code: 1y0z) [@pone.0074477-Bitto1] in the more stringent structure network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}, not shown with arrows) suggests that the protein product of AT3G21360 may be an alpha ketoglutarate-dependent TauD. This inference is further supported by the exclusive partitioning of this protein with PDB codes "3v15" and "10s7" (both confirmed TauDs) in the more stringent domain-based network ([Figure 4B](#pone-0074477-g004){ref-type="fig"}). Similarly, the uncharacterized *Escherichia coli* protein yeaR (PDB code: 3bb6) shares an edge exclusively with the *Vibrio fisheri* tellurite resistance protein B (PDB code: 3dl3) in both the more stringent structure network ([Figure 3B](#pone-0074477-g003){ref-type="fig"}) and the less stringent domain-based network ([Figure 4A](#pone-0074477-g004){ref-type="fig"}) (not shown with arrows in either network).

JmjC Proteins {#s2c}
-------------

The Jumonji C domain-containing proteins are a subfamily of the Fe(II)/2-oxoglutarate-dependent oxygenases. Sequence similarity between JmjC and cupin metalloenzyme domains allowed the prediction of active-site residues in JmjC domains and provided the insights that guided the experimental determination of the reactions catalyzed by these enzymes [@pone.0074477-Clissold1]. The JmjC domain-containing factor inhibiting hypoxia (FIH) ([Figures 1I](#pone-0074477-g001){ref-type="fig"} and [2H](#pone-0074477-g002){ref-type="fig"}) was identified as an asparaginyl hydroxylase that transcriptionally regulates the activity of hypoxia-inducible factor (HIF) [@pone.0074477-Lando1]. JmjC domain-containing proteins were proposed to function as histone demethylases in regulating chromatin structure [@pone.0074477-Trewick1] and many of the JmjC domain-containing enzymes have been shown to comprise the largest class of histone demethylases that catalyze lysine demethylation of mono-, di-, and trimethylated lysine residues through the formation of a hemiaminal intermediate that yields the demethylated product and formaldehyde [@pone.0074477-Klose1], [@pone.0074477-Loenarz1]. The lysine demethylases are represented by lysine-specific demethylase 6B in [Figures 1J](#pone-0074477-g001){ref-type="fig"}, [2I, and 2J](#pone-0074477-g002){ref-type="fig"}. Additionally, JmjC domain-containing enzymes have been identified that have RNA hydroxylase activity [@pone.0074477-Iyer1], [@pone.0074477-Noma1].

At the time of this writing, the cupin clan of the Pfam database contains structures of 13 unique JmjC domain-containing proteins ([Table 1](#pone-0074477-t001){ref-type="table"}). Most (12) of these are given the Pfam family designation 'JmjC'. These 12 structures form an exclusive cluster together in [Figure 3B](#pone-0074477-g003){ref-type="fig"}. The other JmjC domain-containing protein structure (PDB code: 3ld8) is designated 'cupin 8' in Pfam along with an additional three proteins that have structural information available ([Table 1](#pone-0074477-t001){ref-type="table"}). These four structures form an exclusive cluster in [Figure 3B](#pone-0074477-g003){ref-type="fig"} (not labeled under the AlkB-containing cluster). Finally, the three proteins with structures available designated 'cupin 4' in Pfam have been included in [Table 1](#pone-0074477-t001){ref-type="table"} for comparison. In [Figure 3B](#pone-0074477-g003){ref-type="fig"} these structures exist as a pair of nodes and as a single node. In the more stringent sequence similarity network ([Figure 6B](#pone-0074477-g006){ref-type="fig"}), six JmjC domain-containing protein structures and the 'cupin 8' structures cluster together (PDB codes: 2yu2, 3k3n, 3n9l, 3pu8, 3u78, 3ld8, 4 aap, 3a16, 2y0i), while five JmjC domain-containing protein structures form an exclusive cluster (PDB codes: 2gp3, 2w2i, 2 xml, 3 dxu, 3 opt). The three 'cupin 4' structures partition into an exclusive triad. These networks capture a current snapshot of relationships within this subfamily and can be used to update relationships and guide experimental design as new structures become available.

10.1371/journal.pone.0074477.t001

###### JmjC protein clusters.

![](pone.0074477.t001){#pone-0074477-t001-1}

  Protein                                                                   Species             PDB code   Pfam family   Metal ion(s)              Ref.
  --------------------------------------------------------------- ---------------------------- ---------- ------------- -------------- -----------------------------
  Lysine-specific demethylase 6B                                         *Homo sapiens*          2 xue        JmjC          Fe, Zn      [@pone.0074477-Kruidenier1]
  JmjC domain-containing histone demethylation protein 1^a^              *Homo sapiens*           2yu2        JmjC            Fe       
  Lysine-specific demethylase 6A                                         *Homo sapiens*          3 avr        JmjC          Ni, Zn       [@pone.0074477-Sengoku1]
  catalytic domain of PHD finger protein 8                               *Homo sapiens*           3k3n        JmjC            Fe            [@pone.0074477-Yu1]
  Lysine-specific demethylase 7                                     *Caenorhabditis elegans*      3n9l        JmjC          Fe, Zn         [@pone.0074477-Yang2]
  PHD finger protein 2                                                   *Homo sapiens*           3pu8        JmjC            Fe          [@pone.0074477-Horton1]
  Lysine-specific demethylase 7                                          *Homo sapiens*           3u78        JmjC            Ni         [@pone.0074477-Upadhyay1]
  Jumonji domain-containing protein 2A                                   *Homo sapiens*           2gp3        JmjC          Fe, Ni         [@pone.0074477-Chen1]
  Lysine-specific demethylase 4E                                         *Homo sapiens*           2w2i        JmjC            Ni       
  JMJD2C catalytic domain                                                *Homo sapiens*          2 xml        JmjC          Ni, Zn     
  JmjC domain-containing histone demethylation protein 3D                *Homo sapiens*          3 dxu        JmjC            Fe       
  DNA damage-responsive transcriptional repressor RPH1             *Saccharomyces cerevisiae*    3 opt        JmjC            Ni          [@pone.0074477-Chang1]
  Bifunctional arginine demethylase and lysyl-hydroxylase JMJD6          *Homo sapiens*           3ld8       cupin 8        Fe, Hg         [@pone.0074477-Hong1]
  Lysine-specific demethylase 8                                          *Homo sapiens*          4 aap       cupin 8          Ni       
  JmjC domain-containing protein C2orf60                                 *Homo sapiens*           3al6       cupin 8          Ni           [@pone.0074477-Kato1]
  factor inhibiting HIF1                                            *Caenorhabditis elegans*      2y0i       cupin 8          Fe           [@pone.0074477-Yang1]
  Putative asparaginyl hydroxylase                                       *Homo sapiens*          1 vrb       cupin 4          Fe       
  Lysine-specific demethylase NO66                                       *Homo sapiens*          4 diq       cupin 4          Ni       
  MYC-induced nuclear antigen                                            *Homo sapiens*          2 xdv       cupin 4      Cd, Ni, Mn   

RmlC Epimerases {#s2d}
---------------

The RmlC epimerases in [Figure 3A](#pone-0074477-g003){ref-type="fig"} cluster together but share edges with other protein structures. At a more stringent threshold ([Figure 3B](#pone-0074477-g003){ref-type="fig"}), however, these same 10 structures ([Table 2](#pone-0074477-t002){ref-type="table"}) cluster independently from other structures. The epimerases are monocupins. Members of this group do not bind a metal ion and are represented in [Figure 1D](#pone-0074477-g001){ref-type="fig"} by NovW, a 4-keto-6-deoxy sugar epimerase (PDB code: 2c0z) [@pone.0074477-Jakimowicz1]. This grouping is in agreement with a published structure-based phylogenetic analysis of the cupin superfamily generated using a structure dissimilarity matrix through pairwise structure-based alignment of 52 cupin proteins [@pone.0074477-Agarwal1]. Our analysis includes 2 additional proteins whose structures were solved after the phylogenetic analysis was published. The network in [Figure 3B](#pone-0074477-g003){ref-type="fig"} clearly groups the functionally similar RmlC epimerases together as does the phylogenetic analysis, providing further validation that PSNs recapitulate much of the information present in phylogenetic trees [@pone.0074477-Atkinson1]. When the network is constructed of isolated domains, the 10 monocupin RmlC epimerase domains form the same cluster as in the whole protein network ([Figure 4B](#pone-0074477-g004){ref-type="fig"}). Furthermore, a sequence-based network clusters 9 of the 10 epimerases together only when edges between nodes are drawn if the E-value is better than of 1E-6.0 ([Figure 6B](#pone-0074477-g006){ref-type="fig"}). In this case the enzyme from *Aneurinibacillus thermoaerophilus* is excluded from the cluster.

10.1371/journal.pone.0074477.t002

###### The RmlC epimerase cluster.

![](pone.0074477.t002){#pone-0074477-t002-2}

  Protein                                                         Species                  PDB code   \# of chains   Metal ion                             Ref.
  ----------------------------------------------- --------------------------------------- ---------- -------------- ----------- ----------------------------------------------------------
  dTDP-4-keto-6-deoxy-d-hexulose 3, 5-epimerase    *Methanobaterium thermoautotrophicum*     1ep0          1           none                    [@pone.0074477-Christendat1]
  4-keto-6-deoxy sugar epimerase (NovW)                  *Streptomyces spheroides*           2c0z          1           none      [@pone.0074477-Jakimowicz2], [@pone.0074477-Jakimowicz3]
  dTDP-6-deoxy-D-xylo-4-hexulose 3,5-epimerase           *Pseudomonas aeruginosa*           1 rtv          1           none                       [@pone.0074477-Dong1]
  dTDP-4-dehydrorhamnose 3,5 epimerase                 *Mycobacterium tuberculosis*         1 upi          1           none        [@pone.0074477-Dong1], [@pone.0074477-Kantardjieff1]
  dTDP-4-keto-6-deoxy-glucose-5 epimerase EvaD          *Amycolatopsis orientalis*          1 ofn          2           none                      [@pone.0074477-Merkel1]
  dTDP-4-dehydrorhamnose 3,5 epimerase                     *Sulfolobus tokodaii*            1 wlt          2           none     
  dTDP-4-dehydrorhamnose 3,5 epimerase                     *Streptococcus suis*             1 nxm          2           none                       [@pone.0074477-Dong2]
  dTDP-4-dehydrorhamnose 3,5 epimerase                   *Salmonella typhimurium*           1 dzr          2           none                      [@pone.0074477-Giraud1]
  dTDP-4-dehydrorhamnose 3,5 epimerase                     *Bacillus anthracis*             3 ryk          2           none     
  dTDP-6-deoxy-3,4-keto-hexulose isomerase          *Aneurinibacillus thermoaerophilus*      2pa7          2           none                       [@pone.0074477-Davis1]

Conclusions {#s2e}
-----------

Protein similarity networks (structure and sequence) of the cupin superfamily recapitulate and complement phylogenetic studies. Structure- and sequence-based clustering of cupin members reflects functional clustering. Networks based only on cupin domains and networks based on the whole proteins provide complementary information. Domain-clustering supports phylogenetic conclusions that the N- and C- terminal domains of bicupin proteins evolved independently. Interestingly, although many functionally similar enzymatic cupin members bind the same active site metal ion, the structure and sequence clustering does not correlate with the identity of the bound metal. PSNs are expected to be a valuable tool for directing experimental work and for predicting the functions of uncharacterized members of the cupin superfamily.

Methods {#s3}
=======

Dataset Curation {#s3a}
----------------

All structures associated with the PFAM cupin clan (CL0029) were downloaded 31 March 2013. There were 945 chains from 484 different structures (many structures had multiple chains). However, many structures were that of the same protein, a problem that could cause an undue distortion of network topology. Duplicates were manually removed to avoid this issue, reducing the set of proteins to 183 unique structures`. `The structure used to represent duplicate structures was selected arbitrarily. Initially the structures of individual chains of whole protein structures were compared to each other. This approach, however, also resulted in the unequal representation of the 183 structures in the network. Therefore we treated the overall quaternary organization of proteins with multiple chains as a single structure. Biological, ligand, and domain information were obtained using the RESTful web service interface provided by the RCSB. Only the biologically-significant transition metals were considered when painting the networks by bound metal. The Taxonomy Database from the NCBI was used to classify species into their respective domains and phyla. For the sequence similarity network, UniProt was used to map the PDB IDs used in the structure similarity network to their respective sequences in the UniProtKB database.

Building the PSNs {#s3b}
-----------------

Pythoscape was integral in the construction of the network. Pythoscape (1) imported sequences and structures into a database, (2) deployed TM-align and BLASTp for edge calculations, and (3) exported the completed database as a Cytoscape network file. Because TM-Align edge scores are directional (i.e., the score for the edge from A to B might be different than the score from B to A), edges were filtered based on the average of the two scores. Cytoscape 3.0 was used for network visualization in the as well as edge filtering. The Organic Layout was used as described in Atkinson, Babbitt et al. To construct the structure network composed only of domains, a script was written for PyMOL to extract and save the parts of the PDB files that had a domain as defined as by the RCSB. Domains outside of the PFAM cupin clan were manually removed from the network. Networks were visualized using two thresholds in order to illustrate effects edge stringency had on certain clusters. TM-scores of 0.53 and 0.65 were used as thresholds for the structure networks, and E-values of 1×10^−3.5^ and 1×10^−6.0^ were used as thresholds for the sequence networks. These values were selected based on overall visual appeal of the resulting layout in Cytoscape.

Supporting Information {#s4}
======================

###### 

**Contains:** Table S1. Iron-containing structures used in the networks. Table S2. Nickel-containing structures used in the networks.Table S3. Zinc-containing structures used in the networks. Table S4. Manganese-containing structures used in the networks. Table S5. Structures used in the networks which contain copper, mercury, and cadmium, respectively. Table S6. Structures used in the networks that contain multiple metals. Table S7. Structures used in the networks that contain no transition metal.

(DOCX)

###### 

Click here for additional data file.

The authors wish to thank Dr. Patrick Frantom for introducing us to the usefulness of protein similarity networks and for his generous help in getting us started on this project.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: RU EWM. Performed the experiments: RU EWM. Analyzed the data: RU EWM. Contributed reagents/materials/analysis tools: EWM. Wrote the paper: RU EWM.
