Mechanistic complexities of sulfite oxidase: An enzyme with multiple domains, subunits, and cofactors.

Sulfite oxidase (SO) deficiency, an inherited disease that causes severe neonatal neurological problems and early death, arises from defects in the biosynthesis of the molybdenum cofactor (Moco) (general sulfite oxidase deficiency) or from inborn errors in the SUOX gene for SO (isolated sulfite oxidase deficiency, ISOD). The X-ray structure of the highly homologous homonuclear dimeric chicken sulfite oxidase (cSO) provides a template for locating ISOD mutation sites in human sulfite oxidase (hSO). Catalysis occurs within an individual subunit of hSO, but mutations that disrupt the hSO dimer are pathological. The catalytic cycle of SO involves five metal oxidation states (MoVI, MoV, MoIV, FeIII, FeII), two intramolecular electron transfer (IET) steps, and couples a two-electron oxygen atom transfer reaction at the Mo center with two one-electron transfers from the integral b-type heme to exogenous cytochrome c, the physiological oxidant. Several ISOD examples are analyzed using steady-state, stopped-flow, and laser flash photolysis kinetics and physical measurements of recombinant variants of hSO and native cSO. In the structure of cSO, Mo…Fe = 32 Å, much too long for efficient IET through the protein. Interdomain motion that brings the Mo and heme centers closer together to facilitate IET is supported indirectly by decreasing the length of the interdomain tether, by changes in the charges of surface residues of the Mo and heme domains, as well as by preliminary molecular dynamics calculations. However, direct dynamic measurements of interdomain motion are in their infancy.

{Moco}n, (n = 0-8): A general formalism for describing the highly covalent molybdenum cofactor of sulfite oxidase and related Mo enzymes.

Over 50 molybdenum enzymes in three distinct families (sulfite oxidase, xanthine oxidase, DMSO reductase) are known, and representative X-ray crystal structures are available for all families. Structural analogues that replicate the coordination about the Mo atom in the absence of surrounding protein have been synthesized and characterized. The properties of metal complexes of non-innocent dithiolene ligands and their oxidized counter parts, dithiones, are summarized. Pulsed electron paramagnetic resonance (EPR) spectroscopy of the 33S-labeled molybdenum domain of catalytically active bioengineered sulfite oxidase has clearly demonstrated delocalization of electron density from MoV to the dithiolene component of the molybdenum cofactor (Moco) of the enzyme. Moco is highly covalent and has three redox active components: the Mo atom; the dithiolene; and the pterin. In principle, Moco can have a total of eight redox states, making it one of the most redox rich cofactors in biology. The {Moco}n formalism, developed here, gives the total number of electrons (n) associated with a particular redox state of Moco. This flexible notation eliminates the need to assign a specific oxidation state to each of the three components of Moco and allows for internal redistribution of electrons among the components upon substrate binding, changes in the surrounding network of hydrogen bonds, conformational changes, and catalysis. An unexpected result is that sulfite oxidase (an oxotransferase) is predicted to utilize the {Moco}4-6 electron distributions during catalysis, whereas xanthine oxidase (a hydroxylase) is predicted to utilize the {Moco}6-8 electron distributions during catalysis.

Metal-Dithiolene Bonding Contributions to Pyranopterin Molybdenum Enzyme Reactivity.

Here we highlight past work on metal-dithiolene interactions and how the unique electronic structure of the metal-dithiolene unit contributes to both the oxidative and reductive half reactions in pyranopterin molybdenum and tungsten enzymes. The metallodithiolene electronic structures detailed here were interrogated using multiple ground and excited state spectroscopic probes on the enzymes and their small molecule analogs. The spectroscopic results have been interpreted in the context of bonding and spectroscopic calculations, and the pseudo-Jahn-Teller effect. The dithiolene is a unique ligand with respect to its redox active nature, electronic synergy with the pyranopterin component of the molybdenum cofactor, and the ability to undergo chelate ring distortions that control covalency, reduction potential, and reactivity in pyranopterin molybdenum and tungsten enzymes.

Vibrational Control of Covalency Effects Related to the Active Sites of Molybdenum Enzymes.

A multitechnique spectroscopic and theoretical study of the Cp2M(benzenedithiolato) (M = Ti, V, Mo; Cp = η5-C5H5) series provides deep insight into dithiolene electronic structure contributions to electron transfer reactivity and reduction potential modulation in pyranopterin molybdenum enzymes. This work explains the magnitude of the dithiolene folding distortion and the concomitant changes in metal-ligand covalency that are sensitive to electronic structure changes as a function of d-electron occupancy in the redox orbital. It is shown that the large fold angle differences correlate with covalency, and the fold angle distortion is due to a pseudo-Jahn-Teller (PJT) effect. The PJT effect in these and related transition metal dithiolene systems arises from the small energy differences between metal and sulfur valence molecular orbitals, which uniquely poise these systems for dramatic geometric and electronic structure changes as the oxidation state changes. Herein, we have used a combination of resonance Raman, magnetic circular dichroism, electron paramagnetic resonance, and UV photoelectron spectroscopies to explore the electronic states involved in the vibronic coupling mechanism. Comparison between the UV photoelectron spectroscopy (UPS) of the d2 M = Mo complex and the resonance Raman spectra of the d1 M = V complex reveals the power of this combined spectroscopic approach. Here, we observe that the UPS spectrum of Cp2Mo(bdt) contains an intriguing vibronic progession that is dominated by a "missing-mode" that is composed of PJT-active distortions. We discuss the relationship of the PJT distortions to facile electron transfer in molybdenum enzymes.

Consensus structures of the Mo(v) sites of sulfite-oxidizing enzymes derived from variable frequency pulsed EPR spectroscopy, isotopic labelling and DFT calculations.

Sulfite-oxidizing enzymes from eukaryotes and prokaryotes have five-coordinate distorted square-pyramidal coordination about the molybdenum atom. The paramagnetic Mo(v) state is easily generated, and over the years four distinct CW EPR spectra have been identified, depending upon enzyme source and the reaction conditions, namely high and low pH (hpH and lpH), phosphate inhibited (Pi) and sulfite (or blocked). Extensive studies of these paramagnetic forms of sulfite-oxidizing enzymes using variable frequency pulsed electron spin echo (ESE) spectroscopy, isotopic labeling and density functional theory (DFT) calculations have led to the consensus structures that are described here. Errors in some of the previously proposed structures are corrected.

Sulfite-oxidizing enzymes.

Sulfite-oxidizing enzymes (SOEs) are molybdenum enzymes that exist in almost all forms of life where they carry out important functions in protecting cells and organisms against sulfite-induced damage. Due to their nearly ubiquitous presence in living cells, these enzymes can be assumed to be evolutionarily ancient, and this is reflected in the fact that the basic domain architecture and fold structure of all sulfite-oxidizing enzymes studied so far are similar. The Mo centers of all SOEs have five-coordinate square pyramidal coordination geometry, which incorporates a pyranopterin dithiolene cofactor. However, significant differences exist in the quaternary structure of the enzymes, as well as in the kinetic properties and the nature of the electron acceptors used. In addition, some SOEs also contain an integral heme group that participates in the overall catalytic cycle. Catalytic turnover involves the paramagnetic Mo(V) oxidation state, and EPR spectroscopy, especially high-resolution pulsed EPR spectroscopy, provides detailed information about the molecular and electronic structure of the Mo center and the Mo-based sulfite oxidation reaction.

Kinetic results for mutations of conserved residues H304 and R309 of human sulfite oxidase point to mechanistic complexities.

Several point mutations in the gene of human sulfite oxidase (hSO) result in isolated sulfite oxidase deficiency, an inherited metabolic disorder. Three conserved residues (H304, R309, K322) are hydrogen bonded to the phosphate group of the molybdenum cofactor, and the R309H and K322R mutations are responsible for isolated sulfite oxidase deficiency. The kinetic effects of the K322R mutation have been previously reported (Rajapakshe et al., Chem. Biodiversity, 2012, 9, 1621-1634); here we investigate several mutants of H304 and R309 by steady-state kinetics, laser flash photolysis studies of intramolecular electron transfer (IET), and spectroelectrochemistry. An unexpected result is that all of the mutants show decreased rates of IET but increased steady-state rates of catalysis. However, in all cases the rate of IET is greater than the overall turnover rate, showing that IET is not the rate determining step for any of the mutations.

Pulsed electron paramagnetic resonance spectroscopy of (33)S-labeled molybdenum cofactor in catalytically active bioengineered sulfite oxidase.

Molybdenum enzymes contain at least one pyranopterin dithiolate (molybdopterin, MPT) moiety that coordinates Mo through two dithiolate (dithiolene) sulfur atoms. For sulfite oxidase (SO), hyperfine interactions (hfi) and nuclear quadrupole interactions (nqi) of magnetic nuclei (I ≠ 0) near the Mo(V) (d(1)) center have been measured using high-resolution pulsed electron paramagnetic resonance (EPR) methods and interpreted with the help of density functional theory (DFT) calculations. These have provided important insights about the active site structure and the reaction mechanism of the enzyme. However, it has not been possible to use EPR to probe the dithiolene sulfurs directly since naturally abundant (32)S has no nuclear spin (I = 0). Here we describe direct incorporation of (33)S (I = 3/2), the only stable magnetic sulfur isotope, into MPT using controlled in vitro synthesis with purified proteins. The electron spin echo envelope modulation (ESEEM) spectra from (33)S-labeled MPT in this catalytically active SO variant are dominated by the "interdoublet" transition arising from the strong nuclear quadrupole interaction, as also occurs for the (33)S-labeled exchangeable equatorial sulfite ligand [ Klein, E. L., et al. Inorg. Chem. 2012 , 51 , 1408 - 1418 ]. The estimated experimental hfi and nqi parameters for (33)S (aiso = 3 MHz and e(2)Qq/h = 25 MHz) are in good agreement with those predicted by DFT. In addition, the DFT calculations show that the two (33)S atoms are indistinguishable by EPR and reveal a strong intermixing between their out-of-plane pz orbitals and the dxy orbital of Mo(V).

Probing the role of a conserved salt bridge in the intramolecular electron transfer kinetics of human sulfite oxidase.

Sulfite oxidase (SO) is a vital metabolic enzyme that catalyzes the oxidation of toxic sulfite to sulfate. The proposed mechanism of this molybdenum cofactor dependent enzyme involves two one-electron intramolecular electron transfer (IET) steps from the molybdenum center to the iron of the b 5-type heme and two one-electron intermolecular electron transfer steps from the heme to cytochrome c. This work focuses on how the electrostatic interaction between two conserved amino acid residues, R472 and D342, in human SO (hSO) affects catalysis. The hSO variants R472M, R472Q, R472K, R472D, and D342K were created to probe the effect of the position of the salt bridge charges, along with the interaction between these two residues. With the exception of R472K, these variants all showed a significant decrease in their IET rate constants, k et, relative to wild-type hSO, indicating that the salt bridge between residues 472 and 342 is important for rapid IET. Surprisingly, however, except for R472K and R472D, all of the variants show k cat values higher than their corresponding k et values. The turnover number for R472D is about the same as k et, which suggests that the change in this variant is rate-limiting in catalysis. Direct spectroelectrochemical determination of the Fe(III/II) reduction potentials of the heme and calculation of the Mo(VI/V) potentials revealed that all of the variants affected the redox potentials of both metal centers, probably due to changes in their environments. Thus, the position of the positive charge of R472 and that of the negative charge of D342 are both important in hSO, and changing either the position or the nature of these charges perturbs IET and catalysis.

Applications of pulsed EPR spectroscopy to structural studies of sulfite oxidizing enzymes().

Sulfite oxidizing enzymes (SOEs), including sulfite oxidase (SO) and bacterial sulfite dehydrogenase (SDH), catalyze the oxidation of sulfite (SO(3) (2-)) to sulfate (SO(4) (2-)). The active sites of SO and SDH are nearly identical, each having a 5-coordinate, pseudo-square-pyramidal Mo with an axial oxo ligand and three equatorial sulfur donor atoms. One sulfur is from a conserved Cys residue and two are from a pyranopterindithiolene (molybdopterin, MPT) cofactor. The identity of the remaining equatorial ligand, which is solvent-exposed, varies during the catalytic cycle. Numerous in vitro studies, particularly those involving electron paramagnetic resonance (EPR) spectroscopy of the Mo(V) states of SOEs, have shown that the identity and orientation of this exchangeable equatorial ligand depends on the buffer pH, the presence and concentration of certain anions in the buffer, as well as specific point mutations in the protein. Until very recently, however, EPR has not been a practical technique for directly probing specific structures in which the solvent-exposed, exchangeable ligand is an O, OH(-), H(2)O, SO(3) (2-), or SO(4) (2-) group, because the primary O and S isotopes ((16)O and (32)S) are magnetically silent (I = 0). This review focuses on the recent advances in the use of isotopic labeling, variable-frequency high resolution pulsed EPR spectroscopy, synthetic model compounds, and DFT calculations to elucidate the roles of various anions, point mutations, and steric factors in the formation, stabilization, and transformation of SOE active site structures.

Effects of mutating aromatic surface residues of the heme domain of human sulfite oxidase on its heme midpoint potential, intramolecular electron transfer, and steady-state kinetics.

Human sulfite oxidase (hSO), an essential molybdoheme enzyme, catalyzes the oxidation of toxic sulfite to sulfate. The proposed catalytic cycle includes two, one-electron intramolecular electron transfers (IET) between the molybdenum (Mo) and the heme domains. Rapid IET rates are ascribed to conformational changes that bring the two domains into close proximity to one another. Previous studies of hSO have focused on the roles of conserved residues near the Mo active site and on the tether that links the two domains. Here four aromatic surface residues on the heme domain (phenylalanine 57 (F57), phenylalanine 79 (F79), tyrosine 83 (Y83), and histidine 90 (H90)) have been mutated, and their involvement in IET rates, the heme midpoint potential, and the catalytic activity of hSO have been investigated using laser flash photolysis, spectroelectrochemistry, and steady-state kinetics, respectively. The results indicate that the size and hydrophobicity of F57 play an important role in modulating the heme potential and that F57 also affects the IET rates. The data also suggest that important interactions of H90 with a heme propionate group destabilize the Fe(III) state of the heme. The positive charge on H90 at pH ≤ 7.0 may decrease the electrostatic interaction between the Mo and heme domains, thereby decreasing the IET rates of wt hSO at low pH. Lastly, mutations of F79 and Y83, which are located on the surface of the heme domain, but not in direct contact with the heme or the propionate groups, have little effect on either IET or the heme potential.

Identity of the exchangeable sulfur-containing ligand at the Mo(V) center of R160Q human sulfite oxidase.

In our previous study of the fatal R160Q mutant of human sulfite oxidase (hSO) at low pH (Astashkin et al. J. Am. Chem. Soc.2008, 130, 8471-8480), a new Mo(V) species, denoted "species 1", was observed at low pH values. Species 1 was ascribed to a six-coordinate Mo(V) center with an exchangeable terminal oxo ligand and an equatorial sulfate group on the basis of pulsed EPR spectroscopy and (33)S and (17)O labeling. Here we report new results for species 1 of R160Q, based on substitution of the sulfur-containing ligand by a phosphate group, pulsed EPR spectroscopy in K(a)- and W-bands, and extensive density functional theory (DFT) calculations applied to large, more realistic molecular models of the enzyme active site. The combined results unambiguously show that species 1 has an equatorial sulfite as the only exchangeable ligand. The two types of (17)O signals that are observed arise from the coordinated and remote oxygen atoms of the sulfite ligand. A typical five-coordinate Mo(V) site is compatible with the observed and calculated EPR parameters.

Determination of the distance between the Mo(V) and Fe(III) heme centers of wild type human sulfite oxidase by pulsed EPR spectroscopy.

Intramolecular electron transfer (IET) between the molybdenum and heme centers of vertebrate sulfite oxidase (SO) is proposed to be a key step in the catalytic cycle of the enzyme. However, the X-ray crystallographic distance between these centers, R(MoFe) = 32.3 Å, appears to be too long for the rapid IET rates observed in liquid solution. The Mo and heme domains are linked by a flexible tether, and it has been proposed that dynamic interdomain motion brings the two metal centers closer together and thereby facilitates rapid IET. To date, there have been no direct distance measurements for SO in solution that would support or contradict this model. In this work, pulsed electron-electron double resonance (ELDOR) and relaxation induced dipolar modulation enhancement (RIDME) techniques were used to obtain information about R(MoFe) in the Mo(V)Fe(III) state of wild type recombinant human SO in frozen glassy solution. Surprisingly, the data obtained suggest a fixed structure with R(MoFe) = 32 Å, similar to that determined by X-ray crystallography for chicken SO, although the orientation of the R(MoFe) radius-vector with respect to the heme center was found to be somewhat different. The implications of these findings for the flexible tether model are discussed.

Intramolecular electron transfer in sulfite-oxidizing enzymes: probing the role of aromatic amino acids.

Sulfite oxidase (SO) is a molybdoheme enzyme that is important in sulfur catabolism, and mutations in the active site region are known to cause SO deficiency disorder in humans. This investigation probes the effects that mutating aromatic residues (Y273, W338, and H337) in the molybdenum-containing domain of human SO have on both the intramolecular electron transfer (IET) rate between the molybdenum and iron centers using laser flash photolysis and on catalytic turnover via steady-state kinetic analysis. The W338 and H337 mutants show large decreases in their IET rate constants (k (ET)) relative to the wild-type values, suggesting the importance of these residues for rapid IET. In contrast, these mutants are catalytically competent and exhibit higher k (cat) values than their corresponding k (ET), implying that these two processes involve different conformational states of the protein. Redox potential investigations using spectroelectrochemistry revealed that these aromatic residues close to the molybdenum center affect the potential of the presumably distant heme center in the resting state (as shown by the crystal structure of chicken SO), suggesting that the heme may be interacting with these residues during IET and/or catalytic turnover. These combined results suggest that in solution human SO may adopt different conformations for IET and for catalysis in the presence of the substrate. For IET the H337/W338 surface residues may serve as an alternative-docking site for the heme domain. The similarities between the mutant and wild-type EPR spectra indicate that the active site geometry around the Mo(V) center is not changed by the mutations studied here.

Metal-sulfur valence orbital interaction energies in metal-dithiolene complexes: determination of charge and overlap interaction energies by comparison of core and valence ionization energy shifts.

The electronic interactions between metals and dithiolenes are important in the biological processes of many metalloenzymes as well as in diverse chemical and material applications. Of special note is the ability of the dithiolene ligand to support metal centers in multiple coordination environments and oxidation states. To better understand the nature of metal-dithiolene electronic interactions, new capabilities in gas-phase core photoelectron spectroscopy for molecules with high sublimation temperatures have been developed and applied to a series of molecules of the type Cp(2)M(bdt) (Cp = η(5)-cyclopentadienyl, M = Ti, V, Mo, and bdt = benzenedithiolato). Comparison of the gas-phase core and valence ionization energy shifts provides a unique quantitative energy measure of valence orbital overlap interactions between the metal and the sulfur orbitals that is separated from the effects of charge redistribution. The results explain the large amount of sulfur character in the redox-active orbitals and the 'leveling' of oxidation state energies in metal-dithiolene systems. The experimentally determined orbital interaction energies reveal a previously unidentified overlap interaction of the predominantly sulfur HOMO of the bdt ligand with filled π orbitals of the Cp ligands, suggesting that direct dithiolene interactions with other ligands bound to the metal could be significant for other metal-dithiolene systems in chemistry and biology.

Structural studies of the molybdenum center of mitochondrial amidoxime reducing component (mARC) by pulsed EPR spectroscopy and 17O-labeling.

Mitochondrial amidoxime reducing components (mARC-1 and mARC-2) represent a novel group of Mo-containing enzymes in eukaryotes. These proteins form the catalytic part of a three-component enzyme complex known to be responsible for the reductive activation of several N-hydroxylated prodrugs. No X-ray crystal structures are available for these enzymes as yet. A previous biochemical investigation [Wahl, B., et al. (2010) J. Biol. Chem., 285, 37847-37859 ] has revealed that two of the Mo coordination positions are occupied by sulfur atoms from a pyranopterindithiolate (molybdopterin, MPT) cofactor. In this work, we have used continuous wave and pulsed electron paramagnetic resonance (EPR) spectroscopy and density functional theoretical (DFT) calculations to determine the nature of remaining ligands in the Mo(V) state of the active site of mARC-2. Experiments with samples in D(2)O have identified the exchangeable equatorial ligand as a hydroxyl group. Experiments on samples in H(2)(17)O-enriched buffer have shown the presence of a slowly exchangeable axial oxo ligand. Comparison of the experimental (1)H and (17)O hyperfine interactions with those calculated using DFT has shown that the remaining nonexchangeable equatorial ligand is, most likely, protein-derived and that the possibility of an equatorial oxo ligand can be excluded.

Implications for the mechanism of sulfite oxidizing enzymes from pulsed EPR spectroscopy and DFT calculations for "difficult" nuclei.

The catalytic mechanisms of sulfite oxidizing enzymes (SOEs) have been investigated by multi-frequency pulsed EPR measurements of "difficult" magnetic nuclei (35.37Cl, 33S, 17O) associated with the Mo(v) center. Extensive DFT calculations have been used to relate the experimental magnetic resonance parameters of these nuclei to specific active site structures. This combined spectroscopic and computational approach has provided new insights concerning the structure/function relationships of the active sites of SOEs, including: (i) the exchange of oxo ligands; (ii) the nature of the blocked forms; and (iii) the role of Cl- in low pH forms.

Effects of large-scale amino acid substitution in the polypeptide tether connecting the heme and molybdenum domains on catalysis in human sulfite oxidase.

Sulfite oxidase (SO) is a molybdenum-cofactor-dependent enzyme that catalyzes the oxidation of sulfite to sulfate, the final step in the catabolism of the sulfur-containing amino acids, cysteine and methionine. The catalytic mechanism of vertebrate SO involves intramolecular electron transfer (IET) from molybdenum to the integral b-type heme of SO and then to exogenous cytochrome c. However, the crystal structure of chicken sulfite oxidase (CSO) has shown that there is a 32 Å distance between the Fe and Mo atoms of the respective heme and molybdenum domains, which are connected by a flexible polypeptide tether. This distance is too long to be consistent with the measured IET rates. Previous studies have shown that IET is viscosity dependent (Feng et al., Biochemistry, 2002, 41, 5816) and also dependent upon the flexibility and length of the tether (Johnson-Winters et al., Biochemistry, 2010, 49, 1290). Since IET in CSO is more rapid than in human sulfite oxidase (HSO) (Feng et al., Biochemistry, 2003, 42, 12235) the tether sequence of HSO has been mutated into that of CSO, and the resultant chimeric HSO enzyme investigated by laser flash photolysis and steady-state kinetics in order to study the specificity of the tether sequence of SO on the kinetic properties. Surprisingly, the IET kinetics of the chimeric HSO protein with the CSO tether sequence are slower than wildtype HSO. This observation raises the possibility that the composition of the non-conserved tether sequence of animal SOs may be optimized for individual species.