Effect of exchange of the cysteine molybdenum ligand with selenocysteine on the structure and function of the active site in human sulfite oxidase.

Sulfite oxidase (SO) is an essential molybdoenzyme for humans, catalyzing the final step in the degradation of sulfur-containing amino acids and lipids, which is the oxidation of sulfite to sulfate. The catalytic site of SO consists of a molybdenum ion bound to the dithiolene sulfurs of one molybdopterin (MPT) molecule, carrying two oxygen ligands, and is further coordinated by the thiol sulfur of a conserved cysteine residue. We have exchanged four non-active site cysteines in the molybdenum cofactor (Moco) binding domain of human SO (SOMD) with serine using site-directed mutagenesis. This facilitated the specific replacement of the active site Cys207 with selenocysteine during protein expression in Escherichia coli. The sulfite oxidizing activity (kcat/KM) of SeSOMD4Ser was increased at least 1.5-fold, and the pH optimum was shifted to a more acidic value compared to those of SOMD4Ser and SOMD4Cys(wt). X-ray absorption spectroscopy revealed a Mo(VI)-Se bond length of 2.51 Å, likely caused by the specific binding of Sec207 to the molybdenum, and otherwise rather similar square-pyramidal S/Se(Cys)O2Mo(VI)S2(MPT) site structures in the three constructs. The low-pH form of the Mo(V) electron paramagnetic resonance (EPR) signal of SeSOMD4Ser was altered compared to those of SOMD4Ser and SOMD4Cys(wt), with g1 in particular shifted to a lower magnetic field, due to the Se ligation at the molybdenum. In contrast, the Mo(V) EPR signal of the high-pH form was unchanged. The substantially stronger effect of substituting selenocysteine for cysteine at low pH as compared to high pH is most likely due to the decreased covalency of the Mo-Se bond.

Structure-based alteration of substrate specificity and catalytic activity of sulfite oxidase from sulfite oxidation to nitrate reduction.

Eukaryotic sulfite oxidase is a dimeric protein that contains the molybdenum cofactor and catalyzes the metabolically essential conversion of sulfite to sulfate as the terminal step in the metabolism of cysteine and methionine. Nitrate reductase is an evolutionarily related molybdoprotein in lower organisms that is essential for growth on nitrate. In this study, we describe human and chicken sulfite oxidase variants in which the active site has been modified to alter substrate specificity and activity from sulfite oxidation to nitrate reduction. On the basis of sequence alignments and the known crystal structure of chicken sulfite oxidase, two residues are conserved in nitrate reductases that align with residues in the active site of sulfite oxidase. On the basis of the crystal structure of yeast nitrate reductase, both positions were mutated in human sulfite oxidase and chicken sulfite oxidase. The resulting double-mutant variants demonstrated a marked decrease in sulfite oxidase activity but gained nitrate reductase activity. An additional methionine residue in the active site was proposed to be important in nitrate catalysis, and therefore, the triple variant was also produced. The nitrate reducing ability of the human sulfite oxidase triple mutant was nearly 3-fold greater than that of the double mutant. To obtain detailed structural data for the active site of these variants, we introduced the analogous mutations into chicken sulfite oxidase to perform crystallographic analysis. The crystal structures of the Mo domains of the double and triple mutants were determined to 2.4 and 2.1 Å resolution, respectively.

Nature of halide binding to the molybdenum site of sulfite oxidase.

Valuable information on the active sites of molybdenum enzymes has been provided from both Mo(V) electron paramagnetic resonance (EPR) spectroscopy and X-ray absorption spectroscopy (XAS). One of three major categories of Mo(V) EPR signals from the molybdenum enzyme sulfite oxidase is the low-pH signal, which forms in the presence of chloride. Two alternative structures for this species have been proposed, one in which the chloride is coordinated directly to Mo and a second in which chloride is held in the arginine-rich basic pocket some 5 Šfrom Mo. Here we present an independent assessment of the structure of this species by using XAS of the analogous bromide and iodide complexes. We show that there is no evidence of direct Mo-I coordination, and that the data are consistent with a structure in which the halide is bound at ∼5 Šfrom Mo.

The History of the Discovery of the Molybdenum Cofactor and Novel Aspects of its Biosynthesis in Bacteria.

Biosynthesis of the molybdenum cofactor in bacteria is described with a detailed analysis of each individual reaction leading to the formation of stable intermediates during the synthesis of molybdopterin from GTP. As a starting point, the discovery of molybdopterin and the elucidation of its structure through the study of stable degradation products are described. Subsequent to molybdopterin synthesis, the molybdenum atom is added to the molybdopterin dithiolene group to form the molybdenum cofactor. This cofactor is either inserted directly into specific molybdoenzymes or is further modified by the addition of nucleotides to the molybdopterin phosphate group or the replacement of ligands at the molybdenum center.

The structures of the C185S and C185A mutants of sulfite oxidase reveal rearrangement of the active site.

Sulfite oxidase (SO) catalyzes the physiologically critical conversion of sulfite to sulfate. Enzymatic activity is dependent on the presence of the metal molybdenum complexed with a pyranopterin-dithiolene cofactor termed molybdopterin. Comparison of the amino acid sequences of SOs from a variety of sources has identified a single conserved Cys residue essential for catalytic activity. The crystal structure of chicken liver sulfite oxidase indicated that this residue, Cys185 in chicken SO, coordinates the Mo atom in the active site. To improve our understanding of the role of this residue in the catalytic mechanism of sulfite oxidase, serine and alanine variants at position 185 of recombinant chicken SO were generated. Spectroscopic and kinetic studies indicate that neither variant is capable of sulfite oxidation. The crystal structure of the C185S variant was determined to 1.9 A resolution and to 2.4 A resolution in the presence of sulfite, and the C185A variant to 2.8 A resolution. The structures of the C185S and C185A variants revealed that neither the Ser or Ala side chains appeared to closely interact with the Mo atom and that a third oxo group replaced the usual cysteine sulfur ligand at the Mo center, confirming earlier extended X-ray absorption fine structure spectroscopy (EXAFS) work on the human C207S mutant. An unexpected result was that in the C185S variant, in the absence of sulfite, the active site residue Tyr322 became disordered as did the loop region flanking it. In the C185S variant crystallized in the presence of sulfite, the Tyr322 residue relocalized to the active site. The C185A variant structure also indicated the presence of a third oxygen ligand; however, Tyr322 remained in the active site. EXAFS studies of the Mo coordination environment indicate the Mo atom is in the oxidized Mo(VI) state in both the C185S and C185A variants of chicken SO and show the expected trioxodithiolene active site. Density functional theory calculations of the trioxo form of the cofactor reasonably reproducd the Mo horizontal lineO distances of the complex; however, the calculated Mo-S distances were slightly longer than either crystallographic or EXAFS measurements. Taken together, these results indicate that the active sites of the C185S and C185A variants are essentially catalytically inactive, the crystal structures of C185S and C185A variants contain a fully oxidized, trioxo form of the cofactor, and Tyr322 can undergo a conformational change that is relevant to the reaction mechanism. Additional DFT calculations demonstrated that such methods can reasonably reproduce the geometry and bond lengths of the active site.

Structural studies of the molybdenum center of the pathogenic R160Q mutant of human sulfite oxidase by pulsed EPR spectroscopy and 17O and 33S labeling.

Electron paramagnetic resonance (EPR) investigation of the Mo(V) center of the pathogenic R160Q mutant of human sulfite oxidase (hSO) confirms the presence of three distinct species whose relative abundances depend upon pH. Species 1 is exclusively present at pH < or = 6, and remains in significant amounts even at pH 8. Variable-frequency electron spin echo envelope modulation (ESEEM) studies of this species prepared with (33)S-labeled sulfite clearly show the presence of coordinated sulfate, as has previously been found for the "blocked" form of Arabidopsis thaliana at low pH (Astashkin, A. V.; Johnson-Winters, K.; Klein, E. L.; Byrne, R. S.; Hille, R.; Raitsimring, A. M.; Enemark, J. H. J. Am. Chem. Soc. 2007, 129, 14800). The ESEEM spectra of Species 1 prepared in (17)O-enriched water show both strongly and weakly magnetically coupled (17)O atoms that can be assigned to an equatorial sulfate ligand and the axial oxo ligand, respectively. The nuclear quadrupole interaction (nqi) of the axial oxo ligand is substantially stronger than those found for other oxo-Mo(V) centers studied previously. Additionally, pulsed electron-nuclear double resonance (ENDOR) measurements reveal a nearby weakly coupled exchangeable proton. The structure for Species 1 proposed from the pulsed EPR results using isotopic labeling is a six-coordinate Mo(V) center with an equatorial sulfate ligand that is hydrogen bonded to an exchangeable proton. Six-coordination is supported by the (17)O nqi parameters for the axial oxo group of the model compound, (dttd)Mo(17)O((17)Otms), where H2dttd = 2,3:8,9-dibenzo-1,4,7,10-tetrathiadecane; tms = trimethylsilyl. Reduction of R160Q to Mo(V) with Ti(III) gives primarily Species 2, another low pH form, whereas reduction with sulfite at higher pH values gives a mixture of Species 1 and 2, as well as the "primary" high pH form of wild-type SO. The occurrence of significant amounts of the "sulfate-blocked" form of R160Q (Species 1) at physiological pH suggests that this species may be a contributing factor to the lethality of this mutation.

Studies of the Mo(V) Center of the Y343F Mutant of Human Sulfite Oxidase by Variable Frequency Pulsed EPR Spectroscopy.

The Mo(V) forms of the Tyr343Phe (Y343F) mutant of human sulfite oxidase (SO) have been investigated by continuous wave (CW) and variable frequency pulsed EPR spectroscopies as a function of pH. The CW EPR spectrum recorded at low pH (∼6.9) has g-values similar to those known for the low-pH form of the native vertebrate SO (original lpH form); however, unlike the spectrum of original lpH SO, it does not show any hyperfine splittings from a nearby exchangeable proton. The detailed electron spin echo (ESE) envelope modulation (ESEEM) and pulsed electron-nuclear double resonance (ENDOR) experiments also did not reveal any nearby protons that could belong to an exchangeable ligand at the molybdenum center. These results suggest that under low-pH conditions the active site of Y343F SO is in the "blocked" form, with the Mo(V) center coordinated by sulfate. With increasing pH the EPR signal from the "blocked" form decreases, while a signal similar to that of the original lpH form appears and becomes the dominant signal at pH>9. In addition, both the CW EPR and ESE-detected field sweep spectra reveal a considerable contribution from a signal similar to that usually detected for the high-pH form of native vertebrate SO (original hpH form). The nearby exchangeable protons in both of the component forms observed at high pH were studied by the ESEEM spectroscopy. These results indicate that the Y343F mutation increases the apparent pK(a) of the transition from the lpH to hpH forms by ∼2 pH units.

MoV electron paramagnetic resonance of sulfite oxidase revisited: the low-pH chloride signal.

Valuable information on the active sites of molybdenum enzymes has been provided by Mo(V) electron paramagnetic resonance (EPR) spectroscopy. In recent years, multiple resonance techniques have been extensively used to examine details of the active-site structure, but basic continuous-wave (CW) EPR has not been re-evaluated in several decades. Here, we present a re-examination of the CW EPR spectroscopy of the sulfite oxidase low-pH chloride species and provide evidence for direct coordination of molybdenum by chloride.

Structure of the molybdenum site of Escherichia coli trimethylamine N-oxide reductase.

We report a structural characterization of the molybdenum site of recombinant Escherichia coli trimethylamine N-oxide (TMAO) reductase using X-ray absorption spectroscopy. The enzyme active site shows considerable similarity to that of dimethyl sulfoxide (DMSO) reductase, in that, like DMSO reductase, the TMAO reductase active site can exist in multiple forms. Examination of the published crystal structure of TMAO oxidase from Shewanella massilia indicates that the postulated Mo coordination structure is chemically impossible. The presence of multiple active site structures provides a potential explanation for the anomalous features reported from the crystal structure.

Electrocatalytically functional multilayer assembly of sulfite oxidase and cytochrome c.

An electrocatalytically functional multilayer has been designed using two proteins, cytochrome c and sulfite oxidase, and a polyelectrolyte (polyaniline sulfonate). The two proteins were co-immobilized on the surface of a gold electrode in alternating layers by electrostatic interactions using the layer-by-layer technique. The formation of this fully electro-active multilayer is characterized by quartz crystal microbalance and electrochemical experiments. The electro-catalytic characterization of the device containing up to 12 layers is based on generation of an oxidation current after sulfite addition. Besides the electron-transfer mechanism, the role of the different components in the electron-transport chain is clarified. Kinetic data were extracted to characterize the multilayer function. This artificial multilayer assembly is expected to be useful in the biosensor and biofuel cell development.

Crystal structure of a molybdopterin synthase-precursor Z complex: insight into its sulfur transfer mechanism and its role in molybdenum cofactor deficiency.

In almost all biological life forms, molybdenum and tungsten are coordinated by molybdopterin (MPT), a tricyclic pyranopterin containing a cis-dithiolene group. Together, the metal and the pterin moiety form the redox reactive molybdenum cofactor (Moco). Mutations in patients with deficiencies in Moco biosynthesis usually occur in the enzymes catalyzing the first and second steps of biosynthesis, leading to the formation of precursor Z and MPT, respectively. The second step is catalyzed by the heterotetrameric MPT synthase protein consisting of two large (MoaE) and two small (MoaD) subunits with the MoaD subunits located at opposite ends of a central MoaE dimer. Previous studies have determined that the conversion of the sulfur- and metal-free precursor Z to MPT by MPT synthase involves the transfer of sulfur atoms from a C-terminal MoaD thiocarboxylate to the C-1' and C-2' positions of precursor Z. Here, we present the crystal structures of non-thiocarboxylated MPT synthase from Staphylococcus aureus in its apo form and in complex with precursor Z. A comparison of the two structures reveals conformational changes in a loop that participates in interactions with precursor Z. In the complex, precursor Z is bound by strictly conserved residues in a pocket at the MoaE dimer interface in close proximity of the C-terminal glycine of MoaD. Biochemical evidence indicates that the first dithiolene sulfur is added at the C-2' position.

Modified active site coordination in a clinical mutant of sulfite oxidase.

The molybdenum site of the Arginine 160 --> Glutamine clinical mutant of the physiologically vital enzyme sulfite oxidase has been investigated by a combination of X-ray absorption spectroscopy and density functional theory calculations. We conclude that the mutant enzyme has a six-coordinate pseudo-octahedral active site with coordination of Glutamine Oepsilon to molybdenum. This contrasts with the wild-type enzyme which is five-coordinate with approximately square-based pyramidal geometry. This difference in the structure of the molybdenum site explains many of the properties of the mutant enzyme which have previously been reported.

Interaction of product analogues with the active site of rhodobacter sphaeroides dimethyl sulfoxide reductase.

We report a structural characterization using X-ray absorption spectroscopy of Rhodobacter sphaeroides dimethyl sulfoxide (DMSO) reductase reduced with trimethylarsine and show that this is structurally analogous to the physiologically relevant dimethyl sulfide reduced DMSO reductase. Our data unambiguously indicate that these species should be regarded as formal MoIV species and indicate a classical coordination complex of trimethylarsine oxide, with no special structural distortions. The similarity of the trimethylarsine and dimethyl sulfide complexes suggests, in turn, that the dimethyl sulfide reduced enzyme possesses a classical coordination of DMSO with no special elongation of the S-O bond, as previously suggested.

Mutational analysis of Escherichia coli MoeA: two functional activities map to the active site cleft.

The molybdenum cofactor is ubiquitous in nature, and the pathway for Moco biosynthesis is conserved in all three domains of life. Recent work has helped to illuminate one of the most enigmatic steps in Moco biosynthesis, ligation of metal to molybdopterin (the organic component of the cofactor) to form the active cofactor. In Escherichia coli, the MoeA protein mediates ligation of Mo to molybdopterin while the MogA protein enhances this process in an ATP-dependent manner. The X-ray crystal structures for both proteins have been previously described as well as two essential MogA residues, Asp49 and Asp82. Here we describe a detailed mutational analysis of the MoeA protein. Variants of conserved residues at the putative active site of MoeA were analyzed for a loss of function in two different, previously described assays, one employing moeA- crude extracts and the other utilizing a defined system. Oddly, no correlation was observed between the activity in the two assays. In fact, our results showed a general trend toward an inverse relationship between the activity in each assay. Moco binding studies indicated a strong correlation between a variant's ability to bind Moco and its activity in the purified component assay. Crystal structures of the functionally characterized MoeA variants revealed no major structural changes, indicating that the functional differences observed are not due to disruption of the protein structure. On the basis of these results, two different functional areas were assigned to regions at or near the MoeA active site cleft.

Role of the C-terminal Gly-Gly motif of Escherichia coli MoaD, a molybdenum cofactor biosynthesis protein with a ubiquitin fold.

In Escherichia coli, the MoaD protein plays a central role in the conversion of precursor Z to molybdopterin (MPT) during molybdenum cofactor biosynthesis. MoaD has a fold similar to that of ubiquitin and contains a highly conserved C-terminal Gly-Gly motif, which in its active form contains a transferrable sulfur in the form of a thiocarboxylate group. During MPT biosynthesis, MoaD cycles between two different heterotetrameric complexes, one with MoaE to form MPT synthase and the other with MoeB, a protein similar to E1 in the ubiquitin pathway, to regenerate its transferrable sulfur. To determine the specific roles of each of the two terminal Gly residues with regard to the MoaD cycle, variants at the penultimate (Gly80) or terminal (Gly81) residues of both MoaD and thiocarboxylated MoaD were created. These variants were analyzed to determine their effects on complex formation with MoaE and MoeB, formation of the MoaD-acyl-adenylate complex, transfer of sulfur to precursor Z to form MPT, and total cofactor biosynthesis. The combined results show that while conservative substitutions at Gly80 had little effect on any of the processes that were examined, the terminal MoaD residue (Gly81) is important for transfer of sulfur to precursor Z and essential for formation of the MoaD-AMP complex. These results further our understanding of the mechanistic similarities of the MoaD-MoeB reaction to that of the ubiquitin-E1 system.

The G473D mutation impairs dimerization and catalysis in human sulfite oxidase.

Among the mutations identified in patients with isolated sulfite oxidase deficiency, the G473D variant is of particular interest since sedimentation analysis reveals that this variant is a monomer, and the importance of the wild-type dimeric state of mammalian sulfite oxidase is not yet well understood. Analysis of recombinant G473D sulfite oxidase indicated that it is severely impaired both in the ability to bind sulfite and in catalysis, with a second-order rate constant 5 orders of magnitude lower than that of the wild type. To elucidate the specific reasons for the severe effects seen in the G473D variant, several other variants were created, including G473A, G473W, and the double mutant R212A/G473D. Despite the inability to form a stable dimer, the G473W variant had 5-fold higher activity than G473D and nearly wild-type activity at pH 7.0 when ferricyanide was the electron acceptor. In contrast, the R212A/G473D variant demonstrated some ability to oligomerize but had undetectable activity. The G473A variant retained the ability to dimerize and had steady-state activity that was comparable to that of the wild type. Furthermore, stopped-flow analysis of the reductive half-reaction of this variant yielded a rate constant nearly 3 times higher than that of the wild type. Examination of the secondary structures of the variants by CD spectroscopy indicated significant random-coil formation in G473D, G473W, and R212A/G473D. These results demonstrate that both the charge and the large size of an Asp residue in this position contribute to the severe effects seen in a patient with the G473D mutation, by causing partial misfolding and monomerization of sulfite oxidase and attenuating both substrate binding and catalytic efficiency during the reaction cycle.

High-resolution EXAFS of the active site of human sulfite oxidase: comparison with density functional theory and X-ray crystallographic results.

Much of our knowledge about molybdenum enzymes has originated from EXAFS spectroscopy. This technique provides excellent bond-length accuracy but has only limited bond-length resolution. We have used EXAFS spectroscopy with an extended data range in an attempt to improve bond-length resolution for the molybdenum enzyme sulfite oxidase. The Mo site of sulfite oxidase has two oxygen and three Mo-S ligands (two from cofactor dithiolene plus a cysteine). For the oxidized (Mo(VI)) enzyme, we find that the three Mo-S bond lengths are very similar (within 0.05 A) at 2.41 A, as are the Mo=O ligands at 1.72 A. Density functional theory shows that this is consistent with the proposed active-site structure. The reduced (Mo(IV)) enzyme shows two Mo-S bond lengths at 2.35 A and one at 2.41 A (assigned to cofactor dithiolene and cysteine, respectively, from DFT), together with one Mo=O at 1.72 A and one Mo-OH(2) at 2.30 A.

The pathogenic human sulfite oxidase mutants G473D and A208D are defective in intramolecular electron transfer.

Mutations G473D and A208D were identified in patients with isolated sulfite oxidase (SO) deficiency, and the equivalent amino acids (G451 and A186, respectively) have been localized to the vicinity of the molybdopterin active site in the X-ray structure of chicken SO [Kisker, C., Schindelin, H., Pacheco, A., Wehbi, W., Garrett, R. M., Rajagopalan, K. V., Enemark, J. H., and Rees, D. C. (1997) Cell 91, 973-983]. To assess the effects of these mutations in human SO, steady-state kinetic studies of enzyme turnover and laser flash photolysis measurements of intramolecular electron transfer (IET) rate constants between the reduced heme [Fe(II)] and Mo(VI) centers were carried out in the recombinant G473D, G473A, G473W, G473D/R212A, and A208D human SO mutants. In the G473D and A208D mutants, the IET rate constants at pH 6.0 are decreased by 3 orders of magnitude relative to that of the wild type. Steady-state kinetic measurements indicate that the IET process is the rate-limiting step in the catalytic cycle of these two mutants. Thus, the large decreases in the IET rate constants and the kcat values, and the large increases in the Km(sulfite) values, rationalize the fatal impact of these mutations. Far-UV CD spectra of G473D indicate that the protein backbone conformation is remarkably changed, and the sedimentation equilibrium indicates that the protein is monomeric. Furthermore, EPR studies also suggest that the active site structure of the Mo(V) form of A208D is different from that of the wild type. In contrast, similar studies on G473A show that it is dimeric, that its Mo(V) active site structure is similar to that of the wild type, and that its IET rate constant is only 2.6-fold smaller than that of the wild type. IET in G473W is severely impaired, and no IET is observed for G473D/R212A. In chicken SO, the equivalent residues (G451 and A186) are both buried inside the protein. Thus, for human SO, the mutations to charged residues at the equivalent sites most likely cause crucial global or localized structural changes, and expose an alternative docking site that may compete with the Mo domain for docking of the heme, thereby retarding IET and efficient catalytic turnover of the sulfite oxidation reaction.

Structural insights into sulfite oxidase deficiency.

Sulfite oxidase deficiency is a lethal genetic disease that results from defects either in the genes encoding proteins involved in molybdenum cofactor biosynthesis or in the sulfite oxidase gene itself. Several point mutations in the sulfite oxidase gene have been identified from patients suffering from this disease worldwide. Although detailed biochemical analyses have been carried out on these mutations, no structural data could be obtained because of problems in crystallizing recombinant human and rat sulfite oxidases and the failure to clone the chicken sulfite oxidase gene. We synthesized the gene for chicken sulfite oxidase de novo, working backward from the amino acid sequence of the native chicken liver enzyme by PCR amplification of a series of 72 overlapping primers. The recombinant protein displayed the characteristic absorption spectrum of sulfite oxidase and exhibited steady state and rapid kinetic parameters comparable with those of the tissue-derived enzyme. We solved the crystal structures of the wild type and the sulfite oxidase deficiency-causing R138Q (R160Q in humans) variant of recombinant chicken sulfite oxidase in the resting and sulfate-bound forms. Significant alterations in the substrate-binding pocket were detected in the structure of the mutant, and a comparison between the wild type and mutant protein revealed that the active site residue Arg-450 adopts different conformations in the presence and absence of bound sulfate. The size of the binding pocket is thereby considerably reduced, and its position relative to the cofactor is shifted, causing an increase in the distance of the sulfur atom of the bound sulfate to the molybdenum.