Leucine 41 is a gate for water entry in the reduction of Clostridium pasteurianum rubredoxin

Biological electron transfer is an efficient process even though the distances between the redox moieties are often quite large. It is therefore of great interest to gain an understanding of the physical basis of the rates and driving forces of these reactions. The structural relaxation of the protein that occurs upon change in redox state gives rise to the reorganizational energy, which is important in the rates and the driving forces of the proteins involved. To determine the structural relaxation in a redox protein, we have developed methods to hold a redox protein in its final oxidation state during crystallization while maintaining the same pH and salt conditions of the crystallization of the protein in its initial oxidation state. Based on 1.5 A resolution crystal structures and molecular dynamics simulations of oxidized and reduced rubredoxins (Rd) from Clostridium pasteurianum (Cp), the structural rearrangements upon reduction suggest specific mechanisms by which electron transfer reactions of rubredoxin should be facilitated. First, expansion of the [Fe-S] cluster and concomitant contraction of the NH...S hydrogen bonds lead to greater electrostatic stabilization of the extra negative charge. Second, a gating mechanism caused by the conformational change of Leucine 41, a nonpolar side chain, allows transient penetration of water molecules, which greatly increases the polarity of the redox site environment and also provides a source of protons. Our method of producing crystals of Cp Rd from a reducing solution leads to a distribution of water molecules not observed in the crystal structure of the reduced Rd from Pyrococcus furiosus. How general this correlation is among redox proteins must be determined in future work. The combination of our high-resolution crystal structures and molecular dynamics simulations provides a molecular picture of the structural rearrangement that occurs upon reduction in Cp rubredoxin.     

The role of backbone stability near Ala44 in the high reduction potential class of rubredoxins

Rubredoxins may be separated into high and low reduction potential classes, with reduction potentials differing by approximately 50 mV. Our previous work showed that a local shift in the polar backbone due to an A(44) versus V(44) side-chain size causes this reduction potential difference. However, this work also indicated that in the low potential Clostridium pasteurianum (Cp) rubredoxin, a V(44) --> A(44) mutation causes larger local backbone flexibility, because the V(44) side-chain present in the wild-type (wt) is no longer present to interlock with neighboring residues to stabilize the subsequent G(45). Since Pyrococcus furiosus (Pf) and other high potential rubredoxins generally have a P(45), it was presumed that a G(45) --> P(45) mutation might stabilize a V(44) --> A(44) mutation in Cp rubredoxin. Here crystal structure analysis, energy minimization, and molecular dynamics (MD) were performed for wt V(44)G(45), single mutant A(44)G(45) and double mutant A(44)P(45) Cp, and for wt A(44)P(45) Pf rubredoxins. The local structural, dynamical, and electrostatic properties of Cp gradually approach wt Pf in the order wt Cp to single to double mutant because of greater sequence similarity, as expected. The double mutant A(44)P(45) Cp exhibits increased backbone stability near residue 44 and thus enhances the probability that the backbone dipoles point toward the redox site, which favors an increase in the electrostatic contribution to the reduction potential. It appears that the electrostatic potential of residue 44 and the solvent accessibility to the redox are both determinants for the reduction potentials of homologous rubredoxins. Overall, these results indicate that an A(44) in a rubredoxin may require a P(45) for backbone stability whereas a V(44) can accommodate a G(45), since the valine side-chain can interlock with its neighbors.     

Prediction of reduction potential changes in rubredoxin: a molecular mechanics approach

Predicting the effects of mutation on the reduction potential of proteins is crucial in understanding how reduction potentials are modulated by the protein environment. Previously, we proposed that an alanine vs. a valine at residue 44 leads to a 50-mV difference in reduction potential found in homologous rubredoxins because of a shift in the polar backbone relative to the iron site due to the different side-chain sizes. Here, the aim is to determine the effects of mutations to glycine, isoleucine, and leucine at residue 44 on the structure and reduction potential of rubredoxin, and if the effects are proportional to side-chain size. Crystal structure analysis, molecular mechanics simulations, and experimental reduction potentials of wild-type and mutant Clostridium pasteurianum rubredoxin, along with sequence analysis of homologous rubredoxins, indicate that the backbone position relative to the redox site as well as solvent penetration near the redox site are both structural determinants of the reduction potential, although not proportionally to side-chain size. Thus, protein interactions are too complex to be predicted by simple relationships, indicating the utility of molecular mechanics methods in understanding them.     

Redox properties of mesophilic and hyperthermophilic rubredoxins as a function of pressure and temperature.

The formal equilibrium reduction potentials of recombinant electron transport protein, rubredoxin (MW = 7500 Da), from both the mesophilic Clostridium pasteurianum (Topt = 37 degrees C) and hyperthermophilic Pyrococcus furiosus (Topt = 95 degrees C) were recorded as a function of pressure and temperature. Measurements were made utilizing a specially designed stainless steel electrochemical cell that easily maintains pressures between 1 and 600 atm and a temperature-controlled cell that maintains temperatures between 4 and 100 degrees C. The reduction potential of P. furiosus rubredoxin was determined to be 31 mV at 25 degrees C and 1 atm, -93 mV at 95 degrees C and 1 atm, and 44 mV at 25 degrees C and 400 atm. Thus, the reduction potential of P. furiosus rubredoxin obtained under standard conditions is likely to be dramatically different from the reduction potential obtained under its normal operating conditions. Thermodynamic parameters associated with electron transfer were determined for both rubredoxins (for C. pasteurianum, DeltaV degrees = -27 mL/mol, DeltaS degrees = -36 cal K-1 mol-1, and DeltaH degrees = -10 kcal/mol, and for P. furiosus, DeltaV degrees = -31 mL/mol, DeltaS degrees = -41 cal K-1 mol-1, and DeltaH degrees = -13 kcal/mol) from its pressure- and temperature-reduction potential profiles. The thermodynamic parameters for electron transfer (DeltaV degrees, DeltaS degrees, and DeltaH degrees ) for both proteins were very similar, which is not surprising considering their structural similarities and sequence homology. Despite the fact that these two proteins exhibit dramatic differences in thermostability, it appears that structural changes that confer dramatic differences in thermostability do not significantly alter electron transfer reactivity. The experimental changes in reduction potential as a function of pressure and temperature were simulated using a continuum dielectric electrostatic model (DELPHI). A reasonable estimate of the protein dielectric constant (epsilonprotein) of 6 for both rubredoxins was determined from these simulations. A discussion is presented regarding the analysis of electrostatic interaction energies of biomolecules through pressure- and temperature-controlled electrochemical studies.     

Influence of protein flexibility on the redox potential of rubredoxin: Energy minimization studies

A theoretical investigation of the protein contribution to the redox potential of the iron–sulfur protein rubredoxin is presented. Structures of the oxidized and reduced forms of the protein were obtained by energy minimizing the oxidized crystal structure of Clostridium pasteurianum rubredoxin with appropriate charges and parameters. By including 102 crystal waters, structures close to the original crystal structure were obtained (rms difference of 1.16 Å), even with extensive minimization, thus allowing accurate calculations of comparative energies. Our calculations indicate an energy change of about –60 kcal/mol (2.58 eV) in the protein alone upon reduction. This energy change was due to both the change in charge of the redox site and the subsequent relaxation of the protein. An energy minimization procedure for the relaxation gives rms differences between the oxidized and reduced states of about 0.2 Å. The changes were small and occurred in both the backbone and sidechain mainly near the Fe–S center but contributed about – 16 kcal/mol (0.69 eV) to the total protein contribution. Although the neglect of certain effects such as electronic polarization may make the relaxation energies calculated an upper limit, the results indicate that protein relaxation contributes substantially to the redox potential. © 1993 Wiley-Liss, Inc.     

Structural origins of redox potentials in Fe-S proteins: electrostatic potentials of crystal structures.

Redox potentials often differ dramatically for homologous proteins that have identical redox centers. For two types of iron-sulfur proteins, the rubredoxins and the high-potential iron-sulfur proteins (HiPIPs), no structural explanations for these differences have been found. We calculated the classical electrostatic potential at the redox site using static crystal structures of four rubredoxins and four HiPIPs to identify important structural determinants of their redox potentials. The contributions from just the backbone and polar side chains are shown to explain major features of the experimental redox potentials. For instance, in the rubredoxins, the presence of Val 44 versus Ala 44 causes a backbone shift that explains a approximately 50 mV lower redox potential in one of the four rubredoxins. This result is consistent with experimental redox potentials of five additional rubredoxins with known sequence. Also, we attribute the unusually lower redox potentials of two of the HiPIPs studied to a less positive electrostatic environment around their redox sites. Finally, molecular dynamics simulations of solvent around static rubredoxin crystal structures indicate that water alone is a major factor in dampening the contribution of charged side chains, in accord with experiments showing that mutations of surface charges produce relatively little effect on redox potentials.     

A novel parameterization scheme for energy equations and its use to calculate the structure of protein molecules

A novel scheme for the parameterization of a type of “potential energy” function for protein molecules is introduced. The function is parameterized based on the known conformations of previously determined protein structures and their sequence similarity to a molecule whose conformation is to be calculated. Once parameterized, minima of the potential energy function can be located using a version of simulated annealing which has been previously shown to locate global and near-global minima with the given functional form. As a test problem, the potential was parameterized based on the known structures of the rubredoxins from Desulfovibrio vulgaris, Desulfovibrio desulfuricans, and Clostridium pasteurianum, which vary from 45 to 54 amino acids in length, and the sequence alignments of these molecules with the rubredoxin sequence from Desulfovibrio gigas. Since the Desulfovibrio gigas rubredeoxin conformation has also been determined, it is possible to check the accuracy of the results. Ten simulated-annealing runs from random starting conformations were performed. Seven of the 10 resultant conformations have an all-C_α rms deviation from the crystallographically determined conformation of less than 1.7 Å. For five of the structures, the rms deviation is less than 0.8 Å. Four of the structures have conformations which are virtually identical to each other except for the position of the carboxy-terminal residue. This is also the conformation which is achieved if the determined crystal structure is minimized with the same potential. The all-C_α rms difference between the crystal and minimized crystal structures is 0.6 Å. It is further observed that the “energies” of the structures according to the potential function exhibit a strong correlation with rms deviation from the native structure. The conformations of the individual model structures and the computational aspects of the modeling procedure are discussed. © 1993 Wiley-Liss, Inc.     

The Molecular Determinants of the Increased Reduction Potential of the Rubredoxin Domain of Rubrerythrin Relative to Rubredoxin

Based on the crystal structures, three possible sequence determinants have been suggested as the cause of a 285 mV increase in reduction potential of the rubredoxin domain of rubrerythrin over rubredoxin by modulating the polar environment around the redox site. Here, electrostatic calculations of crystal structures of rubredoxin and rubrerythrin and molecular dynamics simulations of rubredoxin wild-type and mutants are used to elucidate the contributions to the increased reduction potential. Asn¹⁶⁰ and His¹⁷⁹ in rubrerythrin versus valines in rubredoxins are predicted to be the major contributors, as the polar side chains contribute significantly to the electrostatic potential in the redox site region. The mutant simulations show both side chains rotating on a nanosecond timescale between two conformations with different electrostatic contributions. Reduction also causes a change in the reduction energy that is consistent with a linear response due to the interesting mechanism of shifting the relative populations of the two conformations. In addition to this, a simulation of a triple mutant indicates the side-chain rotations are approximately anticorrelated so whereas one is in the high potential conformation, the other is in the low potential conformation. However, Ala¹⁷⁶ in rubrerythrin versus a leucine in rubredoxin is not predicted to be a large contributor, because the solvent accessibility increases only slightly in mutant simulations and because it is buried in the interface of the rubrerythrin homodimer.     

Tetrathiolate coordination of germanium(IV) in a protein active site

The tetracysteine metal coordination site of the rubredoxins from Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) are shown to stably bind the inorganic Ge(IV) ion. This is the first characterized coordination complex of tetravalent germanium with a biological macromolecule. Zn(II), Ga(III) and Ge(IV) substitution yields differential NMR chemical shifts for the 1H and 15N amide resonances throughout much of the protein structure. The differential shifts for the six backbone amides that hydrogen bond to the metal-coordinated sulfurs indicate that the pseudo 2-fold symmetry of the active site is more closely maintained in the hyperthermophile Pf rubredoxin than in its mesophile Cp homolog. These three metal substitutions form an isoelectronic series of small diamagnetic proteins for which reference structures are known to 1A resolution. These series provide a promising system to analyze theoretical predictions of the effects of differential charge distribution on chemical shifts from both proximal and long range interactions.     

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).     

Isolation and characterization of rubrerythrin, a non-heme iron protein from Desulfovibrio vulgaris that contains rubredoxin centers and a hemerythrin-like binuclear iron cluster

A new non-heme iron protein from the periplasmic fraction of Desulfovibrio vulgaris (Hildenbourough NCIB 8303) has been purified to homogeneity, and its amino acid composition, molecular weight, redox potential, iron content, and optical, EPR, and M?ssbauer spectroscopic properties have been determined. This new protein is composed of two identical subunits with subunit molecular weight of 21,900 and contains four iron atoms per molecule. The as-purified oxidized protein exhibits an optical spectrum with absorption maxima at 492, 365, and 280 nm, and its EPR spectrum shows resonances at g = 4.3 and 9.4, characteristic of oxidized rubredoxin. The M?ssbauer data indicate the presence of approximately equal amounts of two types of iron; we named them the Rd-like and the Hr-like iron due to their similarity to the iron centers of rubredoxins (Rds) and hemerythrins (Hrs), respectively. For the Rd-like iron, the measured fine and hyperfine parameters (D = 1.5 cm-1, E/D = 0.26, delta EQ = -0.55 mm/s, delta = 0.27 mm/s, Axx/gn beta n = -16.5 T, Ayy/gn beta n = -15.6 T, and Azz/gn beta n = -17.0 T) are almost identical with those obtained for the rubredoxin from Clostridium pasteurianum. Redox-titration studies monitored by EPR, however, showed that these Rd-like centers have a midpoint redox potential of +230 +/- 10 mV, approximately 250 mV more positive than those reported for rubredoxins. Another unusual feature of this protein is the presence of the Hr-like iron atoms.(ABSTRACT TRUNCATED AT 250 WORDS)     

Phylogenetic studies of two rubredoxins from sulfate reducing bacteria

Summary The sequences of two rubredoxins isolated from the sulfate reducing bacteria:Desulfovibrio vulgaris andDesulfovibrio gigas have been elucidated. They have similar sequences but many more differences occur than would be expected from two bacteria of the same genus. Of the 52 sites, only 37 are occupied by identical residues. The primary structures are compared with those of the anaerobic bacteria rubredoxins ofClostridium pasteurianum, Micrococcus aerogenes, Pseudomonas oleovorans andPeptostreptococcus elsdenii: only 12 identities are found, mostly in the two clusters that contain two iron-bound cysteines each. A phylogenetic tree based on the primary structures is presented and possible relations with plant and bacterial ferredoxins are discussed. A secondary and tertiary structure, stereochemically compatible with the sequence data, is proposed.To whom reprint requests should be addressed     

Reduced temperature dependence of collective conformational opening in a hyperthermophile rubredoxin.

Spatially localized differences in the conformational dynamics of the rubredoxins from the hyperthermophile Pyrococcus furiosus (Pf) and the mesophile Clostridium pasteurianum (Cp) are monitored via amide exchange measurements. As shown previously for the hyperthermophile protein, nearly all backbone amides of the Cp rubredoxin exhibit EX(2) hydrogen exchange kinetics with conformational opening rates of >1 s(-)(1). Significantly slower amide exchange is observed for Pf rubredoxin in the region surrounding the metal site and the proximal end of the three-stranded beta-sheet, while for the rest of the structure, the exchange rates at 23 degrees C are similar for both proteins. For the multiple-turn region comprising residues 14-32 in both rubredoxins, the uniformity of both the exchange rate constants and the values of the activation energy at the slowly exchanging sites is consistent with a model of solvent exposure via a subglobal cooperative conformational opening. In contrast to the common expectation of increased rigidity in the hyperthermophile proteins, below room temperature Pf rubredoxin exhibits a larger apparent flexibility in this multiple-turn region. The smaller enthalpy for the conformational opening process of this region in Pf rubredoxin reflects the much weaker temperature dependence of the underlying conformational equilibrium in the hyperthermophile protein.     

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.     

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.     

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.     

Crystallographic investigation of the role of aspartate 95 in the modulation of the redox potentials of Desulfovibrio vulgaris flavodoxin

The side chain of aspartate 95 in flavodoxin from Desulfovibrio vulgaris provides the closest negative charge to N(1) of the bound FMN in the protein. Site-directed mutagenesis was used to substitute alanine, asparagine, or glutamate for this amino acid to assess the effect of this charge on the semiquinone/hydroquinone redox potential (E(1)) of the FMN cofactor. The D95A mutation shifts the E(1) redox potential positively by 16 mV, while a negative shift of 23 mV occurs in the oxidized/semiquinone midpoint redox potential (E(2)). The crystal structures of the oxidized and semiquinone forms of this mutant are similar to the corresponding states of the wild-type protein. In contrast to the wild-type protein, a further change in structure occurs in the D95A mutant in the hydroquinone form. The side chain of Y98 flips into an energetically more favorable edge-to-face interaction with the bound FMN. Analysis of the structural changes in the D95A mutant, taking into account electrostatic interactions at the FMN binding site, suggests that the pi-pi electrostatic repulsions have only a minor contribution to the very low E(1) redox potential of the FMN cofactor when bound to apoflavodoxin. Substitution of D95 with glutamate causes only a slight perturbation of the two one-electron redox potentials of the FMN cofactor. The structure of the D95E mutant reveals a large movement of the 60-loop (residues 60-64) away from the flavin in the oxidized structure. Reduction of this mutant to the hydroquinone causes the conformation of the 60-loop to revert back to that occurring in the structures of the wild-type protein. The crystal structures of the D95E mutant imply that electrostatic repulsion between a carboxylate on the side chain at position 95 and the phenol ring of Y98 prevents rotation of the Y98 side chain to a more energetically favorable conformation as occurs in the D95A mutant. Replacement of D95 with asparagine has no effect on E(2) but causes E(1) to change by 45 mV. The D95N mutant failed to crystallize. The K(d) values of the protein FMN complex in all three oxidation-reduction states differ from those of the wild-type complexes. Molecular modeling showed that the conformational energy of the protein changes with the redox state, in qualitative agreement with the observed changes in K(d), and allowed the electrostatic interactions between the FMN and the surrounding groups on the protein to be quantified.     

Thermostabilization of proteins by diglycerol phosphate, a new compatible solute from the hyperthermophile Archaeoglobus fulgidus

Diglycerol phosphate accumulates under salt stress in the archaeon Archaeoglobus fulgidus (L. O. Martins, R. Huber, H. Huber, K. O. Stetter, M. S. da Costa, and H. Santos, Appl. Environ. Microbiol. 63:896-902, 1997). This solute was purified after extraction from the cell biomass. In addition, the optically active and the optically inactive (racemic) forms of the compound were synthesized, and the ability of the solute to act as a protecting agent against heating was tested on several proteins derived from mesophilic or hyperthermophilic sources. Diglycerol phosphate exerted a considerable stabilizing effect against heat inactivation of rabbit muscle lactate dehydrogenase, baker's yeast alcohol dehydrogenase, and Thermococcus litoralis glutamate dehydrogenase. Highly homologous and structurally well-characterized rubredoxins from Desulfovibrio gigas, Desulfovibrio desulfuricans (ATCC 27774), and Clostridium pasteurianum were also examined for their thermal stabilities in the presence or absence of diglycerol phosphate, glycerol, and inorganic phosphate. These proteins showed different intrinsic thermostabilities, with half-lives in the range of 30 to 100 min. Diglycerol phosphate exerted a strong protecting effect, with approximately a fourfold increase in the half-lives for the loss of the visible spectra of D. gigas and C. pasteurianum rubredoxins. In contrast, the stability of D. desulfuricans rubredoxin was not affected. These different behaviors are discussed in the light of the known structural features of rubredoxins. The data show that diglycerol phosphate is a potentially useful protein stabilizer in biotechnological applications.     

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.     

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.