A Density Functional Study of EPR Parameters for Vanadyl Complexes Containing Schiff Base Ligands

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.

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.

Transport of the anti-diabetic VO2+ complexes formed by pyrone derivatives in the blood serum

The biotransformation in the blood serum of the two anti-diabetic agents [VO(ema)(2)] - or BEOV - and [VO(koj)(2)] formed by ethylmaltol (Hema) and kojic acid (Hkoj) was studied with EPR spectroscopy, pH-potentiometry and DFT calculations. For comparison, the behavior of the systems with tropolone (Htrop) was also analyzed. The interaction of [VO(ema)(2)] and [VO(koj)(2)] with the most important bioligands of the serum, lactic (Hlact) and citric acid (H(3)citr), human serum transferrin (hTf), human serum albumin (HSA) and immunoglobulin G (IgG) was examined and discussed. Among the several mixed species observed, cis-VO(carrier)(2)(hTf), cis-VO(carrier)(2)(HSA) and cis-VO(carrier)(2)(IgG), where carrier is ethylmaltolate or kojate, with a His-N of the protein coordinated in the equatorial position, are plausible candidates for the transport processes of the drug toward the target organs. The values of the logβ are in the range 19.6-19.8 for the species formed by ethylmaltol and 17.4-17.6 for those formed by kojic acid. The formation of such species was confirmed through pH-titrations of the model systems VO(2+)/carrier/1-MeIm and VO(2+)/carrier/Ac-his, where 1-MeIm and Ac-his are 1-methylimidazole and N-acetylhistamine, and DFT calculations of (51)V A(z) of the model species cis-[VO(carrier)(2)(1-MeIm)] and cis-[VO(carrier)(2)(Ac-his)]. The values of the stability constants for the mixed species observed were used to predict the biodistribution of VO(2+) ion between the blood serum components for concentrations of 1, 10 and 50 μM.

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.

Density functional study of EPR parameters and spin-density distribution of azurin and other blue copper proteins

Modern density functional methods have been used to study spin-density distribution, g tensors, as well as Cu and ligand hyperfine tensors for azurin models, for two more blue copper proteins plastocyanin and stellacyanin, and for small model complexes. The aim was to establish a consistent computational protocol that provides a realistic description of the EPR parameters as probes of the spin-density distribution between metal and coordinated ligands in copper proteins. In agreement with earlier conclusions for plastocyanin, hybrid functionals with appreciable exact-exchange admixtures, roughly around 50%, provide the best overall agreement with all parameters. Then the bulk of the spin density is almost equally shared by the copper atom and the sulfur atom of the equatorial cysteine ligand, and the best values are obtained for copper, histidine nitrogen, and cysteine beta-proton hyperfine couplings, as well as for g(parallel). Spin-orbit effects on the EPR parameters may be appreciable and have to be treated carefully to obtain agreement with experiment. Most notably, spin-orbit effects on the (65)Cu hyperfine coupling tensors in blue copper sites are unusually large compared to more regularly coordinated Cu(II) complexes with similar spin density on copper. In addition to the often emphasized high covalency of the Cu-S(Cys) bond, the characteristically small A(parallel) component of blue copper proteins is shown to derive to a large part from a near-cancellation between negative first-order (Fermi contact and dipolar) and unusually large positive second-order (spin-orbital) contributions. The large spin-orbit effects relate to the distorted tetrahedral structures. Square planar dithiolene complexes with similar spin density on copper exhibit much more negative A(parallel) values, as the cancellation between nonrelativistic and spin-orbit contributions is less complete. Calculations on a selenocysteine-substituted variant of azurin have provided further insight into the relations between bonding and EPR parameters.

Comparative Density-Functional Study of the Electron Paramagnetic Resonance Parameters of Amavadin

Electronic g tensors and hyperfine coupling tensors have been calculated for amavadin, an unusual eight-coordinate vanadium(IV) complex isolated from Amanita muscaria mushrooms. Different density-functional methods have been compared, ranging from local via gradient-corrected to hybrid functionals with a variable Hartree-Fock exchange admixture. For both electron paramagnetic resonance (EPR) properties, hybrid functionals with an appreciable exact-exchange admixture provide the closest agreement with experimental data. Second-order spin-orbit corrections provide non-negligible contributions to the 51V hyperfine tensor. The orientation of g and A tensors relative to each other also depends on spin-orbit corrections to the A tensor. A rationalization for the close resemblance of the EPR parameters of amavadin to those of the structurally rather different vanadyl complexes is provided, based on the nature of the relevant frontier orbitals.

Carboxylate Binding in Copper Histidine Complexes in Solution and in Zeolite Y: X- and W-band Pulsed EPR/ENDOR Combined with DFT Calculations

The complexes of copper with histidine exhibit a wide variety of coordination modes in aqueous solution. This stems from the three potential coordination sites of the histidine molecule and the existence of mono- and bis-complexes. The present work concentrates on the determination of the carboxylate binding mode, via the (13)C hyperfine coupling of the carboxyl, in a number of copper complexes in frozen solutions. These are then used as references for the determination of the coordination mode of two zeolite encapsulated complexes. The (13)C hyperfine coupling (sign and magnitude) was determined by a variety of advanced pulsed EPR and electron-nuclear double resonance (ENDOR) techniques carried out at conventional and high magnetic fields. These showed that while the carboxyl (13)C isotropic hyperfine coupling of an equatorially coordinated carboxylate is negative with a magnitude of 3-4 MHz, that of a free carboxylate is small ( approximately 1 MHz) and positive. To rationalize the experimentally determined ligand hyperfine couplings ((1)H and (13)C) and further understand their dependence on the coordination mode and degree of protonation, density functional theory (DFT) calculations were carried out on a number of model complexes, representing the different Cu-histidine complexes studied experimentally. The exchange-correlation functional used for the calculation of the EPR parameters was B3LYP with triple-zeta plus polarization (TZP) quality basis sets. While the polarization agreement between the magnitudes of the calculated and experimental values varied among the various nuclei, sometimes exhibiting deviations of up to 40%, an excellent agreement was found for the sign prediction. This shows the unique advantage of combining high field ENDOR, by which the sign of the hyperfine can often be determined, with DFT predictions for structure determination.