The interstitial atom of the nitrogenase FeMo-cofactor: ENDOR and ESEEM evidence that it is not a nitrogen

X-ray crystallographic study of the nitrogenase MoFe protein revealed electron density from an atom (denoted X) inside the active-site metal cluster, the [MoFe7S9:homocitrate] FeMo-cofactor. The electron density associated with X is consistent with a single N, O, or C atom. We now have tested whether X is an N or not by comparing the Q-band ENDOR and ESEEM signals from resting-state (S = 3/2) MoFe protein and NMF-extracted FeMo-co from bacteria grown with either 14N or 15N as the exclusive N source. All of the 14N or 15N signals associated with the protein are lost upon extraction of the FeMo-co. We interpret this as strong evidence that X is not an N.     

The interstitial atom of the nitrogenase FeMo-cofactor: ENDOR and ESEEM show it is not an exchangeable nitrogen

A recent high-resolution X-ray crystallographic study (1.16 A) of the Azotobacter vinelandii nitrogenase MoFe protein revealed a previously undetected electron density associated with the active site FeMo-cofactor. The density is located inside the cluster at the center of the "trigonal prism" of six irons and is assigned to a species "X". The identity of species X was not resolved, although the electron density is consistent with a single N, O, or C atom. One proposal is that X is an N atom that derives from and exchanges with N from N2 during catalysis. In the present study, we have examined this possibility by employing 14N and 15N isotopes of N2 along with ENDOR and ESEEM spectroscopies. The WT MoFe protein and alpha-359Arg-->Lys and alpha-381Phe-->Leu variants were allowed to turn over in the presence of 14N2 or 15N2, and then were examined as resting enzymes by ENDOR and ESEEM at X- and Q-bands to look for all 14N and 15N signals coupled to the electron spin of the FeMo-cofactor and to determine if any exchanged during turnover. We have found five peaks in Q-band pulsed ENDOR spectra that appear to arise not only from previously reported N1/N2, which give rise to the ESEEM, but also from one or two additional coupled nitrogens. None of the ENDOR and ESEEM signals vanish or are altered by catalytic turnover with 15N2, and no new 15N signal is detected, leading to the conclusion that if species X is a nitrogen atom, it does not exchange during dinitrogen reduction.     

57Fe ENDOR spectroscopy and 'electron inventory' analysis of the nitrogenase E4 intermediate suggest the metal-ion core of FeMo-cofactor cycles through only one redox couple

N(2) binds to the active-site metal cluster in the nitrogenase MoFe protein, the FeMo-cofactor ([7Fe-9S-Mo-homocitrate-X]; FeMo-co) only after the MoFe protein has accumulated three or four electrons/protons (E(3) or E(4) states), with the E(4) state being optimally activated. Here we study the FeMo-co (57)Fe atoms of E(4) trapped with the α-70(Val→Ile) MoFe protein variant through use of advanced ENDOR methods: 'random-hop' Davies pulsed 35 GHz ENDOR; difference triple resonance; the recently developed Pulse-Endor-SaTuration and REcovery (PESTRE) protocol for determining hyperfine-coupling signs; and Raw-DATA (RD)-PESTRE, a PESTRE variant that gives a continuous sign readout over a selected radiofrequency range. These methods have allowed experimental determination of the signed isotropic (57)Fe hyperfine couplings for five of the seven iron sites of the reductively activated E(4) FeMo-co, and given the magnitude of the coupling for a sixth. When supplemented by the use of sum-rules developed to describe electron-spin coupling in FeS proteins, these (57)Fe measurements yield both the magnitude and signs of the isotropic couplings for the complete set of seven Fe sites of FeMo-co in E(4). In light of the previous findings that FeMo-co of E(4) binds two hydrides in the form of (Fe-(μ-H(-))-Fe) fragments, and that molybdenum has not become reduced, an 'electron inventory' analysis assigns the formal redox level of FeMo-co metal ions in E(4) to that of the resting state (M(N)), with the four accumulated electrons residing on the two Fe-bound hydrides. Comparisons with earlier (57)Fe ENDOR studies and electron inventory analyses of the bio-organometallic intermediate formed during the reduction of alkynes and the CO-inhibited forms of nitrogenase (hi-CO and lo-CO) inspire the conjecture that throughout the eight-electron reduction of N(2) plus 2H(+) to two NH(3) plus H(2), the inorganic core of FeMo-co cycles through only a single redox couple connecting two formal redox levels: those associated with the resting state, M(N), and with the one-electron reduced state, M(R). We further note that this conjecture might apply to other complex FeS enzymes.     

Analysing and understanding the active site by IR spectroscopy

IR spectroscopy is a technique particularly adapted for understanding the mechanism of catalytic reactions, being able to probe the surface mechanisms at the molecular level. In this critical review the main advances in the field are presented, both under the aspects of the in situ and operando approaches. A broad view of the most authoritative literature of the domain is given, based largely on the experience built up at the LCS laboratory in the last decades. After having presented the general methodology to observe a potential active site directly or by probe molecule adsorption, several examples illustrate the qualitative and quantitative analysis of the physical-chemical properties of the surface entities. The last part of the review is dedicated to the discrimination of the role of the active site and its links with the catalytic steps; the hot problem of the reaction intermediates and their visibility via spectroscopic techniques is critically addressed (138 references).     

Demonstration by 2H ENDOR spectroscopy that myo-inositol binds via an alkoxide bridge to the mixed-valent diiron center of myo-inositol oxygenase

myo-Inositol oxygenase (MIOX) is a non-heme diiron oxygenase that cleaves cyclohexane-(1,2,3,5/4,6-hexa)-ol (myo-inositol, MI) to d-glucuronate. Here, we use 2H ENDOR spectroscopy to demonstrate that MI binds to the diiron(II/III) cofactor of MIOX via an alkoxide bridge, most likely involving O1. Analysis shows that MI adopts a symmetrical geometry in which the O-C-2H plane of the bridge is approximately orthogonal to the Fe-O-Fe plane.     

Proton positions in the Mn(2+) binding site of concanavalin A as determined by single-crystal high-field ENDOR spectroscopy

High-field (95 GHz) pulsed EPR and electron-nuclear double resonance (ENDOR) techniques have been used for the first time to determine coordinates of ligand protons of a high-spin metal center in a protein single crystal. The protein concanavalin A contains a Mn(2+) ion which is coordinated to two water molecules, a histidine residue, and three carboxylates. Single crystals of concanavalin A were grown in H(2)O and in D(2)O to distinguish the exchangeable water protons from the nonexchangeable protons of the imidazole group. Distinct EPR transitions were selected by performing the ENDOR measurements at different magnetic fields within the EPR spectrum. This selection, combined with the large thermal polarization achieved at 4.5 K and a magnetic field of approximately 3.4 T allowed us to assign the ENDOR signals to their respective M(S) manifolds, thus providing the signs of the hyperfine couplings. Rotation patterns were acquired in the ac and ab crystallographic planes. Two distinct crystallographic sites were identified in each plane, and the hyperfine tensors of two of the imidazole protons and the four water protons were determined by simulations of the rotation patterns. All protons have axially symmetric hyperfine tensors and, by applying the point-dipole approximation, the positions of the various protons relative to the Mn(2+) ion were determined. Likewise, the water protons involved in H-bonding to neighboring residues were identified using the published, ultrahigh-resolution X-ray crystallographic coordinates of the protein (Deacon et al. J. Chem. Soc., Faraday Trans. 1997, 93(24), 4305-4312).     

ENDOR/HYSCORE studies of the common intermediate trapped during nitrogenase reduction of N2H2, CH3N2H, and N2H4 support an alternating reaction pathway for N2 reduction

Enzymatic N(2) reduction proceeds along a reaction pathway composed of a sequence of intermediate states generated as a dinitrogen bound to the active-site iron-molybdenum cofactor (FeMo-co) of the nitrogenase MoFe protein undergoes six steps of hydrogenation (e(-)/H(+) delivery). There are two competing proposals for the reaction pathway, and they invoke different intermediates. In the 'Distal' (D) pathway, a single N of N(2) is hydrogenated in three steps until the first NH(3) is liberated, and then the remaining nitrido-N is hydrogenated three more times to yield the second NH(3). In the 'Alternating' (A) pathway, the two N's instead are hydrogenated alternately, with a hydrazine-bound intermediate formed after four steps of hydrogenation and the first NH(3) liberated only during the fifth step. A recent combination of X/Q-band EPR and (15)N, (1,2)H ENDOR measurements suggested that states trapped during turnover of the α-70(Ala)/α-195(Gln) MoFe protein with diazene or hydrazine as substrate correspond to a common intermediate (here denoted I) in which FeMo-co binds a substrate-derived [N(x)H(y)] moiety, and measurements reported here show that turnover with methyldiazene generates the same intermediate. In the present report we describe X/Q-band EPR and (14/15)N, (1,2)H ENDOR/HYSCORE/ESEEM measurements that characterize the N-atom(s) and proton(s) associated with this moiety. The experiments establish that turnover with N(2)H(2), CH(3)N(2)H, and N(2)H(4) in fact generates a common intermediate, I, and show that the N-N bond of substrate has been cleaved in I. Analysis of this finding leads us to conclude that nitrogenase reduces N(2)H(2), CH(3)N(2)H, and N(2)H(4) via a common A reaction pathway, and that the same is true for N(2) itself, with Fe ion(s) providing the site of reaction.     

Unambiguous determination of the ionization state of a glycoside hydrolase active site lysine by 1H-15N heteronuclear correlation spectroscopy

We have investigated the lysine side chain amines in the 34 kDa catalytic domain from Cellulomonas fimi beta-(1,4)-glycosidase Cex (or CfXyn10A) using 1H-detected 15N heteronuclear correlation NMR spectroscopy. Signals from the 1Hzeta ( approximately 8 ppm) and 15Nzeta ( approximately 35 ppm) of Lys302 in the unmodified enzyme and Lys47 in a trapped cellobiosyl-enzyme intermediate were detected in a 1H-15N HMQC spectrum (pH 6.5 and 30 degrees C). The amine of Lys302 forms a buried ion pair, and that of Lys47 is hydrogen bonded to the cellobioside. Both lysines are positively charged, as unambiguously demonstrated by the splitting of their 15Nzeta signals into quartets (|1JNH| approximately 75 Hz) in a 1H-15N HSQC spectrum recorded without 1H decoupling during 15N evolution. Qualitative insights into the dynamic properties of these lysines are also provided by the deviations of their quartet intensity ratios from that of approximately 3:1:1:3 expected for a highly mobile amine. On the basis of the observed ratios of approximately 1:1:1:1 for the quartet of Lys302 and approximately 0.5:1:1:0.5 for Lys47, the amine of the latter active site residue is most rigidly positioned. Signals from at least 8 and 10 additional positively charged, mobile amines in Cex were observed at 10 degrees C and pH 6.5 and 5.6, respectively. By using conditions of reduced temperature, slightly acidic pH, and low general base concentrations, as well as water flipback pulses to minimize the effects of hydrogen exchange, 1H-15N correlation experiments provide a sensitive route to directly investigate the charge states and dynamic properties of the N-terminal and side chain amines in proteins and protein complexes.     

Carbon monoxide as an intrinsic ligand to iron in the active site of the iron-sulfur-cluster-free hydrogenase H2-forming methylenetetrahydromethanopterin dehydrogenase as revealed by infrared spectroscopy

The iron-sulfur-cluster-free hydrogenase Hmd (H(2)-forming methylenetetrahydromethanopterin dehydrogenase) from methanogenic archaea has recently been found to contain one iron associated tightly with an extractable cofactor of yet unknown structure. We report here that Hmd contains intrinsic CO bound to the Fe. Chemical analysis of Hmd revealed the presence of 2.4 +/- 0.2 mol of CO/mol of iron. Fourier transform infrared spectra of the native enzyme showed two bands of almost equal intensity at 2011 and 1944 cm(-)(1), interpreted as the stretching frequencies of two CO molecules bound to the same iron in an angle of 90 degrees . We also report on the effect of extrinsic (12)CO, (13)CO, (12)CN(-), and (13)CN(-) on the IR spectrum of Hmd.     

Ligand-bound S = 1/2 FeMo-cofactor of nitrogenase: hyperfine interaction analysis and implication for the central ligand X identity

Broken symmetry density functional theory (BS-DFT) has been used to address the hyperfine parameters of the single atom ligand X, proposed to be coordinated by six iron ions in the center of the paramagnetic FeMo-cofactor (FeMoco) of nitrogenase. Using the X = N alternative, we recently found that any hyperfine signal from X would be small (calculated A(iso)(X = (14)N) = 0.3 MHz) due to both structural and electronic symmetry properties of the [Mo-7Fe-9S- X] FeMoco core in its resting S = 3/2 state. Here, we extend our BS-DFT approach to the 2e(-) reduced S = 1/2 FeMoco state. Alternative substrates coordinated to this FeMoco state effectively perturb the electronic and/or structural symmetry properties of the cofactor. Using an example of an allyl alcohol (H2C=CH-CH2-OH) product ligand, we consider three different binding modes at single iron site and three different BS-DFT spin state structures and show that this binding would enhance the key hyperfine signal A(iso)(X) by at least 1 order of magnitude (3.8 < or = A(iso)(X = (14)N) < or = 14.7 MHz), and this result should not depend strongly on the exact identity of X (nitrogen, carbon, or oxygen). The interstitial atom, when the nucleus has a nonzero magnetic moment, should therefore be observable by ESR methods for some ligand-bound FeMoco states. In addition, our results illustrate structural details and likely spin-coupling patterns for models for early intermediates in the catalytic cycle.     

Studying how different terminal groups change the motion of H2NSO2C6H4CONH(EG)3R when bound to the active site of human carbonic anhydrase II

Molecular dynamics simulations have been used to explore the motions of series of ligands containing coupled benzenesulfonamide and oligoethylene glycol moieties (H2NSO2C6H4CONH(CH2-CH2OCH2CH2OCH2CH2)R+; R+ = NH3+, NHCOCH2NH3+, NHCOCH(CH2Ph)NH3+) bound at the active site of human carbonic anhydrase II (HCAII; E.C. 4.2.1.1). These complexes have been examined previously by X-ray crystallography; the locations of the terminal groups of these ligands were not defined in the crystal structures. These stimulations, carried out in the presence of water, provide dynamic information about the motion of the bound ligand that supplements the quasistatic information from crystallography. Our results suggested that the Gly and Phe groups of these ligands interacted weakly with the protein adjacent to the active site. Quantitative estimates of energies of binding did not correlate usefully with observed free energies of binding, but in the absence of information about entropies, it is not possible to tell if the lack of correlation between calculated energies and observed free energies represents inaccuracies in the energies, or a compensation between enthalpies and entropies. When the terminal Phe group was placed near a previously identified hydrophobic patch in the active site (Phe20 and Pro202) the average conformation of the ligand inferred from this simulation was inconsistent with that from the crystal structure; this result illustrates the problems of misleading local minima in these types of simulations.     

Insights into the structure of the active site of the O2-tolerant membrane bound [NiFe] hydrogenase of R. eutropha H16 by molecular modelling

Structural models for the Ni-B state of the wild-type and C81S protein variant of the membrane-bound [NiFe] hydrogenase from Ralstonia eutropha H16 were derived by applying the homology model technique combined with molecular simulations and a hybrid quantum mechanical/molecular mechanical approach. The active site structure was assessed by comparing calculated and experimental IR spectra, confirming the view that the active site structure is very similar to those of anaerobic standard hydrogenases. In addition, the data suggest the presence of a water molecule in the second coordination sphere of the active centre.     

Catalytic behavior of a polynuclear Mg-Mo complex and nitrogenase active site (FeMoco) isolated from the enzyme in reactions with C2H2, N2, and CO: A comparative study

In order to identify common and distinctive features in the catalytic behavior of natural and artificial nitrogen-fixation clusters, the kinetics of the catalytic reduction of C₂H₂ in the presence of Mg-Mo-cluster (1) was investigated and compared with the kinetics of acetylene reduction catalyzed by the cluster FeMoco (2) isolated from the enzyme nitrogenase we studied previously. The reactions were conducted in the presence of Zn/Hg and Eu/Hg as reducing agents and PhSH and C₆F₅SH as proton donors, i.e., under the same conditions as had been used in the case of 2. Both polynuclear Mg-Mo-complex and the europium amalgam-reduced FeMoco have multiple interdependent binding sites for substrates and/or inhibitors. Carbon monoxide inhibits the acetylene reduction much less efficiently in systems with cluster 1 than in systems with cluster 2, although the type of inhibition is mixed in both systems: CO binds to multiple sites of the cluster and affects both C₂H₂ complexation to the reduced cluster and decomposition of the catalyst-substrate complex to give the products. Unlike isolated FeMoco, the Mg-Mo-cluster efficiently catalyzes the reduction of molecular nitrogen. The reaction is greatly inhibited by acetylene, while no inhibiting effect of N₂ is observed in acetylene reduction, as was found earlier for a system with the natural cluster as the catalyst. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 766–774, May, 2006.     

Variational quantum mechanical and active database approaches to the rotational-vibrational spectroscopy of ketene, H2CCO

A variational quantum mechanical protocol is presented for the computation of rovibrational energy levels of semirigid molecules using discrete variable representation of the Eckart-Watson Hamiltonian, a complete, "exact" inclusion of the potential energy surface, and selection of a vibrational subspace. Molecular symmetry is exploited via a symmetry-adapted Lanczos algorithm. Besides symmetry labels, zeroth-order rigid-rotor and harmonic-oscillator quantum numbers are employed to characterize the computed rovibrational states. Using the computational molecular spectroscopy algorithm presented, a large number of rovibrational states, up to J = 50, of the ground electronic state of the parent isotopologue of ketene, H(2) (12)C=(12)C=(16)O, were computed and characterized. Based on 12 references, altogether 3982 measured and assigned rovibrational transitions of H(2) (12)C=(12)C=(16)O have been collected, from which 3194 were validated. These transitions form two spectroscopic networks (SN). The ortho and the para SNs contain 2489 and 705 validated transitions and 1251 and 471 validated energy levels, respectively. The computed energy levels are compared with energy levels obtained, up to J = 41, via an inversion protocol based on this collection of validated measured rovibrational transitions. The accurate inverted energy levels allow new assignments to be proposed. Some regularities and irregularities in the rovibrational spectrum of ketene are elucidated.     

Insights into the mechanism of electrocatalytic hydrogen evolution mediated by Fe2(S2C3H6)(CO)6: the simplest functional model of the Fe-hydrogenase active site

The di-iron complex Fe2(S2C3H6)(CO)6 (a), one of the simplest functional models of the Fe-hydrogenases active site, is able to electrocatalyze proton reduction. In the present study, the H2 evolving path catalyzed by a has been characterized using density functional theory. It is showed that, in the early stages of the catalytic cycle, a neutral mu-H adduct is formed; monoelectron reduction and subsequent protonation can give rise to a diprotonated neutral species (a-muH-SH), which is characterized by a mu-H group, a protonated sulfur atom, and a CO group bridging the two iron centers, in agreement with experimental IR data indicating the formation of a long-lived mu7-CO species. H2 release from a-muH-SH, and its less stable isomer a-H2 is kinetically unfavorable, while the corresponding monoanionic compounds (a-muH-SH- and a-H2-) are more reactive in terms of dihydrogen evolution, in agreement with experimental data. The key species involved in electrocatalysis have structural features different from the hypothetical intermediates recently proposed to be involved in the enzymatic process, an observation that is possibly correlated with the reduced catalytic efficiency of the biomimetic di-iron assembly.     

Structures of E. coli peptide deformylase bound to formate: insight into the preference for Fe2+ over Zn2+ as the active site metal

E. coli peptide deformylase (PDF) catalyzes the deformylation of nascent polypeptides generated during protein synthesis. While PDF was originally thought to be a zinc enzyme, subsequent studies revealed that the active site metal is iron. In an attempt to understand this unusual metal preference, high-resolution structures of Fe-, Co-, and Zn-PDF were determined in complex with its deformylation product, formate. In all three structures, the formate ion binds the metal and forms hydrogen-bonding interactions with the backbone nitrogen of Leu91, the amide side chain of Gln50, and the carboxylate side chain of Glu133. One key difference, however, is how the formate binds the metal. In Fe-PDF and Co-PDF, formate binds in a bidentate fashion, while in Zn-PDF, it binds in a monodentate fashion. Importantly, these structural results provide the first clues into the origins of PDF's metal-dependent activity differences. On the basis of these structures, we propose that the basis for the higher activity of Fe-PDF stems from the better ability of iron to bind and activate the tetrahedral transition state required for cleavage of the N-terminal formyl group.     

Ring current effects in the active site of medium-chain Acyl-CoA dehydrogenase revealed by NMR spectroscopy

Medium-chain acyl-CoA dehydrogenase (MCAD) catalyzes the flavin-dependent oxidation of fatty acyl-CoAs to the corresponding trans-2-enoyl-CoAs. The interaction of hexadienoyl-CoA (HD-CoA), a product analogue, with recombinant pig MCAD (pMCAD) has been studied using (13)C NMR and (1)H-(13)C HSQC spectroscopy. Upon binding to oxidized pMCAD, the chemical shifts of the C1, C2, and C3 HD carbons are shifted upfield by 12.8, 2.1, and 13.8 ppm, respectively. In addition, the (1)H chemical shift of the C3-H is also shifted upfield by 1.31 ppm while the chemical shift of the C4 HD-CoA carbon is unchanged upon binding. These changes in chemical shift are unexpected given the results of previous Raman studies which revealed that the C3=C2-C1=O HD enone fragment is polarized upon binding to MCAD such that the electron density at the C3 and C1 carbons is reduced, not increased (Pellet et al. Biochemistry 2000, 39, 13982-13992). To investigate the apparent discrepancy between the NMR and Raman data for HD-CoA bound to MCAD, (13)C NMR spectra have been obtained for HD-CoA bound to enoyl-CoA hydratase, an enzyme system that has also previously been studied using Raman spectroscopy. Significantly, binding to enoyl-CoA hydratase causes the chemical shifts of the C1 and C3 HD carbons to move downfield by 4.8 and 5.6 ppm, respectively, while the C2 resonance moves upfield by 2.2 ppm, in close agreement with the alterations in electron density at these carbons predicted from Raman spectroscopy (Bell, A. F.; Wu, J.; Feng, Y.; Tonge, P. J. Biochemistry 2001, 40, 1725-33). The large increase in shielding experienced by the C1 and C3 HD carbons in the HD-CoA/MCAD complex is proposed to arise from the ring current field from the isoalloxazine portion of the flavin cofactor. The flavin ring current, which is only present when the enzyme is placed in an external magnetic field, also explains the differences in (13)C NMR chemical shifts for acetoacetyl-CoA when bound as an enolate to MCAD and enoyl-CoA hydratase and is used to rationalize the observation that the line widths of the C1 and C3 resonances are narrower when the ligands are bound to MCAD than when they are free in the protein solution.     

Quantum chemical calculations of the NHA bound nitric oxide synthase active site: O2 binding and implications for the catalytic mechanism

Density functional theory has been employed to model the binding of the intermediate substrate NHA, by nitric oxide synthases. In particular, the orientation and interactions of possibly catalytically important substrate hydrogens, with and without molecular oxygen bound to the active site heme group, are considered. Without O(2), three possible conformers have been found, with the energetically most favored structure being that in which both protons of the -NHOH moiety of NHA are directed toward the heme group. With oxygen bound, four different structures were found. The energetically lowest structure is again found to have both hydrogens of the -NHOH group pointing toward the heme group, thus forming hydrogen bonds between -NH- and the terminal oxygen, and between -OH and the inner oxygen of the heme-O(2) group. In addition, unprotonated structures of the substrate bound to the active site are considered and the proton affinity calculated.     

On the electrochemistry of diiron dithiolate complexes related to the active site of the [FeFe]H2ase

A few recent electrochemical studies of diiron models of the iron-only hydrogenases' active site are summarized. Emphasis is put on the reduction mechanisms of hexacarbonyl complexes and on the different mechanisms of proton reduction that may operate depending on the nature of the complex and the strength of the acid. An attempt is made to discuss the thermodynamic and kinetic limitations of proton reduction processes supported by these compounds.