A c-type hexaheme nitrite reductase (NiR) isolated from nitrate-grown cells of Desulfovibrio desulfuricans (Dd) ATCC 27774 catalyses the six-electron reduction of nitrite to ammonia. Previous electrochemical studies demonstrated that a simple electrocatalytic mechanism can be applied to this system (Moreno, C., Costa, C., Moura, I., LeGall, J., Liu, M. Y., Payne, W. J., Van Dijk, C. and Moura, J. J. G. (1992) Eur.J.Biochem. 212, 79-86). Its substrate specificity, availability and stability under ambient conditions makes this enzymatic system a promising candidate for use in a biosensor device. An electrochemical study of gel-immobilized Dd NiR on a glassy carbon electrode revealed both enzymatic activity and amperometric response to nitrite. In this study it was observed that the catalytic current density is a function of the nitrite concentration in solution and follows a characteristic Michaelis-Menten-type substrate dependence. Such a biosensor device (NiR-electrode) bears the option to be used for analytical determination of nitrite in complex media.

The di-haem cytochrome c peroxidase of Paracoccus denitrificans is a calcium binding dimer of 37.5 kDa subunits. It is responsible for reduction of H(2)O(2) to H(2)O with oxidation of cytochrome c(550) and is isolated in a fully oxidised state (inactive) in which one haem (centre I) is in a high-spin/low-spin equilibrium and high potential and the other (centre II) is low-spin and low potential. The enzyme undergoes direct electron transfer (without the need for mediators) with a 4,4'-dithiodipyridine-modified gold electrode and the response of both haem groups can be observed. By combination of the cyclic and pulse voltammetric data with the established spectroscopic information, it was demonstrated that entry of one electron to the high potential haem leads (in a mechanism involving strong haem-haem interactions) to a complex change of spin states and redox potentials of both haems in order to attain a "ready state" for binding, reduction and cleavage of the hydrogen peroxide. In the absence of endogenous calcium, haem communication can be completely disconnected and is recovered only when Ca(2+) is added, an essential step for the formation of the peroxidatic site. The intricate electrochemical behaviour of this enzyme was interpreted as a mechanism involving, both reduction and oxidation of the high potential haem, an interfacial electron transfer coupled to a homogenous chemical reaction (EC mechanism). We discuss two different models for the sequence of events leading to the appearance of the active pentacoordinated peroxidatic haem.

Desulfovibrio gigas hydrogenase (EC 1.12.2.1) is a complex enzyme containing one nickel, one 3Fe, and two [Fe4S4] clusters (Teixeira, M., Moura, I., Xavier, A. V., Der Vartanian, D. V., LeGall, J., Peck, H. D., Jr., Huynh, B. H., and Moura, J. J. G. (1983) Eur. J. Biochem. 130, 481-484). This hydrogenase belongs to a class of enzymes that are inactive "as isolated" (the so-called "oxygen-stable hydrogenases") and must go through an activation process in order to express full activity. The state of characterization of the active centers of the enzyme as isolated prompted us to do a detailed analysis of the redox patterns, activation profile, and catalytic redox cycle of the enzyme in the presence of either the natural substrate (H2) or chemical reductants. The effect of natural cofactors, as cytochrome C3, was also studied. Special focus was given to the intermediate redox species generated during the catalytic cycle of the enzyme and to the midpoint redox potentials associated. The available information is discussed in terms of a "working hypothesis" for the mechanism of the [NiFe] hydrogenases from sulfate reducing organisms in the context of activation process and catalytic cycle.

The multicopper enzyme nitrous oxide reductase (N 2OR) catalyzes the final step of denitrification, the two-electron reduction of N 2O to N 2. This enzyme is a functional homodimer containing two different multicopper sites: CuA and CuZ. CuA is a binuclear copper site that transfers electrons to the tetranuclear copper sulfide CuZ, the catalytic site. In this study, Pseudomonas nautica cytochrome c 552 was identified as the physiological electron donor. The kinetic data show differences when physiological and artificial electron donors are compared [cytochrome vs methylviologen (MV)]. In the presence of cytochrome c 552, the reaction rate is dependent on the ET reaction and independent of the N 2O concentration. With MV, electron donation is faster than substrate reduction. From the study of cytochrome c 552 concentration dependence, we estimate the following kinetic parameters: K m c 552 = 50.2 +/- 9.0 muM and V max c 552 = 1.8 +/- 0.6 units/mg. The N 2O concentration dependence indicates a K mN 2 O of 14.0 +/- 2.9 muM using MV as the electron donor. The pH effect on the kinetic parameters is different when MV or cytochrome c 552 is used as the electron donor (p K a = 6.6 or 8.3, respectively). The kinetic study also revealed the hydrophobic nature of the interaction, and direct electron transfer studies showed that CuA is the center that receives electrons from the physiological electron donor. The formation of the electron transfer complex was observed by (1)H NMR protein-protein titrations and was modeled with a molecular docking program (BiGGER). The proposed docked complexes corroborated the ET studies giving a large number of solutions in which cytochrome c 552 is placed near a hydrophobic patch located around the CuA center.

Identifying redox partners and the interaction surfaces is crucial for fully understanding electron flow in a respiratory chain. In this study, we focused on the interaction of nitrous oxide reductase (N(2)OR), which catalyzes the final step in bacterial denitrification, with its physiological electron donor, either a c-type cytochrome or a type 1 copper protein. The comparison between the interaction of N(2)OR from three different microorganisms, Pseudomonas nautica, Paracoccus denitrificans, and Achromobacter cycloclastes, with their physiological electron donors was performed through the analysis of the primary sequence alignment, electrostatic surface, and molecular docking simulations, using the bimolecular complex generation with global evaluation and ranking algorithm. The docking results were analyzed taking into account the experimental data, since the interaction is suggested to have either a hydrophobic nature, in the case of P. nautica N(2)OR, or an electrostatic nature, in the case of P. denitrificans N(2)OR and A. cycloclastes N(2)OR. A set of well-conserved residues on the N(2)OR surface were identified as being part of the electron transfer pathway from the redox partner to N(2)OR (Ala495, Asp519, Val524, His566 and Leu568 numbered according to the P. nautica N(2)OR sequence). Moreover, we built a model for Wolinella succinogenes N(2)OR, an enzyme that has an additional c-type-heme-containing domain. The structures of the N(2)OR domain and the c-type-heme-containing domain were modeled and the full-length structure was obtained by molecular docking simulation of these two domains. The orientation of the c-type-heme-containing domain relative to the N(2)OR domain is similar to that found in the other electron transfer complexes.

According to the model proposed in previous papers [Pettigrew, G. W., Prazeres, S., Costa, C., Palma, N., Krippahl, L., and Moura, J. J. (1999) The structure of an electron-transfer complex containing a cytochrome c and a peroxidase, J. Biol. Chem. 274, 11383-11389; Pettigrew, G. W., Goodhew, C. F., Cooper, A., Nutley, M., Jumel, K., and Harding, S. E. (2003) Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans, Biochemistry 42, 2046-2055], cytochrome c peroxidase of Paracoccus denitrificans can accommodate horse cytochrome c and Paracoccus cytochrome c(550) at different sites on its molecular surface. Here we use (1)H NMR spectroscopy, analytical ultracentrifugation, molecular docking simulation, and microcalorimetry to investigate whether these small cytochromes can be accommodated simultaneously in the formation of a ternary complex. The pattern of perturbation of heme methyl and methionine methyl resonances in binary and ternary solutions shows that a ternary complex can be formed, and this is confirmed by the increase in the sedimentation coefficient upon addition of horse cytochrome c to a solution in which cytochrome c(550) fully occupies its binding site on cytochrome c peroxidase. Docking experiments in which favored binary solutions of cytochrome c(550) bound to cytochrome c peroxidase act as targets for horse cytochrome c and the reciprocal experiments in which favored binary solutions of horse cytochrome c bound to cytochrome c peroxidase act as targets for cytochrome c(550) show that the enzyme can accommodate both cytochromes at the same time on adjacent sites. Microcalorimetric titrations are difficult to interpret but are consistent with a weakened binding of horse cytochrome c to a binary complex of cytochrome c peroxidase and cytochrome c(550) and binding of cytochrome c(550) to the cytochrome c peroxidase that is affected little by the presence of horse cytochrome c in the other site. The presence of a substantial capture surface for small cytochromes on the cytochrome c peroxidase has implications for rate enhancement mechanisms which ensure that the two electrons required for re-reduction of the enzyme after reaction with hydrogen peroxide are delivered efficiently.

A hypothetical model of the complex formed between the iron-sulfur protein rubredoxin and the tetraheme cytochrome c3 from the sulfate-reducing bacteria Desulfovibrio vulgaris (Hildenborough) has been proposed utilizing computer graphic modeling, computational methods and NMR spectroscopy. The proposed complex appears feasible on the basis of complementary electrostatic interaction and steric factors and is consistent with the data from NMR experiments. In this model, the non-heme iron atom of rubredoxin is in close proximity to heme 1 of cytochrome c3. The complex is stabilized by charge-pair interactions and hydrogen bonds. This complex is compared to the flavodoxin-cytochrome c3 complex previously proposed [Stewart, D. E., LeGall, J., Moura, I., Moura, J. J. G., Peck, H. D. Jr, Xavier, A. V., Weiner, P. K. & Wampler, J. E. (1988) Biochemistry 27, 2444-2450] and new NMR data shows that both proteins interact with the same heme group of the cytochrome as postulated.

Metalloenzymes control enzymatic activity by changing the characteristics of the metal centers where catalysis takes place. The conversion between inactive and active states can be tuned by altering the coordination number of the metal site, and in some cases by an associated conformational change. These processes will be illustrated using heme proteins (cytochrome c nitrite reductase, cytochrome c peroxidase and cytochrome cd1 nitrite reductase), non-heme proteins (superoxide reductase and [NiFe]-hydrogenase), and copper proteins (nitrite and nitrous oxide reductases) as examples. These examples catalyze electron transfer reactions that include atom transfer, abstraction and insertion.

Electron transfer proteins and redox enzymes containing paramagnetic redox centers with different relaxation rates are widespread in nature. Despite both the long distances and chemical paths connecting these centers, they can present weak magnetic couplings produced by spin-spin interactions such as dipolar and isotropic exchange. We present here a theoretical model based on the Bloch-Wangsness-Redfield theory to analyze the dependence with temperature of EPR spectra of interacting pairs of spin 1/2 centers having different relaxation rates, as is the case of the molybdenum-containing enzyme aldehyde oxidoreductase from Desulfovibrio gigas. We analyze the changes of the EPR spectra of the slow relaxing center (Mo(V)) induced by the faster relaxing center (FeS center). At high temperatures, when the relaxation time T(1) of the fast relaxing center is very short, the magnetic coupling between centers is averaged to zero. Conversely, at low temperatures when T(1) is longer, no modulation of the coupling between metal centers can be detected.

The periplasmic hydrogenase containing equivalent amounts of nickel and selenium plus non-heme iron [NiFeSe) hydrogenase) has been purified from cells of the sulfate reducing bacterium Desulfovibrio baculatus (DSM 1748) grown on a lactate/sulfate medium containing natural Se isotopes and the nuclear isotope, 77Se. Both the 77Se-enriched and unenriched hydrogenases were shown to be free of other hydrogenases and characterized with regard to their Se contents. EPR studies of the reduced nickel signal generated by redox titrations of the enriched and unenriched (NiFeSe) hydrogenases demonstrated that the gx = 2.23 and gy = 2.17 resonances are appreciably broadened by the spin of the 77Se nucleus (I = 1/2). This observation demonstrates unambiguously that the unpaired electron is shared by the Ni and Se atoms and that Se serves as a ligand to the nickel redox center of the (NiFeSe) hydrogenase.

The periplasmic hydrogenase of Desulfovibrio vulgaris (Hildenbourough NCIB 8303) belongs to the category of [Fe] hydrogenase which contains only iron-sulfur clusters as its prosthetic groups. Amino acid analyses were performed on the purified D. vulgaris hydrogenase. The amino acid composition obtained compared very well with the result derived from the nucleotide sequence of the structural gene (Voordouw, G., Brenner, S. (1985) Eur. J. Biochem. 148, 515-520). Detailed EPR reductive titration studies on the D. vulgaris hydrogenase were performed to characterize the metal centers in this hydrogenase. In addition to the three previously observed EPR signals (namely, the "isotropic" 2.02 signal, the rhombic 2.10 signal, and the complex signal of the reduced enzyme), a rhombic signal with resonances at the g-values of 2.06, 1.96, and 1.89 (the rhombic 2.06 signal) was detected when the samples were poised at potentials between 0 and -250 mV (with respect to normal hydrogen electrode). The midpoint redox potentials for each of the four EPR-active species were determined, and the characteristics of each EPR signal are described. Both the rhombic 2.10 and 2.06 signals exhibit spectral properties that are distinct from a ferredoxin-type [4Fe-4S] cluster and are proposed to originate from the same H2-binding center but in two different conformations. The complex signal of the reduced hydrogenase has been shown to represent two spin-spin interacting ferredoxin-type [4Fe-4S]1+ clusters (Grande, H. J., Dunham, W. R., Averill, B., Van Dijk, C., and Sands, R. H. (1983) Eur. J. Biochem. 136, 201-207). The titration data indicated a strong cooperative effect between these two clusters during their reduction. In an effort to accurately estimate the number of iron atoms/molecule of hydrogenase, plasma emission and chemical methods were used to determine the iron contents in the samples; and four different methods, including amino acid analysis, were used for protein determination. The resulting iron stoichiometries were found to be method-dependent and vary over a wide range (+/- 20%). The uncertainties involved in the determination of iron stoichiometry are discussed.

Small electron-transfer proteins such as flavodoxin (16 kDa) and the tetraheme cytochrome c3 (13 kDa) have been used to mimic, in vitro, part of the complex electron-transfer chain operating between substrate electron donors and respiratory electron acceptors, in sulfate-reducing bacteria (Desulfovibrio species). The nature and properties of the complex formed between these proteins are revealed by 1H-NMR and molecular modeling approaches. Our previous study with the Desulfovibrio vulgaris proteins [Moura, I., Moura, J.J. G., Santos, M.H., & Xavier, A. V. (1980) Cienc. Biol. (Portugal) 5, 195-197; Stewart, D.E. LeGall, J., Moura, I., Moura, J. J. G., Peck, H.D. Jr., Xavier, A. V., Weiner, P. K., & Wampler, J.E. (1988) Biochemistry 27, 2444-2450] indicated that the complex between cytochrome c3 and flavodoxin could be monitored by changes in the NMR signals of the heme methyl groups of the cytochrome and that the electrostatic surface charge (Coulomb's law) on the two proteins favored interaction between one unique heme of the cytochrome with flavodoxin. If the interaction is indeed driven by the electrostatic complementarity between the acidic flavodoxin and a unique positive region of the cytochrome c3, other homologous proteins from these two families of proteins might be expected to interact similarly. In this study, three homologous Desulfovibrio cytochromes c3 were used, which show a remarkable variation in their individual isoelectric points (ranging from 5.5 to 9.5). On the basis of data obtained from protein-protein titrations followed at specific proton NMR signals (i.e., heme methyl resonances), a binding model for this complex has been developed with evaluation of stoichiometry and binding constants. This binding model involves one site on the cytochromes c3 and two sites on the flavodoxin, with formation of a ternary complex at saturation. In order to understand the potential chemical form of the binding model, a structural model for the hypothetical ternary complex, formed between one molecule of Desulfovibrio salexigens flavodoxin and two molecules of cytochrome c3, is proposed. These molecular models of the complexes were constructed on the basis of complementarity of Coulombic electrostatic surface potentials, using the available X-ray structures of the isolated proteins and, when required, model structures (D. salexigens flavodoxin and Desulfovibrio desulfuricans ATCC 27774 cytochrome c3) predicted by homology modeling.

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.

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 Mossbauer 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 Mossbauer 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 Mossbauer data also revealed two [4Fe-4S]2+ clusters iun each molecule of hydrogenase II. The EPR and Mossbauer data of hydrogenase I were found to be identical to those of hydrogenase II, indicating that both enzymes have common metallic centers.

Ni and Se x-ray absorption spectroscopic studies of the [NiFeSe]hydrogenases from Desulfovibrio baculatus are described. The Ni site geometry is pseudo-octahedral with a coordinating ligand composition of 3-4 (N,O) at 2.06 A, 1-2 (S,Cl) at 2.17 A, and 1 Se at 2.44 A. The Se coordination environment consists of 1 C at 2.0 A and a heavy scatterer M (M = Ni or Fe) at approximately 2.4 A. These results are interpreted in terms of a selenocysteine residue coordinated to the Ni site. The possible role of the Ni-Se site in the catalytic activation of H2 is discussed.

Two cyclen-derived Gd probes, [Gd-DOTAM](3+) and [Gd-DOTP](5-) (DOTAM = 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetamide; DOTP = 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetrakis(methylenephosphonate)), were assessed as paramagnetic relaxation enhancement (PRE)-inducing probes for characterization of protein-protein interactions. Two proteins, Desulfovibrio gigas rubredoxin and Desulfovibrio gigas cytochrome c(3), were used as model partners. In a (1)H NMR titration it was shown that [Gd-DOTP](5-) binds to cytochrome c(3) near heme IV, causing pronounced PREs, characterized by line width broadenings of the heme methyl resonances at ratios as low as 0.08. A K(d) of 23 +/- 1 muM was calculated based on chemical shift perturbation of selected heme methyl resonances belonging to three different heme groups, caused by allosteric effects upon [Gd-DOTP](5-) binding to cytochrome c(3) at a molar ratio of 2. The other probe, [Gd-DOTAM](3+), caused PREs on a well-defined patch near the metal center of rubredoxin (especially the patch constituted by residues D19-G23 and W37-S45, which broaden beyond detection). This effect was partially reversed for some resonances (C6-Y11, in particular) when cytochrome c(3) was added to this system. Both probes were successful in causing reversible PREs at the partner binding site, thus showing to be good probes to identify partners' binding sites and since the interaction is reversible to structurally characterize protein complexes by better defining the complex interface.

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.

Mossbauer and EPR spectroscopy were used to characterize the heme prosthetic groups of the nitrite reductase isolated from Desulfovibrio desulfuricans (ATCC 27774), which is a membrane-bound multiheme cytochrome capable of catalyzing the 6-electron reduction of nitrite to ammonia. At pH 7.6, the as-isolated enzyme exhibited a complex EPR spectrum consisting of a low-spin ferric heme signal at g = 2.96, 2.28, and 1.50 plus several broad resonances indicative of spin-spin interactions among the heme groups. EPR redox titration studies revealed yet another low-spin ferric heme signal at g = 3.2 and 2.14 (the third g value was undetected) and the presence of a high-spin ferric heme. Mossbauer measurements demonstrated further that this enzyme contained six distinct heme groups: one high-spin (S = 5/2) and five low-spin (S = 1/2) ferric hemes. Characteristic hyperfine parameters for all six hemes were obtained through a detailed analysis of the Mossbauer spectra. D. desulfuricans nitrite reductase can be reduced by chemical reductants, such as dithionite or reduced methyl viologen, or by hydrogenase under hydrogen atmosphere. Addition of nitrite to the fully reduced enzyme reoxidized all five low-spin hemes to their ferric states. The high-spin heme, however, was found to complex NO, suggesting that the high-spin heme could be the substrate binding site and that NO could be an intermediate present in an enzyme-bound form.

This article aims to study hydrogen production/consumption in Desulfovibrio (D.) desulfuricans strain New Jersey, a sulfate reducer isolated from a medium undergoing active biocorrosion and to compare its hydrogen metabolism with two other Desulfovibrio species, D. gigas and D. vulgaris Hildenborough. Hydrogen production was followed during the growth of these three bacterial species under different growth conditions: no limitation of sulfate and lactate, sulfate limitation, lactate limitation, pyruvate/sulfate medium and in the presence of molybdate. Hydrogen production/consumption by D. desulfuricans shows a behavior similar to that of D. gigas but a different one from that of D. vulgaris, which produces higher quantities of hydrogen on lactate/sulfate medium. The three species are able to increase the hydrogen production when the sulfate became limiting. Moreover, in a pyruvate/sulfate medium hydrogen production was lower than on lactate/sulfate medium. Hydrogen production by D. desulfuricans in presence of molybdate is extremely high. Hydrogenases are key enzymes on production/consumption of hydrogen in sulfate reducing organisms. The specific activity, number and cellular localization of hydrogenases vary within the three Desulfovibrio species used in this work, which could explain the differences observed on hydrogen utilization.

A hypothetical model of the flavodoxin-tetraheme cytochrome c3 electron-transfer complex from the sulfate-reducing bacterium Desulfovibrio vulgaris has been constructed by using interactive computer graphics based on electrostatic potential field calculations and previous NMR experiments. Features of the proposed complex are (1) van der Waals contact between the flavin mononucleotide prosthetic group of flavodoxin and one heme of the cytochrome, (2) unique complementarity of electrostatic fields between the region surrounding this heme and the region surrounding the exposed portion of the flavin mononucleotide group of flavodoxin, and (3) no steric interferences between the two polypeptide chains in the complex. This complex is consistent with all structural and spectroscopic data available.