Rhodanese domains are ubiquitous structural modules occurring in the three major evolutionary phyla. They are found as tandem repeats, with the C‐terminal domain hosting the properly structured active‐site Cys residue, as single domain proteins or in combination with distinct protein domains. An increasing number of reports indicate that rhodanese modules are versatile sulfur carriers that have adapted their function to fulfill the need for reactive sulfane sulfur in distinct metabolic and regulatory pathways. Recent investigations have shown that rhodanese domains are also structurally related to the catalytic subunit of Cdc25 phosphatase enzymes and that the two enzyme families are likely to share a common evolutionary origin. In this review, the rhodanese/Cdc25 phosphatase superfamily is analyzed. Although the identification of their biological substrates has thus far proven elusive, the emerging picture points to a role for the amino‐acid composition of the active‐site loop in substrate recognition/specificity. Furthermore, the frequently observed association of catalytically inactive rhodanese modules with other protein domains suggests a distinct regulatory role for these inactive domains, possibly in connection with signaling.
Rhodaneses (thiosulfate:cyanide sulfurtransferases or TSTs; E.C. 126.96.36.199) are widespread enzymes that catalyze the transfer of a sulfane sulfur atom from thiosulfate to cyanide in vitro (Figure 1). 3‐mercaptopyruvate:cyanide sulfurtransferases (MSTs; E.C. 188.8.131.52) catalyze a similar reaction, using 3‐mercaptopyruvate as the sulfur donor. MSTs display clear sequence homology with rhodaneses, with which they share up to 66% identical residues (Nagahara et al., 1995). The most well characterized sufurtransferase is bovine liver rhodanese, which has been the subject of numerous functional investigations (e.g. Westley et al., 1983; Berni et al., 1991; Luo and Horowitz, 1994; Nandi et al., 2000). The catalysis occurs via a double displacement mechanism involving the transient formation of a persulfide‐containing intermediate (Rhod‐S), in which the transferring sulfur is bound to the invariant catalytic Cys residue (Figure 1).
The crystal structures of bovine (Rhobov) and Azotobacter vinelandii (RhdA) rhodaneses are known (Ploegman et al., 1978; Bordo et al., 2000). These two proteins, which exhibit 22% sequence identity, display very similar three‐dimensional conformations. They are each composed of two identically folded domains (Figure 2), called hereafter rhodanese domains or modules, each about 120 amino acids long, which display weak sequence similarity to one another (13 and 21% identical residues, respectively). Due to this peculiar tandem domain architecture, rhodanese has been considered the prototype of divergent evolution from a common ancestor protein which, after gene duplication and under the constraint of tertiary structure conservation, led to the almost complete obliteration of sequence similarity between the two domains (Ploegman et al., 1978). Only the C‐terminal rhodanese domain (referred to as the catalytic domain) hosts the catalytic Cys residue, which is the first residue of a six‐amino‐acid active‐site loop that folds in a cradle‐like structure defining the enzyme catalytic pocket; in the N‐terminal domain the same position is occupied by Asp (Figure 2). The recent biochemical and structural characterization of GlpE, a TST from Escherichia coli composed of a single rhodanese domain (Ray et al., 2000; Spallarossa et al., 2001), indicates that the catalytically inactive N‐terminal domain found in TSTs and MSTs is not essential for catalysis. Furthermore, the structure of the catalytic domain of two human Cdc25 phosphatases, Cdc25A and Cdc25B, two key enzymes involved in cell cycle control, has revealed that these proteins display a rhodanese‐like three‐dimensional fold and maintain the same location of the active‐site Cys residue (Figure 2; Fauman et al., 1998; Reynolds et al., 1999). In spite of the very weak sequence similarity (17% identical residues), sensitive profile analysis indicated that the two enzyme families stem from a common ancestor (Hofmann et al., 1998). The most relevant structural difference between rhodanese and Cdc25 enzymes is the length of the active‐site loop, which in Cdc25 proteins is formed by seven residues instead of the six in sulfurtransferases; this results in a wider catalytic pocket that can accommodate a phosphorous atom, which has a van der Waals radius slightly larger than a sulfur atom (1.9 versus 1.85 Å).
The biological role of sulfurtransferases is still largely debated, since the identification of their in vivo substrates has thus far proven unsuccessful, and it appears unlikely that these substrates are those identified by in vitro reactions (Nandi et al., 2000; Ray et al., 2000). Proposed functions include cyanide detoxification (Sorbo, 1957), formation of prosthetic groups in iron–sulfur proteins (Pagani et al., 1984), maintenance of the sulfane pool (Westley, 1989), selenium metabolism (Ogasawara et al., 2001) and thiamin biosynthesis (Palenchar et al., 2000). The involvement of rhodanese enzymes in distinct biological functions is also suggested by the wide variability in the amino acids of the active‐site loop. In fact, the amino‐acid side chains extending from this structural element define the ridge of the catalytic pocket and are expected, therefore, to play a key role in substrate recognition and catalytic activity.
In this review, the information on rhodanese‐like homology domains contained in the protein database is analyzed together with the available structural and biological information to provide a general view of the distinct structural/functional contexts in which these structural domains occur in living organisms.
Rhodanese homology domains encoded in sequenced genomes
The complete sequencing of more than 80 genomes has shown that proteins bearing weak but significant sequence similarity to known rhodaneses or Cdc25 phosphatase enzymes are widely distributed among living organisms. Sensitive homology search programs based on profile analysis or Hidden Markov Models have detected rhodanese homology domains in about 500 genes in the three major evolutionary phyla (see, for example, SMART, http://smart.embl‐heidelberg.de; DART, http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml; Hofmann et al., 1998). Often, distinct proteins containing rhodanese domains are encoded by the same genome. For example, in the human genome, 47 such instances are observed, whereas in the E. coli K12 genome, seven rhodanese‐like proteins are found (see, for example, the COG database; Tatusov et al., 2000). From a structural viewpoint, rhodanese‐like proteins are either composed of two rhodanese domains, with the C‐terminal domain displaying the putative catalytic Cys as observed in Rhobov and RhdA, or composed of a single catalytic rhodanese domain, as found in GlpE. Rhodanese domains, either catalytic or inactive (i.e. where the active‐site Cys is replaced by another residue), are also found associated with other protein domains such as MAPK‐phosphatases (Keyse and Ginsburg, 1993; Fauman et al., 1998; Hofmann et al., 1998) or ThiI, an E. coli enzyme involved in thiamin and thiouridine biosynthesis (Palenchar et al., 2000).
A multiple alignment of the rhodanese homology domain, containing a representative subset of the entries detected in the genomic database, is provided by the SMART resource (Schultz et al., 1998). The deduced neighbor‐joining (N‐J) tree provides a graphical representation of the mutual similarity among distinct amino‐acid sequences (Figure 3). As expected from the wide amino‐acid sequence divergence found within this family, bootstrap analysis indicates that the detailed tree structure might be unreliable, especially in its deep branching. Nevertheless, the tree is useful for identifying protein subfamilies that may relate to distinct biological functions, as discussed below.
Catalytically active rhodanese and Cdc25 phosphatase domains
The N‐ and C‐terminal domains of TST enzymes form two distinct subfamilies (Figure 3). Each subfamily also contains previously characterized MST sequences from Rattus norvegicus, Arabidopsis thaliana and E. coli (Nagahara et al., 1995; Nagahara and Nishino, 1996; Papenbrock and Schmidt, 2000; Colnaghi et al., 2001). In TSTs and MSTs, the putative catalytic Cys is invariably found in the C‐terminal domain, whereas in the N‐terminal domain, this residue is often replaced by Asp or Gly. In three‐dimensional structures of RhdA and Rhobov, the side chain of this Asp residue is involved in multiple hydrogen bonds with the main chain amide groups of the active‐site loop (Ploegman et al., 1978; Bordo et al., 2000), which makes this side chain unavailable for substrate or ion binding. From an evolutionary viewpoint, as the clusters of both N‐ and C‐terminal domains of TST/MST subfamilies contain eukaryotic, bacterial and archaeal sequences (Figure 3), the gene duplication event that gave rise to the tandem domain architecture found in MSTs and TSTs must have occurred very early during evolution, probably before the diversification of living organisms into the three major evolutionary phyla. Notably, active‐site loop amino‐acid sequences are distinct in TSTs and in MSTs, which display motifs CRXGX[R/T] and CG[S/T]GVT, respectively (square brackets indicate alternative residues, ‘X’ any amino acid; Figure 4). The active‐site loop of known TSTs always contains two basic residues (i.e. if the last loop residue is Thr, a basic residue is observed at one of the X positions), whereas no charged residues are observed in biochemically characterized MSTs; this may relate to the distinct ionic charge of the respective in vitro substrates, thiosulfate (2−) and 3‐mercaptopyruvate (1−).
Cdc25 phosphatases are composed of variable N‐terminal domains, which are poorly characterized but are suspected to serve regulatory functions (Draetta and Eckstein, 1997; Forrest and Gabrielli, 2001), and conserved catalytic, rhodanese‐like, C‐terminal domains (Figure 2; Fauman et al., 1998; Reynolds et al., 1999). In the N‐J tree, these catalytic domains form a distinct clade (Figure 3). Similarly to MSTs/TSTs, the active‐site amino‐acid sequences of Cdc25 enzymes display a conserved motif (CE[Y/F]SXXR; Figure 4), which is consistent with the role of this structural element in substrate recognition.
An elongated seven‐amino‐acid active‐site loop, with the putative catalytic Cys at the first position, is also found in the yeast arsenic‐resistance protein Acr2 and in the closely related yeast protein Yg4E (Figure 3; Bobrowicz et al., 1997; Mukhopadhyay et al., 2000). Although both proteins lack the regulatory N‐terminal domain observed in Cdc25 phosphatase enzymes, the N‐J tree suggests a close evolutionary relationship with this family of enzymes. A seven‐amino‐acid active‐site loop bearing the Cys at the first position is also found in a small rhodanese subfamily, formed by eukaryotic proteins from Mycobacterium tuberculosis and Streptomyces peucetis. These rhodanese‐like domains are found in association with the helix–turn–helix motif of the arsenic resistance operon repressor ArsR (Figures 3 and 4). This structural association supports the hypothesis that these rhodanese modules are able to bind arsenic, and it is tempting to speculate that they may act as sensory domains modulating ArsR activity.
Taken together, the available biochemical data support the hypothesis that rhodanese domains displaying a seven‐amino‐acid active‐site loop are able to bind substrates containing phosphorous or the chemically similar arsenic, whereas rhodanese‐like domains displaying a six‐amino‐acid loop with Cys at the first position interact with substrates containing reactive sulfur or, in some cases, selenium (Ogasawara et al., 2001). Considering the very distant relationship between bacterial ArsR proteins and eukaryotic Cdc25/Acr2 proteins, it is likely that these two subfamilies result from two independent mutational events, each resulting from a single amino‐acid insertion in the active‐site loop of sulfur‐binding rhodanese modules.
Proteins composed of a single rhodanese‐like module, displaying the putative catalytic Cys and a six‐amino‐acid active‐site loop, are widespread, especially among bacteria (Figure 3). These domains display largely divergent amino‐acid sequences, being a cluster of uncharacterized proteins similar to E. coli YceA and are the only example of bacterial and eukaryotic proteins forming a well‐defined subfamily of closely related sequences (Figure 4). Besides GlpE from E. coli, only the periplasmic sulfide dehydrogenase (Sud) from Wolinella succinogenes has been described as a (polysulfide:cyanide) sulfur‐ transferase (Krafft et al., 1995; Klimmek et al., 1998). At the genetic level, single rhodanese domain proteins have been associated with specific stress conditions. Examples include the Drosophila melanogaster heat shock protein 67B2, the Vibrio cholerae shock protein Q9KN65 or the E. coli phage shock protein GspE. Single rhodanese domains are also associated with the process of leaf senescence in A. thaliana, Nicotiana tabacum and Raphanus sativus (Sen1, Ntdin and Din1, respectively; Azumi and Watanabe, 1991; Oh et al., 1996).
Putatively catalytic rhodanese domains also frequently associate with proteins involved in thiamin, thiouridine or molybdopterin biosynthesis (Figure 3). In particular, a rhodanese domain is found at the C‐terminal end of several ThiF/MoeB proteins (see, for example, the DART resource). In these cases the active‐site loop consensus sequence is C[R/K]XGXR and C[R/K]XGXD in bacteria and eukaryota, respectively (Figure 4). ThiF, which is involved in the synthesis of thiamin in bacteria, is deemed responsible for the adenylation of ThiS. Therefore, the C‐terminal carboxy‐adenylate ThiS might interact with a sulfur donor to proceed in the biosynthesis of thiazole, which is an intermediate compound in thiamin biosynthesis (Taylor et al., 1998; Xi et al., 2001). In this context, the role of the rhodanese module might be to provide the sulfane sulfur required by ThiS.
In bacteria and archaea, a catalytic rhodanese module is also recurrently found in association with ThiI, another enzyme participating in thiouridine and thiamin biosynthesis (Figure 3). The sulfurtransferase activity of this rhodanese domain has been established in vitro using ThiI from E. coli (Palenchar et al., 2000). Also, in this case, the significance of sulfur‐delivering proteins in the synthesis of thiouridine and thiamin has been recognized (Begley et al., 1999). In particular, the sulfur transfer between IscS and ThiI has recently been shown and the essential role of the rhodanese Cys residue established (Kambampati and Lauhon, 1999; Mueller et al., 2001).
Catalytically inactive rhodanese domains
Rhodanese domains lacking the active‐site Cys, which is often replaced by acidic or non‐polar residues, are found as N‐terminal domains in about one‐third of the known dual‐specific phosphatases (DSP) and in several ubiquitinating enzymes (UBP), including yeast UBP4, UBP5, UBP7, human UBP8 and the murine deubiquitinating enzyme UbpY (Fauman et al., 1998; Hofmann et al., 1998). In the N‐J tree, these inactive rhodanese domains form a distinct group of sequences that appear to form a clade (Figure 3). The role of the inactive rhodanese domains, either those found in TST and MST enzymes or those in association with DSP or UBP, is still to be elucidated. Although in bovine rhodanese the N‐terminal domain is deemed to be involved in substrate recognition mediated by a loop folding over the C‐terminal domain in the vicinity of the active site (Ploegman et al., 1979), the absence of a structurally equivalent loop in A. vinelandii RhdA indicates that substrate recognition in sulfurtransferases, if mediated by the N‐terminal domain at all, must rely upon more complex mechanisms involving non‐contiguous surface patches. In DSP and UBP, the N‐terminal inactive rhodanese domain might play a regulatory role, perhaps in connection with signaling. In this regard, the multiple sequence alignment of the rhodanese‐like domains associated to DSP and UBP (see, for example, SMART) highlights the conservation of two residues, Tyr37 and His41 (RhdA numbering; Bordo et al., 2000), which occupy neighboring positions at the protein surface (Figure 2). This points to the possibility that either of these two amino acids could be a target for tyrosine or histidine kinases. Obviously, further experimental evidence is needed to prove this hypothesis.
The complete sequencing of a variety of genomes has shown that distinct proteins containing rhodanese‐like modules are often encoded by the same genome and now enables the identification of distinct subfamilies of rhodanese‐like proteins, possibly associated with different biological functions (paralogs). Each subfamily, displayed in the N‐J family tree as a cluster of closely related sequences, harbors conserved motifs in the amino‐acid sequence of the (putative) active‐site loop. However, despite the increasing amount of experimental information, there remains the open question of the identity of the true biological substrates, either proteins or small molecules, that interact with rhodanese domains. The information now available in the genomic databases should be very helpful in designing appropriate experiments to answer these questions.
We wish to thank Toby Gibson for critical reading of the manuscript and Jose Castresana for many helpful suggestions. This research was supported in part by the MIUR grant ‘Sulfotransferasi procariotiche’ and CNR strategic project ‘Molecular genetics’.
- Copyright © 2002 European Molecular Biology Organization