Sequence determination of reduction potentials by cysteinyl hydrogen bonds and peptide pipoles in [4Fe-4S] ferredoxins

Redox Potentials of Iron-Sulfur Clusters in Type I Photosynthetic Reaction Centers

The electron transfer pathways in type I photosynthetic reaction centers, such as photosystem I (PSI) and reaction centers from green sulfur bacteria (GsbRC), are terminated by two Fe4S4 clusters, FA and FB. The protein structures are the basis of understanding how the protein electrostatic environment interacts with the Fe4S4 clusters and facilitates electron transfer. Using the protein structures, we calculated the redox potential (Em) values for FA and FB in PSI and GsbRC, solving the linear Poisson-Boltzmann equation. The FA-to-FB electron transfer is energetically downhill in the cyanobacterial PSI structure, while it is isoenergetic in the plant PSI structure. The discrepancy arises from differences in the electrostatic influences of conserved residues, including PsaC-Lys51 and PsaC-Arg52, located near FA. The FA-to-FB electron transfer is slightly downhill in the GsbRC structure. Em(FA) and Em(FB) exhibit similar levels upon isolation of the membrane-extrinsic PsaC and PscB subunits from the PSI and GsbRC reaction centers, respectively. The binding of the membrane-extrinsic subunit at the heterodimeric/homodimeric reaction center plays a key role in tuning Em(FA) and Em(FB).

Computational engineering of previously crystallized pyruvate formate-lyase activating enzyme reveals insights into SAM binding and reductive cleavage

Radical S-adenosyl-l-methionine (SAM) enzymes are ubiquitous in nature and carry out a broad variety of difficult chemical transformations initiated by hydrogen atom abstraction. Although numerous radical SAM (RS) enzymes have been structurally characterized, many prove recalcitrant to crystallization needed for atomic-level structure determination using X-ray crystallography, and even those that have been crystallized for an initial study can be difficult to recrystallize for further structural work. We present here a method for computationally engineering previously observed crystallographic contacts and employ it to obtain more reproducible crystallization of the RS enzyme pyruvate formate-lyase activating enzyme (PFL-AE). We show that the computationally engineered variant binds a typical RS [4Fe-4S]2+/+ cluster that binds SAM, with electron paramagnetic resonance properties indistinguishable from the native PFL-AE. The variant also retains the typical PFL-AE catalytic activity, as evidenced by the characteristic glycyl radical electron paramagnetic resonance signal observed upon incubation of the PFL-AE variant with reducing agent, SAM, and PFL. The PFL-AE variant was also crystallized in the [4Fe-4S]2+ state with SAM bound, providing a new high-resolution structure of the SAM complex in the absence of substrate. Finally, by incubating such a crystal in a solution of sodium dithionite, the reductive cleavage of SAM is triggered, providing us with a structure in which the SAM cleavage products 5'-deoxyadenosine and methionine are bound in the active site. We propose that the methods described herein may be useful in the structural characterization of other difficult-to-resolve proteins.

The Di-Iron Protein YtfE Is a Nitric Oxide-Generating Nitrite Reductase Involved in the Management of Nitrosative Stress

Previously characterized nitrite reductases fall into three classes: siroheme-containing enzymes (NirBD), cytochrome c hemoproteins (NrfA and NirS), and copper-containing enzymes (NirK). We show here that the di-iron protein YtfE represents a physiologically relevant new class of nitrite reductases. Several functions have been previously proposed for YtfE, including donating iron for the repair of iron-sulfur clusters that have been damaged by nitrosative stress, releasing nitric oxide (NO) from nitrosylated iron, and reducing NO to nitrous oxide (N2O). Here, in vivo reporter assays confirmed that Escherichia coli YtfE increased cytoplasmic NO production from nitrite. Spectroscopic and mass spectrometric investigations revealed that the di-iron site of YtfE exists in a mixture of forms, including nitrosylated and nitrite-bound, when isolated from nitrite-supplemented, but not nitrate-supplemented, cultures. Addition of nitrite to di-ferrous YtfE resulted in nitrosylated YtfE and the release of NO. Kinetics of nitrite reduction were dependent on the nature of the reductant; the lowest Km, measured for the di-ferrous form, was ∼90 μM, well within the intracellular nitrite concentration range. The vicinal di-cysteine motif, located in the N-terminal domain of YtfE, was shown to function in the delivery of electrons to the di-iron center. Notably, YtfE exhibited very low NO reductase activity and was only able to act as an iron donor for reconstitution of apo-ferredoxin under conditions that damaged its di-iron center. Thus, YtfE is a high-affinity, low-capacity nitrite reductase that we propose functions to relieve nitrosative stress by acting in combination with the co-regulated NO-consuming enzymes Hmp and Hcp.

A unified cost-effective method for the construction of reliable potential energy surfaces for H2S and H2O clusters

A DFT-based methodology has been used to construct the potential energy surface of H2S clusters up to pentamers. Geometrical parameters and energetics show very good agreement with the existing experimental and high-level theoretical results. Distinct stable conformers of three dimers, six trimers, eleven tetramers and twenty-three pentamers have been identified. Both S-HS H-bond and SS interactions are identified in dimers, trimers and pentamers, while no SS interactions could be found in any of the 11 tetramer conformers. The binding energies of the most stable dimer, trimer, tetramer and pentamer are -1.66, -5.21, -8.57 and -12.54 kcal mol-1, respectively. The PES has been found to be exceedingly flat and the energy gap between the most and the least stable conformers was found to be only 0.09, 2.13, 1.65 and 1.13 kcal mol-1, from the dimer to the pentamer, respectively. The proposed method has also been used for water clusters up to the pentamer. The results obtained were found to agree closely with the existing results. Only one conformer was found for the water dimer, whereas four, five and fifteen conformers were obtained for the trimer, tetramer and pentamer, respectively. Atoms in molecular calculations were found to corroborate with the geometric and energetic results for both clusters.

Observation of an Unusually Large IR Red-Shift in an Unconventional S-H···S Hydrogen-Bond

The S-H···S non-covalent interaction is generally known as an extremely unconventional weak hydrogen-bond in the literature. The present gas-phase spectroscopic investigation shows that the S-H···S hydrogen-bond can be as strong as any conventional hydrogen-bond in terms of the IR red-shift in the stretching frequency of the hydrogen-bond donor group. Herein, the strength of the S-H···S hydrogen-bond has been determined by measuring the red-shift (∼150 cm^-1) of the S-H stretching frequency in a model complex of 2-chlorothiophenol and dimethyl sulfide using isolated gas-phase IR spectroscopy coupled with quantum chemistry calculations. The observation of an unusually large IR red-shift in the S-H···S hydrogen-bond is explained in terms of the presence of a significant amount of charge-transfer interactions in addition to the usual electrostatic interactions. The existence of ∼750 S-H···S interactions between the cysteine and methionine residues in 642 protein structures determined from an extensive Protein Data Bank analysis also indicates that this interaction is important for the structures of proteins.

Cooperative nature of the sulfur centered hydrogen bond: investigation of (H2S)n (n = 2-4) clusters using an affordable yet accurate level of theory

Existing studies have shown that appreciably high level quantum chemical calculations are required to reproduce experimental energetic and geometric features of a H₂S dimer. This condition severely restricts any practical possibility of obtaining reliable results for higher order H₂S clusters. We have shown here that the binding energies calculated at the CCSD(T)/CBS level with counterpoise corrected geometries calculated at the MP2/aug-cc-pV(Q+d)Z level of theory excellently match with the experimental results for the H₂S dimer. Subsequently, the above mentioned levels of theory were used for trimers and tetramers. (H₂S)_n (n = 2-4) clusters were found to show cooperative strengthening of S-HS hydrogen bonds, which is clearly evident from the evolution of binding energies and hydrogen bond lengths, with increasing cluster size. Localized molecular orbital energy decomposition analyses have been carried out to understand how the contributions of various energy components modulate with the size of the clusters and what are their relative contributions towards the overall stabilization of the clusters. Natural bond orbital and atoms in molecules analyses were also carried out in order to look into the evolution of the electronic charge transfer and electron density topology with cluster size.

A designed heme-[4Fe-4S] metalloenzyme catalyzes sulfite reduction like the native enzyme

Multielectron redox reactions often require multicofactor metalloenzymes to facilitate coupled electron and proton movement, but it is challenging to design artificial enzymes to catalyze these important reactions, owing to their structural and functional complexity. We report a designed heteronuclear heme-[4Fe-4S] cofactor in cytochrome c peroxidase as a structural and functional model of the enzyme sulfite reductase. The initial model exhibits spectroscopic and ligand-binding properties of the native enzyme, and sulfite reduction activity was improved-through rational tuning of the secondary sphere interactions around the [4Fe-4S] and the substrate-binding sites-to be close to that of the native enzyme. By offering insight into the requirements for a demanding six-electron, seven-proton reaction that has so far eluded synthetic catalysts, this study provides strategies for designing highly functional multicofactor artificial enzymes.

Web-based computational chemistry education with CHARMMing III: Reduction potentials of electron transfer proteins

A module for fast determination of reduction potentials, E°, of redox-active proteins has been implemented in the CHARMM INterface and Graphics (CHARMMing) web portal (www.charmming.org). The free energy of reduction, which is proportional to E°, is composed of an intrinsic contribution due to the redox site and an environmental contribution due to the protein and solvent. Here, the intrinsic contribution is selected from a library of pre-calculated density functional theory values for each type of redox site and redox couple, while the environmental contribution is calculated from a crystal structure of the protein using Poisson-Boltzmann continuum electrostatics. An accompanying lesson demonstrates a calculation of E°. In this lesson, an ionizable residue in a [4Fe-4S]-protein that causes a pH-dependent E° is identified, and the E° of a mutant that would test the identification is predicted. This demonstration is valuable to both computational chemistry students and researchers interested in predicting sequence determinants of E° for mutagenesis.

Identifying sequence determinants of reduction potentials of metalloproteins

The reduction potential of an electron transfer protein is one of its most important functional characteristics. Although the type of redox site and the protein fold are the major determinants of the reduction potential of a redox-active protein, its amino acid sequence may tune the reduction potential as well. Thus, homologous proteins can often be divided into different classes, with each class characterized by a biological function and a reduction potential. Site-specific mutagenesis of the sequence determinants of the differences in the reduction potential between classes should change the reduction potential of a protein in one class to that of the other class. Here, a procedure is presented that combines energetic and bioinformatic analysis of homologous proteins to identify sequence determinants that are also good candidates for site-specific mutations, using the [4Fe-4S] ferredoxins and the [4Fe-4S] high-potential iron-sulfur proteins as examples. This procedure is designed to guide site-specific mutations or more computationally expensive studies, such as molecular dynamics simulations. To make the procedure more accessible to the general scientific community, it is being implemented into CHARMMing, a Web-based portal, with a library of density functional theory results for the redox site that are used in the setting up of Poisson-Boltzmann continuum electrostatics calculations for the protein energetics.

Calculating standard reduction potentials of [4Fe-4S] proteins

The oxidation-reduction potentials of electron transfer proteins determine the driving forces for their electron transfer reactions. Although the type of redox site determines the intrinsic energy required to add or remove an electron, the electrostatic interaction energy between the redox site and its surrounding environment can greatly shift the redox potentials. Here, a method for calculating the reduction potential versus the standard hydrogen electrode, E°, of a metalloprotein using a combination of density functional theory and continuum electrostatics is presented. This work focuses on the methodology for the continuum electrostatics calculations, including various factors that may affect the accuracy. The calculations are demonstrated using crystal structures of six homologous HiPIPs, which give E° that are in excellent agreement with experimental results.

Characterizing the effects of the protein environment on the reduction potentials of metalloproteins

The reduction potentials of electron transfer proteins are critically determined by the degree of burial of the redox site within the protein and the degree of permanent polarization of the polypeptide around the redox site. Although continuum electrostatics calculations of protein structures can predict the net effect of these factors, quantifying each individual contribution is a difficult task. Here, the burial of the redox site is characterized by a dielectric radius R(p) (a Born-type radius for the protein), the polarization of the polypeptide is characterized by an electret potential ϕ(p) (the average electrostatic potential at the metal atoms), and an electret-dielectric spheres (EDS) model of the entire protein is then defined in terms of R(p) and ϕ(p). The EDS model shows that for a protein with a redox site of charge Q, the dielectric response free energy is a function of Q(2), while the electret energy is a function of Q. In addition, R(p) and ϕ(p) are shown to be characteristics of the fold of a protein and are predictive of the most likely redox couple for redox sites that undergo different redox couples.

Fold versus sequence effects on the driving force for protein-mediated electron transfer

Electron transport chains composed of electron transfer reactions mainly between proteins provide fast efficient flow of energy in a variety of metabolic pathways. Reduction potentials are essential characteristics of the proteins because they determine the driving forces for the electron transfers. As both polar and charged groups from the backbone and side chains define the electrostatic environment, both the fold and the sequence will contribute. However, although the role of a specific sequence may be determined by experimental mutagenesis studies of reduction potentials, understanding the role of the fold by experiment is much more difficult. Here, continuum electrostatics and density functional theory calculations are used to analyze reduction potentials in [4Fe-4S] proteins. A key feature is that multiple homologous proteins in three different folds are compared: six high potential iron-sulfur proteins, four bacterial ferredoxins, and four nitrogenase iron proteins. Calculated absolute reduction potentials are shown to be in quantitative agreement with electrochemical reduction potentials. Calculations further demonstrate that the contribution of the backbone is larger than that of the side chains and is consistent for homologous proteins but differs for nonhomologous proteins, indicating that the fold is the major protein factor determining the reduction potential, whereas the specific amino acid sequence tunes the reduction potential for a given fold. Moreover, the fold contribution is determined mainly by the proximity of the redox site to the protein surface and the orientation of the dipoles of backbone near the redox site.

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(160) and His(179) 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(176) 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.

Redox potential of the Rieske iron-sulfur protein quantum-chemical and electrostatic study

Quantum-chemical study of structures, energies, and effective partial charge distribution for several models of the Rieske protein redox center is performed in terms of the B3LYP density functional method in combination with the broken symmetry approach using three different atomic basis sets. The structure of the redox complex optimized in vacuum differs markedly from that inside the protein. This means that the protein matrix imposes some stress on the active site resulting in distortion of its structure. The redox potentials calculated for the real active site structure are in a substantially better agreement with the experiment than those calculated for the idealized structure. This shows an important role of the active site distortion in tuning its redox potential. The reference absolute electrode potential of the standard hydrogen electrode is used that accounts for the correction caused by the water surface potential. Electrostatic calculations are performed in the framework of the polarizable solute model. Two dielectric permittivities of the protein are employed: the optical permittivity for calculation of the intraprotein electric field, and the static permittivity for calculation of the dielectric response energy. Only this approach results in a reasonable agreement of the calculated and experimental redox potentials.

Cleavage of [4Fe-4S]-type clusters: breaking the symmetry

The cleavage of [4Fe-4S]-type clusters is thought to be important in proteins such as Fe-S scaffold proteins and nitrogenase. However, most [4Fe-4S](2+) clusters in proteins have two antiferromagnetically coupled high-spin layers in which a minority spin is delocalized in each layer, thus forming a symmetric Fe(2.5+)-Fe(2.5+) pair, and how cleavage occurs between the irons is puzzling because of the shared electron. Previously, we proposed a novel mechanism for the fission of a [4Fe-4S] core into two [2Fe-2S] cores in which the minority spin localizes on one iron, thus breaking the symmetry and creating a transition state with two Fe(3+)-Fe(2+) pairs. Cleavage first through the weak Fe(2+)-S bonds lowers the activation energy. Here, we propose a test of this mechanism: break the symmetry of the cluster by changing the ligands to promote spin localization, which should enhance reactivity. The cleavage reactions for the homoligand [Fe(4)S(4)L(4)](2-) (L = SCH(3), Cl, H) and heteroligand [Fe(4)S(4)(SCH(3))(2)L(2)](2-) (L = Cl, H) clusters in the gas phase were examined via broken-symmetry density functional theory calculations. In the heteroligand clusters, the minority spin localized on the iron coordinated by the weaker electron-donor ligand, and the reaction energy and activation barrier of the cleavage were lowered, which is in accord with our proposed mechanism and consistent with photoelectron spectroscopy and collision-induced dissociation experiments. These studies suggest that proteins requiring facile fission of their [4Fe-4S] cluster in their biological function might have spin-localized [4Fe-4S] clusters.

Geometric characteristics of hydrogen bonds involving sulfur atoms in proteins

Sulfur atoms have been known to participate in hydrogen bonds (H-bonds) and these sulfur-containing H-bonds (SCHBs) are suggested to play important roles in certain biological processes. This study aims to comprehensively characterize all the SCHBs in 500 high-resolution protein structures (< or =1.8 A). We categorized SCHBs into six types according to donor/acceptor behaviors and used explicit hydrogen approach to distinguish SCHBs from those of nonhydrogen bonding interactions. It is revealed that sulfur atom is a very poor H-bond acceptor, but a moderately good H-bond donor. In alpha-helix, considerable SCHBs were found between the sulphydryl group of cysteine residue i and the carbonyl oxygen of residue i-4, and these SCHBs exert effects in stabilizing helices. Although for other SCHBs, they possess no specific secondary structural preference, their geometric characteristics in proteins and in free small compounds are significantly distinct, indicating the protein SCHBs are geometrically distorted. Interestingly, sulfur atom in the disulfide bond tends to form bifurcated H-bond whereas in cysteine-cysteine pairs prefer to form dual H-bond. These special H-bonds remarkably boost the interaction between H-bond donor and acceptor. By oxidation/reduction manner, the mutual transformation between the dual H-bonds and disulfide bonds for cysteine-cysteine pairs can accurately adjust the structural stability and biological function of proteins in different environments. Furthermore, few loose H-bonds were observed to form between the sulphydryl groups and aromatic rings, and in these cases the donor H is almost over against the rim rather than the center of the aromatic ring.

Limited Effectiveness of Antioxidants in the Protection of Yeast Defective in Antioxidant Proteins

Efficacy of several antioxidants in the protection of the yeast Saccharomyces cerevisiae mutants deficient in CuZnSOD and deficient in glutaredoxin 5 to growth restriction induced by oxidants was studied. Ascorbate and glutathione protected the Deltasod1 and Deltagrx5 mutants against the effects of t-butyl hydroperoxide and cumene hydroperoxide, Deltasod1 mutants against oxytetracycline and Deltagrx5 mutants against menadione and 2,2'-azobis-(2-amidinopropane). However, Tempol, Trolox and melatonin were much less effective, showing prooxidative effects and, at high concentrations, hampering the growth of the mutants in the absence of exogenous oxidants. These results point to a complication of cellular effects of antioxidants by their prooxidative effects and to the usefulness of cellular tests to evaluate the biological effectiveness of antioxidants.

Insight into environmental effects on bonding and redox properties of [4Fe-4S] clusters in proteins

The large differences in redox potentials between the HiPIPs and ferredoxins are generally attributed to hydrogen bonds and electrostatic effects from the protein and solvent. Recent ligand K-edge X-ray absorption studies by Solomon and co-workers show that the Fe-S covalencies of [4Fe-4S] clusters in the two proteins differ considerably apparently because of hydrogen bonds from water, indicating electronic effects may be important. However, combined density function theory (DFT) and photoelectron spectroscopy studies by our group and Wang and co-workers indicate that hydrogen bonds tune the potential of [4Fe-4S] clusters by mainly electrostatics. The DFT studies here rationalize both results, namely that the observed change in the Fe-S covalency is due to differences in ligand conformation between the two proteins rather than hydrogen bonds. Moreover, the ligand conformation affects the calculated potentials by approximately 100 mV and, thus, is a heretofore unconsidered means of tuning the potential.

Semi-continuum electrostatic calculations of redox potentials in photosystem I

The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill.

Redox-dependent Structural Reorganization in Putidaredoxin, a Vertebrate-type [2Fe-2S] Ferredoxin from Pseudomonas putida

Putidaredoxin (Pdx), a vertebrate-type [2Fe-2S] ferredoxin from Pseudomonas putida, transfers electrons from NADH-putidaredoxin reductase to cytochrome P450cam. Pdx exhibits redox-dependent binding affinities for P450cam and is thought to play an effector role in the monooxygenase reaction catalyzed by this hemoprotein. To understand how the reduced form of Pdx is stabilized and how reduction of the [2Fe-2S] cluster affects molecular properties of the iron-sulfur protein, crystal structures of reduced C73S and C73S/C85S Pdx were solved to 1.45 angstroms and 1.84 angstroms resolution, respectively, and compared to the corresponding 2.0 angstroms and 2.03 angstroms X-ray models of the oxidized mutants. To prevent photoreduction, the latter models were determined using in-house radiation source and the X-ray dose received by Pdx crystals was significantly decreased. Structural analysis showed that in reduced Pdx the Cys45-Ala46 peptide bond flip initiates readjustment of hydrogen bonding interactions between the [2Fe-2S] cluster, the Sgamma atoms of the cysteinyl ligands, and the backbone amide nitrogen atoms that results in tightening of the Cys39-Cys48 metal cluster binding loop around the prosthetic group and shifting of the metal center toward the Cys45-Thr47 peptide. From the metal center binding loop, the redox changes are transmitted to the linked Ile32-Asp38 peptide triggering structural rearrangement between the Tyr33-Asp34, Ser7-Asp9 and Pro102-Asp103 fragments of Pdx. The newly established hydrogen bonding interactions between Ser7, Asp9, Tyr33, Asp34, and Pro102, in turn, not only stabilize the tightened conformation of the [2Fe-2S] cluster binding loop but also assist in formation of a specific structural patch on the surface of Pdx that can be recognized by P450cam. This redox-linked change in surface properties is likely to be responsible for different binding affinity of oxidized and reduced Pdx to the hemoprotein.

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.