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

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.

Robust Production, Crystallization, Structure Determination, and Analysis of [Fe-S] Proteins: Uncovering Control of Electron Shuttling and Gating in the Respiratory Metabolism of Molybdopterin Guanine Dinucleotide Enzymes

[Fe-S] clusters are essential cofactors in all domains of life. They play many biological roles due to their unique abilities for electron transfer and conformational control. Yet, producing and analyzing Fe-S proteins can be difficult and even misleading if not done anaerobically. Due to unique redox properties of [Fe-S] clusters and their oxygen sensitivity, they pose multiple challenges and can lose enzymatic activity or cause their component proteins to be structurally disordered due to [Fe-S] cluster oxidation and loss in air. Here we highlight tested protocols and strategies enabling efficient and stable [Fe-S] protein production, purification, crystallization, X-ray diffraction data collection, and structure determination. From multiple high-resolution anaerobic crystal structures, we furthermore analyze exemplary data defining [Fe-S] clusters, substrate entry, and product exit for the functional oxidation states of type II molybdo-bis(molybdopterin guanine dinucleotide) (Mo-bisMGD) enzymes. Notably, these enzymes perform electron shuttling between quinone pools and specific substrates to catalyze respiratory metabolism. The identified structure-activity relationships for this enzyme class have broad implications germane to perchlorate environments on Earth and Mars extending to an alternative mechanism underlying metabolic origins for the evolution of the oxygen atmosphere. Integrated structural analyses of type II Mo-bisMGD enzymes unveil novel distinctive shared molecular mechanisms for dynamic control of substrate entry and product release gated by hydrophobic residues. Collective findings support a prototypic model for type II Mo-bisMGD enzymes including insights for a fundamental molecular mechanistic understanding of selectivity and regulation by a conformationally gated channel with general implications for [Fe-S] cluster respiratory enzymes.

Density Functional Theory Calculation of Bonding and Charge Parameters for Molecular Dynamics Studies on [FeFe] Hydrogenases

We have developed and tested molecular mechanics parameters for [FeS] clusters found in known [FeFe] hydrogenases. Bond stretching, angle bending, dihedral and improper torsion parameters for models of the oxidized and reduced catalytic H-cluster, [4Fe4S](+,2+)Cys4, [4Fe4S](+,2+)Cys3His, and [2Fe2S](+,2+)Cys4, were calculated solely from Kohn-Sham density functional theory and Natural Population Analysis. Circumsphere analysis of the cubane clusters in the energy-minimized structure of the full Clostridium pasteurianum hydrogenase I showed the resulting metallocluster structures to be similar to known cubane structures. All clusters were additionally stable in molecular dynamics simulations over the course of 1.0 ns in the fully oxidized and fully reduced enzyme models. Normal modes calculated by quasiharmonic analysis from the dynamics data show unexpected couplings among internal coordinate motions, which may reflect the effects of the protein structure on metallocluster dynamics.

Structural principles for computational and de novo design of 4Fe-4S metalloproteins

Iron-sulfur centers in metalloproteins can access multiple oxidation states over a broad range of potentials, allowing them to participate in a variety of electron transfer reactions and serving as catalysts for high-energy redox processes. The nitrogenase FeMoCO cluster converts di-nitrogen to ammonia in an eight-electron transfer step. The 2(Fe4S4) containing bacterial ferredoxin is an evolutionarily ancient metalloprotein fold and is thought to be a primordial progenitor of extant oxidoreductases. Controlling chemical transformations mediated by iron-sulfur centers such as nitrogen fixation, hydrogen production as well as electron transfer reactions involved in photosynthesis are of tremendous importance for sustainable chemistry and energy production initiatives. As such, there is significant interest in the design of iron-sulfur proteins as minimal models to gain fundamental understanding of complex natural systems and as lead-molecules for industrial and energy applications. Herein, we discuss salient structural characteristics of natural iron-sulfur proteins and how they guide principles for design. Model structures of past designs are analyzed in the context of these principles and potential directions for enhanced designs are presented, and new areas of iron-sulfur protein design are proposed. This article is part of a Special issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, protein networks, edited by Ronald L. Koder and J.L Ross Anderson.

Structural Analysis of Cubane-Type Iron Clusters

The generalized cluster type [M₄(μ₃-Q)₄L _n ] ^x contains the cubane-type [M₄Q₄] ^z core unit that can approach, but typically deviates from, perfect T_d symmetry. The geometric properties of this structure have been analyzed with reference to T_d symmetry by a new protocol. Using coordinates of M and Q atoms, expressions have been derived for interatomic separations, bond angles, and volumes of tetrahedral core units (M₄, Q₄) and the total [M₄Q₄] core (as a tetracapped M₄ tetrahedron). Values for structural parameters have been calculated from observed average values for a given cluster type. Comparison of calculated and observed values measures the extent of deviation of a given parameter from that required in an exact tetrahedral structure. The procedure has been applied to the structures of over 130 clusters containing [Fe₄Q₄] (Q = S^2-, Se^2-, Te^2-, [NPR₃]^-, [NR]^2-) units, of which synthetic and biological sulfide-bridged clusters constitute the largest subset. General structural features and trends in structural parameters are identified and summarized. An extensive database of structural properties (distances, angles, volumes) has been compiled in Supporting Information.

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.

Stabilization of 3:1 site-differentiated cubane-type clusters in the [Fe(4)S(4)](1+) core oxidation state by tertiary phosphine ligation: synthesis, core structural diversity, and S = 1/2 ground states

An extensive series of 3:1 site-differentiated cubane-type clusters [Fe(4)S(4)(PPr(i)(3))(3)L] (L = Cl(-), Br(-), I(-), RO(-), RS(-), RSe(-)) has been prepared in 40-80% yield by two methods: ligand substitution of [Fe(4)S(4)(PPr(i)(3))(4)](1+) in tetrahydrofuran (THF)/acetonitrile by reaction with monoanions, and reductive cleavage of ligand substrates (RSSR, RSeSeR, I(2)) by the all-ferrous clusters [Fe(8)S(8)(PPr(i)(3))(6)]/[Fe(16)S(16)(PPr(i)(3))(8)] in THF. These neutral clusters are stable and do not undergo ligand redistribution reactions involving charged species in benzene and THF solutions. X-ray structural studies confirm the cubane stereochemistry but with substantial and variable distortions of the [Fe(4)S(4)](1+) core from idealized cubic core geometry. Based on Fe-S bond lengths, seven clusters were found to have compressed tetragonal distortions (4 short and 8 long bonds), and the remaining seven display other types of distortions with different combinations of long, short, and intermediate bond lengths. These results further emphasize the facile deformabililty of this core oxidation state previously observed in [Fe(4)S(4)(SR)(4)](3-) clusters. The Fe(2.25+) mean oxidation state was demonstrated from (57)Fe isomer shifts, and the appearance of two quadrupole doublets arises from the spin-coupled |9/2,4,1/2> state. The S = 1/2 ground state was further supported by electron paramagnetic resonance spectra and magnetic susceptibility data.

Electron spin relaxation enhancement measurements of interspin distances in human, porcine, and Rhodobacter electron transfer flavoprotein-ubiquinone oxidoreductase (ETF-QO)

Electron transfer flavoprotein-ubiquinone oxidoreductase (ETF-QO) is a membrane-bound electron transfer protein that links primary flavoprotein dehydrogenases with the main respiratory chain. Human, porcine, and Rhodobacter sphaeroides ETF-QO each contain a single [4Fe-4S](2+,1+) cluster and one equivalent of FAD, which are diamagnetic in the isolated enzyme and become paramagnetic on reduction with the enzymatic electron donor or with dithionite. The anionic flavin semiquinone can be reduced further to diamagnetic hydroquinone. The redox potentials for the three redox couples are so similar that it is not possible to poise the proteins in a state where both the [4Fe-4S](+) cluster and the flavoquinone are fully in the paramagnetic form. Inversion recovery was used to measure the electron spin-lattice relaxation rates for the [4Fe-4S](+) between 8 and 18K and for semiquinone between 25 and 65K. At higher temperatures the spin-lattice relaxation rates for the [4Fe-4S](+) were calculated from the temperature-dependent contributions to the continuous wave linewidths. Although mixtures of the redox states are present, it was possible to analyze the enhancement of the electron spin relaxation of the FAD semiquinone signal due to dipolar interaction with the more rapidly relaxing [4Fe-4S](+) and obtain point-dipole interspin distances of 18.6+/-1A for the three proteins. The point-dipole distances are within experimental uncertainty of the value calculated based on the crystal structure of porcine ETF-QO when spin delocalization is taken into account. The results demonstrate that electron spin relaxation enhancement can be used to measure distances in redox poised proteins even when several redox states are present.

Substrate recognition, protein dynamics, and iron-sulfur cluster in Pseudomonas aeruginosa adenosine 5'-phosphosulfate reductase

APS reductase catalyzes the first committed step of reductive sulfate assimilation in pathogenic bacteria, including Mycobacterium tuberculosis, and is a promising target for drug development. We report the 2.7 A resolution crystal structure of Pseudomonas aeruginosa APS reductase in the thiosulfonate intermediate form of the catalytic cycle and with substrate bound. The structure, high-resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry, and quantitative kinetic analysis, establish that the two chemically discrete steps of the overall reaction take place at distinct sites on the enzyme, mediated via conformational flexibility of the C-terminal 18 residues. The results address the mechanism by which sulfonucleotide reductases protect the covalent but labile enzyme-intermediate before release of sulfite by the protein cofactor thioredoxin. P. aeruginosa APS reductase contains an [4Fe-4S] cluster that is essential for catalysis. The structure reveals an unusual mode of cluster coordination by tandem cysteine residues and suggests how this arrangement might facilitate conformational change and cluster interaction with the substrate. Assimilatory 3'-phosphoadenosine 5'-phosphosulfate (PAPS) reductases are evolutionarily related, homologous enzymes that catalyze the same overall reaction, but do so in the absence of an [Fe-S] cluster. The APS reductase structure reveals adaptive use of a phosphate-binding loop for recognition of the APS O3' hydroxyl group, or the PAPS 3'-phosphate group.

Structure, redox, pK a, spin. A golden tetrad for understanding metalloenzyme energetics and reaction pathways

After a review of the current status of density functional theory (DFT) for spin-polarized and spin-coupled systems, we focus on the resting states and intermediates of redox-active metalloenzymes and electron transfer proteins, showing how comparisons of DFT-calculated spectroscopic parameters with experiment and evaluation of related energies and geometries provide important information. The topics we examine include (1) models for the active-site structure of methane monooxygenase intermediate Q and ribonucleotide reductase intermediate X; (2) the coupling of electron transfer to proton transfer in manganese superoxide dismutase, with implications for reaction kinetics; (3) redox, pK(a), and electronic structure issues in the Rieske iron-sulfur protein, including their connection to coupled electron/proton transfer, and an analysis of how partial electron delocalization strongly alters the electron paramagnetic resonance spectrum; (4) the connection between protein-induced structural distortion and the electronic structure of oxidized high-potential 4Fe4S proteins with implications for cluster reactivity; (5) an analysis of cluster assembly and central-atom insertion into the FeMo cofactor center of nitrogenase based on DFT structural and redox potential calculations.