Density functional theory studies of model complexes for molybdenum-dependent nitrate reductase active sites

Molybdenum and tungsten complexes as models for the active sites of assimilatory or dissimilatory nitrate reductases (NR) were computed at the CPCM-B98/SDDp//B3LYP/Lanl2DZp* plus zero point energy level of density functional theory. The ligands were chosen on the basis of available experimental protein or small chemical model structures. A water molecule is found to bind to assimilatory NR models [(Me₂C₂S₂)MO(YMe)]⁻ (−11.5 kcal mol⁻¹ for M is Mo, Y is S) and may be replaced by nitrate (−4.5 kcal mol⁻¹) (but a hydroxy group may not). Nature’s choice of M is Mo and Y is S for NR has the largest activation energy for protein-free models (13.3 kcal mol⁻¹) and the least exothermic reaction energy for the nitrate reduction (−14.9 kcal mol⁻¹) compared with M is W and Y is O or Se alternatives. Water binding to dissimilatory NR model complexes [(Me₂C₂S₂)₂M(YR)]⁻ is considerably endothermic (10.3 kcal mol⁻¹); nitrate binding is only slightly so (1.5 kcal mol⁻¹ for RY⁻ is MeS⁻). The exchange of an oxo ligand (assimilatory NR) for a dithiolato ligand (dissimilatory NR model) reduces the exothermicity (−8.6 kcal mol⁻¹ relative to the fivefold-coordinate reduced complex) and raises the barrier for oxygen atom transfer (OAT) in the nitrate complex (19.2 kcal mol⁻¹). Not for the mono but only for the bisdithiolato complexes hydrogen bonding involving the coordinated substrate may significantly lower the OAT barrier as shown by explicitly adding water molecules. Substitution of tungsten for molybdenum generally lowers OAT activation energies and makes nitrate reduction reaction energies more negative. Bidentate carboxylato binding identified in Escherichia coli NarGHI is the preferred binding mode also for an acetato model. However, one dithiolato ligand folds when the Mo^VI center is bare of a good π-donor ligand, e.g., an oxo group. Computations on [(mnt)₂Mo^IV(YR)(PPh₃)]⁻ [mnt is (CN)₂C₂S₂ ²⁻] gave a smaller nitrate reduction activation energy for RY⁻ is Cl⁻, compared with RY⁻ is PhS⁻, although experimentally only the phenyl thiolato complex and not the chloro complex was found to be a functional NR model. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

EPR and redox properties of periplasmic nitrate reductase from Desulfovibrio desulfuricans ATCC 27774

Nitrate reductases are enzymes that catalyze the conversion of nitrate to nitrite. We report here electron paramagnetic resonance (EPR) studies in the periplasmic nitrate reductase isolated from the sulfate-reducing bacteria Desulfovibrio desulfuricans ATCC 27774. This protein, belonging to the dimethyl sulfoxide reductase family of mononuclear Mo-containing enzymes, comprises a single 80-kDa subunit and contains a Mo bis(molybdopterin guanosine dinucleotide) cofactor and a [4Fe–4S] cluster. EPR-monitored redox titrations, carried out with and without nitrate in the potential range from 200 to −500 mV, and EPR studies of the enzyme, in both catalytic and inhibited conditions, reveal distinct types of Mo(V) EPR-active species, which indicates that the Mo site presents high coordination flexibility. These studies show that nitrate modulates the redox properties of the Mo active site, but not those of the [4Fe–4S] center. The possible structures and the role in catalysis of the distinct Mo(V) species detected by EPR are discussed.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Secondary structure analysis of the dissimilatory sulphite reductase in Desulfovibrio desulfuricans

The complete sequences of the dsrA and dsrB genes coding for the alpha- and beta-subunits, respectively, of the sulphite reductase enzyme in Desulfovibrio desulfuricans were determined. Analyses of the amino acid sequences indicated a number of serohaem/Fe4S4 binding consensus sequences whilst predictive secondary structure analysis revealed a similar pattern of alpha-helix and beta-strand structures between the two subunits which was indicative of gene duplication.     

Isolation and preliminary characterization of a soluble nitrate reductase from the sulfate reducing organism Desulfovibrio desulfuricans ATCC 27774

Desulfovibrio desulfuricans ATCC 27774 is a sulfate reducer that can adapt to nitrate respiration, inducing the enzymes required to utilize this alternative metabolic pathway. Nitrite reductase from this organism has been previously isolated and characterized, but no information was available on the enzyme involved in the reduction of nitrate. This is the first report of purification to homogeneity of a nitrate reductase from a sulfate reducing organism, thus completing the enzymatic system required to convert nitrate (through nitrite) to ammonia. D. desulfuricans nitrate reductase is a monomeric (circa 70 kDa) periplasmic enzyme with a specific activity of 5.4 K(m) for nitrate was estimated to be 20 microM. EPR signals due to one [4Fe-4S] cluster and Mo(V) were identified in dithionite reduced samples and in the presence of nitrate.     

The nature of the di-iron site in the bacterioferritin from Desulfovibrio desulfuricans

The first crystal structure of a native di-iron center in an iron-storage protein (bacterio)ferritin is reported. The protein, isolated from the anaerobic bacterium Desulfovibrio desulfuricans, has the unique property of having Fe-coproporphyrin III as its heme cofactor. The three-dimensional structure of this bacterioferritin was determined in three distinct catalytic/redox states by X-ray crystallography (at 1.95, 2.05 and 2.35 A resolution), corresponding to different intermediates of the di-iron ferroxidase site. Conformational changes associated with these intermediates support the idea of a route for iron entry into the protein shell through a pore that passes through the di-iron center. Molecular surface and electrostatic potential calculations also suggest the presence of another ion channel, distant from the channels at the three- and four-fold axes proposed as points of entry for the iron atoms.     

A revised model of the active site of alternative oxidase

The plant mitochondrial protein alternative oxidase catalyses dioxygen dependent ubiquinol oxidation to yield ubiquinone and water. A structure of this protein has previously been proposed based on an assumed structural homology to the di-iron carboxylate family of proteins. However, these authors suggested the protein has a very different topology than the known structures of di-iron carboxylate proteins. We have re-examined this model and based on comparison of recent sequences and structural data on di-iron carboxylate proteins we present a new model of the alternative oxidase which allows prediction of active site residues and a possible membrane binding motif.     

Structural basis of eukaryotic nitrate reduction: crystal structures of the nitrate reductase active site

Nitrate assimilation in autotrophs provides most of the reduced nitrogen on earth. In eukaryotes, reduction of nitrate to nitrite is catalyzed by the molybdenum-containing NAD(P)H:nitrate reductase (NR; EC 1.7.1.1-3). In addition to the molybdenum center, NR contains iron-heme and flavin adenine dinucleotide as redox cofactors involved in an internal electron transport chain from NAD(P)H to nitrate. Recombinant, catalytically active Pichia angusta nitrate-reducing, molybdenum-containing fragment (NR-Mo) was expressed in P. pastoris and purified. Crystal structures for NR-Mo were determined at 1.7 and 2.6 angstroms. These structures revealed a unique slot for binding nitrate in the active site and identified key Arg and Trp residues potentially involved in nitrate binding. Dimeric NR-Mo is similar in overall structure to sulfite oxidases, with significant differences in the active site. Sulfate bound in the active site caused conformational changes, as compared with the unbound enzyme. Four ordered water molecules located in close proximity to Mo define a nitrate binding site, a penta-coordinated reaction intermediate, and product release. Because yeast NAD(P)H:NR is representative of the family of eukaryotic NR, we propose a general mechanism for nitrate reduction catalysis.     

Density functional theory investigation of cocaine water complexes

Twenty cocaine–water complexes were studied using density functional theory (DFT) B3LYP/6-311++G** level to understand their geometries, energies, vibrational frequencies, charge transfer and topological parameters. Among the 20 complexes, 12 are neutral and eight are protonated in the cocaine-water complexes. Based on the interaction energy, the protonated complexes are more stable than the neutral complexes. In both complexes, the most stable structure involves the hydrogen bond with water at nitrogen atom in the tropane ring and C?=?O groups in methyl ester. Carbonyl groups in benzoyl and methyl ester is the most reactive site in both forms and it is responsible for the stability order. The calculated topological results show that the interactions involved in the hydrogen bond are electrostatic dominant. Natural bond orbital (NBO) analysis confirms the presence of hydrogen bond and it supports the stability order. Atoms in molecules (AIM) and NBO analysis confirms the C-H?·?·?·?O hydrogen bonds formed between the cocaine-water complexes are blue shifted in nature.     

Crystallographic study of rubredoxin from the bacterium Desulfovibrio desulfuricans strain 27774

Rubredoxin from Desulfovibrio desulfuricans (strain 27774) has been isolated and crystallized. Preliminary amino acid and crystallographic analyses indicate that this rubredoxin is the smallest rubredoxin isolated so far. The amino acid analysis indicates that the molecule is composed of 45 to 48 residues and contains histidine, which is unusual for rubredoxins from anaerobic bacteria. The X-ray diffraction pattern from these crystals reveals they belong to space group P1 with cell parameters: a = 24.92 A, b = 17.79 A, c = 19.72 A, alpha = 101.0 degrees, beta = 83.4 degrees, gamma = 104.5 degrees. The unit cell volume of 8283 A3 indicates a molecule with molecular weight no greater than 5500 and is consistent with the smaller number of amino acids found in this rubredoxin. The solvent content of this rubredoxin crystal appears to be the lowest observed in crystalline proteins.     

Heme-Cu complexes as oxygen-activating functional models for the active site of cytochrome c oxidase

Tri(2-pyridylmethyl)amineCu complex-linked iron meso-tetraphenylporphyine derivatives were prepared to model the active site of cytochrome c oxidase. Exposure to oxygen converted the reduced forms of the complexes to the corresponding stable mu-peroxo species in spite of the presence of three coordination sites, two on the heme and one on the Cu. The oxy forms were characterized spectroscopically. Kinetic analyses of the oxygenation reactions of the reduced forms suggests that preferential O2 binding occurs at the Cu site over the heme. This mechanism is also supported by examination of the redox potentials of the two metal ions. Since the peroxy complexes of the models exhibit a structure similar to that of the previously reported fully-oxidized form, the relevance of the model chemistry to the enzyme reaction is discussed.     

A new type of metal-binding site in cobalt- and zinc-containing adenylate kinases isolated from sulfate-reducers Desulfovibrio gigas and Desulfovibrio desulfuricans ATCC 27774

Adenylate kinase (AK) mediates the reversible transfer of phosphate groups between the adenylate nucleotides and contributes to the maintenance of their constant cellular level, necessary for energy metabolism and nucleic acid synthesis. The AK were purified from crude extracts of two sulfate-reducing bacteria (SRB), Desulfovibrio (D.) gigas NCIB 9332 and Desulfovibrio desulfuricans ATCC 27774, and biochemically and spectroscopically characterised in the native and fully cobalt- or zinc-substituted forms. These are the first reported adenylate kinases that bind either zinc or cobalt and are related to the subgroup of metal-containing AK found, in most cases, in Gram-positive bacteria. The electronic absorption spectrum is consistent with tetrahedral coordinated cobalt, predominantly via sulfur ligands, and is supported by EPR. The involvement of three cysteines in cobalt or zinc coordination was confirmed by chemical methods. Extended X-ray absorption fine structure (EXAFS) indicate that cobalt or zinc are bound by three cysteine residues and one histidine in the metal-binding site of the “LID” domain. The sequence ¹²⁹Cys-X₅-His-X₁₅-Cys-X₂-Cys of the AK from D. gigas is involved in metal coordination and represents a new type of binding motif that differs from other known zinc-binding sites of AK. Cobalt and zinc play a structural role in stabilizing the LID domain.     

Evidence for nickel and a three-iron center in the hydrogenase of Desulfovibrio desulfuricans

Hydrogenase from Desulfovibrio desulfuricans (ATCC No. 27774) grown in unenriched and in enriched 61Ni and 57Fe media has been purified to apparent homogeneity. Two fractions of enzymes with hydrogenase activity were separated and were termed hydrogenase I and hydrogenase II. they were shown to have similar molecular weights (77,600 for hydrogenase I and 75,500 for hydrogenase II), to be composed of two polypeptide chains, and to contain Ni and non-heme iron. Because of its higher specific activity (152 versus 97) hydrogenase II was selected for EPR and M?ssbauer studies. As isolated, hydrogenase II exhibits an "isotropic" EPR signal at g = 2.02 and a rhombic EPR signal at g = 2.3, 2.2, and 2.0. Isotopic substitution of 61Ni proves that the rhombic signal is due to Ni. Combining the M?ssbauer and EPR data, the isotropic g = 2.02 EPR signal was shown to originate from a 3Fe cluster which may have oxygenous or nitrogenous ligands. In addition, the M?ssbauer data also revealed two [4Fe-4S]2+ clusters iun each molecule of hydrogenase II. The EPR and M?ssbauer data of hydrogenase I were found to be identical to those of hydrogenase II, indicating that both enzymes have common metallic centers.     

Density functional theory study of the mechanism for the formation of glycidyl esters from triglyceride

Glycidyl fatty acid esters (GEs) are by-products of edible oil refinement that have attracted attention globally due to concerns over their possible harmful effects on human health when consumed. It is thus important to improve our understanding of GE formation if we are to suppress GE production during edible oil refinement. In this paper, a first-principles density functional theory study of the formation mechanism of GEs was performed. Triglycerides undergo a self-condensation reaction between two adjacent ester groups to yield GEs and an anhydride as a by-product. This process is energetically unfavorable, having a relatively high activation energy of around 80 kcal/mol, which indicates that GE formation is intrinsically a high-temperature process. Both the thermodynamic and the kinetic energies of the reaction are insensitive to the size of the fatty chain substituents present. If water participates in the self-condensation, the activation barrier is notably decreased by 23.9 kcal/mol, indicating that GE production in the presence of high-temperature water vapor should be more kinetically favorable. Our results suggest that reducing the reaction temperature and avoiding the use of water should suppress GE production during edible oil refinement.     

Access to the active site of periplasmic nitrate reductase: insights from site-directed mutagenesis and zinc inhibition studies

The periplasmic nitrate reductase (NapAB), a member of the DMSO reductase superfamily, catalyzes the first step of the denitrification process in bacteria. In this heterodimer, a di-heme NapB subunit is associated to the catalytic NapA subunit that binds a [4Fe-4S] cluster and a bis(molybdopterin guanine dinucleotide) cofactor. Here, we report the kinetic characterization of purified mutated heterodimers from Rhodobacter sphaeroides. By combining site-directed mutagenesis, redox potentiometry, EPR spectroscopy, and enzymatic characterization, we investigate the catalytic role of two conserved residues (M153 and R392) located in the vicinity of the molybdenum active site. We demonstrate that M153 and R392 are involved in nitrate binding: the Vm measured on the M153A and R392A mutants are similar to that measured on the wild-type enzyme, whereas the Km for nitrate is increased 10-fold and 200-fold, respectively. The use of an alternative enzymatic assay led us to discover that NapAB is uncompetitively inhibited by Zn2+ ions (Ki' = 1 microM). We used this property to further probe the active site access in the mutant enzymes. It is proposed that R392 acts as a filter by preventing a direct reduction of the Mo atom by small reducing molecules and partially protecting the active site against zinc inhibition. In addition, we show that M153 is a key residue mediating this inhibition likely by coordinating Zn2+ ions via its sulfur atom. This residue is not conserved in the DMSO reductase superfamily while it is conserved in the periplasmic nitrate reductase family. Zinc inhibition is therefore likely to be specific and restricted to periplasmic nitrate reductases.     

Mercury methylation by Desulfovibrio desulfuricans ND132 in the presence of polysulfides

The extracellular speciation of mercury may control bacterial uptake and methylation. Mercury-polysulfide complexes have recently been shown to be prevalent in sulfidic waters containing zero-valent sulfur. Despite substantial increases in total dissolved mercury concentration, methylation rates in cultures of Desulfovibrio desulfuricans ND132 equilibrated with cinnabar did not increase in the presence of polysulfides, as expected due to the large size and charged nature of most of the complexes. In natural waters not at saturation with cinnabar, mercury-polysulfide complexes would be expected to shift the speciation of mercury from HgS(0)((aq)) toward charged complexes, thereby decreasing methylation rates.     

IR spectroelectrochemical study of the binding of carbon monoxide to the active site of Desulfovibrio fructosovorans Ni-Fe hydrogenase

The binding of carbon monoxide, a competitive inhibitor of many hydrogenases, to the active site of Desulfovibrio fructosovorans hydrogenase has been studied by infrared spectroscopy in a spectroelectrochemical cell. Direct evidence has been obtained of which redox states of the enzyme can bind extrinsic CO. Redox states A, B and SU do not bind extrinsic CO; only after reductive activation of the hydrogenase can CO bind to the active site. Two states with bound extrinsic CO can be distinguished by FTIR. These two states are in redox equilibrium and are most probably due to different oxidation states of the proximal 4Fe-4S cluster. Vibrational frequencies and theoretical quantum mechanics studies (DFT) of this process preclude the possibility of strong bonding of extrinsic CO to the Fe or Ni atoms of the active site. We propose that CO inhibition is caused by weak interaction of the extrinsic ligand with the Ni atom, blocking electron and proton transfer at the active site. A calculated structure with a weakly bound extrinsic CO at Ni has relative CO frequencies in excellent agreement with the experimental ones.     

Mercury Methylation by Desulfovibrio desulfuricans ND132 in the Presence of Polysulfides

The extracellular speciation of mercury may control bacterial uptake and methylation. Mercury-polysulfide complexes have recently been shown to be prevalent in sulfidic waters containing zero-valent sulfur. Despite substantial increases in total dissolved mercury concentration, methylation rates in cultures of Desulfovibrio desulfuricans ND132 equilibrated with cinnabar did not increase in the presence of polysulfides, as expected due to the large size and charged nature of most of the complexes. In natural waters not at saturation with cinnabar, mercury-polysulfide complexes would be expected to shift the speciation of mercury from HgS⁰_(aq) toward charged complexes, thereby decreasing methylation rates.     

The active site of the [FeFe]-hydrogenase from Desulfovibrio desulfuricans. I. Light sensitivity and magnetic hyperfine interactions as observed by electron paramagnetic resonance

The hydrogen-activating cluster (H cluster) in [FeFe]-hydrogenases consists of two moieties. The [2Fe]_H subcluster is a (L)(CO)(CN)Fe(μ-RS₂)(μ-CO)Fe(CysS)(CO)(CN) centre. The Cys-bound Fe is called Fe1, the other iron Fe2. The Cys-thiol forms a bridge to a [4Fe–4S] cluster, the [4Fe–4S]_H subcluster. We report that electron paramagnetic resonance (EPR) spectra of the ⁵⁷Fe-enriched enzyme from Desulfovibrio desulfuricans in the H_ox–CO state are consistent with a magnetic hyperfine interaction of the unpaired spin with all six Fe atoms of the H cluster. In contrast to the inactive aerobic enzyme, the active enzyme is easily destroyed by light. The [2Fe]_H subcluster in some enzyme molecules loses CO by photolysis, whereupon other molecules firmly bind the released CO to form the H_ox–CO state giving rise to the so-called axial 2.06 EPR signal. Though not destroyed by light, the H_ox–CO state is affected by it. As demonstrated in the accompanying paper [49] two of the intrinsic COs, both bound to Fe2, can be exchanged by extrinsic ¹³CO during illumination at 2 °C. We found that only one of the three ¹³COs, the one at the extrinsic position, gives an EPR-detectable isotropic superhyperfine interaction of 0.6 mT. At 30 K both the inhibiting extrinsic CO bound to Fe2 and one more CO can be photolysed. EPR spectra of the photolysed products are consistent with a 3d ⁷ system of Fe with the formal oxidation state +1. The damaged enzyme shows a light-sensitive g=5 signal which is ascribed to an S=3/2 form of the [2Fe]_H subcluster. The light sensitivity of the enzyme explains the occurrence of the g=5 signal and the axial 2.06 signal in published EPR spectra of nearly all preparations studied thus far.