Isolation and characterization of a rubredoxin and a two-(4Fe-4S) ferredoxin from Thermodesulfobacterium commune

The de novo design of a rubredoxin-like Fe site

A redox center similar to that of rubredoxin was designed into the 56 amino acid immunoglobulin binding B1 domain of Streptococcals protein G. The redox center in rubredoxin contains an iron ion tetrahedrally coordinated by four cysteine residues, [Fe(S-Cys)4](-1),(-2). The design criteria for the target site included taking backbone movements into account, tetrahedral metal-binding, and maintaining the structure and stability of the wild-type protein. The optical absorption spectrum of the Co(II) complex of the metal-binding variant is characteristic of tetrahedral chelation by four cysteine residues. Circular dichroism and nuclear magnetic resonance measurements reveal that the metal-free and Cd(II)-bound forms of the variant are folded correctly and are stable. The Fe(III) complex of the metal-binding mutant reproduces the optical and the electron paramagnetic resonance spectra of oxidized rubredoxin. This demonstrates that the engineered protein chelates Fe(III) in a tetrahedral array, and the resulting center is similar to that of oxidized rubredoxin.

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.

Physical characterization of a totally synthetic rubredoxin

The entire polypeptide of hyperthermophilic Pyrococcus furiosus rubredoxin was synthesized in order to specifically probe structural determinants of protein thermostability. The uv-visible, circular dichroic, electron paramagnetic, and nuclear magnetic resonance spectra, and electrochemical properties, of the native and synthetic proteins were essentially identical. The synthetic protein had a half-life for denaturation of 24 hr at 80 degrees C. The synthetic protein is considerably more thermostable than nonhyperthermophilic rubredoxins, but not as stable as the native protein. Based on the spectroscopic evidence, it appears that the synthetic protein is incorporating iron properly to form holoprotein, but the peptide still may not be folded correctly.

Purification and characterization of an NADH-rubredoxin oxidoreductase involved in the utilization of oxygen by Desulfovibrio gigas

An NADH--rubredoxin oxidoreductase previously isolated from Desulfovibrio gigas [LeGall, J. (1968) Ann. Inst. Pasteur 114, 109-115] has now been fully purified and further characterized. It contains two subunits of 27 kDa and 32 kDa. With two mid-point redox potentials of -295 mV and -325 mV, this FMN- and FAD-containing protein can induce the specific reduction of D. gigas rubredoxin. In contrast, rubredoxins from the other Desulfovibrio species or desulforedoxin from D. gigas show very low reaction rates with the same enzyme. The phylogenetic significance of the narrow specificity of the enzyme toward the rubredoxin from the same organism is discussed. The purified enzyme has NADH oxidase activity with H2O2 as a final product of O2 reduction. The reaction is half-inhibited by 4.2 microM p-chloromercuribenzoate, whereas cyanide and azide are not significant inhibitors in this reaction. The role of this protein as a part of the enzymic equipment that allows the formation of ATP in the presence of oxygen from the degradation of carbon reserves is discussed.

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.

Modeling the structure of Pyrococcus furiosus rubredoxin by homology to other X-ray structures

The three-dimensional structure of rubredoxin from the hyperthermophilic archaebacterium, Pyrococcus furiosus, has been modeled from the X-ray crystal structures of three homologous proteins from Clostridium pasteurianum, Desulfovibrio gigas, and Desulfovibrio vulgaris. All three homology models are similar. When comparing the positions of all heavy atoms and essential hydrogen atoms to the recently solved crystal structure (Day, M. W., et al., 1992, Protein Sci. 1, 1494-1507) of the same protein, the homology model differ from the X-ray structure by 2.09 A root mean square (RMS). The X-ray and the zinc-substituted NMR structures (Blake, P. R., et al., 1992b, Protein Sci. 1, 1508-1521) show a similar level of difference (2.05 A RMS). On average, the homology models are closer to the X-ray structure than to the NMR structures (2.09 vs. 2.42 A RMS).

X-ray crystal structures of the oxidized and reduced forms of the rubredoxin from the marine hyperthermophilic archaebacterium Pyrococcus furiosus

The structures of the oxidized and reduced forms of the rubredoxin from the archaebacterium, Pyrococcus furiosus, an organism that grows optimally at 100 degrees C, have been determined by X-ray crystallography to a resolution of 1.8 A. Crystals of this rubredoxin grow in space group P2(1)2(1)2(1) with room temperature cell dimensions a = 34.6 A, b = 35.5 A, and c = 44.4 A. Initial phases were determined by the method of molecular replacement using the oxidized form of the rubredoxin from the mesophilic eubacterium, Clostridium pasteurianum, as a starting model. The oxidized and reduced models of P. furiosus rubredoxin each contain 414 nonhydrogen protein atoms comprising 53 residues. The model of the oxidized form contains 61 solvent H2O oxygen atoms and has been refined with X-PLOR and TNT to a final R = 0.178 with root mean square (rms) deviations from ideality in bond distances and bond angles of 0.014 A and 2.06 degrees, respectively. The model of the reduced form contains 37 solvent H2O oxygen atoms and has been refined to R = 0.193 with rms deviations from ideality in bond lengths of 0.012 A and in bond angles of 1.95 degrees. The overall structure of P. furiosus rubredoxin is similar to the structures of mesophilic rubredoxins, with the exception of a more extensive hydrogen-bonding network in the beta-sheet region and multiple electrostatic interactions (salt bridge, hydrogen bonds) of the Glu 14 side chain with groups on three other residues (the amino-terminal nitrogen of Ala 1; the indole nitrogen of Trp 3; and the amide nitrogen group of Phe 29). The influence of these and other features upon the thermostability of the P. furiosus protein is discussed.

Determinants of protein hyperthermostability: purification and amino acid sequence of rubredoxin from the hyperthermophilic archaebacterium Pyrococcus furiosus and secondary structure of the zinc adduct by NMR

The purification, amino acid sequence, and two-dimensional 1H NMR results are reported for the rubredoxin (Rd) from the hyperthermophilic archaebacterium Pyrococcus furiosus, an organism that grows optimally at 100 degrees C. The molecular mass (5397 Da), iron content (1.2 +/- 0.2 g-atom of Fe/mol), UV-vis spectrophotometric properties, and amino acid sequence (60% sequence identity with Clostridium pasteurianum Rd) are found to be typical of this class of redox protein. However, P. furiosus Rd is remarkably thermostable, being unaffected after incubation for 24 h at 95 degrees C. One- and two-dimensional 1H nuclear magnetic resonance spectra of the oxidized [Fe(III)Rd] and reduced [Fe(II)Rd] forms of P. furiosus Rd exhibited substantial paramagnetic line broadening, and this precluded detailed 3D structural studies. The apoprotein was not readily amenable to NMR studies due to apparent protein oxidation involving the free cysteine sulfhydryls. However, high-quality NMR spectra were obtained for the Zn-substituted protein, Zn(Rd), enabling detailed NMR signal assignment for all backbone amide and alpha and most side-chain protons. Secondary structural elements were determined from qualitative analysis of 2D Overhauser effect spectra. Residues A1-K6, Y10-E14, and F48-E51 form a three-strand antiparallel beta-sheet, which comprises ca. 30% of the primary sequence. Residues C5-Y10 and C38-A43 form types I and II amide-sulfur tight turns common to iron-sulfur proteins. These structural elements are similar to those observed by X-ray crystallography for native Rd from the mesophile C. pasteurianum. However, the beta-sheet domain in P. furiosus Rd is larger than that in C. pasteurianum Rd and appears to begin at the N-terminal residue. From analysis of the secondary structure, potentially stabilizing electrostatic interactions involving the charged groups of residues Ala(1), Glu(14), and Glu(52) are proposed. These interactions, which are not present in rubredoxins from mesophilic organisms, may prevent the beta-sheet from "unzipping" at elevated temperatures.

Analysis of the transcriptional unit encoding the genes for rubredoxin (rub) and a putative rubredoxin oxidoreductase (rbo) in Desulfovibrio vulgaris Hildenborough

The nucleotide sequence of a 2.0-kilobase-pair EcoRI restriction fragment upstream from the gene (rub, 162 base pairs) encoding rubredoxin from Desulfovibrio vulgaris Hildenborough indicates that it is part of a larger transcriptional unit, containing an additional 378-base-pair open reading frame which terminates 16 nucleotides from the translational start of the rub gene and could encode a polypeptide of 14 kilodaltons (kDa). Northern (RNA) blotting of RNA isolated from both D. vulgaris Hildenborough and Escherichia coli TG2 transformed with plasmid pJK29, which contains both genes on a 1.1-kilobase-pair SalI insert, confirms that the genes for this 14-kDa polypeptide and rubredoxin are present on a single transcript of 680 nucleotides. Strong evidence that the 14-kDa polypeptide is also a redox protein is provided by the fact that its NH2 terminus is homologous to desulforedoxin, which has been isolated from D. gigas as a small dimeric redox protein (36 amino acids per monomer), coordinating two iron atoms. Since rubredoxin is a potential redox partner for the 14-kDa protein, it has been tentatively named rubredoxin oxidoreductase, produced by the rbo gene. Southern blotting indicates that the rbo-rub operon is present in several species and strains of sulfate-reducing bacteria.

Isolation, characterization, and biological activity of the Methanococcus thermolithotrophicus ferredoxin

A ferredoxin has been isolated from the thermophilic methanogen Methanococcus thermolithotrophicus. The native protein was a monomer exhibiting a molecular weight of 7,262, calculated from the amino acid composition. Its absorption spectrum had two maxima at 390 and 283 nm, with an absorbance ratio A390/A283 of 0.79. The absorption at 390 nm (E = 29 mM-1 cm-1) and the content of iron of the protein are in agreement with the presence of two 4Fe-4S clusters in M. thermolithotrophicus ferredoxin. Its amino acid composition showed the presence of eight cysteine residues, which is the required number of cysteines for the binding of two 4Fe-4S clusters. The protein was characterized by the lack of histidine, arginine, and leucine and a high content of valine. It was unusually stable to high temperatures but not to oxygen. The ESR spectrum of the protein in the oxidized state showed a minor signal at g = 2.01, corresponding to an oxidized 3Fe-4S cluster. The protein, which was difficult to reduce with dithionite or reduced mediators, exhibited in its reduced state a spectrum typical of two interacting reduced 4Fe-4S clusters. M. thermolithotrophicus ferredoxin functioned as an electron acceptor for the CO dehydrogenase complex with an extract free of ferredoxin. No reaction was detected with F420 or hydrogenase.

Rubredoxin from Desulfovibrio gigas. A molecular model of the oxidized form at 1.4 A resolution

The crystal structure of rubredoxin from the sulfate-reducing bacterium Desulfovibrio gigas has been determined at 1.4 A resolution (1 A = 0.1 nm) by X-ray diffraction methods; starting with a model of the isostructural rubredoxin from Desulfovibrio vulgaris. Refinement of the molecular model has been carried out by restrained least-squares techniques and Fourier series calculations. The present model includes a formyl at the N-terminal end and 121 possible sites for solvent molecules with full or partial occupancy, which corresponds to the modeling of nearly all the solvent medium. The crystallographic R factor against the data with 10 A greater than d greater than 1.4 A with F greater than 2 sig(F), is 0.136; and R = 0.140 when all the data are considered. The estimated average root-mean-square (r.m.s.) error on the positional parameters is about 0.12 A. The overall structural features of this molecule are close to those of the two highly refined rubredoxins from Clostridium pasteurianum and D. vulgaris. Superposition of these two molecules on the rubredoxin from D. gigas shows in both cases an overall r.m.s. deviation of 0.5 A for the atoms in the main-chain and of 0.4 A for the atoms in the side-chains that make up the hydrophobic core. The iron atom is co-ordinated to four cysteine sulfur atoms forming an almost regular tetrahedron, with Fe-SG distances ranging from 2.27 A to 2.31 A and angles varying from 103 degrees to 115 degrees. The intramolecular hydrogen-bonding pattern is quite comparable to those found in other proteins refined at high resolution. All the polar groups are involved in hydrogen bonds: intramolecular, intermolecular or with solvent molecules. The main structural differences from the other rubredoxins are in the nature and the distribution of some of the charged residues over the molecular surface. The possible influence of several structural factors on the intramolecular and intermolecular electron transfer properties such as the NH...SG bonds, the solvent exposure of the redox center, and the aromatic core is discussed. The conservation, during evolution, of a ring of acidic residues in the proximity of the FeSG4 center suggests that this ring may be implicated in the recognition processes between rubredoxins and their functional partners.