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.

Fuscoredoxin is a unique iron containing protein of yet unknown function originally discovered in the sulfate reducers of the genus Desulfovibrio. It contains two iron-sulfur clusters: a cubane [4Fe-4S] and a mixed oxo- and sulfido-bridged 4Fe cluster of unprecedented structure. The recent determination of the genomic sequence of Escherichia coli (E. coli) has revealed a homologue of fuscoredoxin in this facultative microbe. The presence of this gene in E. coli raises interesting questions regarding the function of fuscoredoxin and whether this gene represents a structural homologue of the better-characterized Desulfovibrio proteins. In order to explore the latter, an overexpression system for the E. coli fuscoredoxin gene was devised. The gene was cloned from genomic DNA by use of the polymerase chain reaction into the expression vector pT7-7 and overexpressed in E. coli BL21(DE3) cells. After two chromatographic steps a good yield of recombinant protein was obtained (approximately 4 mg of pure protein per liter of culture). The purified protein exhibits an optical spectrum characteristic of the homologue from D. desulfuricans, indicating that cofactor assembly was accomplished. Iron analysis indicated that the protein contains circa 8 iron atoms/molecule which were shown by EPR and Mossbauer spectroscopies to be present as two multinuclear clusters, albeit with slightly altered spectroscopic features. A comparison of the primary sequences of fuscoredoxins is presented and differences on cluster coordination modes are discussed on the light of the spectroscopic data.

The production of cytochrome c peroxidase (CCP) from Pseudomonas (Ps.) stutzeri (ATCC 11607) was optimized by adjusting the composition of the growth medium and aeration of the culture. The protein was isolated and characterized biochemically and spectroscopically in the oxidized and mixed valence forms. The activity of Ps. stutzeri CCP was studied using two different ferrocytochromes as electron donors: Ps. stutzeri cytochrome C-551 (the physiological electron donor) and horse heart cytochrome c. These electron donors interact differently with Ps. stutzeri CCP, exhibiting different ionic strength dependence. The CCP from Paracoccus (Pa.) denitrificans was proposed to have two different Ca2+ binding sites: one usually occupied (site I) and the other either empty or partially occupied in the oxidized enzyme (site II). The Ps. stutzeri enzyme was purified in a form with tightly bound Ca2+. The affinity for Ca2+ in the mixed valence enzyme is so high that Ca2+ returns to it from the EGTA which was added to empty the site in the oxidized enzyme. Molecular mass determination by ultracentrifugation and behavior on gel filtration chromatography have revealed that this CCP is isolated as an active dimer, in contrast to the Pa. denitrificans CCP which requires added Ca2+ for formation of the dimer and also for activation of the enzyme. This is consistent with the proposal that Ca2+ in the bacterial peroxidases influences the monomer/dimer equilibrium and the transition to the active form of the enzyme. Additional Ca2+ does affect both the kinetics of oxidation of horse heart cytochrome c (but not cytochrome C-551) and higher aggregation states of the enzyme. This suggests the presence of a superficial Ca2+ binding site of low affinity.

The structural changes in the heme macrocycle and substituents caused by binding of Ca2+ to the diheme cytochrome c peroxidase from Paracoccus pantotrophus were clarified by resonance Raman spectroscopy of the inactive fully oxidized form of the enzyme. The changes in the macrocycle vibrational modes are consistent with a Ca2+-dependent increase in the out-of-plane distortion of the low-potential heme, the proposed peroxidatic heme. Most of the increase in out-of-plane distortion occurs when the high-affinity site I is occupied, but a small further increase in distortion occurs when site II is also occupied by Ca2+ or Mg2+. This increase in the heme distortion explains the red shift in the Soret absorption band that occurs upon Ca2+ binding. Changes also occur in the low-frequency substituent modes of the heme, indicating that a structural change in the covalently attached fingerprint pentapeptide of the LP heme occurs upon Ca2+ binding to site I. These structural changes may lead to loss of the sixth ligand at the peroxidatic heme in the semireduced form of the enzyme and activation.

Nitric Oxide Reductase (NOR) is an integral membrane protein performing the reduction of NO to N(2)O. NOR is composed of two subunits: the large one (NorB) is a bundle of 12 transmembrane helices (TMH). It contains a b type heme and a binuclear iron site, which is believed to be the catalytic site, comprising a heme b and a non-hemic iron. The small subunit (NorC) harbors a cytochrome c and is attached to the membrane through a unique TMH. With the aim to perform structural and functional studies of NOR, we have immunized dromedaries with NOR and produced several antibody fragments of the heavy chain (VHHs, also known as nanobodies). These fragments have been used to develop a faster NOR purification procedure, to proceed to crystallization assays and to analyze the electron transfer of electron donors. BIAcore experiments have revealed that up to three VHHs can bind concomitantly to NOR with affinities in the nanomolar range. This is the first example of the use of VHHs with an integral membrane protein. Our results indicate that VHHs are able to recognize with high affinity distinct epitopes on this class of proteins, and can be used as versatile and valuable tool for purification, functional study and crystallization of integral membrane proteins.

Redox intermediates of D. desulfuricans ATCC 27774 [NiFe] hydrogenase were generated under dihydrogen. Detailed redox titrations, coupled to EPR measurements, give access to the mid-point redox potentials of the iron-sulfur centers and of the Nickel-B signal that represents the ready form of the enzyme. The interaction between the dihydrogen molecule and the nickel centre was probed by the observation of an isotopic effect on the EPR signals detected in turnover conditions, by comparison of the H2O/H2 and D2O/D2-reacted samples.

This communication reports the isolation, purification and characterization of key enzymes involved in dissimilatory sulfate reduction of a sulfate reducing bacterium classified as Desulfovibrio desulfuricans subspecies desulfuricans New Jersey (NCIMB 8313) (Ddd NJ). The chosen strain, originally recovered from a corroding cast iron heat exchanger, was grown in large scale batch cultures. Physico-chemical and spectroscopic studies of the purified enzymes were carried out. These analyses revealed a high degree of similarity between proteins isolated from the DddNJ strain and the homologous proteins obtained from Desulfomicrobium baculatus Norway 4. In view of the results obtained, taxonomic reclassification of Desulfovibrio desulfuricans subspecies desulfuricans New Jersey (NCIMB 8313) into Desulfomicrobium baculatus (New Jersey) is proposed.

We have cloned the gene encoding Desulfovibrio gigas ferredoxin using a photodigoxigenin-labelled probe synthesized with the polymerase chain reaction. The DNA sequence of the gene predicts a polypeptide of 58 residues after removal of the initial formyl methionine (polypeptide M(r) = 6,276). The ferredoxin gene was expressed in aerobically grown E. coli behind the lac promoter of pUC18 resulting in a high level of ferredoxin expression which comprises about 10% of the total cell protein. EPR analysis of recombinant ferredoxin revealed the presence of a [3Fe-4S] cluster which is characteristic of native D. gigas ferredoxin II.

Cytochrome c peroxidase from Paracoccus denitrificans LMD 52.44 was recently identified. The enzyme contains two c-type haems: one is reducible physiologically by cytochrome c550 from the same organism or non-physiologically by ascorbate (high-potential haem) and the other by dithionite (low-potential haem). The enzymatically active form of the peroxidase is the half-reduced enzyme state, in which the high-potential haem is in the iron(II) state and the low-potential haem is in the iron(III) state. It was found that the two haems interact and that the enzyme binds calcium ions near the haem sites which are necessary to promote its activation. In the oxidized form, the high-potential haem is in a high-spin and the low-potential haem is in a low-spin state. The half-reduction of the enzyme with ascorbate-diaminodurol changes the high-potential haem (high-spin) into a low-spin state and the low-potential haem converts from a low- into a high-spin state. This high-spin conversion of the low-potential haem is induced by the presence of calcium ions. These processes of reduction and spin state change can be easily resolved in time by removing the calcium from the enzyme using EDTA, facilitating the observation of the intermediate form by NMR.

Pseudoazurin binds at a single site on cytochrome c peroxidase from Paracoccus pantotrophus with a K(d) of 16.4 microM at 25 degrees C, pH 6.0, in an endothermic reaction that is driven by a large entropy change. Sedimentation velocity experiments confirmed the presence of a single site, although results at higher pseudoazurin concentrations are complicated by the dimerization of the protein. Microcalorimetry, ultracentrifugation, and (1)H NMR spectroscopy studies in which cytochrome c550, pseudoazurin, and cytochrome c peroxidase were all present could be modeled using a competitive binding algorithm. Molecular docking simulation of the binding of pseudoazurin to the peroxidase in combination with the chemical shift perturbation pattern for pseudoazurin in the presence of the peroxidase revealed a group of solutions that were situated close to the electron-transferring heme with Cu-Fe distances of about 14 A. This is consistent with the results of (1)H NMR spectroscopy, which showed that pseudoazurin binds closely enough to the electron-transferring heme of the peroxidase to perturb its set of heme methyl resonances. We conclude that cytochrome c550 and pseudoazurin bind at the same site on the cytochrome c peroxidase and that the pair of electrons required to restore the enzyme to its active state after turnover are delivered one-by-one to the electron-transferring heme.

The nitrite reductase (Nir) isolated from Pseudomonas chlororaphis DSM 50135 is a blue enzyme, with type 1 and type 2 copper centers, as in all copper-containing Nirs described so far. For the first time, a direct determination of the reduction potentials of both copper centers in a Cu-Nir was performed: type 2 copper (T2Cu), 172 mV and type 1 copper (T1Cu), 298 mV at pH 7.6. Although the obtained values seem to be inconsistent with the established electron-transfer mechanism, EPR data indicate that the binding of nitrite to the T2Cu center increases its potential, favoring the electron-transfer process. Analysis of the EPR spectrum of the turnover form of the enzyme also suggests that the electron-transfer process between T1Cu and T2Cu is the fastest of the three redox processes involved in the catalysis: (a) reduction of T1Cu; (b) oxidation of T1Cu by T2Cu; and (c) reoxidation of T2Cu by NO2-. Electrochemical experiments show that azurin from the same organism can donate electrons to this enzyme.

The orange-coloured protein (ORP) from Desulfovibrio gigas is a 12 kDa protein that contains a novel mixed-metal sulfide cluster of the type [S(2)MoS(2)CuS(2)MoS(2)]. Diffracting crystals of the apo form of ORP have been obtained. Data have been collected for the apo form of ORP to 2.25 A resolution in-house and to beyond 2.0 A resolution at ESRF, Grenoble. The crystals belonged to a trigonal space group, with unit-cell parameters a = 43, b = 43, c = 106 A.

Cytochrome c peroxidase was expressed in cells of Pseudomonas nautica strain 617 grown under microaerophilic conditions. The 36.5 kDa dihaemic enzyme was purified to electrophoretic homogeneity in three chromatographic steps. N-terminal sequence comparison showed that the Ps. nautica enzyme exhibits a high similarity with the corresponding proteins from Paracoccus denitrificans and Pseudomonas aeruginosa. UV-visible spectra confirm calcium activation of the enzyme through spin state transition of the peroxidatic haem. Monohaemic cytochrome c(552) from Ps. nautica was identified as the physiological electron donor, with a half-saturating concentration of 122 mu M and allowing a maximal catalytic centre activity of 116 000 min(-1). Using this cytochrome the enzyme retained the same activity even at high ionic strength. There are indications that the interactions between the two redox partners are mainly hydrophobic in nature. (C) 1999 Elsevier Science B.V. All rights reserved.

Below 30 K, oxidized Desulfovibrio gigas hydrogenase presents an intense electron paramagnetic resonance (EPR) signal centered at g = 2.02, typical of an iron-sulfur center. In addition a rhombic EPR signal, attributed to Ni(III) species, is also observed [LeGall, J., Ljungdahl, P., Moura, I., Peck, H.D., Jr, Xavier, A.V., Moura, J.J.G., Teixeira, M., Huynh, B.H., and DerVartanian, D.V. (1982) Biochem. Biophys. Res. Commun. 106, 610-616; and Cammack, R., Patil, D., Aguirre, R., and Hatchikian, E.C., (1982) FEBS Lett. 142, 289-292]. At higher temperatures (77 K) the iron-sulfur EPR signal is broader and all the EPR features of the rhombic nickel signal can easily be observed. We have now obtained additional information concerning the redox properties of these EPR active centers, using an EPR redox titration method in the presence of dye mediators at pH = 8.5. The mid-point potential was determined to be -70 mV for the Fe,S cluster and -220 mV for the Ni center. Intermediate oxidation states were obtained upon partial reduction with either dithionite or hydrogen. Although upon dithionite reduction the centers are reduced in the order of decreasing mid-point reduction potentials, under a hydrogen atmosphere the nickel center reduces preferentially. This suggests a catalytic involvement of the nickel redox center in the binding of hydrogen. Preliminary Mossbauer studies on Desulfovibrio gigas hydrogenase reveal the presence of a paramagnetic 3 Fe center and two 4 Fe centers. The 3 Fe center is responsible for the g = 2.02 EPR signal but the two 4 Fe centers have been so far undetectable by EPR.

The electron-transfer kinetics between three different mediators and the hexahemic enzyme nitrite reductase isolated from Desulfovibrio desulfuricans (ATCC 27774) were investigated by cyclic voltammetry and by chronoamperometry. The mediators, methyl viologen, Desulfovibrio vulgaris (Hildenborough) cytochrome c3 and D. desulfuricans (ATCC 27774) cytochrome c3 differ in structure, redox potential and charge. The reduced form of each mediator exchanged electrons with nitrite reductase. Second-order rate constants, k, were calculated on the basis of the theory for a simple catalytic mechanism and the results, obtained by cyclic voltammetry, were compared with those obtained by chronoamperometry. Values for k are in the range 10(6)-10(8) M-1 s-1 and increase in the direction D. desulfuricans cytochrome c3-->D. vulgaris cytochrome c3-->methyl viologen. An explanation is advanced on the basis of electrostatic interactions and relative orientation between the partners involved. Chronoamperometry (computer controlled) offers advantages over cyclic voltammetry in the determination of homogeneous rate constants (faster, more accurate and better reproducibility). Direct, unmediated electrochemical responses of the hexaheme nitrite reductase were also reported.