Controlled protonation of iron-molybdenum cofactor by nitrogenase: a structural and theoretical analysis

Nitrogen Binding to the FeMo-Cofactor of Nitrogenase

Density functional calculations are presented to unravel the first steps of nitrogen fixation of nitrogenase. The individual steps leading from the resting state to nitrogen binding at the FeMo-cofactor with a central nitrogen ligand are characterized. The calculations indicate that the Fe-Mo cage opens as dinitrogen binds to the cluster. In the resting state, the central cage is overall neutral. Electrons and protons are transferred in an alternating manner. Upon dinitrogen binding, one protonated sulfur bridge is broken. An axial and a bridged binding mode of dinitrogen have been identified. Adsorption at the Mo site has been investigated but appears to be less favorable than binding at Fe sites.

A computer simulation study of the hydrated proton in a synthetic proton channel

Classical molecular dynamics simulations using the multistate empirical valence bond model for aqueous proton transport were performed to characterize the hydration structure of an excess proton inside a leucine-serine synthetic ion channel, LS2. For such a nonuniform pore size ion channel, it is found that the Zundel ion (H(5)O(2)(+)) solvation structure is generally more stable in narrow channel regions than in wider channel regions, which is in agreement with a recent study on idealized hydrophobic proton channels. However, considerable diversity in the relative stability of the Zundel to Eigen cation (H(9)O(4)(+)) was observed. Three of the five wide channel regions, one located at the channel's center and the other two located near the channel mouths, are found to show extraordinary preference for the Eigen solvation structure. This implies that proton hopping is inhibited in these regions and therefore suggests that these regions may behave as barriers in the proton conducting pathway inside the channel. The proton solvation is also greatly influenced by the local molecular environment of the protein. In particular, the polar side chains of the Ser residues, which are intimately involved in the solvation structure, can greatly influence proton solvation. However, no preference of the influence by the various Ser side chains was found; they can either promote or prevent the formation of certain solvation structures.

Structural, spectroscopic, and redox consequences of a central ligand in the FeMoco of nitrogenase: a density functional theoretical study

Broken symmetry density functional and electrostatics calculations have been used to shed light on which of three proposed atoms, C, N, or O, is most likely to be present in the center of the FeMoco, the active site of nitrogenase. At the Mo(4+)4Fe(2+)3Fe(3+) oxidation level, a central N(3-) anion results in (1) calculated Fe-N bond distances that are in very good agreement with the recent high-resolution X-ray data of Einsle et al.; (2) a calculated redox potential of 0.19 eV versus the standard hydrogen electrode (SHE) for FeMoco(oxidized) + e(-) --> FeMoco(resting), in good agreement with the measured value of -0.042 V in Azotobacter vinelandii; and (3) average Mössbauer isomer shift values (IS(av) = 0.48 mm s(-1)) compatible with experiment (IS(av) = 0.40 mm s(-1)). At the more reduced Mo(4+)6Fe(2+)1Fe(3+) level, the calculated geometry around a central N(3-) anion still correlates well with the X-ray data, but the average Mössbauer isomer shift value (IS(av) = 0.54 mm s(-1)) and the redox potential of -2.21 eV show a much poorer agreement with experiment. These calculated structural, spectroscopic, and redox data indicate the most likely iron oxidation state for the resting FeMoco of nitrogenase to be 4Fe(2+)3Fe(3+). At this favored oxidation state, oxygen or carbon coordination leads to (1) Fe-O distances in poor agreement and Fe-C distances in good agreement with experiment and (2) calculated redox potentials of +0.97 eV for O(2-) and -1.31 eV for C(4-). The calculated structural parameters and/or redox data suggest either O(2-) or C(4-) is unlikely as a central anion.

An atomic-level mechanism for molybdenum nitrogenase. Part 1. Reduction of dinitrogen

The properties of the Fe and Mo sites of the iron-molybdenum cofactor of nitrogenase with respect to binding and activation of N(2) have been studied by molecular mechanics calculations on the local protein environment and by density functional theory (DFT) calculations on subsections of the cofactor. The DFT calculations indicate that the homocitrate ligand of the cofactor can become monodentate on reduction, allowing N(2) to bind at Mo. In addition, the neighboring Fe atom plays a crucial role in N(2) reduction by stabilizing the initial reduced N(2) species and by facilitating cleavage of the N-N bond. The various possible isomers for partially reduced N(2) intermediates have been compared by DFT, and a detailed model for the reduction of N(2) is developed based on these results, together with chemical precedents and the available biochemical data for nitrogenase.

Binding modes for the first coupled electron and proton addition to FeMoco of nitrogenase

A combined broken-symmetry density functional and electrostatics approach has been used to model the one-electron reduced and protonated state of the iron-molybdenum cofactor active site of nitrogenase. The active site of the protein contains Fe, Mo, S, N, and O atoms, and many possible sites for protonation have been examined. A novel hydridic proton asymmetrically located in the central cavity created by six Fe sites is most favored from the calculations. Under physiological turnover conditions of low electron flux, the formation of this iron-hydride intermediate may represent a first step towards cofactor liberation of dihydrogen in the absence of dinitrogen.

Interaction of acetylene and cyanide with the resting state of nitrogenase alpha-96-substituted MoFe proteins

The nitrogenase MoFe protein contains the active site metallocluster called FeMo-cofactor [7Fe-9S-Mo-homocitrate] that exhibits an S = 3/2 EPR signal in the resting state. No interaction with FeMo-cofactor is detected when either substrates or inhibitors are incubated with MoFe protein in the resting state. Rather, the detection of such interactions requires the incubation of the MoFe protein together with its obligate electron donor, called the Fe protein, and MgATP under turnover conditions. This indicates that a more reduced state of the MoFe protein is required to accommodate substrate or inhibitor interaction. In the present work, substitution of an arginine residue (alpha-96(Arg)) located next to the active site FeMo-cofactor in the MoFe protein by leucine, glutamine, alanine, or histidine is found to result in MoFe proteins that can interact with acetylene or cyanide in the as-isolated, resting state without the need for the Fe protein, or MgATP. The dithionite-reduced, resting states of the alpha-96(Leu)-, alpha-96(Gln)-, alpha-96(Ala)-, or alpha-96(His)-substituted MoFe proteins show an S = 3/2 EPR signal (g = 4.26, 3.67, 2.00) similar to that assigned to FeMo-cofactor in the wild-type MoFe protein. However, in contrast to the wild-type MoFe protein, the alpha-96-substituted MoFe proteins all exhibit changes in their EPR spectra upon incubation with acetylene or cyanide. The alpha-96(Leu)-substituted MoFe protein was representative of the other alpha-96-substituted MoFe proteins examined. The incubation of acetylene with the alpha-96(Leu) MoFe protein decreased the intensity of the normal FeMo-cofactor signal with the appearance of a new EPR signal having inflections at g = 4.50 and 3.50. Incubation of cyanide with the alpha-96(Leu) MoFe protein also decreased the FeMo-cofactor EPR signal with concomitant appearance of a new EPR signal having an inflection at g = 4.06. The acetylene- and cyanide-dependent EPR signals observed for the alpha-96(Leu)-substituted MoFe protein were found to follow Curie law 1/T dependence, consistent with a ground-state transition as observed for FeMo-cofactor. The microwave power dependence of the EPR signal intensity is shifted to higher power for the acetylene- and cyanide-dependent signals, consistent with a change in the relaxation properties of the spin system of FeMo-cofactor. Finally, the alpha-96(Leu)-substituted MoFe protein incubated with (13)C-labeled cyanide displays a (13)C ENDOR signal with an isotropic hyperfine coupling of 0.42 MHz in Q-band Mims pulsed ENDOR spectra. This indicates the existence of some spin density on the cyanide, and thus suggests that the new component of the cyanide-dependent EPR signals arise from the direct bonding of cyanide to the FeMo-cofactor. These data indicate that both acetylene and cyanide are able to interact with FeMo-cofactor contained within the alpha-96-substituted MoFe proteins in the resting state. These results support a model where effective interaction of substrates or inhibitors with FeMo-cofactor occurs as a consequence of both increased reactivity and accessibility of FeMo-cofactor under turnover conditions. We suggest that, for the wild-type MoFe protein, the alpha-96(Arg) side chain acts as a gatekeeper, moving during turnover in order to permit accessibility of acetylene or cyanide to a specific [4Fe-4S] face of FeMo-cofactor.