Revealing the magnetostructural dynamics of [2Fe-2S] ferredoxins from reduced-dimensionality analysis of antiferromagnetic exchange coupling fluctuations

Effect of H bond removal and changes in the position of the iron-sulphur head domain on the spin-lattice relaxation properties of the [2Fe-2S](2+) Rieske cluster in cytochrome bc(1)

Here, comparative electron spin–lattice relaxation studies of the 2Fe–2S iron–sulphur (Fe–S) cluster embedded in a large membrane protein complex – cytochrome bc₁ – are reported.

Here, comparative electron spin–lattice relaxation studies of the 2Fe–2S iron–sulphur (Fe–S) cluster embedded in a large membrane protein complex – cytochrome bc₁ – are reported. Structural modifications of the local environment alone (mutations S158A and Y160W removing specific H bonds between Fe–S and amino acid side chains) or in combination with changes in global protein conformation (mutations/inhibitors changing the position of the Fe–S binding domain within the protein complex) resulted in different redox potentials as well as g-, g-strain and the relaxation rates (T₁^–1) for the Fe–S cluster. The relaxation rates for T < 25 K were measured directly by inversion recovery, while for T > 60 K they were deduced from simulation of continuous wave EPR spectra of the cluster using a model that included anisotropy of Lorentzian broadening. In all cases, the relaxation rate involved contributions from direct, second-order Raman and Orbach processes, each dominating over different temperature ranges. The analysis of T₁^–1 (T) over the range 5–120 K yielded the values of the Orbach energy (E_Orb), Debye temperature θ_D and Raman process efficiency C_Ram for each variant of the protein. As the Orbach energy was generally higher for mutants S158A and Y160W, compared to wild-type protein (WT), it is suggested that H bond removal influences the geometry leading to increased strength of antiferromagnetic coupling between two Fe ions of the cluster. While θ_D was similar for all variants (∼107 K), the efficiency of the Raman process generally depends on the spin–orbit coupling that is lower for S158A and Y160W mutants, when compared to the WT. However, in several cases C_Ram did not only correlate with spin–orbit coupling but was also influenced by other factors – possibly the modification of protein rigidity and therefore the vibrational modes around the Fe–S cluster that change upon the movement of the iron–sulphur head domain.

Magnetostructural dynamics of Rieske versus ferredoxin iron-sulfur cofactors

The local chemical environment of the [2Fe-2S] cofactor hosted by ferredoxin and Rieske-type proteins is fundamentally different due to the presence of distinct ligands at the two iron centers in the case of Rieske proteins, whereas they are identical in ferredoxins. This renders Rieske [2Fe-2S] cores chemically asymmetric and results in more complex vibrational spectra as compared to ferredoxin. Likewise, one would expect other properties, for instance the dynamics of the magnetic exchange coupling constant J, to be also more complex. Applying ab initio molecular dynamics using our recently introduced spin-constrained two-determinant extended broken symmetry (CEBS) approach to Rieske and ferredoxin model complexes at 300 K, we extract the molecular fluctuations and the resulting magnetostructural cross-correlations involving the antiferromagnetic exchange interaction J(t). This analysis demonstrates that the details of the magnetostructural dynamics are indeed distinctly different for Rieske and ferredoxin cofactors, while the time averages of 〈J〉 are shown to be essentially identical. In particular, the frequency window between about 200 and 350 cm(-1), is a "fingerprint region" that allows one to distinguish chemically asymmetric from symmetric cofactors and thus Rieske proteins from ferredoxins.

The iron-sulfur core in Rieske proteins is not symmetric

At variance with ferredoxins, Rieske-type proteins contain a chemically asymmetric iron-sulfur cluster. Nevertheless, X-ray crystallography apparently finds their [2Fe-2S] cores to be structurally symmetric or very close to symmetric (i.e. the four iron-sulfur bonds in the [2Fe-2S] core are equidistant). Using a combination of advanced density-based approaches, including finite-temperature molecular dynamics to access thermal fluctuations and free-energy profiles, in conjunction with correlated wavefunction-based methods we clearly predict an asymmetric core structure. This reveals a fundamental and intrinsic difference within the iron-sulfur clusters hosted by Rieske proteins and ferredoxins and thus opens up a new dimension for the ongoing efforts in understanding the role of Rieske-type [2Fe-2S] cluster in electron transfer processes that occur in almost all biological systems.

Charge transfer in strongly correlated systems: an exact diagonalization approach to model Hamiltonians

We study charge transfer in bridged di- and triruthenium complexes from a theoretical and computational point of view. Ab initio computations are interpreted from the perspective of a simple empirical Hamiltonian, a chemically specific Mott-Hubbard model of the complexes' π electron systems. This Hamiltonian is coupled to classical harmonic oscillators mimicking a polarizable dielectric environment. The model can be solved without further approximations in a valence bond picture using the method of exact diagonalization and permits the computation of charge transfer reaction rates in the framework of Marcus' theory. In comparison to the exact solution, the Hartree-Fock mean field theory overestimates both the activation barrier and the magnitude of charge-transfer excitations significantly. For triruthenium complexes, we are able to directly access the interruthenium antiferromagnetic coupling strengths.

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.