Density functional calculations for modeling the oxidized states of the active site of nickel-iron hydrogenases. 1. Verification of the method with paramagnetic Ni and Co complexes

ZORA relativistic DFT calculations are presented which aim to reproduce geometric structures and EPR properties of [Ni(mnt)(2)](-) (H(2)mnt = maleonitrildithiol), two other paramagnetic low-spin Ni(III) complexes, and an asymmetric paramagnetic Co(II) complex. The study tests the accuracy of the computational method as a prior step to the modeling of the geometric and electronic structure of the active site of NiFe hydrogenases in its EPR-active oxidized states Ni-A and Ni-B. Systematic deviations from experiment are found for the calculated g-values; relative differences among them are, however, well reproduced. Because no significant improvements have been achieved by using larger basis sets or more sophisticated functionals, g-values may be calculated rather rapidly at the VWN level. This is most important for the modeling of the active site of NiFe hydrogenases because its complexity does not permit calculations at high levels of theory. For [Ni(mnt)(2)](-), excellent agreement between calculated and experimental results is obtained for the (14)N quadrupole coupling, whereas the calculated hyperfine couplings are not always in good agreement with experimental data.     

Hydride-containing models for the active site of the nickel-iron hydrogenases

The [NiFe]-hydrogenase model complex NiFe(pdt)(dppe)(CO)(3) (1) (pdt = 1,3-propanedithiolate) has been efficiently synthesized and found to be robust. This neutral complex sustains protonation to give the first nickel-iron hydride [1H]BF(4). One CO ligand in [1H]BF(4) is readily substituted by organophosphorus ligands to afford the substituted derivatives [HNiFe(pdt)(dppe)(PR(3))(CO)(2)]BF(4), where PR(3) = P(OPh)(3) ([2H]BF(4)); PPh(3) ([3H]BF(4)); PPh(2)Py ([4H]BF(4), where Py = 2-pyridyl). Variable temperature NMR measurements show that the neutral and protonated derivatives are dynamic on the NMR time scale, which partially symmetrizes the phosphine complex. The proposed stereodynamics involve twisting of the Ni(dppe) center, not rotation at the Fe(CO)(2)(PR(3)) center. In MeCN solution, 3, which can be prepared by deprotonation of [3H]BF(4) with NaOMe, is about 10(4) stronger base than is 1. X-ray crystallographic analysis of [3H]BF(4) revealed a highly unsymmetrical bridging hydride, the Fe-H bond being 0.40 ? shorter than the Ni-H distance. Complexes [2H]BF(4), [3H]BF(4), and [4H]BF(4) undergo reductions near -1.46 V vs Fc(0/+). For [2H]BF(4), this reduction process is reversible, and we assign it as a one-electron process. In the presence of trifluoroacetic acid, proton reduction catalysis coincides with this reductive event. The dependence of i(c)/i(p) on the concentration of the acid indicates that H(2) evolution entails protonation of a reduced hydride. For [2H](+), [3H](+), and [4H](+), the acid-independent rate constants are 50-75 s(-1). For [2H](+) and [3H](+), the overpotentials for H(2) evolution are estimated to be 430 mV, whereas the overpotential for the N-protonated pyridinium complex [4H(2)](2+) is estimated to be 260 mV. The mechanism of H(2) evolution is proposed to follow an ECEC sequence, where E and C correspond to one-electron reductions and protonations, respectively. On the basis of their values for its pK(a) and redox potentials, the room temperature values of ΔG(H?) and ΔG(H-) are estimated as respectively as 57 and 79 kcal/mol for [1H](+).     

Density functional calculations of ATP systems. 2. ATP hydrolysis at the active site of actin

The hydrolysis of adenosine 5'-triphosphate (ATP) at the active site of actin has been studied using density functional calculations. The active site is modeled by the triphosphate tail of ATP, an Mg cation, surrounding water molecules, and the nearby protein residues. Four reaction paths have been followed by constraining coordinates that allow phosphate stretching, nucleophilic attack of the catalytic water, and OH(-) formation via water deprotonation. The lowest-energy barrier (21.0 kcal/mol) is obtained for a dissociative reaction where the terminal phosphate breaks on approaching the catalytic water, followed by proton release via a proton wire mechanism. A higher barrier (39.6 kcal/mol) results for an associative reaction path where OH(-) is formed first, with a pentacoordinated phosphorus atom (P-O distances 2.1 A). Stretching the terminal bridging P-O bond results in bond rupture at 2.8 A with an energy barrier of 28.8 kcal/mol. The residues Gln137 and His161 are not important in the reactions, but insight into their roles in vivo has been obtained. The favored coordination of the end products H(2)PO(4)(-) and ADP(3-) includes a hydrogen bond and an O-Mg-O bridge between the phosphates as well as a hydrogen bond between H(2)PO(4)(-) and the Ser14 side chain. The total energy is 2.1 kcal/mol lower than in the initial reactants. Classical simulations of ATP- and ADP.P(i)-actin show few hydrolysis-induced differences in the protein structure, indicating that phosphate migration is necessary for a change in conformation.     

Modeling the active sites of metalloenzymes. 4. Predictions of the unready states of [NiFe] Desulfovibrio gigas hydrogenase from density functional theory

Density functional theory has been used to predict the structures of a variety of active site models for the unready states, Ni-A and Ni-SU, of the [NiFe] hydrogenase from Desulfovibrio gigas. By comparing available experimental results on Ni-A, Ni-SU, and Ni-SI with the computational results on these model complexes, we have been able to identify the most likely formulas and structures for the active sites of Ni-A and Ni-SU. Ni-A is predicted to be a Ni(III)-Fe(II) species with the bridging hydroxo ligand, rather than the expected oxo ligand, while Ni-SU is predicted to be a Ni(II)-Fe(II) species with a water molecule coordinated to the Fe center. Both have one of the terminal S atoms (cysteines) protonated.     

Selenium-bridged diiron hexacarbonyl complexes as biomimetic models for the active site of Fe-Fe hydrogenases

Three N-substituted selenium-bridged diiron complexes [{(mu-SeCH2)2NC6H4R}Fe2(CO)6] (R = 4-NO2, 7; R = H, 8; R = 4-CH3, 9) were firstly prepared as biomimetic models for the Fe-Fe hydrogenases active site. Models could be generated by the convergent reaction of [(mu-HSe)2Fe2(CO)6] (6) with N,N-bis(hydroxymethyl)-4-nitroaniline (1), N,N-bis(hydroxymethyl)aniline (2), and N,N-bis(hydroxymethyl)-4-methylaniline (3) in 46-52% yields. All the new complexes were characterized by IR, 1H and 13C NMR and HRMS spectra and their molecular structures were determined by single-crystal X-ray analysis. The redox properties of and their dithiolate analogues [{(mu-SCH2)2NC6H4R}Fe2(CO)6] (R = 4-NO2, 7s; R = H, 8s; R = 4-CH3, 9s ) were evaluated by cyclic voltammograms. The electrochemical proton reduction by and were investigated in the presence of p-toluenesulfonic acid (HOTs) to evaluate the influence of changing the coordinating S atoms of the bridging ligands to Se atoms on the electrocatalytic activity for proton reduction.     

Active-site models for the nickel-iron hydrogenases: effects of ligands on reactivity and catalytic properties

Described are new derivatives of the type [HNiFe(SR)(2)(diphosphine)(CO)(3)](+), which feature a Ni(diphosphine) group linked to a Fe(CO)(3) group by two bridging thiolate ligands. Previous work had described [HNiFe(pdt)(dppe)(CO)(3)](+) ([1H](+)) and its activity as a catalyst for the reduction of protons (J. Am. Chem. Soc. 2010, 132, 14877). Work described in this paper focuses on the effects on properties of NiFe model complexes of the diphosphine attached to nickel as well as the dithiolate bridge, 1,3-propanedithiolate (pdt) vs 1,2-ethanedithiolate (edt). A new synthetic route to these Ni-Fe dithiolates is described, involving reaction of Ni(SR)(2)(diphosphine) with FeI(2)(CO)(4) followed by in situ reduction with cobaltocene. Evidence is presented that this route proceeds via a metastable μ-iodo derivative. Attempted isolation of such species led to the crystallization of NiFe(Me(2)pdt)(dppe)I(2), which features tetrahedral Fe(II) and square planar Ni(II) centers (H(2)Me(2)pdt = 2,2-dimethylpropanedithiol). The new tricarbonyls prepared in this work are NiFe(pdt)(dcpe)(CO)(3) (2, dcpe = 1,2-bis(dicyclohexylphosphino)ethane), NiFe(edt)(dppe)(CO)(3) (3), and NiFe(edt)(dcpe)(CO)(3) (4). Attempted preparation of a phenylthiolate-bridged complex via the FeI(2)(CO)(4) + Ni(SPh)(2)(dppe) route gave the tetrametallic species [(CO)(2)Fe(SPh)(2)Ni(CO)](2)(μ-dppe)(2). Crystallographic analysis of the edt-dcpe compund [2H]BF(4) and the edt-dppe compound [3H]BF(4) verified their close resemblance. Each features pseudo-octahedral Fe and square pyramidal Ni centers. Starting from [3H]BF(4) we prepared the PPh(3) derivative [HNiFe(edt)(dppe)(PPh(3))(CO)(2)]BF(4) ([5H]BF(4)), which was obtained as a ～2:1 mixture of unsymmetrical and symmetrical isomers. Acid-base measurements indicate that changing from Ni(dppe) (dppe = Ph(2)PCH(2)CH(2)PPh(2)) to Ni(dcpe) decreases the acidity of the cationic hydride complexes by 2.5 pK(a)(PhCN) units, from ～11 to ～13.5 (previous work showed that substitution at Fe leads to more dramatic effects). The redox potentials are more strongly affected by the change from dppe to dcpe, for example the [2](0/+) couple occurs at E(1/2) = -820 for [2](0/+) vs -574 mV (vs Fc(+/0)) for [1](0/+). Changes in the dithiolate do not affect the acidity or the reduction potentials of the hydrides. The acid-independent rate of reduction of CH(2)ClCO(2)H by [2H](+) is about 50 s(-1) (25 °C), twice that of [1H](+). The edt-dppe complex [2H](+) proved to be the most active catalyst, with an acid-independent rate of 300 s(-1).     

High level ab initio and DFT calculations of models of the catalytically active Ni-Fe hydrogenases

Multi-reference M?ller-Plesset calculations of a model of the Ni-SI state of nickel-iron hydrogenase predict a singlet rather than a triplet state for this species, and show that it is better described with a BP86 rather than a B3LYP functional.     

Spectroscopic and crystallographic evidence for the N-protonated FeIFeI azadithiolate complex related to the active site of Fe-only hydrogenases

The complex [{(mu-SCH2)2N(CH2C6H4-2-Br)}Fe2(CO)6] and its N-protonated species, as structural models of the Fe-only hydrogenase active site, were identified spectroscopically and crystallographically, and their molecular structures show the 0.04-0.1 A lengthening of the three N-C bonds and an intramolecular HBr contact (2.82 Angstroms) in the crystalline state of the N-protonated species.     

Density functional theory calculations of dense TiO2 polymorphs: implication for visible-light-responsive photocatalysts

Structural parameters and electronic band gaps of dense TiO(2) polymorphs, i.e., alpha-PbO(2), baddeleyite, fluorite, and cotunnite types of structures, were calculated using a first-principles density functional method with local-density approximation. The ambient phases, i.e., rutile and anatase, with known theoretical and experimental data were used to ensure the validity of the calculations. The fluorite-type TiO(2) turned out to have the narrowest band gap, 1.08 or 2.18 eV after applying a very approximate band gap correction, due to highly symmetrical TiO(8) polyhedra with Ti(3d) and O(2p) orbitals in the most mixed state. Ti with eight coordinated oxygens, as feasible under high pressure or residual stress, may have potential applications as a visible-light-responsive photocatalyst.     

Density functional studies of oxidized and reduced methane monooxygenase. Optimized geometries and exchange coupling of active site clusters

The conflicting protein crystallography data for the oxidized form (MMOH(ox)) of methane monooxygenase present a dilemma regarding the identity of the solvent-derived bridging ligands within the active site: do they comprise a diiron unit bridged by 1H2O and 1OH(-) as postulated for Methylococcus capsulatus or 2OH(-) ligands as suggested for Methylosinus trichosporium? Using models derived explicitly from the M. capsulatus and M. trichosporium protein data, spin-unrestricted density functional methods have been used to study two structurally characterized forms of the hydroxylase component of methane monooxygenase. The active site geometries of the oxidized (MMOH(ox)) and two-electron-reduced (MMOH(red)) states have been geometry optimized using several quantum cluster models which take into account the antiferromagnetic (AF) and ferromagnetic (F) coupling of electron spins. Trends in cluster geometries, energetics, and Heisenberg J values have been evaluated. For the majority of models, calculated geometries are in good agreement with the X-ray analyses and appear relatively insensitive to the F or AF alignment of electron spins on adjacent Fe sites. Discrepancies between calculation and experiment appear in the orientation of the coordinated His and Glu amino acid side chains for both MMOH(ox) and MMOH(red) and also in unexpected intramolecular proton transfer in the MMOH(ox) cluster models. There is additional dispersion between (and among) calculated and experimental Fe(3+)-OH(-) distances with relevance to the correct protonation state of the solvent-derived ligands. In an accompanying paper (Lovell, T.; Li, J.; Noodleman, L. Inorg. Chem. 2001, 40, 5267), a comparison of the related energetics of the active site models examined herein is further evaluated in the full protein and solvent environment.     

Density functional calculations of ATP systems. 1. Crystalline ATP hydrates and related molecules

Adenosine 5'-triphosphate (ATP) is an essential energy carrier in mammalian and other cells, and its hydrolysis to the diphosphate (ADP) in the presence of metal cations (e.g., Mg(2+) or Ca(2+)) is one of the most prevalent biochemical reactions. We describe here density functional (DF) calculations on closely related systems and compare the results with other calculations and available experimental data: Na(H2O)n +, Mg(H2O)n 2+, and Ca(H2O)n 2+ clusters (n = 1, 4-7), the crystalline pyrophosphates Mg(2)P(2)O(7).6H2O and alpha-CaNa(2)P(2)O(7).4H2O, and crystalline Na(2)ATP.3H2O. The last of these comprises asymmetric units of ATP dimers (monomers A and B) in a double-protonated state H(2)(ATP)(2-). The calculated structures agree well with available measurements and provide additional information, including the location of the H atoms. Analysis of the dipole moments of individual ATP monomers and their dimers shows that the crystal comprises blocks of opposing dipoles. Replacing one Na+ ion with Mg2+ or Ca2+ results in a significant elongation of the terminal bridging P-O bond. The calculations provide benchmarks for the use of DF methods in ATP systems and are used in the companion paper to study the hydrolysis of ATP at the active site of the protein actin.     

Density functional calculations for predicting structures and vibrational frequencies of aluminum chloride-fluoride complexes

Quantum chemical calculations are used to study AlCl_y−xF_x^3−y (y = 5 or 6, x = 0,…,y) species that can occur in aluminum electrorefining melts. These theoretical studies are included in a wider research program concerning the chemical instabilities in the bulk of molten salts during the refinement process. Stabilization energies, equilibrium geometries and vibrational frequencies of the complexes are calculated using the Delley functional methodology described in Ref. [1] (B. Delley, J. Chem. Phys., 92 (1990) 508). These computational simulations, discussed and compared with the experimental results demonstrate that density functional calculations can be reliably used in the study of complexes existing in molten salts. Quantum chemical calculations are accurate tools for theoretically predicting structures, physical and chemical properties and vibrational frequencies of known entities as well as unknown compounds.     

Fe-S complexes containing five-membered heterocycles: novel models for the active site of hydrogenases with unusual low reduction potential

Three biomimetic 2Fe2S complexes [{(micro-SCH2)2NCH2(2-C4H3O)}](Fe2(CO)6), [{(micro-SCH2)2 NCH2(2-C4H3S)}](Fe2(CO)6) and [{(micro-SCH2)2NCH2(5-Br-2-C4H2S)}Fe2(CO)6] were prepared as models for the active site of Fe-only hydrogenase by the convergent process from [(micro-S2)Fe2(CO)6] and N,N-bis(hydromethyl)-2-furan and thiophene. The structures of these complexes were identified spectroscopically and crystallographically. The electrochemical behavior of the complexes and was unique as they showed catalytic proton reduction with a low reduction potential at -1.13 and -1.09 V vs Fc/Fc+, respectively, in the presence of HClO4.     

Molecular modeling of ceftriaxone activation in the active sites of penicillin-binding proteins 2

Russian Chemical Bulletin - The enzyme—substrate complexes of penicillin-binding proteins PBP2 from FA19, 35/02, and H041 strains of Nisseria gonorrhoeae with ceftriaxone were simulated by...     

Synthesis and structural characterization of diiron azadithiolate complexes relevant to the active site of [FeFe] hydrogenases

Two N-functionally substituted diiron azadithiolate complexes, [(µ-SCH₂)₂NCH₂CH₂OC(O)C₆H₄I-p]Fe₂(CO)₆ (1) and {[(µ-SCH₂)₂NCH₂CH₂OC(O)C₆H₄I-p]Fe₂(CO)₅Ph₂PCH}₂ (2) as models for the active site of [FeFe] hydrogenases, have been prepared and fully characterized. Complex 1 was prepared by the reaction of [(µ-SCH₂)₂NCH₂CH₂OH]Fe₂(CO)₆ with p-iodobenzoic acid in the presence of 4-dimethylaminopyridine (DMAP) and N,N′-dicyclohexylcarbodiimide (DCC) in 78% yield. Further treatment of 1 with 1 equiv. of Me₃NO?·?2H₂O followed by 0.5 equiv. of trans-1,2-bis(diphenylphosphino)ethylene (dppe) affords 2 in 60% yield. The new complexes 1 and 2 were characterized by IR and ¹H (¹³C, ³¹P) NMR spectroscopic techniques and their molecular structures were confirmed by X-ray diffraction analysis. The molecular structure of 1 has two conformational isomers, in one isomer its N-functional substituent is axial to its bridged nitrogen and in the other isomer its N-functional substituent is equatorial. The crystal structure of 2 revealed that its N-functional substituents are equatorial to its nitrogens and dppe occupies the two apical positions of the square-pyramidal irons.     

Density functional calculations of structures and ionization energies for heavy group V cluster anions

The structure and electronic structure of heavy-group V cluster anions (Sb, Bi) are calculated with density functional methods within the local spin density approximation (LSDA ). The influence of gradient corrections of the exchange and correlation energy is investigated. The calculated vertical and adiabatic ionization energies are in very good agreement with data from photoelectron spectroscopy (PES ) for Sb, whereas the relatively large deviations for Bi can be reduced by the consideration of relativistic effects in a scalar-relativistic manner. Concerning the structures, a strong similarity to the corresponding P clusters was found. In particular, the negatively charged pentamers are planar rings (with similarities to the aromatic [C₅H₅]^? anion) with especially high ionization energies. © 1995 John Wiley & Sons, Inc.     

Density functional theory calculations and exploration of a possible mechanism of N2 reduction by nitrogenase

Density functional theory (DFT) calculations have been performed on the nitrogenase cofactor, FeMoco. Issues that have been addressed concern the nature of M-M interactions and the identity and origin of the central light atom, revealed in a recent crystallographic study of the FeMo protein of nitrogenase (Einsle, O.; et al. Science 2002, 297, 871). Introduction of Se in place of the S atoms in the cofactor and energy minimization results in an optimized structure very similar to that in the native enzyme. The nearly identical, short, lengths of the Fe-Fe distances in the Se and S analogues are interpreted in terms of M-M weak bonding interactions. DFT calculations with O or N as the central atoms in the FeMoco marginally support the assignment of the central atom as N rather than O. The assumption was made that the central atom is the N atom, and steps of a catalytic cycle were calculated starting with either of two possible states for the cofactor and maintaining the same charge throughout (by addition of equal numbers of H(+) and e(-)) between steps. The states were [(Cl)Fe(II)(6)Fe(III)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)](2-), [I-N-3H](2-), and [(Cl)Fe(II)(4)Fe(III)(3)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)], [I-N-3H](0) (Gl = deprotonated glycol; Im = imidazole). These are the triply protonated ENDOR/ESEEM [I-N](5-) and M?ssbauer [I-N](3-) models, respectively. The proposed mechanism explores the possibilities that (a) redox-induced distortions facilitate insertion of N(2) and derivative substrates into the Fe(6) central unit of the cofactor, (b) the central atom in the cofactor is an exchangeable nitrogen, and (c) the individual steps are related by H(+)/e(-) additions (and reduction of substrate) or aquation/dehydration (and distortion of the Fe(6) center). The Delta E's associated with the individual steps of the proposed mechanism are small and either positive or negative. The largest positive Delta E is +121 kJ/mol. The largest negative Delta E of -333 kJ/mol is for the FeMoco with a N(3-) in the center (the isolated form) and an intermediate in the proposed mechanism.