Rubredoxin from Pyrococcus furiosus

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Rubredoxin Protein Scaffolds Sourced from Diverse Environmental Niches as an Artificial Hydrogenase Platform

Nickel-substituted rubredoxin (NiRd) from Desulfovibrio desulfuricans has previously been shown to act as both a structural and functional mimic of the [NiFe] hydrogenase. However, improvements both in turnover frequency and overpotential are needed to rival the native [NiFe] hydrogenase enzymes. Characterization of a library of NiRd mutants with variations in the secondary coordination sphere suggested that protein dynamics played a substantial role in modulating activity. In this work, rubredoxin scaffolds were selected from diverse organisms to study the effects of distal sequence variation on catalytic activity. It was found that though electrochemical catalytic activity was only slightly impacted across the series, the Rd sequence from a psychrophilic organism exhibited substantially higher levels of solution-phase hydrogen production. Additionally, Eyring analyses suggest that catalytic activation properties relate to the growth temperature of the parent organism, implying that the general correlation between the parent organism environment and catalytic activity often seen in naturally occurring enzymes may also be observed in artificial enzymes. Selecting protein scaffolds from hosts that inhabit diverse environments, particularly low-temperature environments, represents an alternative approach for engineering artificial metalloenzymes.

A natural fusion of flavodiiron, rubredoxin, and rubredoxin oxidoreductase domains is a self-sufficient water-forming oxidase of Trichomonas vaginalis

Microaerophilic pathogens such as Giardia lamblia, Entamoeba histolytica, and Trichomonas vaginalis have robust oxygen consumption systems to detoxify oxygen and maintain intracellular redox balance. This oxygen consumption results from H2O-forming NADH oxidase (NOX) activity of two distinct flavin-containing systems: H2O-forming NOXes and multicomponent flavodiiron proteins (FDPs). Neither system is membrane bound, and both recycle NADH into oxidized NAD+ while simultaneously removing O2 from the local environment. However, little is known about the specific contributions of these systems in T. vaginalis. In this study, we use bioinformatics and biochemical analyses to show that T. vaginalis lacks a NOX-like enzyme and instead harbors three paralogous genes (FDPF1-3), each encoding a natural fusion product between the N-terminal FDP, central rubredoxin (Rb), and C-terminal NADH:Rb oxidoreductase domains. Unlike a "stand-alone" FDP that lacks Rb and oxidoreductase domains, this natural fusion protein with fully populated flavin redox centers directly accepts reducing equivalents of NADH to catalyze the four-electron reduction of oxygen to water within a single polypeptide with an extremely high turnover. Furthermore, using single-particle cryo-EM, we present structural insights into the spatial organization of the FDP core within this multidomain fusion protein. Together, these results contribute to our understanding of systems that allow protozoan parasites to maintain optimal redox balance and survive transient exposure to oxic conditions.

Small Proteins in Archaea, a Mainly Unexplored World

In recent years, increasing numbers of small proteins have moved into the focus of science. Small proteins have been identified and characterized in all three domains of life, but the majority remains functionally uncharacterized, lack secondary structure, and exhibit limited evolutionary conservation. While quite a few have already been described for bacteria and eukaryotic organisms, the amount of known and functionally analyzed archaeal small proteins is still very limited. In this review, we compile the current state of research, show strategies for systematic approaches for global identification of small archaeal proteins, and address selected functionally characterized examples. Besides, we document exemplarily for one archaeon the tool development and optimization to identify small proteins using genome-wide approaches.

The Pyrococcus furiosus ironome is dominated by [Fe4S4]2+ clusters or thioferrate-like iron depending on the availability of elemental sulfur

Pyrococcus furiosus is a hyperthermophilic anaerobic archaeon whose metabolism depends on whether elemental sulfur is (+S⁰) or is not (-S⁰) included in growth medium. Under +S⁰ conditions, expression of respiratory hydrogenase declines while respiratory membrane-bound sulfane reductase and the putative iron-storage protein IssA increase. Our objective was to investigate the iron content of WT and ΔIssA cells under these growth conditions using Mössbauer spectroscopy. WT-S⁰ cells contained ∼1 mM Fe, with ∼85% present as two spectroscopically distinct forms of S = 0 [Fe₄S₄]²⁺ clusters; the remainder was mainly high-spin Fe^II. WT+S⁰ cells contained 5 to 9 mM Fe, with 75 to 90% present as magnetically ordered thioferrate-like (TFL) iron nanoparticles. TFL iron was similar to chemically defined thioferrates; both consisted of Fe^III ions coordinated by an S₄ environment, and both exhibited strong coupling between particles causing high applied fields to have little spectral effect. At high temperatures with magnetic hyperfine interactions abolished, TFL iron exhibited two doublets overlapping those of [Fe₄S₄]²⁺ clusters in -S⁰ cells. This coincidence arose because of similar coordination environments of TFL iron and cluster iron. The TFL structure was more heterogeneous in the presence of IssA. Presented data suggest that IssA may coordinate insoluble iron sulfides as TFL iron, formed as a byproduct of anaerobic sulfur respiration under high iron conditions, which thereby reduces its toxicity to the cell. This was the first Mössbauer characterization of the ironome of an archaeon, and it illustrates differences relative to the iron content of better-studied bacteria such as Escherichia coli.

Streamlined purification and characterization of Pyrococcus furiosus rubredoxins with different N-terminal modifications by reversed-phase HPLC

Rubredoxins (Rds), like those from Pyrococcus furious (Pf), have largely been found to be expressed in Escherichia coli (E. coli) as a mixture of different N-terminal forms, which may affect the properties of the protein. The typical procedures for the purification of Rds are cumbersome and usually with low yield. We present herein a streamlined purification strategy based on the reversed-phase high performance liquid chromatography (RP-HPLC), which offers high yield and high resolution after simply one-step purification following pre-treatment. We also show that RP-HPLC can be a valuable tool to gain information related to the thermal decomposition pathway of Pf-Rds.

Catalysis and Electron Transfer in De Novo Designed Helical Scaffolds

The relationship between protein structure and function is one of the greatest puzzles within biochemistry. De novo metalloprotein design is a way to wipe the board clean and determine what is required to build in function from the ground up in an unrelated structure. This Review focuses on protein design efforts to create de novo metalloproteins within alpha-helical scaffolds. Examples of successful designs include those with carbonic anhydrase or nitrite reductase activity by incorporating a ZnHis3 or CuHis3 site, or that recapitulate the spectroscopic properties of unique electron-transfer sites in cupredoxins (CuHis2 Cys) or rubredoxins (FeCys4 ). This work showcases the versatility of alpha helices as scaffolds for metalloprotein design and the progress that is possible through careful rational design. Our studies cover the invariance of carbonic anhydrase activity with different site positions and scaffolds, refinement of our cupredoxin models, and enhancement of nitrite reductase activity up to 1000-fold.

Characterization of the metal status of natively purified alpha-synuclein from human blood, brain tissue, or recombinant sources using size exclusion ICP-MS reveals no significant binding of Cu, Fe or Zn

The binding of Cu, Fe or Zn to alpha-synuclein has been implicated in neurodegenerative disease, such as Parkinson's.

Development of a Rubredoxin-Type Center Embedded in a de Dovo-Designed Three-Helix Bundle

Protein design is a powerful tool for interrogating the basic requirements for the function of a metal site in a way that allows for the selective incorporation of elements that are important for function. Rubredoxins are small electron transfer proteins with a reduction potential centered near 0 mV (vs normal hydrogen electrode). All previous attempts to design a rubredoxin site have focused on incorporating the canonical CXXC motifs in addition to reproducing the peptide fold or using flexible loop regions to define the morphology of the site. We have produced a rubredoxin site in an utterly different fold, a three-helix bundle. The spectra of this construct mimic the ultraviolet-visible, Mössbauer, electron paramagnetic resonance, and magnetic circular dichroism spectra of native rubredoxin. Furthermore, the measured reduction potential suggests that this rubredoxin analogue could function similarly. Thus, we have shown that an α-helical scaffold sustains a rubredoxin site that can cycle with the desired potential between the Fe(II) and Fe(III) states and reproduces the spectroscopic characteristics of this electron transport protein without requiring the classic rubredoxin protein fold.

Cluster-Dependent Charge-Transfer Dynamics in Iron-Sulfur Proteins

Photoinduced charge-transfer dynamics and the influence of cluster size on the dynamics were investigated using five iron-sulfur clusters: the 1Fe-4S cluster in Pyrococcus furiosus rubredoxin, the 2Fe-2S cluster in Pseudomonas putida putidaredoxin, the 4Fe-4S cluster in nitrogenase iron protein, and the 8Fe-7S P-cluster and the 7Fe-9S-1Mo FeMo cofactor in nitrogenase MoFe protein. Laser excitation promotes the iron-sulfur clusters to excited electronic states that relax to lower states. The electronic relaxation lifetimes of the 1Fe-4S, 8Fe-7S, and 7Fe-9S-1Mo clusters are on the picosecond time scale, although the dynamics of the MoFe protein is a mixture of the dynamics of the latter two clusters. The lifetimes of the 2Fe-2S and 4Fe-4S clusters, however, extend to several nanoseconds. A competition between reorganization energies and the density of electronic states (thus electronic coupling between states) mediates the charge-transfer lifetimes, with the 2Fe-2S cluster of Pdx and the 4Fe-4S cluster of Fe protein lying at the optimum leading to them having significantly longer lifetimes. Their long lifetimes make them the optimal candidates for long-range electron transfer and as external photosensitizers for other photoactivated chemical reactions like solar hydrogen production. Potential electron-transfer and hole-transfer pathways that possibly facilitate these charge transfers are proposed.

Enhanced heat tolerance in transgenic silkworm via overexpression of Pyrococcus furiosus superoxide reductase

Heat shock causes a serious harm to organisms by accelerating the production of reactive oxygen species (ROS). Pyrococcus furiosus superoxide reductase (PfSOR) is an enzyme that efficiently detoxifies ROS. In order to generate a silkworm strain with high heat tolerance for sericulture, we synthesized an artificial DNA sequence encoding PfSOR based on the codon bias of Bombyx mori. PfSOR was successfully overexpressed in transgenic silkworm (named A4SOR) and BmE cells, as determined by RT-PCR and western blot analyses. An SOR activity assay confirmed that the expressed enzyme was functional in A4SOR. After exposure to a temperature of 35 °C for 44 h, the mortality rate was about 30% lower in transgenic A4SOR than in non-transgenic silkworms. Moreover, transgene expression had no apparent effect on economic characteristics of silkworms. The heat tolerance of silkworm was thus enhanced by expressing an archaeal SOR; this can be useful for sericulture in regions where the average temperature exceeds the optimal environmental temperature for B. mori of 25 °C.

Macromolecular structure phasing by neutron anomalous diffraction

In this report we show for the first time that neutron anomalous dispersion can be used in a practical manner to determine experimental phases of a protein crystal structure, providing a new tool for structural biologists. The approach is demonstrated through the use of a state-of-the-art monochromatic neutron diffractometer at the Institut Laue-Langevin (ILL) in combination with crystals of perdeuterated protein that minimise the level of hydrogen incoherent scattering and enhance the visibility of the anomalous signal. The protein used was rubredoxin in which cadmium replaced the iron at the iron-sulphur site. While this study was carried out using a steady-state neutron beam source, the results will be of major interest for capabilities at existing and emerging spallation neutron sources where time-of-flight instruments provide inherent energy discrimination. In particular this capability may be expected to offer unique opportunities to a rapidly developing structural biology community where there is increasing interest in the identification of protonation states, protein/water interactions and protein-ligand interactions - all of which are of central importance to a wide range of fundamental and applied areas in the biosciences.

Synchrotron-based Nickel Mössbauer Spectroscopy

We used a novel experimental setup to conduct the first synchrotron-based (61)Ni Mössbauer spectroscopy measurements in the energy domain on Ni coordination complexes and metalloproteins. A representative set of samples was chosen to demonstrate the potential of this approach. (61)NiCr2O4 was examined as a case with strong Zeeman splittings. Simulations of the spectra yielded an internal magnetic field of 44.6 T, consistent with previous work by the traditional (61)Ni Mössbauer approach with a radioactive source. A linear Ni amido complex, (61)Ni{N(SiMe3)Dipp}2, where Dipp = C6H3-2,6-(i)Pr2, was chosen as a sample with an "extreme" geometry and large quadrupole splitting. Finally, to demonstrate the feasibility of metalloprotein studies using synchrotron-based (61)Ni Mössbauer spectroscopy, we examined the spectra of (61)Ni-substituted rubredoxin in reduced and oxidized forms, along with [Et4N]2[(61)Ni(SPh)4] as a model compound. For each of the above samples, a reasonable spectrum could be obtained in ∼1 d. Given that there is still room for considerable improvement in experimental sensitivity, synchrotron-based (61)Ni Mössbauer spectroscopy appears to be a promising alternative to measurements with radioactive sources.

Nickel-Substituted Rubredoxin as a Minimal Enzyme Model for Hydrogenase

A simple, functional mimic of [NiFe] hydrogenases based on a nickel-substituted rubredoxin (NiRd) protein is reported. NiRd is capable of light-initiated and solution-phase hydrogen production and demonstrates high electrocatalytic activity using protein film voltammetry. The catalytic voltammograms are modeled using analytical expressions developed for hydrogenase enzymes, revealing maximum turnover frequencies of approximately 20-100 s(-1) at 4 °C with an overpotential of 540 mV. These rates are directly comparable to those observed for [NiFe] hydrogenases under similar conditions. Like the native enzymes, the proton reduction activity of NiRd is strongly inhibited by carbon monoxide. This engineered rubredoxin-based enzyme is chemically and thermally robust, easily accessible, and highly tunable. These results have implications for understanding the enzymatic mechanisms of native hydrogenases, and, using NiRd as a scaffold, it will be possible to optimize this catalyst for application in sustainable fuel generation.

Nuclear resonance vibrational spectroscopy (NRVS) of rubredoxin and MoFe protein crystals

We have applied ⁵⁷Fe nuclear resonance vibrational spectroscopy (NRVS) for the first time to study the dynamics of Fe centers in Fe-S protein crystals, including oxidized wild type rubredoxin crystals from Pyrococcus furiosus, and the MoFe protein of nitrogenase from Azotobacter vinelandii. Thanks to the NRVS selection rule, selectively probed vibrational modes have been observed in both oriented rubredoxin and MoFe protein crystals. The NRVS work was complemented by extended X-ray absorption fine structure spectroscopy (EXAFS) measurements on oxidized wild type rubredoxin crystals from Pyrococcus furiosus. The EXAFS spectra revealed the Fe-S bond length difference in oxidized Pf Rd protein, which is qualitatively consistent with the X-ray crystal structure.

Rapid visualization of hydrogen positions in protein neutron crystallographic structures

Neutron crystallography is a powerful technique for experimental visualization of the positions of light atoms, including hydrogen and its isotope deuterium. In recent years, structural biologists have shown increasing interest in the technique as it uniquely complements X-ray crystallographic data by revealing the positions of D atoms in macromolecules. With this regained interest, access to macromolecular neutron crystallography beamlines is becoming a limiting step. In this report, it is shown that a rapid data-collection strategy can be a valuable alternative to longer data-collection times in appropriate cases. Comparison of perdeuterated rubredoxin structures refined against neutron data sets collected over hours and up to 5 d shows that rapid neutron data collection in just 14 h is sufficient to provide the positions of 269 D atoms without ambiguity.

A computational framework for proteome-wide pursuit and prediction of metalloproteins using ICP-MS and MS/MS data

Background

Metal-containing proteins comprise a diverse and sizable category within the proteomes of organisms, ranging from proteins that use metals to catalyze reactions to proteins in which metals play key structural roles. Unfortunately, reliably predicting that a protein will contain a specific metal from its amino acid sequence is not currently possible. We recently developed a generally-applicable experimental technique for finding metalloproteins on a genome-wide scale. Applying this metal-directed protein purification approach (ICP-MS and MS/MS based) to the prototypical microbe Pyrococcus furiosus conclusively demonstrated the extent and diversity of the uncharacterized portion of microbial metalloproteomes since a majority of the observed metal peaks could not be assigned to known or predicted metalloproteins. However, even using this technique, it is not technically feasible to purify to homogeneity all metalloproteins in an organism. In order to address these limitations and complement the metal-directed protein purification, we developed a computational infrastructure and statistical methodology to aid in the pursuit and identification of novel metalloproteins.

Results

We demonstrate that our methodology enables predictions of metal-protein interactions using an experimental data set derived from a chromatography fractionation experiment in which 870 proteins and 10 metals were measured over 2,589 fractions. For each of the 10 metals, cobalt, iron, manganese, molybdenum, nickel, lead, tungsten, uranium, vanadium, and zinc, clusters of proteins frequently occurring in metal peaks (of a specific metal) within the fractionation space were defined. This resulted in predictions that there are from 5 undiscovered vanadium- to 13 undiscovered cobalt-containing proteins in Pyrococcus furiosus. Molybdenum and nickel were chosen for additional assessment producing lists of genes predicted to encode metalloproteins or metalloprotein subunits, 22 for nickel including seven from known nickel-proteins, and 20 for molybdenum including two from known molybdo-proteins. The uncharacterized proteins are prime candidates for metal-based purification or recombinant approaches to validate these predictions.

Conclusions

We conclude that the largely uncharacterized extent of native metalloproteomes can be revealed through analysis of the co-occurrence of metals and proteins across a fractionation space. This can significantly impact our understanding of metallobiochemistry, disease mechanisms, and metal toxicity, with implications for bioremediation, medicine and other fields.

Microbial metalloproteomes are largely uncharacterized

Metal ion cofactors afford proteins virtually unlimited catalytic potential, enable electron transfer reactions and have a great impact on protein stability. Consequently, metalloproteins have key roles in most biological processes, including respiration (iron and copper), photosynthesis (manganese) and drug metabolism (iron). Yet, predicting from genome sequence the numbers and types of metal an organism assimilates from its environment or uses in its metalloproteome is currently impossible because metal coordination sites are diverse and poorly recognized. We present here a robust, metal-based approach to determine all metals an organism assimilates and identify its metalloproteins on a genome-wide scale. This shifts the focus from classical protein-based purification to metal-based identification and purification by liquid chromatography, high-throughput tandem mass spectrometry (HT-MS/MS) and inductively coupled plasma mass spectrometry (ICP-MS) to characterize cytoplasmic metalloproteins from an exemplary microorganism (Pyrococcus furiosus). Of 343 metal peaks in chromatography fractions, 158 did not match any predicted metalloprotein. Unassigned peaks included metals known to be used (cobalt, iron, nickel, tungsten and zinc; 83 peaks) plus metals the organism was not thought to assimilate (lead, manganese, molybdenum, uranium and vanadium; 75 peaks). Purification of eight of 158 unexpected metal peaks yielded four novel nickel- and molybdenum-containing proteins, whereas four purified proteins contained sub-stoichiometric amounts of misincorporated lead and uranium. Analyses of two additional microorganisms (Escherichia coli and Sulfolobus solfataricus) revealed species-specific assimilation of yet more unexpected metals. Metalloproteomes are therefore much more extensive and diverse than previously recognized, and promise to provide key insights for cell biology, microbial growth and toxicity mechanisms.

Unambiguous determination of H-atom positions: comparing results from neutron and high-resolution X-ray crystallography

The locations of H atoms in biological structures can be difficult to determine using X-ray diffraction methods. Neutron diffraction offers a relatively greater scattering magnitude from H and D atoms. Here, 1.65 A resolution neutron diffraction studies of fully perdeuterated and selectively CH(3)-protonated perdeuterated crystals of Pyrococcus furiosus rubredoxin (D-rubredoxin and HD-rubredoxin, respectively) at room temperature (RT) are described, as well as 1.1 A resolution X-ray diffraction studies of the same protein at both RT and 100 K. The two techniques are quantitatively compared in terms of their power to directly provide atomic positions for D atoms and analyze the role played by atomic thermal motion by computing the sigma level at the D-atom coordinate in simulated-annealing composite D-OMIT maps. It is shown that 1.65 A resolution RT neutron data for perdeuterated rubredoxin are approximately 8 times more likely overall to provide high-confidence positions for D atoms than 1.1 A resolution X-ray data at 100 K or RT. At or above the 1.0sigma level, the joint X-ray/neutron (XN) structures define 342/378 (90%) and 291/365 (80%) of the D-atom positions for D-rubredoxin and HD-rubredoxin, respectively. The X-ray-only 1.1 A resolution 100 K structures determine only 19/388 (5%) and 8/388 (2%) of the D-atom positions above the 1.0sigma level for D-rubredoxin and HD-rubredoxin, respectively. Furthermore, the improved model obtained from joint XN refinement yielded improved electron-density maps, permitting the location of more D atoms than electron-density maps from models refined against X-ray data only.

Expression of Pyrococcus furiosus superoxide reductase in Arabidopsis enhances heat tolerance

Plants produce reactive oxygen species (ROS) in response to environmental stresses sending signaling cues, which, if uncontrolled, result in cell death. Like other aerobic organisms, plants have ROS-scavenging enzymes, such as superoxide dismutase (SOD), which removes superoxide anion radical (O(2)(-)) and prevents the production and buildup of toxic free radicals. However, increasing the expression of cytosolic SODs is complex, and increasing their production in vivo has proven to be challenging. To avoid problems with endogenous regulation of gene expression, we expressed a gene from the archaeal hyperthermophile Pyrococcus furiosus that reduces O(2)(-). P. furiosus uses superoxide reductase (SOR) rather than SOD to remove superoxide. SOR is a thermostable enzyme that reduces O(2)(-) in a one-electron reduction without producing oxygen. We show that P. furiosus SOR can be produced as a functional enzyme in planta and that plants producing SOR have enhanced tolerance to heat, light, and chemically induced ROS. Stress tolerance in the SOR-producing plants correlates positively with a delayed increase in ROS-sensitive transcripts and a decrease in ascorbate peroxidase activity. The SOR plants provide a good model system to study the impact of cytosolic ROS on downstream signaling in plant growth and development. Furthermore, this work demonstrates that this synthetic approach for reducing cytosolic ROS holds promise as a means for improving stress tolerance in crop plants.

Robust, high-throughput solution structural analyses by small angle X-ray scattering (SAXS)

We present an efficient pipeline enabling high-throughput analysis of protein structure in solution with small angle X-ray scattering (SAXS). Our SAXS pipeline combines automated sample handling of microliter volumes, temperature and anaerobic control, rapid data collection, data analysis, and couples structural analysis with automated archiving. We subjected 50 representative proteins, mostly from Pyrococcus furiosus, to this pipeline, revealing that 30 were multimeric structures in solution. SAXS analysis allowed us to distinguish aggregated and unfolded proteins, define global structural parameters and oligomeric states for most samples, identify shapes and similar structures for 25 unknown structures, and determine envelopes for 41 proteins. We believe that high throughput SAXS is an enabling technology that may change the way that structural genomics research is done.

Characterization of an NADH oxidase of the flavin-dependent disulfide reductase family from Methanocaldococcus jannaschii

Methanocaldococcus jannaschii, a deeply rooted hyperthermophilic anaerobic methanarchaeon from a deep-sea hydrothermal vent, carries an NADH oxidase (Nox) homologue (MJ0649). According to the characteristics described here, MJ0649 represents an unusual member within group 3 of the flavin-dependent disulfide reductase (FDR) family. This FDR group comprises Nox, NADH peroxidases (Npx) and coenzyme A disulfide reductases (CoADRs); each carries a Cys residue that forms Cys-sulfenic acid during catalysis. A sequence analysis identified MJ0649 as a CoADR homologue. However, recombinant MJ0649 (rMJNox), expressed in Escherichia coli and purified to homogeneity an 86 kDa homodimer with 0.27 mol FAD (mol subunit)−1, showed Nox but not CoADR activity. Incubation with FAD increased FAD content to 1 mol (mol subunit)−1 and improved NADH oxidase activity 3.4-fold. The FAD-incubated enzyme was characterized further. The optimum pH and temperature were ≥10 and ≥95 °C, respectively. At pH 7 and 83 °C, apparent K m values for NADH and O2 were 3 μM and 1.9 mM, respectively, and the specific activity at 1.4 mM O2 was 60 μmol min−1 mg−1; 62 % of NADH-derived reducing equivalents were recovered as H2O2 and the rest probably generated H2O. rMjNox had poor NADPH oxidase, NADH peroxidase and superoxide formation activities. It reduced ferricyanide, plumbagin and 5,5′-dithiobis(2-nitrobenzoic acid), but not disulfide coenzyme A and disulfide coenzyme M. Due to a high K m, O2 is not a physiologically relevant substrate for MJ0649; its true substrate remains unknown.

Identifying stabilizing key residues in proteins using interresidue interaction energy matrix

We are proposing an interresidue interaction energy map (IEM)--a new tool for protein structure analysis and protein bioinformatics. This approach employs the sum of pair-wise interaction energies of a particular residue as a measure of its structural importance. We will show that the IEM can serve as a means for identifying key residues responsible for the stability of a protein. Our method can be compared with the interresidue contact map but has the advantage of weighting the contacts by the stabilization energy content which they bring to the protein structure. For the theoretical adjustment of the proposed method, we chose the Trp-cage mini protein as a model system to compare a spectrum of computational methods ranging from the ab initio MP2 level through the DFT method to empirical force-field methods. The IEM method correctly identifies Tryptophane 6 as the key residue in the Trp-cage. The other residues with the highest stabilizing contributions correspond to the structurally important positions in the protein. We have further tested our method on the Trp2Cage miniprotein--a P12W mutant of the Trp-cage and on two proteins from the rubredoxin family that differ in their thermostability. Our method correctly identified the thermodynamically more stable variants in both cases and therefore can also be used as a tool for the relative measurement of protein stability. Finally, we will point out the important role played by dispersion energy, which contributes significantly to the total stabilization energy and whose role in aromatic pairs is clearly dominant. Surprisingly, the dispersion energy plays an even more important role in the interaction of prolines with aromatic systems.

Preliminary neutron crystallographic analysis of selectively CH3-protonated deuterated rubredoxin from Pyrococcus furiosus

Neutron crystallography is used to locate H atoms in biological materials and can distinguish between negatively scattering hydrogen-substituted and positively scattering deuterium-substituted positions in isomorphous neutron structures. Recently, Hauptman & Langs (2003; Acta Cryst. A59, 250-254) have shown that neutron diffraction data can be used to solve macromolecular structures by direct methods and that solution is aided by the presence of negatively scattering H atoms in the structure. Selective-labeling protocols allow the design and production of H/D-labeled macromolecular structures in which the ratio of H to D atoms can be precisely controlled. Methyl selective-labeling protocols were applied to introduce (1H-delta methyl)-leucine and (1H-gamma methyl)-valine into deuterated rubredoxin from Pyrococcus furiosus (PfRd). Here, the production, crystallization and preliminary neutron analysis of a selectively CH3-protonated deuterated PfRd sample, which provided a high-quality neutron data set that extended to 1.75 A resolution using the new LADI-III instrument at the Institut Laue-Langevin, are reported. Preliminary analysis of neutron density maps allows unambiguous assignment of the positions of H atoms at the methyl groups of the valine and leucine residues in the otherwise deuterated rubredoxin structure.

Dispersion interactions govern the strong thermal stability of a protein

Rubredoxin from the hyperthermophile Pyrococcus furiosus (Pf Rd) is an extremely thermostable protein, which makes it an attractive subject of protein folding and stability studies. A fundamental question arises as to what the reason for such extreme stability is and how it can be elucidated from a complex set of interatomic interactions. We addressed this issue first theoretically through a computational analysis of the hydrophobic core of the protein and its mutants, including the interactions taking place inside the core. Here we show that a single mutation of one of phenylalanine's residues inside the protein's hydrophobic core results in a dramatic decrease in its thermal stability. The calculated unfolding Gibbs energy as well as the stabilization energy differences between a few core residues follows the same trend as the melting temperature of protein variants determined experimentally by microcalorimetry measurements. NMR spectroscopy experiments have shown that the only part of the protein affected by mutation is the reasonably rearranged hydrophobic core. It is hence concluded that stabilization energies, which are dominated by London dispersion, represent the main source of stability of this protein.

Biochemical and structural characterization of a novel family of cystathionine beta-synthase domain proteins fused to a Zn ribbon-like domain

We have identified a novel family of proteins, in which the N-terminal cystathionine beta-synthase (CBS) domain is fused to the C-terminal Zn ribbon domain. Four proteins were overexpressed in Escherichia coli and purified: TA0289 from Thermoplasma acidophilum, TV1335 from Thermoplasma volcanium, PF1953 from Pyrococcus furiosus, and PH0267 from Pyrococcus horikoshii. The purified proteins had a red/purple color in solution and an absorption spectrum typical of rubredoxins (Rds). Metal analysis of purified proteins revealed the presence of several metals, with iron and zinc being the most abundant metals (2-67% of iron and 12-74% of zinc). Crystal structures of both mercury- and iron-bound TA0289 (1.5-2.0 A resolution) revealed a dimeric protein whose intersubunit contacts are formed exclusively by the alpha-helices of two cystathionine beta-synthase subdomains, whereas the C-terminal domain has a classical Zn ribbon planar architecture. All proteins were reversibly reduced by chemical reductants (ascorbate or dithionite) or by the general Rd reductase NorW from E. coli in the presence of NADH. Reduced TA0289 was found to be capable of transferring electrons to cytochrome C from horse heart. Likewise, the purified Zn ribbon protein KTI11 from Saccharomyces cerevisiae had a purple color in solution and an Rd-like absorption spectrum, contained both iron and zinc, and was reduced by the Rd reductase NorW from E. coli. Thus, recombinant Zn ribbon domains from archaea and yeast demonstrate an Rd-like electron carrier activity in vitro. We suggest that, in vivo, some Zn ribbon domains might also bind iron and therefore possess an electron carrier activity, adding another physiological role to this large family of important proteins.

Observation of terahertz vibrations in Pyrococcus furiosus rubredoxin via impulsive coherent vibrational spectroscopy and nuclear resonance vibrational spectroscopy – interpretation by molecular mechanics

We have used impulsive coherent vibrational spectroscopy (ICVS) to study the Fe(S-Cys)(4) site in oxidized rubredoxin (Rd) from Pyrococcus furiosus (Pf). In this experiment, a 15 fs visible laser pulse is used to coherently pump the sample to an excited electronic state, and a second <10 fs pulse is used to probe the change in transmission as a function of the time delay. PfRd was observed to relax to the ground state by a single exponential decay with time constants of approximately 255-275 fs. Superimposed on this relaxation are oscillations caused by coherent excitation of vibrational modes in both excited and ground electronic states. Fourier transformation reveals the frequencies of these modes. The strongest ICV mode with 570 nm excitation is the symmetric Fe-S stretching mode near 310 cm(-1), compared to 313 cm(-1) in the low temperature resonance Raman. If the rubredoxin is pumped at 520 nm, a set of strong bands occurs between 20 and 110 cm(-1). Finally, there is a mode at approximately 500 cm(-1) which is similar to features near 508 cm(-1) in blue Cu proteins that have been attributed to excited state vibrations. Normal mode analysis using 488 protein atoms and 558 waters gave calculated spectra that are in good agreement with previous nuclear resonance vibrational spectra (NRVS) results. The lowest frequency normal modes are identified as collective motions of the entire protein or large segments of polypeptide. Motion in these modes may affect the polar environment of the redox site and thus tune the electron transfer functions in rubredoxins.

Production of a thermostable archaeal superoxide reductase in plant cells

Pyrococcus furiosus superoxide reductase (SOR) is a thermostable archaeal enzyme that reduces superoxide without producing oxygen. When produced as a fusion protein with the green fluorescent protein in plant cells, P. furiosus SOR is located in the cytosol and nucleus. The recombinant SOR enzyme retains its function and heat stability when assayed in vitro. Importantly, expressing SOR in plant cells enhances their survival at high temperature indicating that it functions in vivo. The archaeal SOR provides a novel mechanism to reduce superoxide and demonstrates the potential for using archaeal genes to alter eukaryotic metabolism.

In vitro reconstitution of an NADPH-dependent superoxide reduction pathway from Pyrococcus furiosus

A scheme for the detoxification of superoxide in Pyrococcus furiosus has been previously proposed in which superoxide reductase (SOR) reduces (rather than dismutates) superoxide to hydrogen peroxide by using electrons from reduced rubredoxin (Rd). Rd is reduced with electrons from NAD(P)H by the enzyme NAD(P)H:rubredoxin oxidoreductase (NROR). The goal of the present work was to reconstitute this pathway in vitro using recombinant enzymes. While recombinant forms of SOR and Rd are available, the gene encoding P. furiosus NROR (PF1197) was found to be exceedingly toxic to Escherichia coli, and an active recombinant form (rNROR) was obtained via a fusion protein expression system, which produced an inactive form of NROR until cleavage. This allowed the complete pathway from NAD(P)H to the reduction of SOR via NROR and Rd to be reconstituted in vitro using recombinant proteins. rNROR is a 39.9-kDa protein whose sequence contains both flavin adenine dinucleotide (FAD)- and NAD(P)H-binding motifs, and it shares significant similarity with known and putative Rd-dependent oxidoreductases from several anaerobic bacteria, both mesophilic and hyperthermophilic. FAD was shown to be essential for activity in reconstitution assays and could not be replaced by flavin mononucleotide (FMN). The bound FAD has a midpoint potential of -173 mV at 23 degrees C (-193 mV at 80 degrees C). Like native NROR, the recombinant enzyme catalyzed the NADPH-dependent reduction of rubredoxin both at high (80 degrees C) and low (23 degrees C) temperatures, consistent with its proposed role in the superoxide reduction pathway. This is the first demonstration of in vitro superoxide reduction to hydrogen peroxide using NAD(P)H as the electron donor in an SOR-mediated pathway.

Increased peptide deformylase activity for N-formylmethionine processing of proteins overexpressed in Escherichia coli: application to homogeneous rubredoxin production

Deformylation of the initiator N-formylmethionine does not always proceed to completion for proteins overexpressed in Escherichia coli. To overcome this limitation, the def gene encoding the Escherichia coli peptide deformylase was cloned into the plysS plasmid under the tetracycline (Tc) promoter control. The efficiency of this constitutive level of peptide deformylase expression was demonstrated for the case of the rubredoxins from both mesophilic and hyperthermophilic organisms which normally retain a majority of their N-formyl terminal form. Indicating the potential structural/functional significance of residual formylation, the presence of a highly solvent exposed N-formyl group in rubredoxin is discernable in the amide NMR chemical shifts for the active site metal-coordinating cysteines more than 21A away.

Rubrerythrin from the hyperthermophilic archaeon Pyrococcus furiosus is a rubredoxin-dependent, iron-containing peroxidase

Rubrerythrin was purified by multistep chromatography under anaerobic, reducing conditions from the hyperthermophilic archaeon Pyrococcus furiosus. It is a homodimer with a molecular mass of 39.2 kDa and contains 2.9 +/- 0.2 iron atoms per subunit. The purified protein had peroxidase activity at 85 degrees C using hydrogen peroxide with reduced P. furiosus rubredoxin as the electron donor. The specific activity was 36 micromol of rubredoxin oxidized/min/mg with apparent K(m) values of 35 and 70 microM for hydrogen peroxide and rubredoxin, respectively. When rubrerythrin was combined with rubredoxin and P. furiosus NADH:rubredoxin oxidoreductase, the complete system used NADH as the electron donor to reduce hydrogen peroxide with a specific activity of 7.0 micromol of H(2)O(2) reduced/min/mg of rubrerythrin at 85 degrees C. Strangely, as-purified (reduced) rubrerythrin precipitated when oxidized by either hydrogen peroxide, air, or ferricyanide. The gene (PF1283) encoding rubrerythrin was expressed in Escherichia coli grown in medium with various metal contents. The purified recombinant proteins each contained approximately three metal atoms/subunit, ranging from 0.4 Fe plus 2.2 Zn to 1.9 Fe plus 1.2 Zn, where the metal content of the protein depended on the metal content of the E. coli growth medium. The peroxidase activities of the recombinant forms were proportional to the iron content. P. furiosus rubrerythrin is the first to be characterized from a hyperthermophile or from an archaeon, and the results are the first demonstration that this protein functions in an NADH-dependent, hydrogen peroxide:rubredoxin oxidoreductase system. Rubrerythrin is proposed to play a role in the recently defined anaerobic detoxification pathway for reactive oxygen species.

Neutron crystallographic study on rubredoxin from Pyrococcus furiosus by BIX-3, a single-crystal diffractometer for biomacromolecules

The structure of a partially deuterated rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus, an organism that grows optimally at 100 degrees C, was determined by using the neutron single-crystal diffractometer dedicated for biological macromolecules (BIX-3) at the JRR-3M reactor of the Japan Atomic Energy Research Institute. Data were collected at room temperature up to a resolution of 1.5 A, and the completeness factor of the data set was 81.9%. The model contains 306 H and 50 D atoms. A total of 37 hydration water molecules were identified, with 15 having all three atoms fully located and the remaining D2O molecules partially defined. The model has been refined to final agreement factors of R = 18.6% and Rfree = 21.7%. Several orientations of the O-D bonds of side chains, whose assignments from x-ray data were previously ambiguous, were clearly visible in the neutron structure. Although most backbone N-H bonds had undergone some degree of H/D exchange throughout the rubredoxin molecule, 5 H atom positions still had distinctly negative (H) peaks. The neutron Fourier maps clearly showed the details of an extensive set of H bonds involving the ND3+ terminus that may contribute to the unusual thermostability of this molecule.

Hydration in proteins observed by high-resolution neutron crystallography

It is well known that water molecules surrounding a protein play important roles in maintaining its structural stability. Water molecules are known to participate in several physiological processes through the formation of hydrogen bonds. However, the hydration structures of most proteins are not known well at an atomic level at present because X-ray protein crystallography has difficulties to localize hydrogen atoms. In contrast, neutron crystallography has no problem in determining the position of hydrogens with high accuracy.1 In this article, the hydration structures of three proteins are described- myoglobin, wild-type rubredoxin, and a mutant rubredoxin-the structures of which were solved at 1.5- or 1.6-A resolution by neutron structure determination. These hydration patterns show fascinating features and the water molecules adopt a variety of shapes in the neutron Fourier maps, revealing details of intermolecular hydrogen bond formation and dynamics of hydration. Our results further show that there are strong relationships between these shapes and the water environments.

A solution NMR molecular model for the aspartate-ligated, cubane cluster containing ferredoxin from the hyperthermophilic archeaon Pyrococcus furiosus

A solution molecular model for the conformationally dynamically heterogeneous Pyrococcus furiosus ferredoxin with an intact disulfide bond has been constructed on the basis of reported (1)H NMR spectral parameters using distance geometry and simulated annealing protocols. Conventional long-mixing time NOESY and H-bonding constraints have been augmented by previously reported short-mixing time NOESY, steady-state NOE, and cluster paramagnetism-induced relaxation. The family of 15 structures with inconsequential violations exhibited low rms deviations for backbone atoms for the overwhelming majority of the residues, including the cluster ligating loop with the unprecedented ligated Asp14. Larger rms deviations were observed across the disulfide bond, but closer inspection revealed that the 15 structures can be factored into 10 substructures exhibiting an "S" or right-handed disulfide orientation and 5 exhibiting an "R" or left-handed disulfide orientation. The remainder of the structure is indistinguishable for the two disulfide orientations but confirms stabilizing extensions of secondary structural elements in the lengthening of the long helix and both the lengthening and incorporation of a third strand into the beta-sheet involving the termini, with these extensions interacting strongly in a modular fashion through the rings of Tyr46 and Trp2. These extensions of stabilizing interactions in Pyrococcus furiosus Fd, however, lead to strong destabilization of the disulfide bond and destabilization of the highly conserved first and last beta-turns in the sequence. It is concluded that the structural alternations in Pyrococcus Fd relative to other hyperthermostable Fds are not to increase thermostability but to place "stress" on the disulfide bond and render it more reducible. The possible physiological implications of this unique reducible disulfide bond are discussed.