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.

Iron and zinc porphyrin linked MoO(dithiolene) complexes in relevance to electron transfer between Mo-cofactor and cytochrome b5 in sulfite oxidase.

Oxo-molybdenum (dithiolene) complexes covalently linked individually to iron and zinc porphyrin have been synthesized to show an electron transfer between the two metal centres in relevance to electron transfer from Mo-cofactor to cytochrome b5 domains in the oxidative half of the catalytic cycle of native sulfite oxidase. This association has been investigated by electrochemical, EPR measurement and X-ray absorbance spectroscopy techniques.

{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.

Functional mononuclear molybdenum enzymes: challenges and triumphs in molecular cloning, expression, and isolation.

Mononuclear molybdenum enzymes catalyze a variety of reactions that are essential in the cycling of nitrogen, carbon, arsenic, and sulfur. For decades, the structure and function of these crucial enzymes have been investigated to develop a fundamental knowledge for this vast family of enzymes and the chemistries they carry out. Therefore, obtaining abundant quantities of active enzyme is necessary for exploring this family's biochemical capability. This mini-review summarizes the methods for overexpressing mononuclear molybdenum enzymes in the context of the challenges encountered in the process. Effective methods for molybdenum cofactor synthesis and incorporation, optimization of expression conditions, improving isolation of active vs. inactive enzyme, incorporation of additional prosthetic groups, and inclusion of redox enzyme maturation protein chaperones are discussed in relation to the current molybdenum enzyme literature. This article summarizes the heterologous and homologous expression studies providing underlying patterns and potential future directions.

Hot Carriers and Photothermal Effects of Monolayer MoOx for Promoting Sulfite Oxidase Mimetic Activity.

Local surface plasmon resonance (LSPR)-enhanced catalysis has brought a substantial amount of opportunities across various disciplines such as photocatalysis, photodetection, and photothermal therapeutics. Plasmon-induced photothermal and hot carriers effects have also been utilized to activate the enzyme-like reactions. Compared with natural enzymes, the relatively low catalytic performance of nanozymes severely hampered the potential applications in the field of biomedicine. For these issues mentioned above, herein, we demonstrate a highly efficient sulfite oxidase (SuOx) mimetic performance of plasmonic monolayer MoOx (ML-MoOx) upon LSPR excitation. We also established that the considerable photothermal effect and the injection of hot carriers induced by LSPR are responsible for promoting the SuOx activity of ML-MoOx. The high transient local temperature on the surface of ML-MoOx generated by the photothermal effect facilitates to impact the reaction velocity and feed the SuOx-like activity, while the generation of hot carriers which are suggested as predominant effects catalyzes the oxidation of sulfite to sulfate through significantly decreasing the activation energy for the SuOx-like reaction. These investigations present a contribution to the basic understanding of plasmon-enhanced enzyme-like reaction and provided an insight into the optimization of the SuOx mimetic performance of nanomaterials.

Impaired mitochondrial maturation of sulfite oxidase in a patient with severe sulfite oxidase deficiency.

Sulfite oxidase (SO) is encoded by the nuclear SUOX gene and catalyzes the final step in cysteine catabolism thereby oxidizing sulfite to sulfate. Oxidation of sulfite is dependent on two cofactors within SO, a heme and the molybdenum cofactor (Moco), the latter forming the catalytic site of sulfite oxidation. SO localizes to the intermembrane space of mitochondria where both-pre-SO processing and cofactor insertion-are essential steps during SO maturation. Isolated SO deficiency (iSOD) is a rare inborn error of metabolism caused by mutations in the SUOX gene that lead to non-functional SO. ISOD is characterized by rapidly progressive neurodegeneration and death in early infancy. We diagnosed an iSOD patient with homozygous mutation of SUOX at c.1084G>A replacing Gly362 to serine. To understand the mechanism of disease, we expressed patient-derived G362S SO in Escherichia coli and surprisingly found full catalytic activity, while in patient fibroblasts no SO activity was detected, suggesting differences between bacterial and human expression. Moco reconstitution of apo-G362S SO was found to be approximately 90-fold reduced in comparison to apo-WT SO in vitro. In line, levels of SO-bound Moco in cells overexpressing G362S SO were significantly reduced compared to cells expressing WT SO providing evidence for compromised maturation of G362S SO in cellulo. Addition of molybdate to culture medium partially rescued impaired Moco binding of G362S SO and restored SO activity in patient fibroblasts. Thus, this study demonstrates the importance of the orchestrated maturation of SO and provides a first case of Moco-responsive iSOD.

Nitrite-dependent nitric oxide synthesis by molybdenum enzymes.

Nitric oxide (NO) is an important gasotransmitter involved in numerous intra- and intercellular signaling events. In addition to the oxidative pathway of NO generation, which includes three NO synthase (NOS) isoforms in mammals, a reductive pathway contributes to NO generation. In this pathway, nitrite is reduced to NO by various metal-containing proteins. Among these, all members of the eukaryotic molybdenum (Mo)-dependent enzyme family were found to be able to reduce nitrite to NO. This Review focuses on the current state of research in the field of Mo-dependent nitrite reduction in eukaryotes. An overview on the five eukaryotic Mo-enzymes is given, and similarities as well as differences in their nitrite reduction mechanisms are presented and discussed in the context of physiological relevance.

QM/MM study of the reaction mechanism of sulfite oxidase.

Sulfite oxidase is a mononuclear molybdenum enzyme that oxidises sulfite to sulfate in many organisms, including man. Three different reaction mechanisms have been suggested, based on experimental and computational studies. Here, we study all three with combined quantum mechanical (QM) and molecular mechanical (QM/MM) methods, including calculations with large basis sets, very large QM regions (803 atoms) and QM/MM free-energy perturbations. Our results show that the enzyme is set up to follow a mechanism in which the sulfur atom of the sulfite substrate reacts directly with the equatorial oxo ligand of the Mo ion, forming a Mo-bound sulfate product, which dissociates in the second step. The first step is rate limiting, with a barrier of 39-49 kJ/mol. The low barrier is obtained by an intricate hydrogen-bond network around the substrate, which is preserved during the reaction. This network favours the deprotonated substrate and disfavours the other two reaction mechanisms. We have studied the reaction with both an oxidised and a reduced form of the molybdopterin ligand and quantum-refinement calculations indicate that it is in the normal reduced tetrahydro form in this protein.

Electrochemical Investigations of Hydrogenases and Other Enzymes That Produce and Use Solar Fuels.

Many enzymes that produce or transform small molecules such as O2, H2, and CO2 embed inorganic cofactors based on transition metals. Their active site, where the chemical reaction occurs, is buried in and protected by the protein matrix, and connected to the solvent in several ways: chains of redox cofactors mediate long-range electron transfer; static or dynamic tunnels guide the substrate, product and inhibitors; amino acids and water molecules transfer protons. The catalytic mechanism of these enzymes is therefore delocalized over the protein and involves many different steps, some of which determine the response of the enzyme under conditions of stress (extreme redox conditions, presence of inhibitors, light), the catalytic rates in the two directions of the reaction and their ratio (the "catalytic bias"). Understanding all the steps in the catalytic cycle, including those that occur on sites of the protein that are remote from the active site, requires a combination of biochemical, structural, spectroscopic, theoretical, and kinetic methods. Here we argue that kinetics should be used to the fullest extent, by extracting quantitative information from the comparison of data and kinetic models and by exploring the combination of experimental kinetics and theoretical chemistry. In studies of these catalytic mechanisms, direct electrochemistry, the technique which we use and contribute to develop, has become unescapable. It simply consists in monitoring the changes in activity of an enzyme that is wired to an electrode by recording an electric current. We have described kinetic models that can be used to make sense of these data and to learn about various aspects of the mechanism that are difficult to probe using more conventional methods: long-range electron transfer, diffusion along gas channels, redox-driven (in)activations, active site chemistry and photoreactivity under conditions of turnover. In this Account, we highlight a few results that illustrate our approach. We describe how electrochemistry can be used to monitor substrate and inhibitor diffusion along the gas channels of hydrogenases and we discuss how the kinetics of intramolecular diffusion relates to global properties such as resistance to oxygen and catalytic bias. The kinetics and/or thermodynamics of intramolecular electron transfer may also affect the catalytic bias, the catalytic potentials on either side of the equilibrium potential, and the overpotentials for catalysis (defined as the difference between the catalytic potentials and the open circuit potential). This is understood by modeling the shape of the steady-state catalytic response of the enzyme. Other determinants of the catalytic rate, such as domain motions, have been probed by examining the transient catalytic response recorded at fast scan rates. Last, we show that combining electrochemical investigations and MD, DFT, and TD-DFT calculations is an original way of probing the reactivity of the H-cluster of hydrogenase, in particular its reactions with CO, O2, and light. This approach contrasts with the usual strategy which aims at stabilizing species that are presumed to be catalytic intermediates, and determining their structure using spectroscopic or structural methods.

Chemical Biology of H2S Signaling through Persulfidation.

Signaling by H2S is proposed to occur via persulfidation, a posttranslational modification of cysteine residues (RSH) to persulfides (RSSH). Persulfidation provides a framework for understanding the physiological and pharmacological effects of H2S. Due to the inherent instability of persulfides, their chemistry is understudied. In this review, we discuss the biologically relevant chemistry of H2S and the enzymatic routes for its production and oxidation. We cover the chemical biology of persulfides and the chemical probes for detecting them. We conclude by discussing the roles ascribed to protein persulfidation in cell signaling pathways.

Bioelectronic Circuit on a 3D Electrode Architecture: Enzymatic Catalysis Interconnected with Photosystem I.

Artificial light-driven signal chains are particularly important for the development of systems converting light into a current, into chemicals or for light-induced sensing. Here, we report on the construction of an all-protein, light-triggered, catalytic circuit based on photosystem I, cytochrome c (cyt c) and human sulfite oxidase (hSOX). The defined assembly of all components using a modular design results in an artificial biohybrid electrode architecture, combining the photophysical features of PSI with the biocatalytic properties of hSOX for advanced light-controlled bioelectronics. The working principle is based on a competitive switch between electron supply from the electrode or by enzymatic substrate conversion.

Does the DFT Self-Interaction Error Affect Energies Calculated in Proteins with Large QM Systems?

We have examined how the self-interaction error in density-functional theory (DFT) calculations affects energies calculated on large systems (600-1000 atoms) involving several charged groups. We employ 18 different quantum mechanical (QM) methods, including Hartree-Fock, as well as pure, hybrid, and range-separated DFT methods. They are used to calculate reaction and activation energies for three different protein models in vacuum, in a point-charge surrounding, or with a continuum-solvent model. We show that pure DFT functionals give rise to a significant delocalization of the charges in charged groups in the protein, typically by ∼0.1 e, as evidenced from the Mulliken charges. This has a clear effect on how the surroundings affect calculated reaction and activation energies, indicating that these methods should be avoided for DFT calculations on large systems. Fortunately, methods such as CAM-B3LYP, BHLYP, and M06-2X give results that agree within a few kilojoules per mole, especially when the calculations are performed in a point-charge surrounding. Therefore, we recommend these methods to estimate the effect of the surroundings with large QM systems (but other QM methods may be used to study the intrinsic reaction and activation energies).

Chemical systems modeling the d1 Mo(V) states of molybdenum enzymes.

This review focuses on the synthesis, properties and electron paramagnetic resonance (EPR), electron spin echo envelope modulation (ESEEM) and electron-nuclear double resonance (ENDOR) spectroscopy of mononuclear d1 oxo- and sulfido-Mo(V) complexes relevant to the understanding of the EPR-active Mo(V) forms of pterin-containing molybdenum enzymes.

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.

Shifting the metallocentric molybdoenzyme paradigm: the importance of pyranopterin coordination.

In this review, we test the hypothesis that pyranopterin coordination plays a critical role in defining substrate reactivities in the four families of mononuclear molybdenum and tungsten enzymes (Mo/W-enzymes). Enzyme families containing a single pyranopterin dithiolene chelate have been demonstrated to have reactivity towards two (sulfite oxidase, SUOX-fold) and five (xanthine dehydrogenase, XDH-fold) types of substrate, whereas the major family of enzymes containing a bis-pyranopterin dithiolene chelate (dimethylsulfoxide reductase, DMSOR-fold) is reactive towards eight types of substrate. A second bis-pyranopterin enzyme (aldehyde oxidoreductase, AOR-fold) family catalyzes a single type of reaction. The diversity of reactions catalyzed by each family correlates with active site variability, and also with the number of pyranopterins and their coordination by the protein. In the case of the AOR-fold enzymes, inflexibility of pyranopterin coordination correlates with their limited substrate specificity (oxidation of aldehydes). In examples of the SUOX-fold and DMSOR-fold enzymes, we observe three types of histidine-containing charge-transfer relays that can: (1) connect the piperazine ring of the pyranopterin to the substrate-binding site (SUOX-fold enzymes); (2) provide inter-pyranopterin communication (DMSOR-fold enzymes); and (3) connect a pyran ring oxygen to deeply buried water molecules (the DMSOR-fold NarGHI-type nitrate reductases). Finally, sequence data mining reveals a number of bacterial species whose predicted proteomes contain large numbers (up to 64) of Mo/W-enzymes, with the DMSOR-fold enzymes being dominant. These analyses also reveal an inverse correlation between Mo/W-enzyme content and pathogenicity.

Sulfite Oxidase Catalyzes Single-Electron Transfer at Molybdenum Domain to Reduce Nitrite to Nitric Oxide.

AIMS: Recent studies suggest that the molybdenum enzymes xanthine oxidase, aldehyde oxidase, and mARC exhibit nitrite reductase activity at low oxygen pressures. However, inhibition studies of xanthine oxidase in humans have failed to block nitrite-dependent changes in blood flow, leading to continued exploration for other candidate nitrite reductases. Another physiologically important molybdenum enzyme­sulfite oxidase (SO)­has not been extensively studied. RESULTS: Using gas-phase nitric oxide (NO) detection and physiological concentrations of nitrite, SO functions as nitrite reductase in the presence of a one-electron donor, exhibiting redox coupling of substrate oxidation and nitrite reduction to form NO. With sulfite, the physiological substrate, SO only facilitates one turnover of nitrite reduction. Studies with recombinant heme and molybdenum domains of SO indicate that nitrite reduction occurs at the molybdenum center via coupled oxidation of Mo(IV) to Mo(V). Reaction rates of nitrite to NO decreased in the presence of a functional heme domain, mediated by steric and redox effects of this domain. Using knockdown of all molybdopterin enzymes and SO in fibroblasts isolated from patients with genetic deficiencies of molybdenum cofactor and SO, respectively, SO was found to significantly contribute to hypoxic nitrite signaling as demonstrated by activation of the canonical NO-sGC-cGMP pathway. INNOVATION: Nitrite binds to and is reduced at the molybdenum site of mammalian SO, which may be allosterically regulated by heme and molybdenum domain interactions, and contributes to the mammalian nitrate-nitrite-NO signaling pathway in human fibroblasts. CONCLUSION: SO is a putative mammalian nitrite reductase, catalyzing nitrite reduction at the Mo(IV) center.

Comparison of the active-site design of molybdenum oxo-transfer enzymes by quantum mechanical calculations.

There are three families of mononuclear molybdenum enzymes that catalyze oxygen atom transfer (OAT) reactions, named after a typical example from each family, viz., dimethyl sulfoxide reductase (DMSOR), sulfite oxidase (SO), and xanthine oxidase (XO). These families differ in the construction of their active sites, with two molybdopterin groups in the DMSOR family, two oxy groups in the SO family, and a sulfido group in the XO family. We have employed density functional theory calculations on cluster models of the active sites to understand the selection of molybdenum ligands in the three enzyme families. Our calculations show that the DMSOR active site has a much stronger oxidative power than the other two sites, owing to the extra molybdopterin ligand. However, the active sites do not seem to have been constructed to make the OAT reaction as exergonic as possible, but instead to keep the reaction free energy close to zero (to avoid excessive loss of energy), thereby making the reoxidation (SO and XO) or rereduction of the active sites (DMSOR) after the OAT reaction facile. We also show that active-site models of the three enzyme families can all catalyze the reduction of DMSO and that the DMSOR model does not give the lowest activation barrier. Likewise, all three models can catalyze the oxidation of sulfite, provided that the Coulombic repulsion between the substrate and the enzyme model can be overcome, but for this harder reaction, the SO model gives the lowest activation barrier, although the differences are not large. However, only the XO model can catalyze the oxidation of xanthine, owing to its sulfido ligand.

A quantum-mechanical study of the reaction mechanism of sulfite oxidase.

The oxidation of sulfite to sulfate by two different models of the active site of sulfite oxidase has been studied. Both protonated and deprotonated substrates were tested. Geometries were optimized with density functional theory (TPSS/def2-SV(P)) and energies were calculated either with hybrid functionals and large basis sets (B3LYP/def2-TZVPD) including corrections for dispersion, solvation, and entropy, or with coupled-cluster theory (LCCSD(T0)) extrapolated toward a complete basis set. Three suggested reaction mechanisms have been compared and the results show that the lowest barriers are obtained for a mechanism where the substrate attacks a Mo-bound oxo ligand, directly forming a Mo-bound sulfate complex, which then dissociates into the products. Such a mechanism is more favorable than mechanisms involving a Mo-sulfite complex with the substrate coordinating either by the S or O atom. The activation energy is dominated by the Coulomb repulsion between the Mo complex and the substrate, which both have a negative charge of -1 or -2.

Sulfur K-edge X-ray absorption spectroscopy and density functional theory calculations on monooxo Mo(IV) and bisoxo Mo(VI) bis-dithiolenes: insights into the mechanism of oxo transfer in sulfite oxidase and its relation to the mechanism of DMSO reductase.

Sulfur K-edge X-ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations have been used to determine the electronic structures of two complexes [Mo(IV)O(bdt)2](2-) and [Mo(VI)O2(bdt)2](2-) (bdt = benzene-1,2-dithiolate(2-)) that relate to the reduced and oxidized forms of sulfite oxidase (SO). These are compared with those of previously studied dimethyl sulfoxide reductase (DMSOr) models. DFT calculations supported by the data are extended to evaluate the reaction coordinate for oxo transfer to a phosphite ester substrate. Three possible transition states are found with the one at lowest energy, stabilized by a P-S interaction, in good agreement with experimental kinetics data. Comparison of both oxo transfer reactions shows that in DMSOr, where the oxo is transferred from the substrate to the metal ion, the oxo transfer induces electron transfer, while in SO, where the oxo transfer is from the metal site to the substrate, the electron transfer initiates oxo transfer. This difference in reactivity is related to the difference in frontier molecular orbitals (FMO) of the metal-oxo and substrate-oxo bonds. Finally, these experimentally related calculations are extended to oxo transfer by sulfite oxidase. The presence of only one dithiolene at the enzyme active site selectively activates the equatorial oxo for transfer, and allows facile structural reorganization during turnover.

Molybdenum trioxide nanoparticles with intrinsic sulfite oxidase activity.

Sulfite oxidase is a mitochondria-located molybdenum-containing enzyme catalyzing the oxidation of sulfite to sulfate in the amino acid and lipid metabolism. Therefore, it plays a major role in detoxification processes, where defects in the enzyme cause a severe infant disease leading to early death with no efficient or cost-effective therapy in sight. Here we report that molybdenum trioxide (MoO3) nanoparticles display an intrinsic biomimetic sulfite oxidase activity under physiological conditions, and, functionalized with a customized bifunctional ligand containing dopamine as anchor group and triphenylphosphonium ion as targeting agent, they selectively target the mitochondria while being highly dispersible in aqueous solutions. Chemically induced sulfite oxidase knockdown cells treated with MoO3 nanoparticles recovered their sulfite oxidase activity in vitro, which makes MoO3 nanoparticles a potential therapeutic for sulfite oxidase deficiency and opens new avenues for cost-effective therapies for gene-induced deficiencies.