The Sixteenth Iron in the Nitrogenase MoFe Protein

Abstract

Another iron in the fire: X-ray anomalous diffraction studies on the nitrogenase MoFe protein show the presence of a mononuclear iron site, designated as Fe16, which was previously identified as either Ca2+ or Mg2+. The position of the absorption edge indicates that this site is in the oxidation state +2. The high sequence conservation of the residues coordinated to Fe16 emphasizes the potential importance of the site in nitrogenase. The nitrogenase molybdenum–iron (MoFe) protein contains two types of metalloclusters, the [8Fe:7S] P-cluster and the [7Fe:9S:Mo:C:homocitrate] FeMo cofactor, which have been the focus of considerable interest to define their roles in the mechanism of biological nitrogen fixation.1 Along with these two metalloclusters, a third metal binding site is present in the MoFe protein that has attracted much less attention. The metal at this site has been assigned to either Ca2+ or Mg2+ in different MoFe protein structures and is positioned at the interface between the two β-subunits: approximately 25 Å and 21 Å away from the P-cluster and FeMo-cofactor, respectively.2 Six ligands coordinate the metal in approximate octahedral geometry: the side-chain oxygen atoms from β-Glu109, β′-Asp353, and β′-Asp357, the backbone carbonyl oxygen of β-Arg108, and two water molecules (β and β′ denote the two distinct β-subunits in an α2β2 MoFe protein tetramer). The assignments as either Ca2+ or Mg2+ were inferred from the value of the electron density at this position, the coordination environment at the metal binding site and/or the presence of Mg2+ in the crystallization solution.2a,2b,2d As the identified species are not redox-active metals, it was speculated that this metal site most likely serves a structural rather than a mechanistic role. However, as shown below, our multiple-wavelength anomalous diffraction (MAD) data around the Fe K-edge unambiguously demonstrate the presence of Fe at this binding site, which we will subsequently refer to as the mononuclear metal binding site (MMB site) and the metal at this site as Fe16. MoFe protein crystals were grown under anaerobic conditions as described previously1b using wild-type protein isolated from Azotobacter vinelandii. Two sets of anomalous diffraction data were collected at 7080 eV and 7130 eV (below and above the Fe K-edge, respectively; Supporting Information, Table S1) to 2.1 Å resolution. As shown in Figure 1, well-defined density is clearly present at the MMB site in the double difference (Δanom7130eV−Δanom7080eV) anomalous density map, and the peak density is about half of that observed at the Fe sites in the FeMo cofactor and P-cluster. In contrast, the density around non-Fe heavy atoms (for example S and Mo) is comparable to the noise level. The double difference anomalous density maps clearly indicate the presence of Fe at this site, and the comparison of the anomalous density to the Fe sites in the two metalloclusters suggests a partial occupancy. The double difference (Δanom7130eV−Δanom7080eV) anomalous density maps (highlighted in cyan) at the MMB site (A), the FeMo cofactor (B), and the P-cluster (C). Fe orange, Mo cyan, S yellow, C gray, N blue, O red. These maps are contoured at the level of 10 σ. To determine the oxidation state and chemical environment of Fe16, site-specific Fe K-edge absorption spectra were extracted from anomalous differences present in a 10° wedge of diffraction data (plus inverse beam) measured from a MoFe protein crystal at 19 different energies across the Fe K-edge (Supporting Information, Table S2). The Δf′′ for heavy atoms (Fe, Mo, and S) in the MoFe protein structure were then refined as described previously.3 As shown in the top panel of Figure 2, the well-defined absorption edge feature of Fe16 confirms unambiguously the presence of an iron at this site, which is consistent with the anomalous density analysis shown above. Significantly, the absorption edge of the Δf′′ spectrum for Fe16 is about 5 eV higher relative to the averaged Δf′′ spectra for the iron atoms in the two metalloclusters. Because the cluster Fe are mainly four-coordinate ferrous irons, the 5 eV shift suggests either a different coordination environment and/or a different oxidation state. Indeed, the edge energy and the fine structure around the absorption edge of Fe16 are similar to that of a six-coordinate ferrous model compound (ferrous sulfate; bottom panel of Figure 2), and differ markedly from the corresponding ferric model compound (ferric sulfate; bottom panel of Figure 2). The average FeO distance for Fe16 is 2.23 Å, which matches the average FeO distances of six-coordinate oxygen-only ferrous compounds (2.16±0.08 Å for 61 ferrous FeO6 compounds and 2.00±0.02 Å for 14 ferric FeO6 compounds with R factor ≤0.05) in the Cambridge Crystallographic Structure Database.4 Combining these observations, we conclude that Fe16 has properties consistent with a six-coordinate ferrous iron. The refined Δf′′ spectra for Fe sites in the MoFe protein in comparison to that of Fe models. Top panel: refined Δf′′ for Fe16 (▵,—) and the averaged Δf′′ for all Fe in the P-cluster (◊,••••) and FeMo cofactor (○,- - - -), respectively. Bottom panel: Δf′′ converted from the X-ray absorption spectra of ferrous sulfate heptahydrate (FeSO4⋅7 H2O, —) and ferric sulfate hydrate (Fe2(SO4)3⋅xH2O, ••••), respectively. No background removal or normalization were applied to the refined Δf′′ spectra for Fe16 and other Fe in the MoFe protein. The absorption edge of Δf′′ curve for Fe16 is about 5 eV higher relative to the metalloclusters; the shift towards higher energy for the FeMo cofactor relative to the P-cluster is consistent with a more reduced state for the Fe in the latter. As a further check, we examined the recent 1 Å resolution structure of the MoFe protein (PDB code: 3U7Q)1f for the feasibility of fitting a partially occupied Fe at the MMB site. Indeed, the refinement results show the electron density at the MMB sites is approximately equivalent to either 0.5 Fe, 0.7 Ca or 1 Mg/Na per site. The shape of the electron density in the MMB-sites, however, indicates it may be described as a superposition of two distinct arrangements (Figure 3 A). Currently, the non-spherical electron density around Fe16 is modeled with two mutually exclusive arrangements that contain either the Fe16 (Figure 3 B) or Na1 (Figure 3 C), offset by 1.2 Å. The electron density also indicates that the side chain of Glu109 has two distinct conformations. Surrounding these peaks are two half-occupied water molecules (designated HOH-2A and HOH-2B) that could be coordinated to Fe16 and Na1, respectively (Figure 3). The average NaO distance of 2.37 Å in the pentahedral coordination environment around Na1 is consistent with that of Na+, but not Mg2+ or H2O binding sites reported in the Cambridge Crystallographic Structure Database and the Protein Data Bank (Supporting Information, Table S3).5 The MO distances around Na1 are also compatible with Ca2+, but the electron density does not account for a half-occupied Ca2+. The refined crystallographic structure at the MMB site in the Av MoFe protein superimposed on the electron density map contoured at 2 σ (A). The observed density of the MMB-site is interpreted as a superposition of partially occupied Fe16-bound (B) and Fe16-free (C) forms, with the latter modeled as a half-occupied Na+ (Na1). Fe orange, O red, N blue, C gray, Na blue spheres. The presence of any transition metals other than Fe and Mo in the MoFe protein has been ruled out by hard X-ray emission spectroscopy (XES) studies and inductively coupled plasma mass spectrometry (ICP-MS) analysis of the metal composition of the MoFe protein. As shown in the Supporting Information, Figure S1, Fe, Zn, Cu and V are the only transition metals with atomic number Z≤30 observed in the X-ray emission spectra of the MoFe protein crystal and the crystallization solution. The Cu and VKα peak intensities in the spectrum of the MoFe protein crystal are present in both the crystallization solution and metal-free water (data not shown), suggesting that they are primarily contributed by the emission from the beamline components. From the ICP-MS analysis, the Fe/Mo ratio averages 15.6±0.3 (Supporting Information, Table S4), which is an excess of 0.6 over the iron present in the FeMo cofactor and P-cluster, but this excess is within the experimental error of ICP-MS analysis on Fe. In addition to Fe and Mo, Mg and Ca are present at 1.7–2.6 and 0.8–1.5 atoms per MoFe-protein tetramer, respectively. The only other metal detected at a level of about 0.1 atom per molecule by ICP-MS analysis is Zn. To quantify the Fe16 occupancy from the MAD data, we compared the edge jump for Fe16 to the average of that observed for all the Fe in the P-cluster and FeMo cofactor, with the assumption that the latter corresponds to unit occupancy. The edge jump here refers to the difference between the Δf′′ values refined above and below the Fe K-edge (at 7130 and 7080 eV, respectively). Eight sets of the MAD data were collected for this study, which were acquired from different MoFe protein crystals and with resolutions ranging from 2.1 Å to 1.39 Å (Supporting Information, Table S5). When all 8 data sets are analyzed, the Fe16 occupancies average 0.5, ranging from 0.1 to 1.1. The ratios of the electron density at the Fe16 site to the average of all the Fe in the metalloclusters average 0.7, with a range of 0.4 to 1.1; the higher value relative to the anomalous density-derived quantity reflects the contribution of all scatterers at this site to the electron density, while only iron contributes to the absorption edge jump. The observed variation could reflect changes in the iron content during protein purification or crystallization; as Fe16 is in an approximate octahedral geometry with O-only ligands, the observed large variation in the edge jumps cannot be explained by the polarization effect.6 Substoichiometric occupancies of metal binding sites are not unprecedented and have been noted in several metalloenzymes, such as the copper-containing particulate methane monooxygenase,7 the tungsten–iron–sulfur enzyme acetylene hydratase,8 acetyl-CoA synthase/carbon monoxide dehydrogenase,9 and in the ferrous form of superoxide reductase,10 as well as for copper-bound metallochaperones.11 A critical question concerns the possible roles of Fe16 at the MMB site. Protein sequence comparisons of nitrogenase homologues indicate that the residues comprising the MMB site may be important (Supporting Information, Table S6); the identities of the three β-subunit residues (Glu109, Asp353, and Asp357) involved in Fe16 binding are highly conserved in almost all the nitrogenase homologues within and across species (Supporting Information, Figure S2). Arg108 is the only Fe16-binding residue not strictly conserved, but as it coordinates through its backbone carbonyl oxygen, it is therefore not restricted to any specific side chain. Intriguingly, one of the Fe16-bound ligands, Asp357, is conserved in the nitrogenase-like, light-independent protochlorophyllide reductase BchNB.12 This residue in BchNB is close to the substrate binding site and has been shown to be involved in proton transfer during the reduction of protochlorophyllide (Supporting Information, Figure S2F).13 The possibility of inferring potential functional roles for Fe16 from the functions of comparable sites in other proteins is precluded by the apparent absence of mononuclear iron sites with similar coordination environments in the protein data bank, as searched with version 10 of the MESPEUS database.14 Indeed, the best match for Fe16 is one of the two iron sites present in diiron centers of rubrerythrins (Supporting Information, Figure S3), although the correspondence is incomplete and there is no evidence for a second iron site in our electron density maps. While the conservation of the residues at the MMB site emphasizes the potential importance of the site in nitrogenase, the precise function has yet to be defined. The identification of Fe16 in the MMB site, particularly at partial occupancy, is a remarkable demonstration of the analytical capability of MAD to site-specifically identify elements, even in complex systems such as nitrogenase that have been extensively studied for over 40 years. In the case of the MoFe protein, Fe16 at full occupancy would account for about 6 % of the Fe, a value not reliably measurable by conventional metal analysis while its presence has potential implications for the interpretation of other biophysical studies. A review of past literature suggests that Fe16 may have been “ghosting” around in spectroscopic investigations of nitrogenase. A minor ferrous “impurity”, varied from 3 % to 8 % of total Fe in nitrogenase, has been reported in Mössbauer studies.15 The characteristics of this ferrous “impurity” are typical for a six-coordinate, high-spin ferrous species with an oxygen and/or nitrogen ligand sphere,15b matching the chemical environment around Fe16. Because a high-spin ferrous iron cannot be easily characterized by conventional electronic paramagnetic resonance spectroscopy studies typically conducted on nitrogenase,16 the presence of this iron species may have been overlooked in past studies. It is important to keep the possible presence of a 16th Fe in nitrogenase in mind when interpreting results from structural and crystallographic and spectroscopic studies on nitrogenase. As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re-organized for online delivery, but are not copy-edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.