The periplasmic hydrogenase of Desulfovibrio vulgaris (Hildenbourough) is an all Fe-containing hydrogenase. It contains two ferredoxin type [4Fe-4S] clusters, termed the F clusters, and a catalytic H cluster. Recent X-ray crystallographic studies on two Fe hydrogenases revealed that the H cluster is composed of two sub-clusters, a [4Fe-4S] cluster ([4Fe-4S]H) and-a binuclear Fe cluster ([2Fe]H), bridged by a cysteine sulfur. The aerobically purified D. vulgaris hydrogenase is stable in air. It is inactive and requires reductive activation. Upon reduction, the enzyme becomes sensitive to O(2) indicating that the reductive activation process is irreversible. Previous EPR investigations showed that upon reoxidation (under argon) the H cluster exhibits a rhombic EPR signal that is not seen in the as-purified enzyme, suggesting a conformational change in association with the reductive activation. For the purpose of gaining more information on the electronic properties of this unique H cluster and to understand further the reductive activation process, variable-temperature and variable-field Mossbauer spectroscopy has been used to characterize the Fe-S clusters in D. vulgaris hydrogenase poised at different redox states generated during a reductive titration, and in the GO-reacted enzyme. The data were successfully decomposed into spectral components corresponding to the F and H clusters,and characteristic parameters describing the electronic and magnetic properties of the F and H clusters were obtained. Consistent with the X-ray crystallographic results, the spectra of the H cluster can be understood as originating from an exchange coupled [4Fe-4S] - [2Fe] system. In particular, detailed analysis of the data reveals that the reductive activation begins with reduction of the [4Fe-4S]H cluster from the 2+ to the If state, followed by transfer of the reducing equivalent from the [4Fe-4S]H subcluster to the binuclear [2Fe]H subcluster. The results also reveal that binding of exogenous CO to the H cluster affects significantly the exchange coupling between the [4Fe-4S]H and the [2Fe]H subclusters. Implication of such a CO binding effect is discussed.

The outcome of O-2 activation at the diiron(II) cluster in the R2 subunit of Escherichia coli (class I) ribonucleotide reductase has been rationally altered from the normal tyrosyl radical (Y122)(1) production to self-hydroxylation of a phenylalanine side-chain by two amino acid substitutions that leave intact the (histidine)(2)-(carboxylate)(4) ligand set characteristic of the diiron-carboxylate family. Iron ligand Asp (D) 84 was replaced with Glu (E), the amino acid found in the cognate position of the structurally similar diiron-carboxylate protein, methane monooxygenase hydroxylase (MMOH). We previously showed that this substitution allows accumulation of a mu -1,2-peroxodiiron(III) intermediate,(2 3) which does not accumulate in the wild-type (wt) protein and is probably a structural homologue of intermediate P (H-peroxo) in O-2 activation by MMOH.(4) In addition, the near-surface residue Trp (W) 48 was replaced with Phe (F), blocking transfer of the ``extra'' electron that occurs in wt R2 during formation of the formally Fe(LII)Fe(IV) cluster X.(5-7) Decay of the mu1,2-peroxodiiron(III) complex in R2-W38F/D84E gives an initial brown product, which contains very little YI22(.) and which converts very slowly (t(1/2) similar to 7 h) upon incubation at 0 degreesC to an intensely purple final product. X-ray crystallographic analysis of the purple product indicates that F208 has undergone epsilon -hydroxylation and the resulting phenol has shifted significantly to become st ligand to Fe2 of the diiron cluster. Resonance Raman (RR) spectra of the purple product generated with O-16(2) or O-18(2) show appropriate isotopic sensitivity in bands assigned to O-phenyl and Fe-O-phenyl vibrational modes, confirming that the oxygen of the Fe(III)-phenolate species is derived from Or. Chemical analysis, experiments involving interception of the hydroxylating intermediate with exogenous reductant, and Mossbauer and EXAFS characterization of the brown and purple species establish that F208 hydroxylation occurs during decay of the peroxo complex and formation of the initial brown product. The slow transition to the purple Fe(LII)-phenolate species is ascribed to a ligand rearrangement in which mu -O2- is lost and the F208-derived phenolate coordinates. The reprogramming to F208 monooxygenase requires both amino acid substitutions, as very little epsilon -hydroxyphenylalanine is formed and pathways leading to Y122(.) formation predominate in both R2-D84E and R2-W48F(2-7).

Ferrochelatase (EC 4.99.1.1) is the terminal enzyme of the haem biosynthetic pathway and catalyses iron chelation into the protoporphyrin IX ring. Glutamate-287 (E287) of murine mature ferrochelatase is a conserved residue in all known sequences of ferrochelatase, is present at the active site of the enzyme, as inferred from the Bacillus subtilis ferrochelatase three-dimensional structure, and is critical for enzyme activity. Substitution of E287 with either glutamine (Q) or alanine (A) yielded variants with lower enzymic activity than that of the wild-type ferrochelatase and with different absorption spectra from the wild-type enzyme. In contrast to the wild-type enzyme, the absorption spectra of the variants indicate that these enzymes, as purified, contain protoporphyrin IX. Identification and quantification of the porphyrin bound to the E287-directed variants indicate that approx. 80% of the total porphyrin corresponds to protoporphyrin IX. Significantly, rapid stopped-flow experiments of the E287A and E287Q Variants demonstrate that reaction with Zn2+ results in the formation of bound Zn-protoporphyrin IX, indicating that the endogenously bound protoporphyrin IX can be used as a substrate. Taken together, these findings suggest that the structural strain imposed by ferrochelatase on the porphyrin substrate as a critical step in the enzyme catalytic mechanism is also accomplished by the E287A and E287Q variants, but without the release of the product. Thus E287 in murine ferrochelatase appears to be critical For the catalytic process by controlling the release of the product.

Mycobacterium tuberculosis KatG is a multifunctional heme enzyme responsible for activation of the antibiotic isoniazid. A KatG(S315T) point mutation is found in >50% of isoniazid-resistant clinical isolates. Since isoniazid activation is thought to involve an oxidation reaction, the redox potential of KatG was determined using cyclic voltammetry, square wave voltammetry, and spectroelectrochemical titrations. Isoniazid activation may proceed via a cytochrome P450-like mechanism. Therefore, the possibility that substrate binding by KatG leads to an increase in the heme redox potential and the possibility that KatG(S315T) confers isoniazid resistance by altering the redox potential were examined. Effects of the heme spin state on the reduction potentials of KatG and KatG(S315T) were also determined. Assessment of the Fe(3+)/Fe(2+) couple gave a midpoint potential of ca. -50 mV for both KatG and KatG(S315T). In contrast to cytochrome P450s, addition of substrate had no significant effect on either the KatG or KatG(S315T) redox potential. Conversion of the heme to a low-spin configuration resulted in a -150 to -200 mV shift of the KatG and KatG(S315T) redox potentials. These results suggest that isoniazid resistance conferred by KatG(S315T) is not mediated through changes in the heme redox potential. The redox potentials of isoniazid were also determined using cyclic and square wave voltammetry, and the results provide evidence that the ferric KatG and KatG(S315T) midpoint potentials are too low to promote isoniazid oxidation without formation of a high-valent enzyme intermediate such as compounds I and II or oxyferrous KatG.

The combination of docking algorithms with NMR data has been developed extensively for the studies of protein-ligand interactions. However, to extend this development for the studies of protein-protein interactions, the intermolecular NOE constraints, which are needed, are more difficult to access. In the present work, we describe a new approach that combines an ab initio docking calculation and the mapping of an interaction site using chemical shift variation analysis. The cytochrome c553-ferredoxin complex is used as a model of numerous electron-transfer complexes. The 15N-labeling of both molecules has been obtained, and the mapping of the interacting site on each partner, respectively, has been done using HSQC experiments. 1H and 15N chemical shift analysis defines the area of both molecules involved in the recognition interface. Models of the complex were generated by an ab initio docking software, the BiGGER program (bimolecular complex generation with global evaluation and ranking). This program generates a population of protein-protein docked geometries ranked by a scoring function, combining relevant stabilization parameters such as geometric complementarity surfaces, electrostatic interactions, desolvation energy, and pairwise affinities of amino acid side chains. We have implemented a new module that includes experimental input (here, NMR mapping of the interacting site) as a filter to select the accurate models. Final structures were energy minimized using the X-PLOR software and then analyzed. The best solution has an interface area (1037.4 A2) falling close to the range of generally observed recognition interfaces, with a distance of 10.0 A between the redox centers.

Nitrous oxide (N20) is a greenhouse gas, the third most significant contributor to global warming. As a key process for N20 elimination from the biosphere, N20 reductases catalyze the two-electron reduction of N20 to N2. These 2 x 65 kDa copper enzymes are thought to contain a CuA electron entry site, similar to that of cytochrome c oxidase, and a CuZ catalytic center. The copper anomalous signal was used to solve the crystal structure of N20 reductase from Pseudomonas nautica by multiwavelength anomalous dispersion, to a resolution of 2.4 A. The structure reveals that the CuZ center belongs to a new type of metal cluster, in which four copper ions are liganded by seven histidine residues. N20 binds to this center via a single copper ion. The remaining copper ions might act as an electron reservoir, assuring a fast electron transfer and avoiding the formation of dead-end products.

A new computationally efficient and automated "soft docking" algorithm is described to assist the prediction of the mode of binding between two proteins, using the three-dimensional structures of the unbound molecules. The method is implemented in a software package called BiGGER (Bimolecular Complex Generation with Global Evaluation and Ranking) and works in two sequential steps: first, the complete 6-dimensional binding spaces of both molecules is systematically searched. A population of candidate protein-protein docked geometries is thus generated and selected on the basis of the geometric complementarity and amino acid pairwise affinities between the two molecular surfaces. Most of the conformational changes observed during protein association are treated in an implicit way and test results are equally satisfactory, regardless of starting from the bound or the unbound forms of known structures of the interacting proteins. In contrast to other methods, the entire molecular surfaces are searched during the simulation, using absolutely no additional information regarding the binding sites. In a second step, an interaction scoring function is used to rank the putative docked structures. The function incorporates interaction terms that are thought to be relevant to the stabilization of protein complexes. These include: geometric complementarity of the surfaces, explicit electrostatic interactions, desolvation energy, and pairwise propensities of the amino acid side chains to contact across the molecular interface. The relative functional contribution of each of these interaction terms to the global scoring function has been empirically adjusted through a neural network optimizer using a learning set of 25 protein-protein complexes of known crystallographic structures. In 22 out of 25 protein-protein complexes tested, near-native docked geometries were found with C(alpha) RMS deviations < or =4.0 A from the experimental structures, of which 14 were found within the 20 top ranking solutions. The program works on widely available personal computers and takes 2 to 8 hours of CPU time to run any of the docking tests herein presented. Finally, the value and limitations of the method for the study of macromolecular interactions, not yet revealed by experimental techniques, are discussed.

Fe-hydrogenase is a 54-kDa iron-sulfur enzyme essential for hydrogen cycling in sulfate-reducing bacteria. The x-ray structure of Desulfovibrio desulfuricans Fe-hydrogenase has recently been solved, but structural information on the recognition of its redox partners is essential to understand the structure-function relationships of the enzyme. In the present work, we have obtained a structural model of the complex of Fe-hydrogenase with its redox partner, the cytochrome c(553), combining docking calculations and NMR experiments. The putative models of the complex demonstrate that the small subunit of the hydrogenase has an important role in the complex formation with the redox partner; 50% of the interacting site on the hydrogenase involves the small subunit. The closest contact between the redox centers is observed between Cys-38, a ligand of the distal cluster of the hydrogenase and Cys-10, a ligand of the heme in the cytochrome. The electron pathway from the distal cluster of the Fe-hydrogenase to the heme of cytochrome c(553) was investigated using the software Greenpath and indicates that the observed cysteine/cysteine contact has an essential role. The spatial arrangement of the residues on the interface of the complex is very similar to that already described in the ferredoxin-cytochrome c(553) complex, which therefore, is a very good model for the interacting domain of the Fe-hydrogenase-cytochrome c(553).

Iron-sulfur clusters with [3Fe-4S] cores are widely distributed in biological systems. In the oxidized state, designated [3Fe-4S](+), these electron-transfer agents have an electronic ground state with S = 1/2, and; they exhibit EPR signals centered at g = 2.01. It has been established by Mossbauer spectroscopy that the three iron sites of the cluster are high-spin Fe3+; and the general properties of the S = 1/2 ground state have been described with the exchange Hamiltonian H-exch = J(12)S(1).S-2 + J(23)S(2).S-3 + J(13)S(1).S-3 Some [3Fe-4S](+) clusters (type 1) have their g-values confined to the range between g = 2.03 and 2.00 while others (type 2) exhibit a continuous distribution of g-values down to g approximate to 1.85. Despite considerable efforts in various laboratories no model has emerged that explains the g-values of type 2 clusters. The 4.2 K spectra of all [3Fe-4S](+) clusters have broad features,which have been simulated in the past by using Fe-57 magnetic hyperfine tensors with anisotropies that are unusually large for high-spin feme sites. It is proposed here that antisymmetric exchange, H-AS = d.(S-1 x S-2 + S-2 x S-3 + S-3 x S-1), is the cause of the g-value shifts in type 2 clusters. We have been able to fit the EPR and Mossbauer spectra of the 3Fe clusters of beef heart aconitase and Desulfovibrio gigas ferredoxin II by using antisymmetric exchange in combination with distributed exchange coupling constants J(12), J(13), and J(23) (J-strain). While antisymmetric exchange is negligible for aconitase (which has a type 1 cluster), fits of the ferredoxin II spectra require \d\ approximate to 0.4 cm(-1). Our studies show that the data of both proteins can lie fit using the same isotropic Fe-57 magnetic hyperfine coupling constant for th three cluster sites, namely a -18.0 MHz for aconitase and a = -18.5 MHz for the D. gigas ferredoxin. The effects of antisymmetric exchange and J-strain on the Mossbauer and EPR spectra are discussed.

The aerobic purification of Pseudomonas nautica 617 nitrous oxide reductase yielded two forms of the enzyme exhibiting different chromatographic behaviors. The protein contains six copper atoms per monomer, arranged in two centers named Cu(A) and Cu(Z). Cu(Z) could be neither oxidized nor further reduced under our experimental conditions, and exhibits a 4-line EPR spectrum (g(x)=2.015, A(x)=1.5 mT, g(y)=2.071, A(y)=2 mT, g(z)=2.138, A(z)=7 mT) and a strong absorption at approximately 640 nm. Cu(A) can be stabilized in a reduced EPR-silent state and in an oxidized state with a typical 7-line EPR spectrum (g(x)=g(y)= 2.021, A(x) = A(y)=0 mT, g(z) = 2.178, A(z)= 4 mT) and absorption bands at 480, 540, and approximately 800 nm. The difference between the two purified forms of nitrous oxide reductase is interpreted as a difference in the oxidation state of the Cu(A) center. In form A, Cu(A) is predominantly oxidized (S = (1)/(2), Cu(1.5+)-Cu(1.5+)), while in form B it is mostly in the one-electron reduced state (S = 0, Cu(1+)-Cu(1+)). In both forms, Cu(Z) remains reduced (S = 1/2). Complete crystallographic data at 2.4 A indicate that Cu(A) is a binuclear site (similar to the site found in cytochrome c oxidase) and Cu(Z) is a novel tetracopper cluster [Brown, K., et al. (2000) Nat. Struct. Biol. (in press)]. The complete amino acid sequence of the enzyme was determined and comparisons made with sequences of other nitrous oxide reductases, emphasizing the coordination of the centers. A 10.3 kDa peptide copurified with both forms of nitrous oxide reductase shows strong homology with proteins of the heat-shock GroES chaperonin family.

Nitrous oxide (N(2)O) is a greenhouse gas, the third most significant contributor to global warming. As a key process for N(2)O elimination from the biosphere, N(2)O reductases catalyze the two-electron reduction of N(2)O to N(2). These 2 x 65 kDa copper enzymes are thought to contain a CuA electron entry site, similar to that of cytochrome c oxidase, and a CuZ catalytic center. The copper anomalous signal was used to solve the crystal structure of N(2)O reductase from Pseudomonas nautica by multiwavelength anomalous dispersion, to a resolution of 2.4 Angstrom. The structure reveals that the CuZ center belongs to a new type of metal cluster, in which four copper ions are liganded by seven histidine residues. N(2)O binds to this center via a single copper ion. The remaining copper ions might act as an electron reservoir, assuring a fast electron transfer and avoiding the formation of dead-end products.

Coordination/detachment of a pendent functionality in the Zn(ii) complex with a macrocyclic ligand L gives rise to on/off switching of exciplex emission, defining an elementary molecular machine whose movements are driven by both pH and light.