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

Characterization of ferredoxins involved in electron transfer pathways for nitrogen fixation implicates differences in electronic structure in tuning 2[4Fe4S] Fd activity

Ferredoxins (Fds) are small proteins which shuttle electrons to pathways like biological nitrogen fixation. Physical properties tune the reactivity of Fds with different pathways, but knowledge on how these properties can be manipulated to engineer new electron transfer pathways is lacking. Recently, we showed that an evolved strain of Rhodopseudomonas palustris uses a new electron transfer pathway for nitrogen fixation. This pathway involves a variant of the primary Fd of nitrogen fixation in R. palustris, Fer1, in which threonine at position 11 is substituted for isoleucine (Fer1T11I). To understand why this substitution in Fer1 enables more efficient electron transfer, we used in vivo and in vitro methods to characterize Fer1 and Fer1T11I. Electrochemical characterization revealed both Fer1 and Fer1T11I have similar redox transitions (-480 mV and - 550 mV), indicating the reduction potential was unaffected despite the proximity of T11 to an iron‑sulfur (FeS) cluster of Fer1. Additionally, disruption of hydrogen bonding around an FeS cluster in Fer1 by substituting threonine with alanine (T11A) or valine (T11V) did not increase nitrogenase activity, indicating that disruption of hydrogen bonding does not explain the difference in activity observed for Fer1T11I. Electron paramagnetic resonance spectroscopy studies revealed key differences in the electronic structure of Fer1 and Fer1T11I, which indicate changes to the high spin states and/or spin-spin coupling between the FeS clusters of Fer1. Our data implicates these electronic structure differences in facilitating electron flow and sets a foundation for further investigations to understand the connection between these properties and intermolecular electron transfer.

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.

Exploiting Molecular Symmetry to Quantitatively Map the Excited-State Landscape of Iron-Sulfur Clusters

Cuboidal [Fe4S4] clusters are ubiquitous cofactors in biological redox chemistry. In the [Fe4S4]1+ state, pairwise spin coupling gives rise to six arrangements of the Fe valences ("valence isomers") among the four Fe centers. Because of the magnetic complexity of these systems, it has been challenging to understand how a protein's active site dictates both the arrangement of the valences in the ground state as well as the population of excited-state valence isomers. Here, we show that the ground-state valence isomer landscape can be simplified from a six-level system in an asymmetric protein environment to a two-level system by studying the problem in synthetic [Fe4S4]1+ clusters with solution C3v symmetry. This simplification allows for the energy differences between valence isomers to be quantified (in some cases with a resolution of <0.1 kcal/mol) by simultaneously fitting the VT NMR and solution magnetic moment data. Using this fitting protocol, we map the excited-state landscape for a range of clusters of the form [(SIMes)3Fe4S4-X/L]n, (SIMes = 1,3-dimesityl-imidazol-4,5-dihydro-2-ylidene; n = 0 for anionic, X-type ligands and n = +1 for neutral, L-type ligands) and find that a single ligand substitution can alter the relative ground-state energies of valence isomers by at least 103 cm-1. On this basis, we suggest that one result of "non-canonical" amino acid ligation in Fe-S proteins is the redistribution of the valence electrons in the manifold of thermally populated excited states.

Characterization of the Geometrical and Electronic Structures of the Active Site and Its Effects on the Surrounding Environment in Reduced High-Potential Iron-Sulfur Proteins Investigated by the Density Functional Theory Approach

The high-potential iron-sulfur protein (HiPIP) is an electron-transporting protein that functions in the photosynthetic electron-transfer system and possesses a cubane-type [4Fe-4S] cluster in the active center. Characterization of the geometrical and electronic structures of the [4Fe-4S] cluster leads to an understanding of the functions in HiPIP, which are expected to be influenced by the environment surrounding the [4Fe-4S] cluster. This work characterized the geometrical and electronic structures of the [4Fe-4S] cluster in the reduced HiPIP and evaluated their effects on the protein environment using the density functional theory (DFT) approach. DFT calculations showed that the structural asymmetry and spin delocalization between iron atoms allowed for the acquisition of a unique stable geometrical and electronic structure in the open-shell singlet. In addition, the formation of an Fe-Fe bond accompanying the spin delocalization was found to depend on the interatomic distance. A comparison of the calculated stable structures with and without consideration of the amino acids around the [4Fe-4S] cluster demonstrated that the surrounding amino acids stabilized the unique geometrical and electronic structure of the [4Fe-4S] cluster in HiPIP.

Revisiting the concept of peptide bond planarity in an iron-sulfur protein by neutron structure analysis

The planarity of the peptide bond is important for the stability and structure formation of proteins. However, substantial distortion of peptide bonds has been reported in several high-resolution structures and computational analyses. To investigate the peptide bond planarity, including hydrogen atoms, we report a 1.2-Å resolution neutron structure of the oxidized form of high-potential iron-sulfur protein. This high-resolution neutron structure shows that the nucleus positions of the amide protons deviate from the peptide plane and shift toward the acceptors. The planarity of the H─N─C═O plane depends strongly on the pyramidalization of the nitrogen atom. Moreover, the orientation of the amide proton of Cys⁷⁵ is different in the reduced and oxidized states, possibly because of the electron storage capacity of the iron-sulfur cluster.

Recombination of 2Fe-2S Ferredoxins Reveals Differences in the Inheritance of Thermostability and Midpoint Potential

Recombination can be used in the laboratory to overcome component limitations in synthetic biology by creating enzymes that exhibit distinct activities and stabilities from native proteins. To investigate how recombination affects the properties of an oxidoreductase that transfers electrons in cells, we created ferredoxin (Fd) chimeras by recombining distantly related cyanobacterial and cyanomyophage Fds (53% identity) that present similar midpoint potentials but distinct thermostabilities. Fd chimeras having a wide range of amino acid substitutions retained the ability to coordinate an iron-sulfur cluster, although their thermostabilities varied with the fraction of residues inherited from each parent. The midpoint potentials of chimeric Fds also varied. However, all of the synthetic Fds exhibited midpoint potentials outside of the parental protein range. Each of the chimeric Fds could also support electron transfer between Fd-NADP reductase and sulfite reductase in Escherichia coli, although the chimeric Fds varied in the expression required for similar levels of cellular electron transfer. These results show how Fds can be diversified through recombination and reveal differences in the inheritance of thermostability and electrochemical properties. Furthermore, they illustrate how electron transfer efficiencies of chimeric Fds can be rapidly evaluated using a synthetic metabolic pathway.

Insight into the reaction mechanism of lipoyl synthase: a QM/MM study

Lipoyl synthase (LipA) catalyses the final step of the biosynthesis of the lipoyl cofactor by insertion of two sulfur atoms at the C6 and C8 atoms of the protein-bound octanoyl substrate. In this reaction, two [4Fe4S] clusters and two molecules of S-adenosyl-L-methionine are used. One of the two FeS clusters is responsible for the generation of a powerful oxidant, the 5'-deoxyadenosyl radical (5'-dA•). The other (the auxiliary cluster) is the source of both sulfur atoms that are inserted into the substrate. In this paper, the spin state of the FeS clusters and the reaction mechanism is investigated by the combined quantum mechanical and molecular mechanics approach. The calculations show that the ground state of the two FeS clusters, both in the [4Fe4S]2+ oxidation state, is a singlet state with antiferromagnetically coupled high-spin Fe ions and that there is quite a large variation of the energies of the various broken-symmetry states, up to 40 kJ/mol. For the two S-insertion reactions, the highest energy barrier is found for the hydrogen-atom abstraction from the octanoyl substrate by 5'-dA•. The formation of 5'-dA• is very facile for LipA, with an energy barrier of 6 kJ/mol for the first S-insertion reaction and without any barrier for the second S-insertion reaction. In addition, the first S ion attack on the C6 radical of octanoyl was found to take place directly by the transfer of the H6 from the substrate to 5'-dA•, whereas for the second S-insertion reaction, a C8 radical intermediate was formed with a rate-limiting barrier of 71 kJ/mol.

Importance of Hydrogen Bonding in Fine Tuning the [2Fe-2S] Cluster Redox Potential of HydC from Thermotoga maritima

Iron-sulfur clusters form one of the largest and most diverse classes of enzyme cofactors in nature. They may serve as structural factors, form electron transfer chains between active sites and external redox partners, or form components of enzyme active sites. Their specific role is a consequence of the cluster type and the surrounding protein environment. The relative effects of these factors are not completely understood, and it is not yet possible to predict the properties of iron-sulfur clusters based on amino acid sequences or rationally tune their properties to generate proteins with new desirable functions. Here, we generate mutations in a [2Fe-2S] cluster protein, the TmHydC subunit of the trimeric [FeFe]-hydrogenase from Thermotoga maritima, to study the factors that affect its redox potential. Saturation mutagenesis of Val131 was used to tune the redox potential over a 135 mV range and revealed that cluster redox potential and electronic properties correlate with amino acid hydrophobicity and the ability to form hydrogen bonds to the cluster. Proline scanning mutagenesis between pairs of ligating cysteines was used to remove backbone amide hydrogen bonds to the cluster and decrease the redox potential by up to 132 mV, without large structural changes in most cases. However, substitution of Gly83 with proline caused a change of HydC to a [4Fe-4S] cluster protein with a redox potential of -526 mV. Together, these results confirm the importance of hydrogen bonding in tuning cluster redox potentials and demonstrate the versatility of iron-sulfur cluster protein folds at binding different types of clusters.

Charge-density analysis of an iron-sulfur protein at an ultra-high resolution of 0.48 Å

The fine structures of proteins, such as the positions of hydrogen atoms, distributions of valence electrons and orientations of bound waters, are critical factors for determining the dynamic and chemical properties of proteins. Such information cannot be obtained by conventional protein X-ray analyses at 3.0-1.5 Å resolution, in which amino acids are fitted into atomically unresolved electron-density maps and refinement calculations are performed under strong restraints. Therefore, we usually supplement the information on hydrogen atoms and valence electrons in proteins with pre-existing common knowledge obtained by chemistry in small molecules. However, even now, computational calculation of such information with quantum chemistry also tends to be difficult, especially for polynuclear metalloproteins. Here we report a charge-density analysis of the high-potential iron-sulfur protein from the thermophilic purple bacterium Thermochromatium tepidum using X-ray data at an ultra-high resolution of 0.48 Å. Residual electron densities in the conventional refinement are assigned as valence electrons in the multipolar refinement. Iron 3d and sulfur 3p electron densities of the Fe4S4 cluster are visualized around the atoms. Such information provides the most detailed view of the valence electrons of the metal complex in the protein. The asymmetry of the iron-sulfur cluster and the protein environment suggests the structural basis of charge storing on electron transfer. Our charge-density analysis reveals many fine features around the metal complex for the first time, and will enable further theoretical and experimental studies of metalloproteins.

Mechanistic insights into the reductive dehydroxylation pathway for the biosynthesis of isoprenoids promoted by the IspH enzyme

Here, we report an integrated quantum mechanics/molecular mechanics (QM/MM) study of the bio-organometallic reaction pathway of the 2H⁺/2e^- reduction of (E)-4-hydroxy-3-methylbut-2-enyl pyrophosphate (HMBPP) into the so called universal terpenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), promoted by the IspH enzyme. Our results support the viability of the bio-organometallic pathway through rotation of the OH group of HMBPP away from the [Fe₄S₄] cluster at the core of the catalytic site, to become engaged in a H-bond with Glu126. This rotation is synchronous with π-coordination of the C2 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 C3 double bond of HMBPP to the apical Fe atom of the [Fe₄S₄] cluster. Dehydroxylation of HMBPP is triggered by a proton transfer from Glu126 to the OH group of HMBPP. The reaction pathway is completed by competitive proton transfer from the terminal phosphate group to the C2 or C4 atom of HMBPP.

Binding of promoter DNA to SoxR protein decreases the reduction potential of the [2Fe-2S] cluster

The [2Fe-2S] transcriptional factor SoxR, a member of the MerR family, functions as a sensor of oxidative stress in Escherichia coli. The transcriptional activity of SoxR is regulated by the reversible oxidation and reduction of [2Fe-2S] clusters. Electrochemistry measurements on DNA-modified electrodes have shown a dramatic shift in the reduction potential of SoxR from -290 to +200 mV with the promoter DNA-bound [ Gorodetsky , A. A. , Dietrich , L. E. P. , Lee , P. E. , Demple , B. , , Newman , D. K. , and Barton , J. K. ( 2008 ) DNA binding shifts the reduction potential of the transcription factor SoxR , Proc. Natl. Acad. Sci. U.S.A. 105 , 3684 - 3689 ]. To determine the change of the SoxR reduction potential using the new condition, the one-electron oxidation-reduction properties of [2Fe-2S] cluster in SoxR were investigated in the absence and presence of the DNA. The [2Fe-2S] cluster of SoxR was completely reduced by nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CRP) in the presence of a NADPH generating system (glucose 6-dehydrogenase and glucose-6 phosphate), indicating that CRP can serve as an NADPH-dependent electron carrier for SoxR. The reduction potential of SoxR was measured from equilibrium data coupled with NADPH and CRP in the presence of electron mediators. The reduction potentials of DNA-bound and DNA-free states of SoxR were -320 and -293 mV versus NHE (normal hydrogen electrode), respectively. These results indicate that DNA binding causes a moderate shift in the reduction potential of SoxR.

With or without light: comparing the reaction mechanism of dark-operative protochlorophyllide oxidoreductase with the energetic requirements of the light-dependent protochlorophyllide oxidoreductase

The addition of two electrons and two protons to the C₁₇=C₁₈ bond in protochlorophyllide is catalyzed by a light-dependent enzyme relying on NADPH as electron donor, and by a light-independent enzyme bearing a (Cys)₃Asp-ligated [4Fe–4S] cluster which is reduced by cytoplasmic electron donors in an ATP-dependent manner and then functions as electron donor to protochlorophyllide. The precise sequence of events occurring at the C₁₇=C₁₈ bond has not, however, been determined experimentally in the dark-operating enzyme. In this paper, we present the computational investigation of the reaction mechanism of this enzyme at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level of theory. The reaction mechanism begins with single-electron reduction of the substrate by the (Cys)₃Asp-ligated [4Fe–4S], yielding a negatively-charged intermediate. Depending on the rate of Fe–S cluster re-reduction, the reaction either proceeds through double protonation of the single-electron-reduced substrate, or by alternating proton/electron transfer. The computed reaction barriers suggest that Fe–S cluster re-reduction should be the rate-limiting stage of the process. Poisson–Boltzmann computations on the full enzyme–substrate complex, followed by Monte Carlo simulations of redox and protonation titrations revealed a hitherto unsuspected pH-dependence of the reaction potential of the Fe–S cluster. Furthermore, the computed distributions of protonation states of the His, Asp and Glu residues were used in conjuntion with single-point ONIOM computations to obtain, for the first time, the influence of all protonation states of an enzyme on the reaction it catalyzes. Despite exaggerating the ease of reduction of the substrate, these computations confirmed the broad features of the reaction mechanism obtained with the medium-sized models, and afforded valuable insights on the influence of the titratable amino acids on each reaction step. Additional comparisons of the energetic features of the reaction intermediates with those of common biochemical redox intermediates suggest a surprisingly simple explanation for the mechanistic differences between the dark-catalyzed and light-dependent enzyme reaction mechanisms.

Assessment of Quantum Mechanical Methods for Copper and Iron Complexes by Photoelectron Spectroscopy

Broken-symmetry density functional theory (BS-DFT) calculations are assessed for redox energetics [Cu(SCH₃)₂]^1–/0, [Cu(NCS)₂]^1–/0, [FeCl₄]^1–/0, and [Fe(SCH₃)₄]^1–/0 against vertical detachment energies (VDE) from valence photoelectron spectroscopy (PES), as a prelude to studies of metalloprotein analogs. The M06 and B3LYP hybrid functionals give VDE that agree with the PES VDE for the Fe complexes, but both underestimate it by ∼400 meV for the Cu complexes; other hybrid functionals give VDEs that are an increasing function of the amount of Hartree–Fock (HF) exchange and so cannot show good agreement for both Cu and Fe complexes. Range-separated (RS) functionals appear to give a better distribution of HF exchange since the negative HOMO energy is approximately equal to the VDEs but also give VDEs dependent on the amount of HF exchange, sometimes leading to ground states with incorrect electron configurations; the LRC-ωPBEh functional reduced to 10% HF exchange at short-range give somewhat better values for both, although still ∼150 meV too low for the Cu complexes and ∼50 meV too high for the Fe complexes. Overall, the results indicate that while HF exchange compensates for self-interaction error in DFT calculations of both Cu and Fe complexes, too much may lead to more sensitivity to nondynamical correlation in the spin-polarized Fe complexes.

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.

Heterometallic [AgFe(3)S (4)] ferredoxin variants: synthesis, characterization, and the first crystal structure of an engineered heterometallic iron-sulfur protein

Heterometallic [AgFe(3)S(4)] iron-sulfur clusters assembled in wild-type Pyrococcus furiosus ferredoxin and two variants, D14C and D14H, are characterized. The crystal structure of the [AgFe(3)S(4)] D14C variant shows that the silver(I) ion is indeed part of the cluster and is coordinated to the thiolate group of residue 14. Cyclic voltammetry shows one redox pair with a reduction potential of +220 mV versus the standard hydrogen electrode which is assigned to the [AgFe(3)S(4)](2+/+) couple. The oxidized form of the [AgFe(3)S(4)] D14C variant is stable in the presence of dioxygen, whereas the oxidized forms of the [AgFe(3)S(4)] wild type and D14H variants convert to the [Fe(3)S(4)] ferredoxin form. The monovalent d (10) silver(I) ion stabilizes the [Fe(3)S(4)](+/0) cluster fragment, as opposed to divalent d (10) metal ions, resulting in more than 0.4 V difference in reduction potentials between the silver(I) and, e.g., zinc(II) heterometallic [MFe(3)S(4)] ferredoxins. The trend in reduction potentials for the variants containing the [AgFe(3)S(4)] cluster is wild type ≤ D14C < D14H and shows the same trend as reported for the variants containing the [Fe(3)S(4)] cluster, but is different from the D14C < D14H < wild type trend reported for the [Fe(4)S(4)] ferredoxin. The similarity in the reduction potential trend for the variants containing the heterometallic [AgFe(3)S(4)] cluster and the [Fe(3)S(4)] cluster can be rationalized in terms of the electrostatic influence of the residue 14 side chains, rather than the dissociation constant of this residue, as is the case for [Fe(4)S(4)] ferredoxins. The trends in reduction potentials are in line with there being no electronic coupling between the silver(I) ion and the Fe(3)S(4) fragment.

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.

Understanding rubredoxin redox sites by density functional theory studies of analogues

Determining the redox energetics of redox site analogues of metalloproteins is essential in unraveling the various contributions to electron transfer properties of these proteins. Since studies of the [4Fe-4S] analogues show that the energies are dependent on the ligand dihedral angles, broken symmetry density functional theory (BS-DFT) with the B3LYP functional and double-ζ basis sets calculations of optimized geometries and electron detachment energies of [1Fe] rubredoxin analogues are compared to crystal structures and gas-phase photoelectron spectroscopy data, respectively, for [Fe(SCH(3))(4)](0/1-/2-), [Fe(S(2)-o-xyl)(2)](0/1-/2-), and Na(+)[Fe(S(2)-o-xyl)(2)](1-/2-) in different conformations. In particular, the study of Na(+)[Fe(S(2)-o-xyl)(2)](1-/2-) is the only direct comparison of calculated and experimental gas phase detachment energies for the 1-/2- couple found in the rubredoxins. These results show that variations in the inner sphere energetics by up to ∼0.4 eV can be caused by differences in the ligand dihedral angles in either or both redox states. Moreover, these results indicate that the protein stabilizes the conformation that favors reduction. In addition, the free energies and reorganization energies of oxidation and reduction as well as electrostatic potential charges are calculated, which can be used as estimates in continuum electrostatic calculations of electron transfer properties of [1Fe] proteins.

Geometric and electrostatic study of the [4Fe-4S] cluster of adenosine-5'-phosphosulfate reductase from broken symmetry density functional calculations and extended X-ray absorption fine structure spectroscopy

Adenosine-5'-phosphosulfate reductase (APSR) is an iron-sulfur protein that catalyzes the reduction of adenosine-5'-phosphosulfate (APS) to sulfite. APSR coordinates to a [4Fe-4S] cluster via a conserved CC-X(~80)-CXXC motif, and the cluster is essential for catalysis. Despite extensive functional, structural, and spectroscopic studies, the exact role of the iron-sulfur cluster in APS reduction remains unknown. To gain an understanding into the role of the cluster, density functional theory (DFT) analysis and extended X-ray fine structure spectroscopy (EXAFS) have been performed to reveal insights into the coordination, geometry, and electrostatics of the [4Fe-4S] cluster. X-ray absorption near-edge structure (XANES) data confirms that the cluster is in the [4Fe-4S](2+) state in both native and substrate-bound APSR while EXAFS data recorded at ~0.1 Å resolution indicates that there is no significant change in the structure of the [4Fe-4S] cluster between the native and substrate-bound forms of the protein. On the other hand, DFT calculations provide an insight into the subtle differences between the geometry of the cluster in the native and APS-bound forms of APSR. A comparison between models with and without the tandem cysteine pair coordination of the cluster suggests a role for the unique coordination in facilitating a compact geometric structure and "fine-tuning" the electronic structure to prevent reduction of the cluster. Further, calculations using models in which residue Lys144 is mutated to Ala confirm the finding that Lys144 serves as a crucial link in the interactions involving the [4Fe-4S] cluster and APS.

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.

Solvation effects on S K-edge XAS spectra of Fe-S proteins: normal and inverse effects on WT and mutant rubredoxin

S K-edge X-ray absorption spectroscopy (XAS) was performed on wild type Cp rubredoxin and its Cys --> Ser mutants in both solution and lyophilized forms. For wild type rubredoxin and for the mutants where an interior cysteine residue (C6 or C39) is substituted by serine, a normal solvent effect is observed, that is, the S covalency increases upon lyophilization. For the mutants where a solvent accessible surface cysteine residue is substituted by serine, the S covalency decreases upon lyophilization which is an inverse solvent effect. Density functional theory (DFT) calculations reproduce these experimental results and show that the normal solvent effect reflects the covalency decrease due to solvent H-bonding to the surface thiolates and that the inverse solvent effect results from the covalency compensation from the interior thiolates. With respect to the Cys --> Ser substitution, the S covalency decreases. Calculations indicate that the stronger bonding interaction of the alkoxide with the Fe relative to that of thiolate increases the energy of the Fe d orbitals and reduces their bonding interaction with the remaining cysteines. The solvent effects support a surface solvent tuning contribution to electron transfer, and the Cys --> Ser result provides an explanation for the change in properties of related iron-sulfur sites with this mutation.