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.

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.

Adenylate kinase (AK; ATP:AMP phosphotransferase; EC 2.7.4.3) is involved in the reversible transfer of the terminal phosphate group from ATP to AMP. AKs contribute to the maintenance of a constant level of cellular adenine nucleotides, which is necessary for the energetic metabolism of the cell. Three metal ions, cobalt, zinc and iron(II), have been reported to be present in AKs from some Gram-negative bacteria. Native zinc-containing AK from Desulfovibrio gigas was purified to homogeneity and crystallized. The crystals diffracted to beyond 1.8 A resolution. Furthermore, cobalt- and iron-containing crystal forms of recombinant AK were also obtained and diffracted to 2.0 and 3.0 A resolution, respectively. Zn(2+)-AK and Fe(2+)-AK crystallized in space group I222 with similar unit-cell parameters, whereas Co(2+)-AK crystallized in space group C2; a monomer was present in the asymmetric unit for both the Zn(2+)-AK and Fe(2+)-AK forms and a dimer was present for the Co(2+)-AK form. The structures of the three metal-bound forms of AK will provide new insights into the role and selectivity of the metal in these enzymes.

Rubredoxin and desulforedoxin both contain an Fe(S-Cys)4 center. However, the spectroscopic properties of the center in desulforedoxin differ from rubredoxin. These differences arise from a distortion of the metal site hypothesized to result from adjacent cysteine residues in the primary sequence of desulforedoxin. Two desulforedoxin mutants were generated in which either a G or P-V were inserted between adjacent cysteines. Both mutants exhibited optical spectra with maxima at 278, 345, 380, 480, and 560 nm while the low temperature X-band EPR spectra indicated highspin Fe3+ ions with large rhombic distortions (E/D = 0.21-0.23). These spectroscopic properties are distinct from wild type desulforedoxin and virtually identical to rubredoxin.

Adenylate kinases (AK) from Gram-negative bacteria are generally devoid of metal ions in their LID domain. However, three metal ions, zinc, cobalt, and iron, have been found in AK from Gram-negative bacteria. Crystal structures of substrate-free AK from Desulfovibrio gigas with three different metal ions (Zn(2+), Zn-AK; Co(2+), Co-AK; and Fe(2+), Fe-AK) bound in its LID domain have been determined by X-ray crystallography to resolutions 1.8, 2.0, and 3.0 A, respectively. The zinc and iron forms of the enzyme were crystallized in space group I222, whereas the cobalt-form crystals were C2. The presence of the metals was confirmed by calculation of anomalous difference maps and by X-ray fluorescence scans. The work presented here is the first report of a structure of a metal-containing AK from a Gram-negative bacterium. The native enzyme was crystallized, and only zinc was detected in the LID domain. Co-AK and Fe-AK were obtained by overexpressing the protein in Escherichia coli. Zn-AK and Fe-AK crystallized as monomers in the asymmetric unit, whereas Co-AK crystallized as a dimer. Nevertheless, all three crystal structures are very similar to each other, with the same LID domain topology, the only change being the presence of the different metal atoms. In the absence of any substrate, the LID domain of all holoforms of AK was present in a fully open conformational state. Normal mode analysis was performed to predict fluctuations of the LID domain along the catalytic pathway.

Membrane electrodes (ME) were constructed using gold, glassy carbon and pyrolytic graphite supports and a dialysis membrane, and used to study the electrochemical behavior of small size electron transfer proteins: monohemic cytochrome c(522) from Pseudomonas nautica and cytochrome c(533) as well as rubredoxin from Desulfovibrio vulgaris. Different electrochemical techniques were used including cyclic voltammetry (CV), square wave voltammetry (SW) and differential pulse voltammetry (DP). A direct electrochemical response was obtained in all cases except with rubredoxin where a facilitator was added to the protein solution entrapped between the membrane and the electrode surface. Formal potentials and heterogeneous charge transfer rate constants were determined from the voltammetric data. The influence of the ionic strength and the pH of the medium on the electrochemical response at the ME were analyzed. The benefits from the use of the ME in protein electrochemistry and its role in modulating the redox behavior are analyzed. A critical comparison is presented with data obtained at non-MEs. Finally, the interactions that must be established between the proteins and the electrode surfaces are discussed, thereby modeling molecular interactions that occur in biological systems. (C) 2002 Elsevier Science B.V. All rights reserved.

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.

The regulation of the properties of Desulfovibrio gigas flavodoxin in AOT/water/iso-octane micellar system was studied. UV-visible spectroscopic studies have shown that photoreduction of flavodoxin in the presence of EDTA leads to hydroquinone formation through the intermediate semiquinone. The [free FMN] - [bound to flavodoxin FMN] equilibrium (and hence, the amount of apoprotein) depends on redox state of FMN and on hydration degree which controls the micellar size. Thus, a new method of reversible cofactor removing under mild conditions (at low hydration degree of micelles) is suggested, accompained by isolation of apo-form of the protein.

The tetrameric form of a Desulfovibrio gigas ferredoxin, named Fd II, mediates electron transfer between cytochrome c3 and sulfite reductase. We have studied two stable oxidation states of this protein with Mossbauer spectroscopy and electron paramagnetic resonance. We found 3 iron atoms/monomer and a spin concentration of 0.9 spins/monomer for the oxidized protein. Taken together, the EPR and Mossbauer data demonstrate conclusively the presence of a spin-coupled structure containing 3 iron atoms and labile sulfur. The Mossbauer data show also that this metal center is structurally similar, if not identical, with the low potential center of a ferredoxin from Azotobacter vinelandii, a novel cluster described recently (Emptage, M.H., Kent, T.A., Huynh, B.H., Rawlings, J., Orme-Johnson, W.H., and Munck, E. (1980) J. Biol. Chem. 255, 1793-1796).

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.

The interaction of reduced rabbit cytochrome b(5) with reduced yeast iso-1 cytochrome c has been studied through the analysis of (1)H-(15)N HSQC spectra, of (15)N longitudinal ( R(1)) and transverse ( R(2)) relaxation rates, and of the solvent exchange rates of protein backbone amides. For the first time, the adduct has been investigated also from the cytochrome c side. The analysis of the NMR data was integrated with docking calculations. The result is that cytochrome b(5) has two negative patches capable of interacting with a single positive surface area of cytochrome c. At low protein concentrations and in equimolar mixture, two different 1:1 adducts are formed. At high concentration and/or with excess cytochrome c, a 2:1 adduct is formed. All the species are in fast exchange on the scale of differences in chemical shift. By comparison with literature data, it appears that the structure of one 1:1 adduct changes with the origin or primary sequence of cytochrome b(5).

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.

The localization of APS reductase and bisulfite reductase in Desulfovibrio gigas, D. vulgaris Hildenborough and D. thermophilus was studied by immunoelectron microscopy. Polyclonal antibodies were raised against the purified enzymes from each strain. Cells fixed with formaldehyde/glutaraldehyde were embedded and ultrathin sections were incubated with antibodies and subsequently labeled with protein A-gold. The bisulfite reductase in all three strains and APS reductase in d. gigas and D. vulgaris were found in the cytoplasm. The labeling of d. thermophilus with APS reductase antibodies resulted in a distribution of gold particles over the cytoplasmic membrane region. The localization of the two enzymes is discussed with respect to the mechanism and energetics of dissimilatory sulfate reduction.

Among the biotargets interacting with vanadium is the calcium pump from the sarcoplasmic reticulum (SR). To this end, initial research efforts were launched with two vanadium(V)-citrate complexes, namely (NH(4))(6)[V(2)O(4)(C(6)H(4)O(7))(2)].6H(2)O and (NH(4))(6)[V(2)O(2)(O(2))(2)(C(6)H(4)O(7))(2)].4H(2)O, potentially capable of interacting with the SR calcium pump by combining kinetic studies with (51)V NMR spectroscopy. Upon dissolution in the reaction medium (concentration range: 4-0.5mM), both vanadium(V):citrate (VC) and peroxovanadium(V):citrate (PVC) complexes are partially converted into vanadate oligomers. A 1mM solution of the PVC complex, containing 184microM of the PVC complex, 94microM oxoperoxovanadium(V) (PV) species, 222microM monomeric (V1), 43microM dimeric (V2) and 53microM tetrameric (V4) species, inhibits Ca(2+) accumulation by 75 %, whereas a solution of the VC complex of the same vanadium concentration, containing 98microM of the VC complex, 263microM monomeric (V1), 64microM dimeric (V2) and 92microM tetrameric (V4) species inhibits the calcium pump activity by 33 %. In contrast, a 1 mM metavanadate solution, containing 460microM monomeric (V1), 90.2microM dimeric (V2) and 80microM tetrameric (V4) species, has no effect on Ca(2+) accumulation. The NMR signals from the VC complex (-548.0ppm), PVC complex (-551.5ppm) and PV (-611.1ppm) are broadened upon SR vesicle addition (2.5mg/ml total protein). The relative order for the half width line broadening of the NMR signals, which reflect the interaction with the protein, was found to be V4>PVC>VC>PV>V2=V1=1, with no effect observed for the V1 and V2 signals. Putting it all together the effects of two vanadium(V)-citrate complexes on the modulation of calcium accumulation and ATP hydrolysis by the SR calcium pump reflected the observed variable reactivity into the nature of key species forming upon dissolution of the title complexes in the reaction media.

We report the kinetic behavior of the enzyme aldehyde oxidoreductase (AOR) from the sulfate reducing bacterium Desulfovibrio gigas (Dg) encapsulated in reverse micelles of sodium bis-(2-ethylhexyl) sulfosuccinate in isooctane using benzaldehyde, octaldehyde, and decylaldehyde as substrates. Dg AOR is a 200-kDa homodimeric protein that catalyzes the conversion of aldehydes to carboxylic acids. Ultrasedimentation analysis of Dg AOR-containing micelles showed the presence of 100-kDa molecular weight species, confirming that the Dg AOR subunits can be dissociated. UV-visible spectra of encapsulated Dg AOR are indistinguishable from the enzyme spectrum in solution, suggesting that both protein fold and metal cofactor are kept intact upon encapsulation. The catalytic constant (k(cat)) profile as a function of the micelle size W(0) (W(0)=[H(2)O]/[AOT]) using benzaldehyde as substrate showed two bell-shaped activity peaks at W(0)=20 and 26. Furthermore, enzymatic activity for octaldehyde and decylaldehyde was detected only in reverse micelles. Like for the benzaldehyde kinetics, two peaks with both similar k(cat) values and W(0) positions were obtained. EPR studies using spin-labeled reverse micelles indicated that octaldehyde and benzaldehyde are intercalated in the micelle membrane. This suggests that, though Dg AOR is found in the cytoplasm of bacterial cells, the enzyme may catalyze the reaction of substrates incorporated into a cell membrane.

After up to 21 days without food, adult male quails (Coturnix coturnix japonica) lost about 45% of the initial body weight (100-150 g). As in naturally fast-adapted and larger birds, three phases were identified during prolonged fasting in quails. Phase I lasted 2-3 days and was characterized by a rapid decrease in the rate of body weight loss and high fat mobilization. Phase II was longer and characterized by a slow and steady decline in the rates of body weight loss and of nitrogen excretion. The third (critical) period was marked by an abrupt increase in the rates of body weight loss and of nitrogen excretion. Despite their small size, the duration of phase II in quails was relatively long, a clear advantage for the study of the relationships between the several metabolic events that occur during this crucial adaptative period. Also, the beginning of phase III could be precisely determined. Changes in blood glucose, plasma FFA and triacylglycerols levels, as well as in liver and carcass lipid content were similar to those found in other species of birds. Therefore, quails seem to be a suitable model to investigate the biochemical mechanisms involved in the metabolic adjustments to prolonged food deprivation in non fasting-adapted birds.

Desulforedoxin and the N-terminus of desulfoferrodoxin share a 36 amino acid domain containing a (Cys-S)(4) metal binding site. Recombinant forms of desulforedoxin, an N-terminal fragment of desulfoferrodoxin, and two desulforedoxin mutant proteins were reconstituted with Fe3+ Cd2+, and Zn2+ and relative metal ion affinities assessed by proton titrations. Protons compete with metal for protein ligands, a process that can be followed by monitoring the optical spectrum of the metal-protein complex as a function of pH. For all polypeptides, Fe3+ bound with the highest affinity, whereas the affinity of Zn2+ was greater than Cd2+ in desulforedoxin and the N-terminal fragment of desulfoferrodoxin, but this order was reversed in desulforedoxin mutant proteins. Metal binding in both mutants was significantly impaired. Furthermore, the Fe3+ complex of both mutants underwent a time-dependent bleaching process which coincided with increased reactivity of cysteine residues to Ellman's reagent and concomitant metal dissociation. It is hypothesized that this results from an autoredox reaction in which Fe3+ is reduced to Fe2+ with attendant oxidation of ligand thiols.

This article describes the method and results of our participation in the Critical Assessment of PRediction of Interactions (CAPRI) experiment, using the protein docking program BiGGER (Bimolecular complex Generation with Global Evaluation and Ranking) (Palma et al., Proteins 2000;39:372-384). Of five target complexes (CAPRI targets 2, 4, 5, 6, and 7), only one was successfully predicted (target 6), but BiGGER generated reasonable models for targets 4, 5, and 7, which could have been identified if additional biochemical information had been available.

Three-dimensional protein structures of the xanthine oxidase family show different solutions for the problem of transferring electrons between the flavin adenine dinucleotide (FAD) group and the molybdenum cofactor. In xanthine oxidase all the cofactors he within domains of the same protein chain, whereas in CO dehydrogenase the Fe-S centres, FAD and Mo cofactors are enclosed in separate chains and the enzyme exists as a stable complex of all three. In aldehyde oxidore-ductase, only Fe-S and Mo co-factors are present in a single protein chain. Flavodoxin is docked to aldehyde oxidoreductase to mimic the flavin component on the intramolecular electron transfer chain of aanthine oxidase and CO dehydrogenase and, remarkably, the main features of the electron-transfer pathway are observed.