g- and A-Tensor Calculations in the Zero-Order Approximation for Relativistic Effects of Ni Complexes and Ni(CO)3H as Model Complexes for the Active Center of [NiFe]-Hydrogenase

Synthesis, Crystal Structure, EPR, and DFT Studies of an Unusually Distorted Vanadium(II) Complex

We report a synthesis and structure of the most highly distorted four-coordinate d3 ion known to date that also serves as the second known example of a bis(biphenolato) transition metal...

H2 Evolution from a Thiolate-Bound Ni(III) Hydride

Terminal Ni^III hydrides are proposed intermediates in proton reduction catalyzed by both molecular electrocatalysts and metalloenzymes, but well-defined examples of paramagnetic nickel hydride complexes are largely limited to bridging hydrides. Herein, we report the synthesis of an S = 1/2, terminally bound thiolate-Ni^III-H complex. This species and its terminal hydride ligand in particular have been thoroughly characterized by vibrational and EPR techniques, including pulse EPR studies. Corresponding DFT calculations suggest appreciable spin leakage onto the thiolate ligand. The hyperfine coupling to the terminal hydride ligand of the thiolate-Ni^III-H species is comparable to that of the hydride ligand proposed for the Ni-C hydrogenase intermediate (Ni^III-H-Fe^II). Upon warming, the featured thiolate-Ni^III-H species undergoes bimolecular reductive elimination of H₂. Associated kinetic studies are discussed and compared with a structurally related Fe^III-H species that has also recently been reported to undergo bimolecular H-H coupling.

A Biochemical Nickel(I) State Supports Nucleophilic Alkyl Addition: A Roadmap for Methyl Reactivity in Acetyl Coenzyme A Synthase

Nickel-containing enzymes such as methyl coenzyme M reductase (MCR) and carbon monoxide dehydrogenase/acetyl coenzyme A synthase (CODH/ACS) play a critical role in global energy conversion reactions, with significant contributions to carbon-centered processes. These enzymes are implied to cycle through a series of nickel-based organometallic intermediates during catalysis, though identification of these intermediates remains challenging. In this work, we have developed and characterized a nickel-containing metalloprotein that models the methyl-bound organometallic intermediates proposed in the native enzymes. Using a nickel(I)-substituted azurin mutant, we demonstrate that alkyl binding occurs via nucleophilic addition of methyl iodide as a methyl donor. The paramagnetic Ni^III-CH₃ species initially generated can be rapidly reduced to a high-spin Ni^II-CH₃ species in the presence of exogenous reducing agent, following a reaction sequence analogous to that proposed for ACS. These two distinct bioorganometallic species have been characterized by optical, EPR, XAS, and MCD spectroscopy, and the overall mechanism describing methyl reactivity with nickel azurin has been quantitatively modeled using global kinetic simulations. A comparison between the nickel azurin protein system and existing ACS model compounds is presented. Ni^III-CH₃ Az is only the second example of two-electron addition of methyl iodide to a Ni^I center to give an isolable species and the first to be formed in a biologically relevant system. These results highlight the divergent reactivity of nickel across the two intermediates, with implications for likely reaction mechanisms and catalytically relevant states in the native ACS enzyme.

A bioorganometallic model for acetyl coenzyme A synthase has been developed. This model protein is able to bind a cationic methyl group via direct addition to the nickel(I) center. The resultant nickel(III)-methyl species has been characterized via optical and electron paramagnetic resonance spectroscopy, and the reduced nickel(II)-methyl state has been characterized using magnetic circular dichroism and X-ray spectroscopy. Implications for further reactivity with CO are gleaned from electronic structure analysis of the nickel-methyl species.

DFT Protocol for EPR Prediction of Paramagnetic Cu(II) Complexes and Application to Protein Binding Sites

With the aim to provide a general protocol to interpret electron paramagnetic resonance (EPR) spectra of paramagnetic copper(II) coordination compounds, density functional theory (DFT) calculations of spin Hamiltonian parameters g and A for fourteen Cu(II) complexes with different charges, donor sets, and geometry were carried out using ORCA software. The performance of eleven functionals was tested, and on the basis of the mean absolute percent deviation (MAPD) and standard deviation (SD), the ranking of the functionals for Az is: B3LYP > B3PW91 ~ B3P86 > PBE0 > CAM-B3LYP > TPSSh > BH and HLYP > B2PLYP > MPW1PW91 > ω-B97x-D >> M06; and for gz is: PBE0 > BH and HLYP > B2PLYP > ω-B97x-D > B3PW91~B3LYP~B3P86 > CAM-B3LYP > TPSSh~MPW1PW91 >> M06. With B3LYP the MAPD with respect to A z exp t l is 8.6% with a SD of 4.2%, while with PBE0 the MAPD with respect to g z exp t l is 2.9% with a SD of 1.1%. The results of the validation confirm the fundamental role of the second order spin-orbit contribution to Az. The computational procedure was applied to predict the values of gz and Az of the adducts formed by Cu(II) with albumin and two fragments of prion protein, 106–126 and 180–193.

Molecular qubits based on potentially nuclear-spin-free nickel ions

Molecular qubits with the longest coherence times thus far are based on nuclear-spin-carrying central ions. These nuclear spins can cause quantum state leakage, which is detrimental to quantum algorithm performance. We present two novel molecular qubits based on potentially nuclear spin-free Ni in the formal oxidation state 3+. (d₂₀-PPh₄)[Ni(mnt)₂] (Ni-mnt, mnt^2- = maleonitrile-1,2-dithiolate) possesses a coherence time of up to 38.7 μs at 7 K. Functionalization of the dithiolate ligand decreases the coherence time by a factor of only four in (HNEt₃)[Ni(dip)₂] (Ni-dip, dip^2- = 3-(diphenylphosphoryl)-methylbenzene-1,2-dithiolate), indicating that monoanionic Ni-dithiolene complexes are promising and robust building blocks for polynuclear molecular qubit gates.

Optimized Basis Sets for Calculation of Electron Paramagnetic Resonance Hyperfine Coupling Constants: aug-cc-pVTZ-J for the 3d Atoms Sc-Zn

The hyperfine coupling tensor of electron paramagnetic resonance (EPR), describing the interaction between an electron and a given nuclei, depends strongly on the electron density at the nucleus. With standard Gaussian-type orbital basis sets (GTOs), employed in most calculations, it is difficult to obtain converged results of the hyperfine coupling tensor, and basis sets with more flexible core regions have therefore been devised. To this class of core property basis sets belong the aug-cc-pVTZ-J basis sets developed for the s- and p-block atoms. Here, we extend the aug-cc-pVTZ-J basis sets to include the 3d elements Sc-Zn. The converged optimal basis sets are throughout the series described by a (25s17p10d3f2g)/[17s10p7d3f2g] contraction scheme, where four tight s-, one tight p-, and one tight d-type function have been added to the original aug-cc-pVTZ basis sets. The basis sets are generally contracted, and molecular orbital coefficients are used as contraction coefficients. By validation studies with different functionals and compounds, it is shown that the values of the contraction coefficient are effectively independent of the compound used in their generation and the exchange-correlation functional employed in the calculation.

Copper(II) complexes of rat amylin fragments

The fragments of rat amylin rIAPP(17-29) (Ac-VRSSNNLGPVLPP-NH(2)), rIAPP(17-22) (Ac-VRSSNN-NH(2)), rIAPP(19-22) (Ac-SSNN-NH(2)) and rIAPP(17-20) (Ac-VRSS-NH(2)) together with the related mutant peptides (Ac-VASS-NH(2) and Ac-VRAA-NH(2)) have been synthesized and their copper(II) complexes studied by potentiometric, UV-Vis, CD and EPR spectroscopic methods. Despite the lack of any common strongly coordinating donor functions some of these fragments are able to bind copper(II) ions in the physiological pH range. The longest fragment rat amylin(17-29) keeps one equivalent copper(II) ion in solution in the whole pH range, while two other peptides Ac-VRSSNN-NH(2) and Ac-SSNN-NH(2) are also able to interact with copper(II) ions in the slightly alkaline pH range. According to the spectral parameters of the complexes, the peptides can be classified into two different categories: (i) the tetrapeptides Ac-VRSS-NH(2), Ac-VASS-NH(2) and Ac-VRAA-NH(2) can interact with copper(II) only under strongly alkaline conditions (pH > 10.0) and the formation of only one species with four amide nitrogen coordination can be detected; (ii) the peptides Ac-VRSSNNLGPVLPP-NH(2), Ac-VRSSNN-NH(2) and Ac-SSNN-NH(2) can form complexes above pH 6.0 with the major stoichiometries [CuH(-2)L], [CuH(-3)L](-) and [CuH(-4)L](2-). These data support that rIAPP(17-29) can interact with copper(II) ions under physiological conditions and the SSNN tetrapeptide fragment can be considered as the shortest sequence responsible for metal binding. Density functional theory (DFT) calculations provide some information on the possible coordination modes of Ac-SSNN-NH(2) towards the copper(II) ion and suggest that for [CuH(-2)L], [CuH(-3)L](-) and [CuH(-4)L](2-), the binding of two, three and four deprotonated amide nitrogens, with NH(-) of the side chain of asparagine as anchoring group, is probable. Moreover, these data reveal that peptides can be effective metal binding ligands even in the absence of anchoring groups, if more polar side chains are present in a specific sequence.

Is the spin-orbit coupling important in the prediction of the 51V hyperfine coupling constants of V(IV) O2+ species? ORCA versus Gaussian performance and biological applications

Density functional theory calculations of the (51)V hyperfine coupling (HFC) tensor A, have been completed for eighteen V(IV)O(2+) complexes with different donor set, electric charge and coordination geometry. A tensor was calculated with ORCA software with several functionals and basis sets taking into account the spin-orbit coupling contribution. The results were compared with those obtained with Gaussian 03 software using the half-and-half functional BHandHLYP and 6-311g(d,p) basis set. The order of accuracy of the functionals in the prediction of A(iso), A(z) and dipolar term A(z,anis) is BHandHLYP > PBE0 >> B3PW > TPSSh >> B3LYP >> BP86 > VWN5 (for A(iso)), BHandHLYP > PBE0 >> B3PW > TPSSh > B3LYP >> BP86 > VWN5 (for A(z)), B3LYP > PBE0 ∼ B3PW ∼ BHandHLYP >> TPSSh > BP86 ∼ VWN5 (for A(z,anis)). The good agreement in the prediction of A(z) with BHandHLYP is due to a compensation between the overestimation of A(iso) and underestimation of A(z,anis) (A(z) = A(iso) + A(z,anis)), whereas among the hybrid functionals PBE0 performs better than the other ones. BHandHLYP functional and Gaussian software are recommended when the V(IV)O(2+) species contains only V-O and/or V-N bonds, whereas PBE0 functional and ORCA software for V(IV)O(2+) complexes with one or more V-S bonds. Finally, the application of these methods to the coordination environment of V(IV)O(2+) ion in V-proteins, like vanadyl-substituted insulin, carbonic anhydrase, collagen and S-adenosylmethionine synthetase, was discussed.

Observation of intermolecular charge transfer in a quasi-one-dimensional molecular alloy system

X-band single-crystal electron paramagnetic resonance (EPR) studies of the molecular alloy [NO 2BzPy][Au 0.57Ni 0.43(mnt) 2] are presented in this paper. At room temperature, EPR spectra show both intense resonance signals (main signals) and weak satellite quartet lines. The characteristics of both intense and weak EPR signals depend on the magnetic field orientation. The main signals arise from two magnetically nonequivalent [Ni(mnt) 2] (-) anions, and their corresponding principal values of the g tensor are ( g x ') 1 = 2.04653, ( g y ') 1 = 2.00096, and ( g z ') 1 = 2.15319 and ( g x ') 2 = 2.04520, ( g y ') 2 = 1.99734, and ( g z ') 2 = 2.15361, respectively. The weak satellite lines, whose patterns strongly depend on the magnetic field direction, can be attributed to the hyperfine coupling of the electron spin with the (197)Au nucleus of the [Au(mnt) 2] (-) species. Density functional theory calculations for the spin and charge distributions of the dimer {[Ni(mnt) 2][Au(mnt) 2]} (2-) indicate that the hyperfine interaction of the electron spin with the (197)Au nuclear spins is caused, in part, by the charge transfer between the [Ni(mnt) 2] (-) and the [Au(mnt) 2] (-) species.

Copper binding sites in the C-terminal domain of mouse prion protein: A hybrid (QM/MM) molecular dynamics study

We present a hybrid QM/MM Car-Parrinello molecular dynamics study of the copper-loaded C-terminal domain of the mouse prion protein. By means of a statistical analysis of copper coordination in known protein structures, we localized the protein regions with the highest propensity for copper ion binding. The identified candidate structures were subsequently refined via QM/MM simulations. Their EPR characteristics were computed to make contact with the experimental data and to probe the sensitivity to structural and chemical changes. Overall best agreement with the experimental EPR data (Van Doorslaer et al., J Phys Chem B 2001; 105: 1631-1639) and the information currently available in the literature is observed for a binding site involving H187. Moreover, a reinterpretation of the experimental proton hyperfine couplings was possible in the light of the present computational findings.

Redox Non-innocence of Thioether Macrocycles: Elucidation of the Electronic Structures of Mononuclear Complexes of Gold(II) and Silver(II)

The mononuclear +2 oxidation state metal complexes [Au([9]aneS(3))(2)](2+) and [Ag([18]aneS(6))](2+) have been synthesized and characterized crystallographically. The crystal structure of the Au(II) species [Au([9]aneS(3))(2)](BF(4))(2) shows a Jahn-Teller tetragonally distorted geometry with Au-S(1) = 2.839(5), Au-S(2) = 2.462(5), and Au-S(3) = 2.452(5) A. The related Ag(II) complex [Ag([18]aneS(6))](ClO(4))(2) has been structurally characterized at both 150 and 30 K and is the first structurally characterized complex of Ag(II) with homoleptic thioether S-coordination. The single-crystal X-ray structure of [Ag([18]aneS(6))](ClO(4))(2) confirms octahedral homoleptic S(6)-thioether coordination. At 150 K, the structure contains two independent Ag(II)-S distances of 2.569(7) and 2.720(6) A. At 30 K, the structure retains two independent Ag(II)-S distances of 2.615(6) and 2.620(6) A, with the complex cation retaining 3-fold symmetry. The electronic structures of [Au([9]aneS(3))(2)](2+) and [Ag([18]aneS(6))](2+) have been probed in depth using multifrequency EPR spectroscopy coupled with DFT calculations. For [Au([9]aneS(3))(2)](2+), the spectra are complex due to large quadrupole coupling to (197)Au. Simulation of the multifrequency spectra gives the principal g values, hyperfine (A) and quadrupole (P) couplings, and furthermore reveals non-co-incidence of the principal axes of the P tensor with respect to the A and g matrices. These results are rationalized in terms of the electronic and geometric structure and reveal that the SOMO has ca. 30% Au 5d(xy)() character, consistent with DFT calculations (27% Au character). For [Ag([18]aneS(6))](2+), detailed EPR spectroscopic analysis confirms that the SOMO has ca. 26% Ag 4d(xy)() character and DFT calculations are consistent with this result (22% Ag character).

Computational Studies of the H-Cluster of Fe-Only Hydrogenases: Geometric, Electronic, and Magnetic Properties and Their Dependence on the [Fe4S4] Cubane

The active sites of Fe-only hydrogenases (FeHases) feature an unusual polynuclear iron-sulfur cluster, known as the H-cluster, that consists of a [Fe4S4] cubane linked to a di-iron subunit (the [2Fe]H component) via a bridging cysteine ligand (SCys). While previous computational studies of FeHases employed H-cluster models that only included the [2Fe]H component, we have utilized density functional theory (DFT), in conjunction with the broken-symmetry (BS) approach, to explore the geometric, electronic, and magnetic properties of the entire H-cluster. These calculations have allowed us to evaluate, for the first time, the influence of the [Fe4S4] cubane on the [2Fe]H component of the H-cluster in its active (Hox) and CO-inhibited (Hox-CO) states, both of which are paramagnetic (S=1/2). Our results reveal that the presence of the cubane tunes both the position and the donor strength of the SCys ligand, which, in turn, modulates the internal geometric and electronic structures of the [2Fe]H subcluster. Importantly, the BS methodology provides an accurate description of the exchange interactions within the H-cluster, permitting insight into the electronic origin of the changes in magnetic properties observed experimentally upon conversion of Hox to Hox-CO. Specifically, while the unpaired spin density in the Hox state is localized on the distal Fe center, in the Hox-CO state, it is delocalized over the [2Fe]H component, such that the proximal Fe center acquires significant spin density (where distal and proximal refer to the positions of the Fe centers relative to the cubane). To validate our H-cluster models on the basis of experimental data, two DFT-based approaches and the semiempirical INDO/S method have been employed to compute electron paramagnetic resonance parameters for the H-cluster states. While most computations yield reasonably accurate g values and ligand hyperfine coupling constants (i.e., A values) for the Hox and Hox-CO states, they fail to reproduce the isotropic 57Fe A tensors found experimentally. Finally, extension of the computational methodology employed successfully for the Hox and Hox-CO states to the metastable Hoxphoto state, generated by irradiation of the Hox-CO state at cryogenic temperatures, has allowed us to discriminate between proposed structural models for this species.

Theoretical Studies on the Magnetic Switching Controlled by Stacking Patterns of Bis(maleonitriledithiolato) Nickelate(III) Dimers

Magnetic switchable maleonitriledithiolate (mnt) complexes were studied by density functional theory. The calculations were performed for anion dimers of [RBzPyR'][Ni(mnt)(2)] (RBzPyR' = derivatives of benzylpyridinium) to elucidate magnetostructural correlations and the nature of the weak intermolecular chemical bonding. The calculated results showed that the spin delocalization, favored by the eclipsed stacking and the shorter interlayer distance, was responsible for the diamagnetic character of [1-benzyl-4-aminopyridinium][Ni(mnt)(2)] at low temperature. The weak antiferromagnetic and ferromagnetic interactions were also reproduced for [1-benzyl-4-aminopyridinium][Ni(mnt)(2)] and [1-(4'-fluorobenzyl)pyridinium][Ni(mnt)(2)] at high temperature, respectively. The natural bond orbital analysis suggested that the cooperative effect of the weak intermolecular bondings may be the intrinsic driving force resulting in the switchable property, which is essentially similar to those in organic radicals exhibiting magnetic bistability. Further investigations with varying interlayer distance d, the extent of slippage (slipping distance r and deviation angle alpha), and rotational angle theta suggested that the extent of slippage played an important role in magnetic interactions. Therefore, the abrupt modulation of the extent of slippage in the [Ni(mnt)(2)](-) complexes by external perturbations provided new possibilities for the design of molecular magnetic switching devices.

Relativistic DFT calculation of the reaction cycle intermediates of [NiFe] hydrogenase: a contribution to understanding the enzymatic mechanism

Structures and spectroscopic observables of the paramagnetic intermediates of the enzymatic reaction cycle of the metalloenzyme [NiFe] hydrogenase were calculated using relativistic density functional theory (DFT) within the zero-order regular approximation (ZORA). By comparing experimental and calculated magnetic resonance parameters (g- and hyperfine tensors) for the states Ni-A, Ni-B, Ni-C, Ni-L, and Ni-CO the details of the atomic composition of these paramagnetic intermediates could be elucidated that are mostly not available from X-ray structure analysis. In general, good agreement between calculated and experimental observables could be obtained. A detailed picture of the changes of the active center during the catalytic cycle was deduced from the obtained structures. Based on these results, a consistent model for the sequence of redox states including protonation steps is proposed which is important for understanding the mechanism of the [NiFe] hydrogenase.

A theoretical study of spin states in Ni-S4 complexes and models of the [NiFe] hydrogenase active site

We have applied density functional theory, using both pure (BP86) and hybrid (B3LYP and B3LYP*) functionals, to investigate structural parameters and reaction energies for nickel(II)-sulfur coordination compounds, as well as for small cluster models of the Ni-SI and Ni-R redox state of [NiFe] hydrogenases. Results obtained investigating experimentally well-characterized complexes show that BP86 is well suited to describe the structural features of this class of compounds. However, the singlet-triplet energy splitting and even the computed ground state are strongly dependent on the applied functional. Results for the cluster models of [NiFe] hydrogenases lead to the conclusion that in the reduced protein structures characterized by X-ray diffraction a hydride bridges the two metal centres. The energy splitting of the singlet and triplet states in Ni-R and Ni-SI models is calculated to be very small and may be overcome at room temperature to allow a spin crossover. Moreover, the relative stability of the Ni-SI and Ni-R structures adopted in the present investigation is fully compatible with their involvement in the reversible heterolytic cleavage of H(2).

Activation of Molecular Hydrogen over a Binuclear Complex with Rh2S2 Core: DFT Calculations and NMR Mechanistic Studies

The dicationic complex [(triphos)Rh(mu-S)(2)Rh(triphos)](2+), 1 (modeled as 1c) [triphos = CH(3)C(CH(2)PPh(2))(3)], is known to activate two dihydrogen molecules and produce the bis(mu-hydrosulfido) product [(triphos)(H)Rh(mu-SH)(2)Rh(H)(triphos)](2+), 2 (modeled as 2b), from which 1 is reversibly obtained. The possible steps of the process have been investigated with DFT calculations. It has been found that each d(6) metal ion in 1c, with local square pyramidal geometry, is able to anchor one H(2) molecule in the side-on coordination. The step is followed by heterolytic splitting of the H-H bond over one adjacent and polarized Rh-S linkage. The process may be completed before the second H(2) molecule is added. Alternatively, both H(2) molecules are trapped by the Rh(2)S(2) core before being split in two distinct steps. Since the ambiguity could not be solved by calculations, (31)P and (1)H NMR experiments, including para-hydrogen techniques, have been performed to identify the actual pathway. In no case is there experimental evidence for any Rh-(eta(2)-H(2)) adduct, probably due to its very short lifetime. Conversely, (1)H NMR analysis of the hydride region indicates only one reaction intermediate which corresponds to the monohydride-mu-hydrosulfide complex [(triphos)Rh(H)(mu-SH)(mu-S)Rh(triphos)](2+) (3) (model 5a). This excludes the second hypothesized pathway. From an energetic viewpoint the computational results support the feasibility of the whole process. In fact, the highest energy for H(2) activation is 8.6 kcal mol(-1), while a larger but still surmountable barrier of 34.6 kcal mol(-1) is in line with the reversibility of the process.

Supramolecular Control over Donor−Acceptor Photoinduced Charge Separation

A novel donor-bridge-acceptor system has been synthesized by covalently linking a p-phenylene vinylene oligomer (OPV) and a perylene diimid (PERY) at opposite ends of a m-phenylene ethynylene oligomer (FOLD) of twelve phenyl rings, containing nonpolar (S)-3,7-dimethyl-1-octanoxy side chains. For comparison, model compounds have been prepared in which either the donor or acceptor is absent. In chloroform, the oligomeric bridge is in a random coil conformation. Upon addition of an apolar solvent (heptane) the oligomeric bridge first folds into a helical stack and subsequently intermolecular self-assembly of the stacks into columnar architectures occurs. Photoexcitation in the random coil conformation, where the interaction between the donor and acceptor chromophores is small, results only in long-range intramolecular energy transfer in which the OPV singlet-excited state is transformed into the PERY singlet-excited state. In the folded conformation of the bridge, donor and acceptor are closer and their enhanced interaction favors the formation the OPV(*)(+)-FOLD-PERY(*)(-) charge-separated state upon photoexcitation. As a result, the extent of photoinduced charge separation depends on the degree of folding of the bridge between donor and acceptor and therefore on the apolar nature of the medium. As a consequence, and contrary to conventional photoinduced charge separation processes, the formation of the OPV(*)(+)-FOLD-PERY(*)(-) charge-separated state is more favored in apolar media.

Spin Distribution and the Location of Protons in Paramagnetic Proteins

▪ Abstract  Two current frontiers in EPR research are high-field ([Formula: see text]) electron paramagnetic resonance (EPR) and high-field electron-nuclear double resonance (ENDOR). This review focuses on recent advances in high-field ENDOR and its applications to the study of proteins containing native paramagnetic sites. It concentrates on two aspects; the first concerns the determination of the location of protons and is related to the site geometry, and the second focuses on the spin density distribution within the site, which is inherent to the electronic structure. Both spin density and proton locations can be derived from ligand hyperfine couplings determined by ENDOR measurements. A brief description of the experimental methods is presented along with a discussion of the advantages and disadvantages of high-field ENDOR compared with conventional X-band (∼ 9.5 GHz) experiments. Specific examples of both protein single crystals and frozen solutions are then presented. These include the determination of the coordinates of water ligand protons in the Mn(II) site of concanavalin A, the detection of hydrogen bonds in a quinone radical in the bacterial photosynthetic reaction center as well as in the tyrosyl radical in ribonuclease reductase, and the study of the spin distribution in copper proteins. The copper proteins discussed are the type I copper of azurin and the binuclear CuA center in a number of proteins. The last part of the review presents a brief discussion of the interpretation of hyperfine couplings using quantum chemical calculations, primarily density functional theory (DFT) methods. Such methods are becoming an integral part of the data analysis tools, as they can facilitate signal assignment and provide the ultimate relation between the experimental hyperfine couplings and the electronic wave function.

Relativistic spin-orbit effects on hyperfine coupling tensors by density-functional theory

A second-order perturbation theory treatment of spin-orbit corrections to hyperfine coupling tensors has been implemented within a density-functional framework. The method uses the all-electron atomic mean-field approximation and/or spin-orbit pseudopotentials in incorporating one- and two-electron spin-orbit interaction within a first-principles framework. Validation of the approach on a set of main-group radicals and transition metal complexes indicates good agreement between all-electron and pseudopotential results for hyperfine coupling constants of the lighter nuclei in the system, except for cases in which scalar relativistic effects become important. The nonrelativistic Fermi contact part of the isotropic hyperfine coupling constants is not always accurately reproduced by the exchange-correlation functionals employed, particularly for the triplet and pi-type doublet radicals in the present work. For this reason, ab initio coupled-cluster singles and doubles with perturbative triples results for the first-order contributions have been combined in the validation calculations with the density-functional results for the second-order spin-orbit contributions. In the cases where spin-orbit corrections are of significant magnitude relative to the nonrelativistic first-order terms, they improve the agreement with experiment. Antisymmetric contributions to the hyperfine tensor arise from the spin-orbit contributions and are discussed for the IO2 radical, whereas rovibrational effects have been evaluated for RhC, NBr, and NI.

Direct Detection of a Hydrogen Ligand in the [NiFe] Center of the Regulatory H2-Sensing Hydrogenase from Ralstonia eutropha in Its Reduced State by HYSCORE and ENDOR Spectroscopy

The regulatory H2-sensing [NiFe] hydrogenase of the beta-proteobacterium Ralstonia eutropha displays an Ni-C "active" state after reduction with H2 that is very similar to the reduced Ni-C state of standard [NiFe] hydrogenases. Pulse electron nuclear double resonance (ENDOR) and four-pulse ESEEM (hyperfine sublevel correlation, HYSCORE) spectroscopy are applied to obtain structural information on this state via detection of the electron-nuclear hyperfine coupling constants. Two proton hyperfine couplings are determined by analysis of ENDOR spectra recorded over the full magnetic field range of the EPR spectrum. These are associated with nonexchangeable protons and belong to the beta-CH(2) protons of a bridging cysteine of the NiFe center. The signals of a third proton exhibit a large anisotropic coupling (Ax = 18.4 MHz, Ay = -10.8 MHz, Az = -18 MHz). They disappear from the 1H region of the ENDOR spectra after exchange of H2O with 2H2O and activation with 2H2 instead of H2 gas. They reappear in the 2H region of the ENDOR and HYSCORE spectra. Based on a comparison with the spectroscopically similar [NiFe] hydrogenase of Desulfovibrio vulgaris Miyazaki F, for which the g-tensor orientation of the Ni-C state with respect to the crystal structure is known (Foerster et al. J. Am. Chem. Soc. 2003, 125, 83-93), an assignment of the 1H hyperfine couplings is proposed. The exchangeable proton resides in a bridging position between the Ni and Fe and is assigned to a formal hydride ion. After illumination at low temperature (T = 10 K), the Ni-L state is formed. For the Ni-L state, the strong hyperfine coupling observed for the exchangeable hydrogen in Ni-C is lost, indicating a cleavage of the metal-hydride bond(s). These experiments give first direct information on the position of hydrogen binding in the active NiFe center of the regulatory hydrogenase. It is proposed that such a binding situation is also present in the active Ni-C state of standard hydrogenases.

Nickel complexes of o-amidochalcogenophenolate(2-)/o-iminochalcogenobenzosemiquinonate(1-) pi-radical: synthesis, structures, electron spin resonance, and x-ray absorption spectroscopic evidence

The preparation of complexes trans-[Ni(-SeC(6)H(4)-o-NH-)(2)](-) (1), cis-[Ni(-TeC(6)H(4)-o-NH-)(2)](-) (2), trans-[Ni(-SC(6)H(4)-o-NH-)(2)](-) (3), and [Ni(-SC(6)H(4)-o-S-)(2)](-) (4) by oxidative addition of 2-aminophenyl dichalcogenides to anionic [Ni(CO)(SePh)(3)](-) proves to be a successful approach in this direction. The cis arrangement of the two tellurium atoms in complex 2 is attributed to the intramolecular Te.Te contact interaction (Te.Te contact distance of 3.455 A). The UV-vis electronic spectra of complexes 1 and 2 exhibit an intense absorption at 936 and 942 nm, respectively, with extinction coefficient epsilon > 10000 L mol(-)(1) cm(-)(1). The observed small g anisotropy, the principal g values at g(1) = 2.036, g(2) = 2.062, and g(3) = 2.120 for 1 and g(1) = 2.021, g(2) = 2.119, and g(3) = 2.250 for 2, respectively, indicates the ligand radical character accompanied by the contribution of the singly occupied d orbital of Ni(III). The X-ray absorption spectra of all four complexes show L(III) peaks at approximately 854.5 and approximately 853.5 eV. This may indicate a variation of contribution of the Ni(II)-Ni(III) valence state. According to the DFT calculation, the unpaired electron of complex 1 and 2 is mainly distributed on the 3d(xz)() orbital of the nickel ion and on the 4p(z)() orbital of selenium (tellurium, 5p(z)()) as well as the 2p(z)() orbital of nitrogen of the ligand. On the basis of X-ray structural data, UV-vis absorption, electron spin resonance, magnetic properties, DFT computation, and X-ray absorption (K- and L-edge) spectroscopy, the monoanionic trans-[Ni(-SeC(6)H(4)-o-NH-)(2)](-) and cis-[Ni(-TeC(6)H(4)-o-NH-)(2)](-) complexes are appositely described as a resonance hybrid form of Ni(III)-bis(o-amidochalcogenophenolato(2-)) and Ni(II)-(o-amidochalcogenophenolato(2-))-(o-iminochalcogenobenzosemiquinonato(1-) pi-radical; i.e., complexes 1 and 2 contain delocalized oxidation levels of the nickel ion and ligands.

Quantum chemical calculations of spectroscopic properties of metalloproteins and model compounds: EPR and Mössbauer properties

Recently developed theoretical methods to predict EPR and Mössbauer parameters open the way for close interactions between theorists and experimentalists to elucidate the geometric and electronic structures of metalloenzymes and model complexes and to obtain insight into their reactive properties. Spectral calculations (g-values, hyperfine couplings, zero-field splittings, isomer shifts and quadrupole splittings) are also a means to validate theoretical approaches and therefore complement the prediction of geometries, reaction energies and transition states.

Single crystal EPR studies of the reduced active site of [NiFe] hydrogenase from Desulfovibrio vulgaris Miyazaki F

In the catalytic cycle of [NiFe] hydrogenase the paramagnetic Ni-C intermediate is of key importance, since it is believed to carry the substrate hydrogen, albeit in a yet unknown geometry. Upon illumination at low temperatures, Ni-C is converted to the so-called Ni-L state with markedly different spectroscopic parameters. It is suspected that Ni-L has lost the "substrate hydrogen". In this work, both paramagnetic states have been generated in single crystals obtained from the [NiFe] hydrogenase from Desulfovibrio vulgaris Miyazaki F. Evaluation of the orientation dependent spectra yielded the magnitudes of the g tensors and their orientations in the crystal axes system for both Ni-C and Ni-L. The g tensors could further be related to the atomic structure by comparison with the X-ray crystallographic structure of the reduced enzyme. Although the g tensor magnitudes of Ni-C and Ni-L are quite different, the orientations of the resulting g tensors are very similar but differ from those obtained earlier for Ni-A and Ni-B (Trofanchuk et al. J. Biol. Inorg. Chem. 2000, 5, 36-44). The g tensors were also calculated by density functional theory (DFT) methods using various structural models of the active site. The calculated g tensor of Ni-C is, concerning magnitudes and orientation, in good agreement with the experimental one for a formal Ni(III) oxidation state with a hydride (H(-)) bridge between the Ni and the Fe atom. Satisfying agreement is obtained for the Ni-L state when a formal Ni(I) oxidation state is assumed for this species with a proton (H(+)) removed from the bridge between the nickel and the iron atom.

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.

Calculation of electronic g-tensors for transition metal complexes using hybrid density functionals and atomic meanfield spin-orbit operators

We report the first implementation of the calculation of electronic g-tensors by density functional methods with hybrid functionals. Spin-orbit coupling is treated by the atomic meanfield approximation. g-Tensors for a set of small main group radicals and for a series of ten 3d and two 4d transition metal complexes have been compared using the local density approximation (VWN functional), the generalized gradient approximation (BP86 functional), as well as B3-type (B3PW91) and BH-type (BHPW91) hybrid functionals. For main group radicals, the effect of exact-exchange mixing is small. In contrast, significant differences between the various functionals arise for transition metal complexes. As has been shown previously, local and in particular gradient-corrected functionals tend to underestimate the "paramagnetic" contributions to the g-tensors in these cases and thereby recover only about 40-50% of the range of experimental g-tensor components. This is improved to ca. 60% by the B3PW91 functional, which also gives slightly reduced standard deviations. The range increases to almost 100% using the half-and-half functional BHPW91. However, the quality of the correlation with experimental data worsens due to a significant overestimate of some intermediate g-tensor values. The worse performance of the BHPW91 functional in these cases is accompanied by spin contamination. Although none of the functionals tested thus appears to be ideal for the treatment of electronic g-tensors in transition metal complexes, the B3PW91 hybrid functional exhibited the overall most satisfactory performance. Apart from the validation of hybrid functionals, some aspects in the treatment of spin-orbit contributions to the g-tensor are discussed.