Density Functional Theory Calculations of the Electron Paramagnetic Resonance Parameters for VO2+ Complexes

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...

Integrated experimental/computational approaches to characterize the systems formed by vanadium with proteins and enzymes

An integrated instrumental/computational approach to characterize metallodrug–protein adducts at the molecular level is reviewed. A series of applications are described, focusing on potential vanadium drugs with a generalization to other metals.

Quantitative prediction of electronic absorption spectra of copper(II)-bioligand systems: Validation and applications

The visible region of the electronic absorption spectra of Cu(II) complexes was studied by time-dependent density functional theory (TD-DFT). The performance of twelve functionals in the prediction of absorption maxima (λ_max) was tested on eleven compounds with different geometry, donors and charge. The ranking of the functionals for λ_max was determined in terms of mean absolute percent deviation (MAPD) and standard deviation (SD) and it is as follows: BHandHLYP > M06 ≫ CAM-B3LYP ≫ MPW1PW91 ~ B1LYP ~ BLYP > HSE06 ~ B3LYP > B3P86 ~ ω-B97x-D ≫ TPSSh ≫ M06-2X (MAPD) and BHandHLYP > M06 ~ HSE06 > ω-B97x-D ~ CAM-B3LYP ~ MPW1PW91 > B1LYP ~ B3LYP > B3P86 > BLYP ≫ TPSSh ≫ M06-2X (SD). With BHandHLYP functional the MAPD is 3.1% and SD is 2.3%, while with M06 the MAPD is 3.7% and SD is 3.7%. The protocol validated in the first step of the study was applied to: i) calculate the number of transitions in the spectra and relate them to the geometry of Cu(II) species; ii) determine the coordination of axial water(s); iii) predict the electronic spectra of the systems where Cu(II) is bound to human serum albumin (HSA) and to the regions 94-97 and 108-112 of prion protein (PrP). The results indicate that the proposed computational protocol allows a successful prediction of the electronic spectra of Cu(II) species and to relate an experimental spectrum to a specific structure.

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.

Speciation in human blood of Metvan, a vanadium based potential anti-tumor drug

The first report on the anti-cancer activity of the compound Metvan, [V^IVO(Me₂phen)₂(SO₄)], where Me₂phen is 4,7-dimethyl-1,10-phenanthroline, dates back to 2001. Although it was immediately identified as one of the most promising multitargeted anti-cancer V compounds, no development on the medical experimentation was carried out. One of the possible reasons is the lack of information on its speciation in aqueous solution and its thermodynamic stability, factors which influence the transport in the blood and the final form which reaches the target organs. To fill this gap, in this work the speciation of Metvan in aqueous solution and human blood was studied by instrumental (EPR, electronic absorption spectroscopy, ESI-MS and ESI-MS/MS), analytical (pH-potentiometry) and computational (DFT) methods. The results suggested that Metvan transforms at physiological pH into the hydrolytic species cis-[VO(Me₂phen)₂(OH)]⁺ and that both citrate and proteins (transferrin and albumin in the blood serum, and hemoglobin in the erythrocytes) form mixed complexes, denoted [VO(Me₂phen)(citrH_-1)]^2- and VO-Me₂phen-Protein with the probable binding of His-N donors. The measurements with erythrocytes suggest that Metvan is able to cross their membrane forming mixed species VO-Me₂phen-Hb. The redox stability in cell culture medium was also examined, showing that ca. 60% is oxidized to V^V after 5 h. Overall, the speciation of Metvan in the blood mainly depends on the V concentration: when it is larger than 50 μM, [VO(Me₂phen)(citrH_-1)]^2- and VO-Me₂phen-Protein are the major species, while for concentrations lower than 10 μM, (VO)(hTf) is formed and Me₂phen is lost. Therefore, it is plausible that the pharmacological activity of Metvan could be due to the synergic action of free Me₂phen, and V^IVO and V^VO/V^VO₂ species.

Biorelevant reactions of the potential anti-tumor agent vanadocene dichloride

The interaction of the potential anti-tumor agent vanadocene dichloride ([Cp2VCl2] or VDC) with some relevant bioligands of the cytosol such as proteins (Hb), amino acids (glycine and histidine), NADH derivatives (NADH, NADPH, NAD(+) and NADP(+)), reductants (GSH and ascorbic acid), phosphates (HPO4(2-), P2O7(4-), cAMP, AMP, ADP and ATP) and carboxylate derivatives (lactate) and its uptake by red blood cells were studied. The results indicated that [Cp2VCl2] transforms at physiological pH into [Cp2V(OH)2] and that only HPO4(2-), P2O7(4-), lactate, ATP and ADP form mixed species with the [Cp2V](2+) moiety replacing the two hydroxide ions. EPR and electronic absorption spectroscopy, agarose gel electrophoresis and spin trapping measurements allow excluding any direct interaction and/or intercalation with DNA and the formation of reactive oxygen species (ROS) in Fenton-like reactions. Uptake experiments by erythrocytes suggested that VDC crosses the membrane and enters inside the cells, whereas 'bare' V(IV) transforms into V(IV)O species with loss of the two cyclopentadienyl rings. This transformation in the cellular environment could be related to the mechanism of action of VDC.

Nonoxido V(IV) Complexes: Prediction of the EPR Spectrum and Electronic Structure of Simple Coordination Compounds and Amavadin

Density functional theory (DFT) calculations of the (51)V hyperfine coupling (HFC) tensor A have been completed for 20 "bare" V(IV) complexes with different donor sets, electric charges, and coordination geometries. Calculations were performed with ORCA and Gaussian software, using functionals BP86, TPSS0, B1LYP, PBE0, B3LYP, B3P, B3PW, O3LYP, BHandHLYP, BHandH, and B2PLYP. Among the basis sets, 6-311g(d,p), 6-311++g(d,p), VTZ, cc-pVTZ, def2-TZVPP, and the "core properties" CP(PPP) were tested. The experimental Aiso and Ai (where i = x or z, depending on the geometry and electronic structure of V(IV) complex) were compared with the values calculated by DFT methods. The results indicated that, based on the mean absolute percentage deviation (MAPD), the best functional to predict Aiso or Ai is the double hybrid B2PLYP. With this functional and the basis set VTZ, it is possible to predict the Aiso and Az of the EPR spectrum of amavadin with deviations of -1.1% and -2.0% from the experimental values. The results allowed us to divide the spectra of nonoxido V(IV) compounds in three types-called "type 1", "type 2", and "type 3", characterized by different composition of the singly occupied molecular orbital (SOMO) and relationship between the values of Ax, Ay, and Az. For "type 1" spectra, Az ≫ Ax ≈ Ay and Az is in the range of (135-155) × 10(-4) cm(-1); for "type 2" spectra, Ax ≈ Ay ≫ Az and Ax ≈ Ay are in the range of (90-120) × 10(-4) cm(-1); and for the intermediate spectra of "type 3", Az > Ay > Ax or Ax > Ay > Az, with Az or Ax values in the range of (120-135) × 10(-4) cm(-1). The electronic structure of the V(IV) species was also discussed, and the results showed that the values of Ax or Az are correlated with the percent contribution of V-dxy orbital in the SOMO. Similarly to V(IV)O species, for amavadin the SOMO is based mainly on the V-dxy orbital, and this accounts for the large experimental value of Az (153 × 10(-4) cm(-1)).

Characterization and biotransformation in the plasma and red blood cells of V(IV)O(2+) complexes formed by ceftriaxone

The coordination mode and geometry in aqueous solution of oxidovanadium(IV) complexes formed by a third-generation cephalosporin, ceftriaxone (H3cef), were studied by spectroscopic (EPR, electron paramagnetic resonance), pH-potentiometric and computational (DFT, density functional theory) methods. The behavior of the model systems containing 6-hydroxy-2-methyl-3-thioxo-3,4-dihydro-1,2,4-triazine-5(2H)-one (H2hmtdt) and 3-benzylthio-6-hydroxy-2-methyl-1,2,4-triazine-5(2H)-one (Hbhmt) was examined for comparison. The stability of the tautomers of ceftriaxone and 6-hydroxy-2-methyl-3-thioxo-3,4-dihydro-1,2,4-triazine-5(2H)-one in the neutral, mono- and bi-anionic form was calculated by DFT methods, both in the gas phase and in aqueous solution, and the electron density on the oxygen atoms of the hydroxytriazinone ring was related to the pKa of the ligands. The data demonstrate that ceftriaxone coordinates V(IV)O(2+) forming mono- and bis-chelated complexes with (Oket, O(-)) donor set and formation of five-membered chelate rings. The geometry of the bis-chelated complex, cis-[VO(Hcef)2(H2O)](2-), is cis-octahedral and this species can deprotonate, around physiological pH, to form the corresponding mono-hydroxido cis-[VO(Hcef)2(OH)](3-). The interaction of cis-[VO(Hcef)2(H2O)](2-) with apo-transferrin (apo-hTf) was studied and the results suggest that V(IV)O(2+) distributes between (VO)apo-hTf/(VO)2apo-hTf and cis-[VO(Hcef)2(H2O)](2-), whereas mixed complexes are not formed for charge and steric effects. The interaction of cis-[VO(Hcef)2(H2O)](2-) with red blood cells shows that ceftriaxone helps V(IV)O(2+) ion to cross the erythrocyte membrane. Inside the cell cis-[VO(Hcef)2(H2O)](2-) decomposes and the same species formed by inorganic V(IV)O(2+) are observed. The relationship between the biotransformation in the plasma and red blood cells and the potential pharmacological activity of V(IV)O(2+) species of ceftriaxone is finally discussed.

Application of DFT methods to the study of the coordination environment of the VO2+ ion in V proteins

Density functional theory (DFT) methods were used to simulate the environment of vanadium in several V proteins, such as vanadyl-substituted carboxypeptidase (sites A and B), vanadyl-substituted chloroplast F(1)-ATPase (CF(1); site 3), the reduced inactive form of vanadium bromoperoxidase (VBrPO; low- and high-pH sites), and vanadyl-substituted imidazole glycerol phosphate dehydratase (IGPD; sites α, β, and γ). Structural, electron paramagnetic resonance, and electron spin echo envelope modulation parameters were calculated and compared with the experimental values. All the simulations were performed in water within the framework of the polarizable continuum model. The angular dependence of [Formula: see text] and [Formula: see text] on the dihedral angle θ between the V=O and N-C bonds and on the angle φ between the V=O and V-N bonds, where N is the coordinated aromatic nitrogen atom, was also found. From the results it emerges that it is possible to model the active site of a vanadium protein through DFT methods and determine its structure through the comparison between the calculated and experimental spectroscopic parameters. The calculations confirm that the donor sets of sites B and A of vanadyl-substituted carboxypeptidase are [[Formula: see text], H(2)O, H(2)O, H(2)O] and [N(His)(||), N(His)(⊥), [Formula: see text], H(2)O], and that the donor set of site 3 of CF(1)-ATPase is [[Formula: see text], OH(Thr), H(2)O, H(2)O, [Formula: see text]]. For VBrPO, the coordination modes [N(His)(||), N(His)(∠), OH(Ser), H(2)O, H(2)O(ax)] for the low-pH site and [N(His)(||), N(His)(∠), OH(Ser), OH(-), H(2)O(ax)] or [N(His)(||), N(His)(∠), [Formula: see text], H(2)O] for the high-pH site, with an imidazole ring of histidine strongly displaced from the equatorial plane, can be proposed. Finally, for sites α, β, and γ of IGPD, the subsequent deprotonation of one, two, and three imidazole rings of histidine and the participation of a carboxylate group of a glutamate residue ([N(His)(||), [Formula: see text], H(2)O, H(2)O], [N(His)(||), N(His)(||), [Formula: see text], H(2)O], and [N(His)(||), N(His)(||), [Formula: see text], OH(-), [Formula: see text]], respectively) seems to be the most plausible hypothesis.

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.

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.

Interaction of VO2+ ion and some insulin-enhancing compounds with immunoglobulin G

Complexation of VO(2+) ion with the most abundant class of human immunoglobulins, immunoglobulin G (IgG), was studied using EPR spectroscopy. Differently from the data in the literature which report no interaction of IgG with vanadium, in the binary system VO(2+)/IgG at least three sites with comparable strength were revealed. These sites, named 1, 2, and 3, seem to be not specific, and the most probable candidates for metal ion coordination are histidine-N, aspartate-O or glutamate-O, and serinate-O or threoninate-O. The mean value for the association constant of (VO)(x)IgG, with x = 3-4, is log β = 10.3 ± 1.0. Examination of the ternary systems formed by VO(2+) with IgG and human serum transferrin (hTf) and human serum albumin (HSA) allows one to find that the order of complexing strength is hTf ≫ HSA ≈ IgG. The behavior of the ternary systems with IgG and one insulin-enhancing agent, like [VO(6-mepic)(2)], cis-[VO(pic)(2)(H(2)O)], [VO(acac)(2)], and [VO(dhp)(2)], where 6-mepic, pic, acac, and dhp indicate the deprotonated forms of 6-methylpicolinic and picolinic acids, acetylacetone, and 1,2-dimethyl-3-hydroxy-4(1H)-pyridinone, is very similar to the corresponding systems with albumin. In particular, at the physiological pH value, VO(6-mepic)(IgG)(OH), cis-VO(pic)(2)(IgG), and cis-VO(dhp)(2)(IgG) are formed. In such species, IgG coordinates nonspecifically VO(2+) through an imidazole-N belonging to a histidine residue exposed on the protein surface. For cis-VO(dhp)(2)(IgG), log β is 25.6 ± 0.6, comparable with that of the analogous species cis-VO(dhp)(2)(HSA) and cis-VO(dhp)(2)(hTf). Finally, with these new values of log β, the predicted percent distribution of an insulin-enhancing VO(2+) agent between the high molecular mass (hTf, HSA, and IgG) and low molecular mass (lactate) components of the blood serum at physiological conditions is calculated.

The equilibrium between the octahedral and square pyramidal form and the influence of an axial ligand on the molecular properties of V(IV)O complexes: a spectroscopic and DFT study

The previously unreported equilibrium in aqueous solution between the V(IV)O square pyramidal and trans octahedral form with an axial water molecule for a number of bidentate ligands with (N,O) and (O,O) donor sets (6-methylpicolinic (6-mepicH) and 6-methyl-2,3-pyridinedicarboxylic (6-me-2,3-pdcH(2)) acids, dipyridin-2-ylmethanol (Hdpmo) and 1,2-dimethyl-3-hydroxy-4(1H)-pyridinone (Hdhp)) has been demonstrated by the combined application of EPR spectroscopy and DFT methods. The EPR spectra suggest that, with increasing ionic strength, the equilibrium is shifted towards the formation of the pentacoordinated species and values of K approximately 4.0 and 7.0 for the systems containing 6-methyl-2,3-pyridinedicarboxylic acid and dipyridin-2-ylmethanol were measured. DFT calculations performed with Gaussian 03 and ORCA software predict the (51)V anisotropic hyperfine coupling constant along the z axis (A(z)), which can be used to demonstrate the presence of an axially bound ligand trans to the V=O bond. The results suggest that an axial donor (charged or not) can lower |A(z)|, in contrast to what was previously believed on the basis of the "additivity rule", and this explains the anomalous behaviour of the V(IV)O complexes formed by N-{2-[(2-pyridylmethylene)amino]phenyl}pyridine-2-carboxamide (Hcapca) and several amidrazone derivatives. The decrease in |A(z)| for the axial binding of a solvent molecule is mainly a result of the reduction of |A(iso)| and this was also observed when the solid [VO(6-methylpicolinato)(2)] was dissolved in DMSO or DMF. The variations in the structural (V=O, V-O and V-N distances, O-V-O and N-V-N angles, and the trigonality index tau) and spectroscopic (|A(z)|, |A(iso)| and nu(V=O)) properties as a function of the axial V-OH(2) distance (R) are also presented. Finally, the electronic structures of the penta- and hexacoordinated complexes are discussed.

Bis- and tris(pyridyl)amine-oxidovanadium complexes: characteristics and insulin-mimetic potential

Two novel vanadium complexes, [V(IV)O(bp-O)(HSO4)] (1) and [V(IV)O(bp-OH)Cl2] x CH3OH (2 x CH3OH), where bp-OH is 2-{[bis(pyrid-2-yl)methyl]amine}methylphenol, were prepared and structurally characterised. EPR spectra of methanol solutions of 2 suggest exchange of Cl- for CH3OH and partial conversion to [VO(bp-OH)(CH3OH)3]2+. Speciation studies on the VO2+-bpOH system in a water/dmso mixture (4:1 v/v) revealed [VO(bp-O)(H2O)n]+ as the dominating species in the pH range 2-7. The insulin-mimetic properties of 1 and 2, [V(IV)O(SO4)tpa] (3), [V(IV)O(pic-trpMe)2] (5) and the new mixed-ligand complexes [V(V)O(pic-trpH)tpa]Cl2 (4Cl2) and [V(V)O(pic-OEt)tpa]Cl2 (6Cl2), tpa = tris(pyrid-2-yl)methylamine, picH-trpH = 2-carboxypyridine-5-(L-tryptophan)carboxamide (picH-trpMe is the respective tryptophanmethyl ester), pic-OEt = 5-carboethoxypyridine-2-carboxylic acid, were evaluated with rat adipocytes, employing two lipolysis assays (release of glycerol and free fatty acids (FFA)), respectively and a lipogenesis assay (incorporation of glucose into lipids). The IC50 values for the inhibition of lipolysis in the FFA assay vary between 0.41 (+/-0.03) (5) and 21.2 (+/-0.6) mM (2), as compared to 0.81 (+/-0.2) mM for VOSO4.

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.

Determining the Geometry and Magnetic Parameters of Fluorinated Radicals by Simulation of Powder ESR Spectra and DFT Calculations: The Case of the Radical RCF2CF2• in Nafion Perfluorinated Ionomers

The ESR spectrum of the chain-end radical RCF2CF2* detected in Nafion perfluorinated membranes exposed to the photo-Fenton reagent was accurately simulated by an automatic fitting procedure, using as input the hyperfine coupling tensors of the two F alpha and two F beta nuclei as well as the corresponding directions of the principal values from density functional theory (DFT) calculations. An accurate fit was obtained only for different orientations of the hyperfine coupling tensors for the two F alpha nuclei, indicating a nonplanar structure about the C alpha radical center. The fitted isotropic hyperfine splittings for the two F beta nuclei in the Nafion radical, 24.9 and 27.5 G, are significantly larger than those for the chain-end radical in Teflon (15 G), implying different radical conformations in the two systems. The excellent fit indicated that the geometry and electronic structure of free radicals can be obtained not only from single-crystal ESR spectroscopy, but also, in certain cases, from powder spectra, by combination with data from DFT calculations. The optimized structures obtained by DFT calculations for the CF3CF2CF2CF2* or CF3OCF2CF2* radicals as models provided additional support for the pyramidal structure determined from the spectral fit. Comparison and analysis of calculated and fitted values for the hyperfine splittings of the two F beta nuclei suggested that the radical detected by ESR in Nafion is ROCF2CF2*, which originates from attack of oxygen radicals on the Nafion side chain. The combination of spectrum fitting and DFT is considered important in terms of understanding the hyperfine splittings from 19F nuclei and the different conformations of fluorinated chain-end-type radicals RCF2CF2* in different systems, and also for elucidating the mechanism of Nafion fragmentation when exposed to oxygen radicals in fuel cell conditions.

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.