Atomic resolution structures of oxidized [4Fe-4S] ferredoxin from Bacillus thermoproteolyticus in two crystal forms: systematic distortion of [4Fe-4S] cluster in the protein

Predicting Iron-Sulfur Cluster Redox Potentials: A Simple Model Derived from Protein Structures

Iron-sulfur (Fe-S) clusters are critical cofactors in metalloproteins, essential for cellular processes such as energy production, DNA repair, enzymatic catalysis, and metabolic regulation. While Fe-S cluster functions are intimately linked to their redox properties, their precise roles in many proteins remain unclear. In this study, we present a regression model based on experimental redox potential (E m ) data, utilizing only two features: the Fe-S cluster's total charge and the Fe atoms' average valence. This model achieves a high correlation with experimental data (R 2 = 0.82) and an average prediction error of 0.12 V. Applying this model across the Protein Data Bank, we predict E m values for all cataloged Fe-S clusters, uncovering redox potential trends across diverse cluster classes. The computed redox potentials showed strong agreement with experimental values, achieving an overall accuracy of 88%. This streamlined, computationally accessible approach enhances the annotation and mechanistic understanding of Fe-S proteins, offering new insights into the redox variability of electron transport proteins. Our model holds promise for advancing studies of metalloprotein function and facilitating the design of bioinspired redox systems.

Parameterization of a Fluctuating Charge Model for Complexes Containing 3d Transition Metals

Metalloproteins widely exist in biology, playing pivotal roles in diverse life processes. Meanwhile, molecular dynamics (MD) simulations based on classical force fields has emerged as an important tool in scientific research. Partial charges are critical parameters within classical force fields and usually derived from quantum mechanical (QM) calculations. However, QM calculations are often time-consuming and prone to basis set dependence. Alternatively, fluctuating charge (FQ) models offer another avenue for partial charge derivation, which has significant speed advantages and can be used for large-scale screening. Building upon our previous work, which introduced an FQ model for zinc-containing complexes, herein we extend this model to include additional 3d transition metals which are important to the life sciences, namely chromium, manganese, iron, cobalt, and nickel. Employing CM5 charges as target for parametrization, our FQ model accurately reproduces partial charges for 3d metal complexes featuring biologically relevant ligands. Furthermore, by using atomic charges derived by our FQ model, MD simulations have been performed. These charges showed excellent performance in simulating proteomic metal sites housing multiple metal ions, specifically, a metalloprotein containing an iron-sulfur cluster and another containing a dimanganese metal site, showcasing comparable performance to those of RESP charges. We anticipate that our study can accelerate the parametrization of atomic charges for metalloproteins featuring 3d transition metals, thereby facilitating simulations of relevant systems.

Ferredoxins: Functions, Evolution, Potential Applications, and Challenges of Subtype Classification

Ferredoxins are proteins found in all biological kingdoms and are involved in essential biological processes including photosynthesis, lipid metabolism, and biogeochemical cycles. Ferredoxins are classified into different groups based on the iron-sulfur (Fe-S) clusters that they contain. A new subtype classification and nomenclature system, based on the spacing between amino acids in the Fe-S binding motif, has been proposed in order to better understand ferredoxins' biological diversity and evolutionary linkage across different organisms. This new classification system has revealed an unparalleled diversity between ferredoxins and has helped identify evolutionarily linked ferredoxins between species. The current review provides the latest insights into ferredoxin functions and evolution, and the new subtype classification, outlining their potential biotechnological applications and the future challenges in streamlining the process.

Two local minima for structures of [4Fe-4S] clusters obtained with density functional theory methods

[4Fe-4S] clusters are essential cofactors in many proteins involved in biological redox-active processes. Density functional theory (DFT) methods are widely used to study these clusters. Previous investigations have indicated that there exist two local minima for these clusters in proteins. We perform a detailed study of these minima in five proteins and two oxidation states, using combined quantum mechanical and molecular mechanical (QM/MM) methods. We show that one local minimum (L state) has longer Fe-Fe distances than the other (S state), and that the L state is more stable for all cases studied. We also show that some DFT methods may only obtain the L state, while others may obtain both states. Our work provides new insights into the structural diversity and stability of [4Fe-4S] clusters in proteins, and highlights the importance of reliable DFT methods and geometry optimization. We recommend r2SCAN for optimizing [4Fe-4S] clusters in proteins, which gives the most accurate structures for the five proteins studied.

Mechanism of Mixed-Valence Fe2.5+···Fe2.5+ Formation in Fe4S4 Clusters in the Ferredoxin Binding Motif

Most low-potential Fe₄S₄ clusters exist in the conserved binding sequence CxxCxxC (CⁿCⁿ⁺³Cⁿ⁺⁶). Fe(II) and Fe(III) at the first (Cⁿ) and third (Cⁿ⁺⁶) cysteine ligand sites form a mixed-valence Fe^2.5+···Fe^2.5+ pair in the reduced Fe(II)₃Fe(III) cluster. Here, we investigate the mechanism of how the conserved protein environment induces mixed-valence pair formation in the Fe₄S₄ clusters, F_X, F_A, and F_B in photosystem I, using a quantum mechanical/molecular mechanical approach. Exchange coupling between Fe sites is predominantly determined by the shape of the Fe₄S₄ cluster, which is stabilized by the preorganized protein electrostatic environment. The backbone NH and CO groups in the conserved CxxCxxC and adjacent helix regions orient along the Fe^Cn···Fe^C(n+6) axis, generating an electric field and stabilizing the Fe^Cn(II)Fe^C(n+6)(III) state in F_A and F_B. The overlap of the d orbitals via -S- (superexchange) is observed for the single Fe^Cn(II)···Fe^C(n+6)(III) pair, leading to the formation of the mixed-valence Fe^2.5+···Fe^2.5+ pair. In contrast, several superexchange Fe(II)···Fe(III) pairs are observed in F_X due to the highly symmetric pair of the CDGPGRGGTC sequences. This is likely the origin of F_X serving as an electron acceptor in the two electron transfer branches.

Benchmark Study of Redox Potential Calculations for Iron-Sulfur Clusters in Proteins

Redox potentials have been calculated for 12 different iron–sulfur sites of 6 different types with 1–4 iron ions. Structures were optimized with combined quantum mechanical and molecular mechanical (QM/MM) methods, and the redox potentials were calculated using the QM/MM energies, single-point QM methods in a continuum solvent or by QM/MM thermodynamic cycle perturbations. We show that the best results are obtained with a large QM system (∼300 atoms, but a smaller QM system, ∼150 atoms, can be used for the QM/MM geometry optimization) and a large value of the dielectric constant (80). For absolute redox potentials, the B3LYP density functional method gives better results than TPSS, and the results are improved with a larger basis set. However, for relative redox potentials, the opposite is true. The results are insensitive to the force field (charges of the surroundings) used for the QM/MM calculations or whether the protein and solvent outside the QM system are relaxed or kept fixed at the crystal structure. With the best approach for relative potentials, mean absolute and maximum deviations of 0.17 and 0.44 V, respectively, are obtained after removing a systematic error of −0.55 V. Such an approach can be used to identify the correct oxidation states involved in a certain redox reaction.

We have studied redox potentials of 12 iron−sulfur sites of 6 types with 1−4 iron ions. Structures were optimized with combined quantum mechanical and molecular mechanical (QM/MM) methods, and the redox potentials were calculated with QM/MM, QM calculations in a continuum solvent or by QM/MM thermodynamic cycle perturbations. The best results are obtained with the second approach using ∼300 atoms in the QM model and a large dielectric constant.

Theoretical Study on Redox Potential Control of Iron-Sulfur Cluster by Hydrogen Bonds: A Possibility of Redox Potential Programming

The effect of hydrogen bonds around the active site of Anabaena [2Fe-2S] ferredoxin (Fd) on a vertical ionization potential of the reduced state (IP(red)) is examined based on the density functional theory (DFT) calculations. The results indicate that a single hydrogen bond increases the relative stability of the reduced state, and shifts IP(red) to a reductive side by 0.31–0.33 eV, regardless of the attached sulfur atoms. In addition, the IP(red) value can be changed by the number of hydrogen bonds around the active site. The results also suggest that the redox potential of [2Fe-2S] Fd is controlled by the number of hydrogen bonds because IP(red) is considered to be a major factor in the redox potential. Furthermore, there is a possibility that the redox potentials of artificial iron-sulfur clusters can be finely controlled by the number of the hydrogen bonds attached to the sulfur atoms of the cluster.

Electron Acceptor-Donor Iron Sites in the Iron-Sulfur Cluster of Photosynthetic Electron-Transfer Pathways

In photosystem I, two electron-transfer pathways via quinones (A_1A and A_1B) are merged at the iron-sulfur Fe₄S₄ cluster F_X into a single pathway toward the other two Fe₄S₄ clusters F_A and F_B. Using a quantum mechanical/molecular mechanical approach, we identify the redox-active Fe sites in the clusters. In F_A and F_B, the Fe site, which does not belong to the CxxCxxCxxxCP motif, serves as an electron acceptor/donor. F_X has two independent electron acceptor Fe sites for A- and B-branch electron transfers, depending on the Asp-B575 protonation state, which causes the A_1A-to-F_X electron transfer to be uphill and the A_1B-to-F_X electron transfer to be downhill. The two asymmetric electron-transfer pathways from A₁ to F_X and the separation of the electron acceptor and donor Fe sites are likely associated with the specific role of F_X in merging the two electron transfer pathways into the single pathway.

Dispensability of the [4Fe-4S] cluster in novel homologues of adenine glycosylase MutY

7,8-Dihydro-8-deoxyguanine (8oG) is one of the most common oxidative lesions in DNA. DNA polymerases misincorporate an adenine across from this lesion. Thus, 8oG is a highly mutagenic lesion responsible for G:C→T:A transversions. MutY is an adenine glycosylase, part of the base excision repair pathway that removes adenines, when mispaired with 8oG or guanine. Its catalytic domain includes a [4Fe-4S] cluster motif coordinated by cysteinyl ligands. When this cluster is absent, MutY activity is depleted and several studies concluded that the [4Fe-4S] cluster motif is an indispensable component for DNA binding, substrate recognition and enzymatic activity. In the present study, we identified 46 MutY homologues that lack the canonical cysteinyl ligands, suggesting an absence of the [4Fe-4S] cluster. A phylogenetic analysis groups these novel MutYs into two different clades. One clade is exclusive of the order Lactobacillales and another clade has a mixed composition of anaerobic and microaerophilic bacteria and species from the protozoan genus Entamoeba. Structural modeling and sequence analysis suggests that the loss of the [4Fe-4S] cluster is compensated by a convergent solution in which bulky amino acids substitute the [4Fe-4S] cluster. We functionally characterized MutYs from Lactobacillus brevis and Entamoeba histolytica as representative members from each clade and found that both enzymes are active adenine glycosylases. Furthermore, chimeric glycosylases, in which the [4Fe-4S] cluster of Escherichia coli MutY is replaced by the corresponding amino acids of LbY and EhY, are also active. Our data indicates that the [4Fe-4S] cluster plays a structural role in MutYs and evidences the existence of alternative functional solutions in nature.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.

Detection of alternative conformations:Shift_plotandAC_predictionprograms

Increased mobility of poorly ordered atoms compared to well ordered ones in unrestrained refinement was adopted as the basis for the development of a method for predicting alternative conformations. Two programs for the analysis of atomic mobility are presented.Shift_plotmay be useful for visual analysis of magnitudes of atomic shifts in refinement.AC_predictionimplements automatic decision-making procedures based on atomic shift magnitudes in unrestrained refinement to highlight the specific residues that are likely to harbour alternative conformations and that hence require checkingviaelectron density maps.

Characterization of [4Fe-4S] cluster vibrations and structure in nitrogenase Fe protein at three oxidation levels via combined NRVS, EXAFS, and DFT analyses

Azotobacter vinelandii nitrogenase Fe protein (Av2) provides a rare opportunity to investigate a [4Fe-4S] cluster at three oxidation levels in the same protein environment. Here, we report the structural and vibrational changes of this cluster upon reduction using a combination of NRVS and EXAFS spectroscopies and DFT calculations. Key to this work is the synergy between these three techniques as each generates highly complementary information and their analytical methodologies are interdependent. Importantly, the spectroscopic samples contained no glassing agents. NRVS and DFT reveal a systematic 10-30 cm(-1) decrease in Fe-S stretching frequencies with each added electron. The "oxidized" [4Fe-4S](2+) state spectrum is consistent with and extends previous resonance Raman spectra. For the "reduced" [4Fe-4S](1+) state in Fe protein, and for any "all-ferrous" [4Fe-4S](0) cluster, these NRVS spectra are the first available vibrational data. NRVS simulations also allow estimation of the vibrational disorder for Fe-S and Fe-Fe distances, constraining the EXAFS analysis and allowing structural disorder to be estimated. For oxidized Av2, EXAFS and DFT indicate nearly equal Fe-Fe distances, while addition of one electron decreases the cluster symmetry. However, addition of the second electron to form the all-ferrous state induces significant structural change. EXAFS data recorded to k = 21 Å(-1) indicates a 1:1 ratio of Fe-Fe interactions at 2.56 Å and 2.75 Å, a result consistent with DFT. Broken symmetry (BS) DFT rationalizes the interplay between redox state and the Fe-S and Fe-Fe distances as predominantly spin-dependent behavior inherent to the [4Fe-4S] cluster and perturbed by the Av2 protein environment.

Dynamics of the [4Fe-4S] cluster in Pyrococcus furiosus D14C ferredoxin via nuclear resonance vibrational and resonance Raman spectroscopies, force field simulations, and density functional theory calculations

We have used (57)Fe nuclear resonance vibrational spectroscopy (NRVS) to study oxidized and reduced forms of the [4Fe-4S] cluster in the D14C variant ferredoxin from Pyrococcus furiosus (Pf D14C Fd). To assist the normal-mode assignments, we conducted NRVS with D14C ferredoxin samples with (36)S substituted into the [4Fe-4S] cluster bridging sulfide positions, and a model compound without ligand side chains, (Ph(4)P)(2)[Fe(4)S(4)Cl(4)]. Several distinct regions of NRVS intensity are identified, ranging from "protein" and torsional modes below 100 cm(-1), through bending and breathing modes near 150 cm(-1), to strong bands from Fe-S stretching modes between 250 and ∼400 cm(-1). The oxidized ferredoxin samples were also investigated by resonance Raman (RR) spectroscopy. We found good agreement between NRVS and RR frequencies, but because of different selection rules, the intensities vary dramatically between the two types of spectra. The (57)Fe partial vibrational densities of states for the oxidized samples were interpreted by normal-mode analysis with optimization of Urey-Bradley force fields for local models of the [4Fe-4S] clusters. Full protein model calculations were also conducted using a supplemented CHARMM force field, and these calculations revealed low-frequency modes that may be relevant to electron transfer with Pf Fd partners. Density functional theory (DFT) calculations complemented these empirical analyses, and DFT was used to estimate the reorganization energy associated with the [Fe(4)S(4)](2+/+) redox cycle. Overall, the NRVS technique demonstrates great promise for the observation and quantitative interpretation of the dynamical properties of Fe-S proteins.

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.

Biomimetic chemistry of iron, nickel, molybdenum, and tungsten in sulfur-ligated protein sites

Biomimetic inorganic chemistry has as its primary goal the synthesis of molecules that approach or achieve the structures, oxidation states, and electronic and reactivity features of native metal-containing sites of variant nuclearity. Comparison of properties of accurate analogues and these sites ideally provides insight into the influence of protein structure and environment on intrinsic properties as represented by the analogue. For polynuclear sites in particular, the goal provides a formidable challenge for, with the exception of iron-sulfur clusters, all such site structures have never been achieved and few have even been closely approximated by chemical synthesis. This account describes the current status of the synthetic analogue approach as applied to the mononuclear sites in certain molybdoenzymes and the polynuclear sites in hydrogenases, nitrogenase, and carbon monoxide dehydrogenases.

Insight into the protein and solvent contributions to the reduction potentials of [4Fe-4S]2+/+ clusters: crystal structures of the Allochromatium vinosum ferredoxin variants C57A and V13G and the homologous Escherichia coli ferredoxin

The crystal structures of the C57A and V13G molecular variants of Allochromatium vinosum 2[4Fe-4S] ferredoxin (AlvinFd) and that of the homologous ferredoxin from Escherichia coli (EcFd) have been determined at 1.05-, 1.48-, and 1.65-A resolution, respectively. The present structures combined with cyclic voltammetry studies establish clear effects of the degree of exposure of the cluster with the lowest reduction potential (cluster I) towards less negative reduction potentials (E degrees ). This is better illustrated by V13G AlvinFd (high exposure, E degrees = -594 mV) and EcFd (low exposure, E degrees = -675 mV). In C57A AlvinFd, the movement of the protein backbone, as a result of replacing the noncoordinating Cys57 by Ala, leads to a +50-mV upshift of the potential of the nearby cluster I, by removal of polar interactions involving the thiolate group and adjustment of the hydrogen-bond network involving the cluster atoms. In addition, the present structures and other previously reported accurate structures of this family of ferredoxins indicate that polar interactions of side chains and water molecules with cluster II sulfur atoms, which are absent in the environment of cluster I, are correlated to the approximately 180-250 mV difference between the reduction potentials of clusters I and II. These findings provide insight into the significant effects of subtle structural differences of the protein and solvent environment around the clusters of [4Fe-4S] ferredoxins on their electrochemical properties.

Structural analysis of metal sites in proteins: non-heme iron sites as a case study

In metalloproteins, the protein environment modulates metal properties to achieve the required goal, which can be protein stabilization or function. The analysis of metal sites at the atomic level of detail provided by protein structures can thus be of benefit in functional and evolutionary studies of proteins. In this work, we propose a structural bioinformatics approach to the study of metalloproteins based on structural templates of metal sites that include the PDB coordinates of protein residues forming the first and the second coordination sphere of the metal. We have applied this approach to non-heme iron sites, which have been analyzed at various levels. Templates of sites located in different protein domains have been compared, showing that similar sites can be found in unrelated proteins as the result of convergent evolution. Templates of sites located in proteins of a large superfamily have been compared, showing possible mechanisms of divergent evolution of proteins to achieve different functions. Furthermore, template comparisons have been used to predict the function of uncharacterized proteins, showing that similarity searches focused on metal sites can be advantageously combined with typical whole-domain comparisons. Structural templates of metal sites, finally, may constitute the basis for a systematic classification of metalloproteins in databases.

X-ray structure of the complete ABC enzyme ABCE1 from Pyrococcus abyssi

The ATP binding cassette enzyme ABCE1 (also known as RNase-L (ribonuclease L) inhibitor, Pixie, and HP68), one of the evolutionary most sequence-conserved enzymes, functions in translation initiation, ribosome biogenesis, and human immunodeficiency virus capsid assembly. However, its structural mechanism and biochemical role in these processes have not been revealed. We determined the crystal structure of Pyrococcus abyssi ABCE1 in complex with Mg(2+) and ADP to 2.8A resolution. ABCE1 consists of four structural domains. Two nucleotide binding domains are arranged in a head-to-tail orientation by a hinge domain, suggesting that these domains undergo the characteristic tweezers-like powerstroke of ABC enzymes. In contrast to all other known ABC enzymes, ABCE1 has a N-terminal iron-sulfur-cluster (FeS) domain. The FeS domain contains two [4Fe-4S] clusters and is structurally highly related to bacterial-type ferredoxins. However, one cluster is coordinated by an unusual CX(4)CX(3/4)C triad. Surprisingly, intimate interactions of the FeS domain with the adenine and ribose binding Y-loop on nucleotide binding domain 1 suggest a linkage between FeS domain function and ATP-induced conformational control of the ABC tandem cassette. The structure substantially expands the functional architecture of ABC enzymes and raises the possibility that ABCE1 is a chemomechanical engine linked to a redox process.

Chemical rescue of a site-modified ligand to a [4Fe–4S] cluster in PsaC, a bacterial-like dicluster ferredoxin bound to Photosystem I

Chemical rescue of site-modified amino acids using externally supplied organic molecules represents a powerful method to investigate structure-function relationships in proteins. Here we provide definitive evidence that aryl and alkyl thiolates, reagents typically used for in vitro iron-sulfur cluster reconstitutions, serve as rescue ligands to a site-specifically modified [4Fe-4S](1+,2+) cluster in PsaC, a bacterial dicluster ferredoxin-like subunit of Photosystem I. PsaC binds two low-potential [4Fe-4S](1+,2+) clusters termed F(A) and F(B). In the C13G/C33S variant of PsaC, glycine has replaced cysteine at position 13 creating a protein that is missing one of the ligating amino acids to iron-sulfur cluster F(B). Using a variety of analytical techniques, including non-heme iron and acid-labile sulfur assays, and EPR, resonance Raman, and Mössbauer spectroscopies, we showed that the C13G/C33S variant of PsaC binds two [4Fe-4S](1+,2+) clusters, despite the absence of one of the biological ligands. (19)F NMR spectroscopy indicated that the external thiolate replaces cysteine 13 as a substitute ligand to the F(B) cluster. The finding that site-modified [4Fe-4S](1+,2+) clusters can be chemically rescued with external thiolates opens new opportunities for modulating their properties in proteins. In particular, it provides a mechanism to attach an additional electron transfer cofactor to the protein via a bound, external ligand.

How sensitive are nanosecond molecular dynamics simulations of proteins to changes in the force field?

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

The structure of the 2[4Fe-4S] ferredoxin from Pseudomonas aeruginosa at 1.32-A resolution: comparison with other high-resolution structures of ferredoxins and contributing structural features to reduction potential values

The structure of the 2[4Fe-4S] ferredoxin (PaFd) from Pseudomonas aeruginosa, which belongs to the Allochromatium vinosum (Alvin) subfamily, has been determined by X-ray crystallography at 1.32-A resolution, which is the highest up to now for a member of this subfamily of Fds. The main structural features of PaFd are similar to those of AlvinFd. However, the significantly higher resolution of the PaFd structure makes possible a reliable comparison with available high-resolution structures of [4Fe-4S]-containing Fds, in an effort to rationalize the unusual electrochemical properties of Alvin-like Fds. Three major factors contributing to the reduction potential values of [4Fe-4S]2+/+ clusters of Fds, namely, the surface accessibility of the clusters, the N-H...S hydrogen-bonding network, and the volume of the cavities hosting the clusters, are extensively discussed. The volume of the cavities is introduced in the present work for the first time, and can in part explain the very negative potential of cluster I of Alvin-like Fds.

Dynamics of an [Fe4S4(SPh)4]2- cluster explored via IR, Raman, and nuclear resonance vibrational spectroscopy (NRVS)-analysis using 36S substitution, DFT calculations, and empirical force fields

We have used four vibrational spectroscopies--FT-IR, FT-Raman, resonance Raman, and 57Fe nuclear resonance vibrational spectroscopy (NRVS)--to study the normal modes of the Fe-S cluster in [(n-Bu)4N]2[Fe4S4(SPh)4]. This [Fe4S4(SR)4]2- complex serves as a model for the clusters in 4Fe ferredoxins and high-potential iron proteins (HiPIPs). The IR spectra exhibited differences above and below the 243 K phase transition. Significant shifts with 36S substitution into the bridging S positions were also observed. The NRVS results were in good agreement with the low temperature data from the conventional spectroscopies. The NRVS spectra were interpreted by normal mode analysis using optimized Urey-Bradley force fields (UBFF) as well as from DFT theory. For the UBFF calculations, the parameters were refined by comparing calculated and observed NRVS frequencies and intensities. The frequency shifts after 36S substitution were used as an additional constraint. A D 2d symmetry Fe4S4S'4 model could explain most of the observed frequencies, but a better match to the observed intensities was obtained when the ligand aromatic rings were included for a D 2d Fe4S4(SPh)4 model. The best results were obtained using the low temperature structure without symmetry constraints. In addition to stretching and bending vibrations, low frequency modes between approximately 50 and 100 cm(-1) were observed. These modes, which have not been seen before, are interpreted as twisting motions with opposing sides of the cube rotating in opposite directions. In contrast with a recent paper on a related Fe4S4 cluster, we find no need to assign a large fraction of the low frequency NRVS intensity to 'rotational lattice modes'. We also reassign the 430 cm(-1) band as primarily an elongation of the thiophenolate ring, with approximately 10% terminal Fe-S stretch character. This study illustrates the benefits of combining NRVS with conventional Raman and IR analysis for characterization of Fe-S centers. DFT theory is shown to provide remarkable agreement with the experimental NRVS data. These results provide a reference point for the analysis of more complex Fe-S clusters in proteins.

Identification of variant molecules of Bacillus thermoproteolyticus ferredoxin: crystal structure reveals bound coenzyme A and an unexpected [3Fe-4S] cluster associated with a canonical [4Fe-4S] ligand motif

During the purification of recombinant Bacillus thermoproteolyticus ferredoxin (BtFd) from Escherichia coli, we have noted that some Fe-S proteins were produced in relatively small amounts compared to the originally identified BtFd carrying a [4Fe-4S] cluster. These variants could be purified into three Fe-S protein components (designated as V-I, V-II, and V-III) by standard chromatography procedures. UV-vis and EPR spectroscopic analyses indicated that each of these variants accommodates a [3Fe-4S] cluster. From mass spectrometric and protein sequence analyses together with native and SDS gel electrophoresis, we established that V-I and V-II contain the polypeptide of BtFd associated with acyl carrier protein (ACP) and with coenzyme A (CoA), respectively, and that V-III is a BtFd dimer linked by a disulfide bond. The crystal structure of the BtFd-CoA complex (V-II) determined at 1.6 A resolution revealed that each of the four complexes in the crystallographic asymmetric unit possesses a [3Fe-4S] cluster that is coordinated by Cys(11), Cys(17), and Cys(61). The polypeptide chain of each complex is superimposable onto that of the original [4Fe-4S] BtFd except for the segment containing Cys(14), the fourth ligand to the [4Fe-4S] cluster of BtFd. In the variant molecules, the side chain of Cys(14) is rotated away to the molecular surface, forming a disulfide bond with the terminal sulfhydryl group of CoA. This covalent modification may have occurred in vivo, thereby preventing the assembly of the [4Fe-4S] cluster as observed previously for Desulfovibrio gigas ferredoxin. Possibilities concerning how the variant molecules are formed in the cell are discussed.

N-H···S Hydrogen Bonds in a Ferredoxin Model

The Fe4S4 complex {(CH3)3NCH2CONH2}2[Fe4S4((tBuS)4] (1) was synthesized to replicate the ferredoxin active site with a subset of its N-H...S hydrogen bonds. The two cationic counterions mimic the polypeptide backbone of ferredoxin (Fd) as amide hydrogen-bond donors to sulfur atoms of the iron-sulfur cluster. X-ray crystallographic data show that the organic sulfur (Sgamma) of one tert-butylthiolate ligand and one inorganic sulfur of the cluster core serve as N-H...S hydrogen-bond acceptors. The cluster core of complex 1 is tetragonally elongated in contrast to that of Fd, which is tetragonally compressed. This is the first observation of an elongated [Fe4S4]2+ cluster core. Additionally, this is the first synthetic Fd model in which N-H...S hydrogen bonding to a cluster has been achieved.

Ligand K-Edge X-ray Absorption Spectroscopy of [Fe4S4]1+,2+,3+ Clusters: Changes in Bonding and Electronic Relaxation upon Redox

Sulfur K-edge X-ray absorption spectroscopy (XAS) is reported for [Fe(4)S(4)](1+,2+,3+) clusters. The results are quantitatively and qualitatively compared with DFT calculations. The change in covalency upon redox in both the [Fe(4)S(4)](1+/2+) (ferredoxin) and the [Fe(4)S(4)](2+/3+) (HiPIP) couple are much larger than that expected from just the change in number of 3d holes. Moreover, the change in the HiPIP couple is higher than that of the ferredoxin couple. These changes in electronic structure are analyzed using DFT calculations in terms of contributions from the nature of the redox active molecular orbital (RAMO) and electronic relaxation. The results indicate that the RAMO of HiPIP has 50% ligand character, and hence, the HiPIP redox couple involves limited electronic relaxation. Alternatively, the RAMO of the ferredoxin couple is metal-based, and the ferredoxin redox couple involves extensive electronic relaxation. The contributions of these RAMO differences to ET processes in the different proteins are discussed.

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.

The circumsphere as a tool to assess distortion in [4Fe-4S] atom clusters

The geometry proposition that "four points not in a plane describe one and only one sphere" provides a novel tool for analyzing protein-induced distortions in [4Fe-4S] clusters. A geometrically perfect reference structure comprises interlaced, regular tetrahedra of Fe, S, and S gamma atoms having T(d) symmetry. Three circumspheres are defined by the three sets of four atoms, the circumcenters of which are unique points within the cluster. The structure is thus re-defined by the positions of the circumcenters in xyz space and the r, theta, phi of each atom on its respective sphere. Analysis of 12 high-resolution structures of protein-bound and small molecule [4Fe-4S](SR)(4) clusters revealed: (a) the circumcenters are generally non-coincident by approximately 0.01 to approximately 0.06 A; (b) the Fe radius, r(Fe), is nominally independent of core oxidation state, having values between 1.66 to 1.69 A, whereas r(S) and r(SG), which have ranges of 2.18-2.24 A and 3.87-3.94 A, respectively, both increase by as much as approximately 3% upon reduction from the 3+ to the 1+ core valence; (c) deviation of some atoms from the theta, phi of a perfect tetrahedron can be large, approximately 10 degrees, and sets of atoms can show patterns of motion on their spheres that result from changes in Fe-S bond lengths. Density functional theory calculations suggest that the [4Fe-4S] core itself requires rather little energy to distort (approximately 2 kcal/mol), whereas significantly more energy is required to distort the Sgamma shell (~4 kcal/mol) to that of cluster I in Clostridium acidurici ferredoxin.