Models of the bis-histidine-coordinated ferricytochromes: Mössbauer and EPR spectroscopic studies of low-spin iron(III) tetrapyrroles of various electronic ground states and axial ligand orientations

57 Fe Mössbauer Spectroscopy Characterization of Electrocatalysts

This work addresses the importance of Mössbauer spectroscopy for the characterization of iron-containing electrocatalysts. The most important aspects of electrocatalysis and Mössbauer spectroscopy are summarized. Next, Fe-N-C catalysts and important conclusions made by this technique on preparation, active site identification and degradation are summarized. Furthermore, recent highlights derived for other iron-containing electrocatalysts are summarized.

A Redox-Induced Spin-State Cascade in a Mixed-Valent Fe3 (μ3 -O) Triangle

One-electron reduction of a pyrazolate-bridged triangular Fe₃ (μ₃ -O) core induces a cascade wherein all three metal centers switch from high-spin Fe³⁺ to low-spin Fe^2.66+ . This hypothesis is supported by spectroscopic data (¹ H-NMR, UV-vis-NIR, infra-red, ⁵⁷ Fe-Mössbauer, EPR), X-ray crystallographic characterization of the cluster in both oxidation states and also density functional theory. The reduction induces substantial contraction in all bond lengths around the metal centers, along with diagnostic shifts in the spectroscopic parameters. This is, to the best of our knowledge, the first example of a one-electron redox event causing concerted change in multiple iron centers.

Magnetic properties of a Kramers doublet. An univocal bridge between experimental results and theoretical predictions

The magnetic response of a Kramers doublet is analyzed in a general case taking into account only the formal properties derived from time reversal operation. It leads to a definition of a matrix G (gyromagnetic matrix) whose expression depends on the chosen reference frame and on the Kramers conjugate basis used to describe the physical system. It is shown that there exists a reference frame and a suitable Kramers conjugate basis that gives a diagonal form for the G-matrix with all non-null elements having the same sign. A detailed procedure for obtaining this canonical expression of G is presented when the electronic structure of the KD is known regardless the level of the used theory. This procedure provides a univocal way to compare the theoretical predictions with the experimental results obtained from a complete set of magnetic experiments. In this way the problems arising from ambiguities in the g-tensor definition are overcome. This procedure is extended to find a spin-Hamiltonian suitable for describing the magnetic behavior of a pair of weakly coupled Kramers systems in the multispin scheme when the interaction between the two moieties as well as the individual Zeeman interaction are small enough as compared with ligand field splitting. Explicit relations between the physical interaction and the parameters of such a spin-Hamiltonian are also obtained.

Determination of the principal g-values of Type I or highly-anisotropic low spin (HALS) ferriheme centers in frozen solutions

Continuous wave (CW) electron paramagnetic resonance (EPR) spectroscopy of highly-anisotropic low spin (HALS) ferric heme centers in frozen solutions is not a very informative approach because usually only one feature is reliably observed in the spectra, that at the maximal principal g-value of, typically, 3.3-3.79. The other two EPR turning points are severely broadened by g-strain and are not easily observed in the first-derivative CW EPR spectra. In this work, we have explored the potential of alternative EPR techniques, the electron spin echo (ESE) field sweep and electron spin transient nutation (TN), for obtaining information about the g-tensors of such systems, using as an example a typical HALS ferric heme center, [Fe(III)((15)N-coproporphyrin)(CN)2]. The analysis of the experimental g-tensor of [Fe(III)((15)N-coproporphyrin)(CN)2](-) has shown that the widths of the underlying energy distributions for this HALS center are comparable to those found for the rhombic bis-imidazole complex. The greater effect on the g-value distributions for HALS centers is determined by near degeneracy of two of the three lower-energy d-orbitals, d(yz) and d(xz), which contain the unpaired electron.

Nuclear inelastic scattering and Mössbauer spectroscopy as local probes for ligand binding modes and electronic properties in proteins: vibrational behavior of a ferriheme center inside a β-barrel protein

In this work, we present a study of the influence of the protein matrix on its ability to tune the binding of small ligands such as NO, cyanide (CN(-)), and histamine to the ferric heme iron center in the NO-storage and -transport protein Nitrophorin 2 (NP2) from the salivary glands of the blood-sucking insect Rhodnius prolixus. Conventional Mössbauer spectroscopy shows a diamagnetic ground state of the NP2-NO complex and Type I and II electronic ground states of the NP2-CN(-) and NP2-histamine complex, respectively. The change in the vibrational signature of the protein upon ligand binding has been monitored by Nuclear Inelastic Scattering (NIS), also called Nuclear Resonant Vibrational Spectroscopy (NRVS). The NIS data thus obtained have also been calculated by quantum mechanical (QM) density functional theory (DFT) coupled with molecular mechanics (MM) methods. The calculations presented here show that the heme ruffling in NP2 is a consequence of the interaction with the protein matrix. Structure optimizations of the heme and its ligands with DFT retain the characteristic saddling and ruffling only if the protein matrix is taken into account. Furthermore, simulations of the NIS data by QM/MM calculations suggest that the pH dependence of the binding of NO, but not of CN(-) and histamine, might be a consequence of the protonation state of the heme carboxyls.

Modulation of function in a minimalist heme-binding membrane protein

De novo designed heme-binding proteins have been used successfully to recapitulate features of natural hemoproteins. This approach has now been extended to membrane-soluble model proteins. Our group designed a functional hemoprotein, ME1, by engineering a bishistidine binding site into a natural membrane protein, glycophorin A (Cordova et al. in J Am Chem Soc 129:512-518, 2007). ME1 binds iron(III) protoporphyrin IX with submicromolar affinity, has a redox potential of -128 mV, and displays peroxidase activity. Here, we show the effect of aromatic residues in modulating the redox potential in the context of a membrane-soluble model system. We designed aromatic interactions with the heme through a single-point mutant, G25F, in which a phenylalanine is designed to dock against the porphyrin ring. This mutation results in roughly tenfold tighter binding to iron(III) protoporphyrin IX (K(d,app) = 6.5 × 10(-8) M), and lowers the redox potential of the cofactor to -172 mV. This work demonstrates that specific design features aimed at controlling the properties of bound cofactors can be introduced in a minimalist membrane hemoprotein model. The ability to modulate the redox potential of cofactors embedded in artificial membrane proteins is crucial for the design of electron transfer chains across membranes in functional photosynthetic devices.

Identification of a thiol/disulfide redox switch in the human BK channel that controls its affinity for heme and CO

Heme is a required prosthetic group in many electron transfer proteins and redox enzymes. The human BK channel, which is a large-conductance Ca(2+) and voltage-activated K(+) channel, is involved in the hypoxic response in the carotid body. The BK channel has been shown to bind and undergo inhibition by heme and activation by CO. Furthermore, evidence suggests that human heme oxygenase-2 (HO2) acts as an oxygen sensor and CO donor that can form a protein complex with the BK channel. Here we describe a thiol/disulfide redox switch in the human BK channel and biochemical experiments of heme, CO, and HO2 binding to a 134-residue region within the cytoplasmic domain of the channel. This region, called the heme binding domain (HBD) forms a linker segment between two Ca(2+)-sensing domains (called RCK1 and RCK2) of the BK channel. The HBD includes a CXXCH motif in which histidine serves as the axial heme ligand and the two cysteine residues can form a reversible thiol/disulfide redox switch that regulates affinity of the HBD for heme. The reduced dithiol state binds heme (K(d) = 210 nm) 14-fold more tightly than the oxidized disulfide state. Furthermore, the HBD is shown to tightly bind CO (K(d) = 50 nm) with the Cys residues in the CXXCH motif regulating affinity of the HBD for CO. This HBD is also shown to interact with heme oxygenase-2. We propose that the thiol/disulfide switch in the HBD is a mechanism by which activity of the BK channel can respond quickly and reversibly to changes in the redox state of the cell, especially as it switches between hypoxic and normoxic conditions.

Low-temperature EPR and Mössbauer spectroscopy of two cytochromes with His-Met axial coordination exhibiting HALS signals

C-type cytochromes with histidine-methionine (His-Met) iron coordination play important roles in electron-transfer reactions and in enzymes. Low-temperature electron paramagnetic resonance (EPR) spectra of low-spin ferric cytochromes c can be divided into two groups, depending on the spread of g values: the normal rhombic ones with small g anisotropy and g(max) below 3.2, and those featuring large g anisotropy with g(max) between 3.3 and 3.8, also denoted as highly axial low spin (HALS) species. Herein we present the detailed magnetic properties of cytochrome c(553) from Bacillus pasteurii (g(max) 3.36) and cytochrome c(552) from Nitrosomonas europaea (g(max) 3.34) over the pH range 6.2 to 8.2. Besides being structurally very similar, cytochrome c(553) shows the presence of a minor rhombic species at pH 6.2 (6 %), whereas cytochrome c(552) has about 25 % rhombic species over pH 7.5. The detailed Mössbauer analysis of cytochrome c(552) confirms the presence of these two low-spin ferric species (HALS and rhombic) together with an 8 % ferrous form with parameters comparable to the horse cytochrome c. Both EPR and Mössbauer data of axial cytochromes c with His-Met iron coordination are consistent with an electronic (d(xy))(2) (d(xz))(2) (d(yz))(1) ground state, which is typical for Type I model hemes.

Distinct reaction pathways followed upon reduction of oxy-heme oxygenase and oxy-myoglobin as characterized by Mössbauer spectroscopy

Activation of O(2) by heme-containing monooxygenases generally commences with the common initial steps of reduction to the ferrous heme and binding of O(2) followed by a one-electron reduction of the O(2)-bound heme. Subsequent steps that generate reactive oxygen intermediates diverge and reflect the effects of protein control on the reaction pathway. In this study, Mössbauer and EPR spectroscopies were used to characterize the electronic states and reaction pathways of reactive oxygen intermediates generated by 77 K radiolytic cryoreduction and subsequent annealing of oxy-heme oxygenase (HO) and oxy-myoglobin (Mb). The results confirm that one-electron reduction of (Fe(II)-O(2))HO is accompanied by protonation of the bound O(2) to generate a low-spin (Fe(III)-O(2)H(-))HO that undergoes self-hydroxylation to form the alpha-meso-hydroxyhemin-HO product. In contrast, one-electron reduction of (Fe(II)-O(2))Mb yields a low-spin (Fe(III)-O(2)(2-))Mb. Protonation of this intermediate generates (Fe(III)-O(2)H(-))Mb, which then decays to a ferryl complex, (Fe(IV)=O(2-))Mb, that exhibits magnetic properties characteristic of the compound II species generated in the reactions of peroxide with heme peroxidases and with Mb. Generation of reactive high-valent states with ferryl species via hydroperoxo intermediates is believed to be the key oxygen-activation steps involved in the catalytic cycles of P450-type monooxygenases. The Mössbauer data presented here provide direct spectroscopic evidence supporting the idea that ferric-hydroperoxo hemes are indeed the precursors of the reactive ferryl intermediates. The fact that a ferryl intermediate does not accumulate in HO underscores the determining role played by protein structure in controlling the reactivity of reaction intermediates.

Pentacoordinate and hexacoordinate ferric hemes in acid medium: EPR, UV–Vis and CD studies of the giant extracellular hemoglobin of Glossoscolex paulistus

The equilibrium complexity involving different axially coordinated hemes is peculiar to hemoglobins. The pH dependence of the spontaneous exchange of ligands in the extracellular hemoglobin from Glossoscolex paulistus was studied using UV-Vis, EPR, and CD spectroscopies. This protein has a complex oligomeric assembly with molecular weight of 3.1 MDa that presents an important cooperative effect. A complex coexistence of different species was observed in almost all pH values, except pH 7.0, where just aquomet species is present. Four new species were formed and coexist with the aquomethemoglobin upon acidification: (i) a "pure" low-spin hemichrome (Type II), also called hemichrome B, with an usual spin state (d(xy))(2)(d(xz),d(yz))(3); (ii) a strong g(max) hemichrome (Type I), also showing an usual spin state (d(xy))(2)(d(xz),d(yz))(3); (iii) a hemichrome with unusual spin state (d(xz),d(yz))(4)(d(xy))(1) (Type III); (iv) and a high-spin pentacoordinate species. CD measurements suggest that the mechanism of species formation could be related with an initial process of acid denaturation. However, it is worth mentioning that based on EPR the aquomet species remains even at acidic pH, indicating that the transitions are not complete. The "pure" low-spin hemichrome presents a parallel orientation of the imidazole ring planes but the strong g(max) hemichrome is a HALS (highly anisotropic low-spin) species indicating a reciprocally perpendicular orientation of the imidazole ring planes. The hemichromes and pentacoordinate formation mechanisms are discussed in detail.

Electron Distribution in Iron Octaethyloxophlorin Complexes. Importance of the Fe(III) Oxophlorin Trianion Form in the Bis-pyridine and Bis-imidazole Complexes

The apportionment of electrons between iron and the porphyrinic macrocycle in complexes of octaethyloxophlorin (H3OEPO) has been a vexing problem. In particular, for (Py)2Fe(OEPO), which is an important intermediate in heme degradation, three resonance structures involving Fe(III), Fe(II), or Fe(I), respectively, have been considered. To clarify this matter, the electronic and geometric structures of (Py)2Fe(III)(OEPO), (Im)2Fe(III)(OEPO).2THF, and (Im)2Fe(III)(OEPO).1.6CHCl3 have been examined by single-crystal X-ray diffraction, measurement of magnetic moments as a function of temperature, and EPR and NMR spectral studies. The results clearly show that both complexes exist in the Fe(III)/oxophlorin trianion form rather than the Fe(II)/oxophlorin radical form previously established for (2,6-xylylNC)(2)Fe(II)(OEPO.). In the solid state from 10 to 300 K, (Py)2Fe(III)(OEPO) exists in the high-spin (S = 5/2) state with the axial ligands in parallel planes, a planar porphyrin, and long axial Fe-N distances. However, in solution it exists predominantly in a low-spin (S = 1/2) form. In contrast, the structures of (Im)2Fe(III)(OEPO).2THF and (Im)2Fe(III)(OEPO).1.6CHCl3 consist of porphyrins with a severe ruffled distortion, axial ligands in nearly perpendicular planes, and relatively short axial Fe-N distances. The crystallographic, magnetic, EPR, and NMR results all indicate that (Im)2Fe(III)(OEPO) exists in the low-spin Fe(III) form in both the solid state and in solution.

Models of the membrane-bound cytochromes: mössbauer spectra of crystalline low-spin ferriheme complexes having axial ligand plane dihedral angles ranging from 0 degree to 90 degrees

Crystalline samples of four low-spin Fe(III) octaalkyltetraphenylporphyrinate and two low-spin Fe(III) tetramesitylporphyrinate complexes, all of which are models of the bis-histidine-coordinated cytochromes of mitochondrial complexes II, III, and IV and chloroplast complex b(6)f, and whose molecular structures and EPR spectra have been reported previously, have been investigated in detail by Mössbauer spectroscopy. The six complexes and the dihedral angles between axial ligand planes of each are [(TMP)Fe(1-MeIm)(2)]ClO(4) (0 degree), paral-[(OMTPP)Fe(1-MeIm)(2)]Cl (19.5 degrees), paral-[(TMP)Fe(5-MeHIm)(2)]ClO(4) (26 degrees, 30 degrees for two molecules in the unit cell whose EPR spectra overlap), [(OETPP)Fe(4-Me(2)NPy)(2)]Cl (70 degrees), perp-[(OETPP)Fe(1-MeIm)(2)]Cl (73 degrees), and perp-[(OMTPP)Fe(1-MeIm)(2)]Cl (90 degrees). Of these, the first three have been shown to exhibit normal rhombic EPR spectra, each with three clearly resolved g-values, while the last three have been shown to exhibit "large g(max)" EPR spectra at 4.2 K. It is found that the hyperfine coupling constants of the complexes are consistent with those reported previously for low-spin ferriheme systems, with the largest-magnitude hyperfine coupling constant, A(zz), being considerably smaller for the "parallel" complexes (400-540 kG) than for the strictly perpendicular complex (902 kG), A(xx) being negative for all six complexes, and A(zz) and A(xx) being of similar magnitude for the "parallel" complexes (for example, for [(TMP)Fe(1-MeIm)(2)]Cl, A(zz) = 400 kG, A(xx) = -400 kG). In all cases, A(yy) is small but difficult to estimate with accuracy. With results for six structurally characterized model systems, we find for the first time qualitative correlations of g(zz), A(zz), and DeltaE(Q) with axial ligand plane dihedral angle Deltavarphi.

NMR and EPR Studies of Low-Spin Fe(III) Complexes of meso-Tetra-(2,6-Disubstituted Phenyl)Porphyrinates Complexed to Imidazoles and Pyridines of Widely Differing Basicities

A series of bis-axially ligated complexes of iron(III) tetramesitylporphyrin, TMPFe(III), tetra-(2,6-dibromophenyl)porphyrin, (2,6-Br2)4TPPFe(III), tetra-(2,6-dichlorophenyl)porphyrin, (2,6-Cl2)4TPPFe(III), tetra-(2,6-difluorophenyl)porphyrin, (2,6-F2)4TPPFe(III), and tetra-(2,6-dimethoxyphenyl)porphyrin, (2,6-(OMe)2)4TPPFe(III), where the axial ligands are 1-methylimidazole, 2-methylimidazole, and a series of nine substituted pyridines ranging in basicity from 4-(dimethylamino)pyridine (pK(a)(PyH(+)) = 9.70) to 3- and 4-cyanopyridine (pKa(PyH+) = 1.45 and 1.1, respectively), have been prepared and characterized by EPR and 1H NMR spectroscopy. The EPR spectra, recorded at 4.2 K, show "large g(max)", rhombic, or axial signals, depending on the iron porphyrinate and axial ligand, with the g(max) value decreasing as the basicity of the pyridine decreases, thus indicating a change in electron configuration from (d(xy))2(d(xz),d(yz)3 to (d(xz),d(yz))4(d(xy))1 through each series at this low temperature. Over the temperature range of the NMR investigations (183-313 K), most of the high-basicity pyridine complexes of all five iron(III) porphyrinates exhibit simple Curie temperature dependence of their pyrrole-H paramagnetic shifts and beta-pyrrole spin densities, rho(C) approximately 0.015-0.017, that are indicative of the S = 1/2 (d(xy))(2)(d(xz),d(yz))(3) electron configuration, while the temperature dependences of the pyrrole-H resonances of the lower-basicity pyridine complexes (pK(a)(PyH(+)) < 6.00) show significant deviations from simple Curie behavior which could be fit to an expanded version of the Curie law using a temperature-dependent fitting program developed in this laboratory that includes consideration of a thermally accessible excited state. In most cases, the ground state of the lower-basicity pyridine complexes is an S = 1/2 state with a mixed (d(xy))2(d(xz),d(yz))3/(d(xz),d(yz))4(d(xy))1 electron configuration, indicating that these two are so close in energy that they cannot be separated by analysis of the NMR shifts; however, for the TMPFe(III) complexes with 3- and 4-CNPy, the ground states were found to be fairly pure (d(xz),d(yz))4(d(xy))1 electron configurations. In all but one case of the intermediate- to low-basicity pyridine complexes of the five iron(III) porphyrinates, the excited state is found to be S = 3/2, with a (d(xz),d(yz))3(d(xy))1(d(z)2)1 electron configuration, lying some 120-680 cm(-1) higher in energy, depending on the particular porphyrinate and axial ligand. Full analysis of the paramagnetic shifts to allow separation of the contact and pseudocontact contributions could be achieved only for the [TMPFe(L)2]+ series of complexes.

Performance of Nonrelativistic and Quasi-Relativistic Hybrid DFT for the Prediction of Electric and Magnetic Hyperfine Parameters in 57Fe Mössbauer Spectra

57Fe electric and magnetic hyperfine parameters were calculated for a series of 10 iron model complexes, covering a wide range of oxidation and spin states. Employing the B3LYP hybrid method, results from nonrelativistic density functional theory (DFT) and quasi-relativistic DFT within the zero-order regular approximation (ZORA) were compared. Electron densities at the iron nuclei were calculated and correlated with experimental isomer shifts. It was shown that the fit parameters do not depend on a specific training set of iron complexes and are, therefore, more universal than might be expected. The nonrelativistic and quasi-relativistic electron densities gave fit parameters of similar quality; the ZORA densities are only shifted by a factor of 1.32, upward in the direction of the four-component Dirac-Fock value. From a correlation of calculated electric field gradients and experimental quadrupole splittings, the value of the 57Fe nuclear quadrupole moment was redetermined to a value of 0.16 barn, in good agreement with other studies. The ZORA approach gave no additional improvement of the calculated quadrupole splittings in comparison to the nonrelativistic approach. The comparison of the calculated and measured 57Fe isotropic hyperfine coupling constants (hfcc's) revealed that both the ZORA approach and the inclusion of spin-orbit contributions lead to better agreement between theory and experiment in comparison to the nonrelativistic results. For all iron complexes with small spin-orbit contributions (high-spin ferric and ferryl systems), a distinct underestimation of the isotropic hfcc's was found. Scaling factors of 1.81 (nonrelativistic DFT) and 1.69 (ZORA) are suggested. The calculated 57Fe isotropic hfcc's of the remaining model systems (low-spin ferric and high-spin ferrous systems) contain 10-50% second-order contributions and were found to be in reasonable agreement with the experimental results. This is assumed to be the consequence of error cancellation because g-tensor calculations for these systems are of poor quality with the existing DFT approaches. Excellent agreement between theory and experiment was found for the 57Fe anisotropic hfcc's. Finally, all of the obtained fit parameters were used for an application study of the [Fe(H2O)6]3+ ion. The calculated spectroscopic data are in good agreement with the Mossbauer and electron paramagnetic resonance results discussed in detail in a forthcoming paper.

De novo design of a D2-symmetrical protein that reproduces the diheme four-helix bundle in cytochrome bc1

An idealized, water-soluble D(2)-symmetric diheme protein is constructed based on a mathematical parametrization of the backbone coordinates of the transmembrane diheme four-helix bundle in cytochrome bc(1). Each heme is coordinated by two His residues from diagonally apposed helices. In the model, the imidazole rings of the His ligands are held in a somewhat unusual perpendicular orientation as found in cytochrome bc(1), which is maintained by a second-shell hydrogen bond to a Thr side chain on a neighboring helix. The resulting peptide is unfolded in the apo state but assembles cooperatively upon binding to heme into a well-folded tetramer. Each tetramer binds two hemes with high affinity at low micromolar concentrations. The equilibrium reduction midpoint potential varies between -76 mV and -124 mV vs SHE in the reducing and oxidizing direction, respectively. The EPR spectrum of the ferric complex indicates the presence of a low-spin species, with a g(max) value of 3.35 comparable to those obtained for hemes b of cytochrome bc(1) (3.79 and 3.44). This provides strong support for the designed perpendicular orientation of the imidazole ligands. Moreover, NMR spectra show that the protein exists in solution in a unique conformation and is amenable to structural studies. This protein may provide a useful scaffold for determining how second-shell ligands affect the redox potential of the heme cofactor.

The Effects of Axial Ligands on Electron Distribution and Spin States in Iron Complexes of Octaethyloxophlorin, Intermediates in Heme Degradation

The results presented here show that the nature of the axial ligand can alter the distribution of electrons between the metal and the porphyrin in complexes where there is an oxygen atom replacing one of the meso protons. The complexes (1-MeIm)(2)Fe(III)(OEPO) and (2,6-xylylNC)(2)Fe(II)(OEPO(*)) (where OEPO is the trianionic octaethyloxophlorin ligand and OEPO(*) is the dianionic octaethyloxophlorin radical) were prepared by addition of an excess of the appropriate axial ligand to a slurry of [Fe(III)(OEPO)](2) in chloroform under anaerobic conditions. The magnetic moment of (2,6-xylylNC)(2)Fe(II)(OEPO(*)) is temperature invariant and consistent with a simple S = (1)/(2) ground state. This complex with an EPR resonance at g = 2.004 may be considered as a model for the free-radical like EPR signal seen when the meso-hydroxylated heme/heme oxygenase complex is treated with carbon monoxide. In contrast, the magnetic moment of (1-MeIm)(2)Fe(III)(OEPO) drops with temperature and indicates a spin-state change from an S = (5)/(2) or an admixed S = (3)/(2),(5)/(2) state at high temperatures (near room temperature) to an S = (1)/(2) state at temperatures below 100 K. X-ray diffraction studies show that each complex crystallizes in centrosymmetric form with the expected six-coordinate geometry. The structure of (1-MeIm)(2)Fe(III)(OEPO) has been determined at 90, 129, and 296 K and shows a gradual and selective lengthening of the Fe-N(axial bond). This behavior is consistent with population of a higher spin state at elevated temperatures.

Axial ligand complexes of the Rhodnius nitrophorins: reduction potentials, binding constants, EPR spectra, and structures of the 4-iodopyrazole and imidazole complexes of NP4

Previously, we utilized 4-iodopyrazole (4IPzH) as a heavy atom derivative for the initial solution of the crystal structure of the nitrophorin from Rhodnius prolixus, NP1, where it was found to bind to the heme with the iodo group disordered in two positions. We have now determined the structure of the 4IPzH complex of NP4 at pH 7.5 and find that the geometry and bond lengths at the iron center are extremely similar to those of the imidazole (ImH) complex of the same protein (structure determined at pH 5.6), except that the G-H loop is not in the closed conformation. 4IPzH binds to the heme of NP4 in an ordered manner, with the iodo substituent pointed toward the opening of the heme pocket, near the surface of the protein. In order to understand the solution chemistry in terms of the relative binding abilities of 4IPzH, ImH, and histamine (Hm, a physiological ligand for the nitrophorins), we have also investigated the equilibrium binding constants and reduction potentials of these three ligand complexes of the four Rhodnius nitrophorins as a function of pH. We have found that, unlike the other Lewis bases, 4IPzH forms less stable complexes with the Fe(III) than the Fe(II) oxidation states of NP1 and NP4, and similar stability for the two oxidation states of NP2 and NP3, suggesting that this ligand is a softer base than ImH or Hm, for both of which the Fe(III) complexes are more stable than those of Fe(II) for all four nitrophorins. Surprisingly, in spite of this and the much lower basicity of 4IPzH than imidazole and histamine, the EPR g-values of all three ligand complexes are very similar.