51V solid-state NMR and density functional theory studies of vanadium environments in V(V)O2 dipicolinic acid complexes

(51)V solid-state NMR and density functional theory (DFT) investigations are reported for a series of pentacoordinate dioxovanadium(V)-dipicolinate [V(V)O(2)-dipicolinate] and heptacoordinate aquahydroxylamidooxovanadium(V)-dipicolinate [V(V)O-dipicolinate] complexes. These compounds are of interest because of their potency as phosphatase inhibitors as well as their insulin enhancing properties and potential for the treatment of diabetes. Experimental solid-state NMR results show that the electric field gradient tensors in the V(V)O(2)-dipicolinate derivatives are affected significantly by substitution on the dipicolinate ring and range from 5.8 to 8.3 MHz. The chemical shift anisotropies show less dramatic variations with respect to the ligand changes and range between -550 and -600 ppm. To gain insights on the origins of the NMR parameters, DFT calculations were conducted for an extensive series of the V(V)O(2)- and V(V)O-dipicolinate complexes. To assess the level of theory required for the accurate calculation of the (51)V NMR parameters, different functionals, basis sets, and structural models were explored in the DFT study. It is shown that the original x-ray crystallographic geometries, including all counterions and solvation water molecules within 5 A of the vanadium, lead to the most accurate results. The choice of the functional and the basis set at a high level of theory has a relatively minor impact on the outcome of the chemical shift anisotropy calculations; however, the use of large basis sets is necessary for accurate calculations of the quadrupole coupling constants for several compounds of the V(V)O(2) series. These studies demonstrate that even though the vanadium compounds under investigations exhibit distorted trigonal bipyramidal coordination geometry, they have a "perfect" trigonal bipyramidal electronic environment. This observation could potentially explain why vanadate and vanadium(V) adducts are often recognized as potent transition state analogs.     

Investigating the vanadium environments in hydroxylamido V(V) dipicolinate complexes using 51V NMR spectroscopy and density functional theory

Using (51)V magic angle spinning solid-state NMR, SSNMR, spectroscopy and quantum chemical DFT calculations we have characterized the chemical shift and quadrupolar coupling parameters of a series of eight hydroxylamido vanadium(V) dipicolinate complexes of the general formula VO(dipic)(ONR1R2)(H2O) where R1 and R2 can be H, CH3, or CH2CH3. This class of vanadium compounds was chosen for investigation because of their seven-coordinate vanadium atom, a geometry for which there is limited (51)V SSNMR data. Furthermore, a systematic series of compounds with different electronic properties are available and allows for the effects of ligand substitution on the NMR parameters to be studied. The quadrupolar coupling constants, C(Q), are small, 3.0-3.9 MHz, but exhibit variations as a function of the ligand substitution. The chemical shift tensors in the solid state are sensitive to changes in both the hydroxylamide substituent and the dipic ligand, a sensitivity which is not observed for isotropic chemical shifts in solution. The chemical shift tensors span approximately 1000 ppm and are nearly axially symmetric. On the basis of DFT calculations of the chemical shift tensors, one of the largest contributors to the magnetic shielding anisotropy is an occupied molecular orbital with significant vanadium d(z)2 character along the V=O bond.     

103Rh NMR chemical shifts in organometallic complexes: a combined experimental and density functional study

Experimental 103Rh NMR chemical shifts of mono- and binuclear rhodium(I) complexes containing s- or as-hydroindacenide and indacenediide bridging ligands with different ancillary ligands (1,5-cyclooctadiene, ethylene, carbonyl) are presented. A protocol, based on density functional theory calculations, was established to determine 103Rh NMR shielding constants in order to rationalise the effects of electronic and structural variations on the spectroscopic signal, and to gain insight into the efficiency of this computational method when applied to organometallic systems. Scalar and spin-orbit relativistic effects based on the ZORA (zeroth order regular approximation) level have been taken into account and discussed. A good agreement was found for model compounds over a wide range of chemical shifts of rhodium (approximately 10,000 ppm). This allowed us to discuss the experimental and calculated delta(103Rh) in larger complexes and to relate it to their electronic structure.     

Peroxovanadate imidazole complexes as catalysts for olefin epoxidation: density functional study of dynamics, 51V NMR chemical shifts, and mechanism

A density functional study of [VO(O(2))(2)(Im)](-) (1, Im = imidazole) is presented, calling special attention to effects of dynamics and solvation on the (51)V chemical shift. According to Car-Parrinello molecular dynamics simulations, rotation of the Im ligand can be fast in the gas phase, but is more hindered in aqueous solution. In the latter, bonding between Im and V is reinforced, and dynamic averaging of GIAO-B3LYP magnetic shieldings affords a gas-to-liquid shift of ca. -100 ppm for delta((51)V). A complete catalytic cycle has been characterized for olefin epoxidation mediated by 1, using H(2)O(2) as oxidant. The rate-determining step is indicated to be initial oxygen atom transfer from 1 to the substrate via a spiro-like transition state. Substituent effects on this barrier are examined, and a significant decrease (by 2-6 kcal/mol) is revealed upon removal of the Im proton or upon complexation with a H-bond acceptor. Implications for the mechanism of the oxidative chemistry of vanadium-dependent haloperoxidases and requirements for prospective biomimetic analogues are discussed.     

Fragment density functional theory calculation of NMR chemical shifts for proteins with implicit solvation

Fragment density functional theory (DFT) calculation of NMR chemical shifts for several proteins (Trp-cage, Pin1 WW domain, the third IgG-binding domain of Protein G (GB3) and human ubiquitin) has been carried out. The present study is based on a recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach but the solvent effects are included by using the PB (Poisson-Boltzmann) model. Our calculated chemical shifts of (1)H and (13)C for these four proteins are in excellent agreement with experimentally measured values and represent clear improvement over that from the gas phase calculation. However, although the inclusion of the solvent effect also improves the computed chemical shifts of (15)N, the results do not agree with experimental values as well as (1)H and (13)C. Our study also demonstrates that AF-QM/MM calculated results accurately reproduce the separation of α-helical and β-sheet chemical shifts for (13)C(α) atoms in proteins, and using the (1)H chemical shift to discriminate the native structure of proteins from decoys is quite remarkable.     

Density functional theory calculation of (29)Si NMR chemical shifts of organosiloxanes

Density functional theory at the B3LYP/6-311++G(d,p) level is applied to calculate the (29)Si NMR chemical shifts of a variety of organosiloxane moieties including monomers or precursors for polymerization and representative segments of organosiloxane polymers or thin films. The calculated shifts of two linear dimethylsiloxane compounds, hexamethylcyclotrisiloxane (D3) and octamethylcyclotetrasiloxane (D4), compare well with their known values, having an average error of 3.4 ppm. The same method is applied to structures believed to occur in organosilicate glass thin films deposited using hot-filament chemical vapor deposition (HFCVD) from D3 and D4. The chemical shift at -15 ppm is identified as a cross-linking Si-Si bond between two strained D groups and has not previously been reported. Retention of the strained ringed structure in HFCVD films deposited from D3 is confirmed. The rings are bonded to the matrix through either Si-O or Si-Si bonds, with the latter only becoming prevalent when higher filament temperatures are employed. The strained ring structure is also observed in films deposited from a precursor with a larger unstrained ring structure, D4. These observations suggest that the known gas-phase conversion pathways of D4 to D3 and dimethylsilanone as well as the methyl abstraction reaction from D3 operate in the HFCVD reaction chemistry.     

51V NMR chemical shifts from quantum-mechanical/molecular-mechanical models of vanadium bromoperoxidase

According to quantum-mechanical/molecular-mechanical (QM/MM) optimizations, the active-site geometries of vanadium-dependent bromoperoxidase (VBPO) and vanadium-dependent chloroperoxidase (VCPO) are very similar. 51V NMR chemical shifts calculated from QM/MM-optimized models of VBPO are critically compared to VCPO and are found to be very similar for the two related proteins. The primary difference between these related structures, the presence of a His411 in VBPO whereas Phe397 is located at that position in VCPO, is studied via analysis of the respective theoretical 51V NMR spectra. The long-range electrostatic effects from more distal residues are also studied to establish their effect. Similar results are obtained for the two active sites of the VBPO homodimer. The experimentally observed shielding of the isotropic 51V NMR chemical shift on going from VCPO to VBPO is somewhat underestimated in the QM/MM models studied. NMR and NQC tensors of both enzymes are predicted to show noticeable differences, suggesting that precise solid-state 51V NMR data, when they become available, can be a sensitive probe for subtle differences in structural details between these enzymes.     

Structure and bonding of vanadium(V) complexes with hydroxyurea in aqueous solution: density functional theory investigations of isomers and intramolecular rearrangements

Extensive geometry optimizations have been performed at the BP86 level of density functional theory, in order to identify the most stable isomer of pentacoordinated [VO(OH)UH2O]+ and [VOU(H2O)2]2+ as well as of hexacoordinated [VO(OH)U(H2O)2]+ and [VOU(H2O)3]2+ complexes (U = hydroxyurea anion). Most of these are conformationally very flexible, with up to 12 isomers within an energy range of 5 kcal/mol. The most stable hexacoordinate forms are characterized by the oxo ligand in trans position to the carbonyl O atom of U. Bulk solvent effects on the relative stabilities, estimated from a polarizable continuum model, are indicated to be small and do not affect the energetic sequence of the isomers significantly. Details of the coordination sphere of the most stable isomers in aqueous solution (coordination number, protonation state) have been studied with Car-Parrinello molecular dynamics simulations. The preferred mechanisms of interconversion between selected [VO(OH)U(H2O)2]+ isomers, according to the DFT computations, involve proton transfers between H2O and OH or between O and OH ligands in the coordination sphere of the metal, assisted by a water molecule from the second hydration sphere.     

C3 vanadium(V) amine triphenolate complexes: vanadium haloperoxidase structural and functional models

The C 3 vanadium(V) amine triphenolate complex 1f has been characterized as a structural and functional model of vanadium haloperoxidases. The complex catalyzes efficiently sulfoxidations at room temperature using hydrogen peroxide as the terminal oxidant, yielding the corresponding sulfoxides in quantitative yields and high selectivities (catalyst loading down to 0.01%, TONs up to 9900, and TOFs up to 8000 h (-1)) as well as bromination of 1,3,5-trimethoxybenzene (catalyst loading down to 0.05%, TONs up to 1260, and TOFs up to 220 h (-1)).     

A simple graphical approach to predict local residue conformation using NMR chemical shifts and density functional theory

The dependency of amino acid chemical shifts on φ and ψ torsion angle is, independently, studied using a five‐residue fragment of ubiquitin and ONIOM(DFT:HF) approach. The variation of absolute deviation of ¹³C^α chemical shifts relative to φ dihedral angle is specifically dependent on secondary structure of protein not on amino acid type and fragment sequence. This dependency is observed neither on any of ¹³C^β, and ¹H^α chemical shifts nor on the variation of absolute deviation of ¹³C^α chemical shifts relative to ψ dihedral angle. The ¹³C^α absolute deviation chemical shifts (ADCC) plots are found as a suitable and simple tool to predict secondary structure of protein with no requirement of highly accurate calculations, priori knowledge of protein structure and structural refinement. Comparison of Full‐DFT and ONIOM(DFT:HF) approaches illustrates that the trend of ¹³C^α ADCC plots are independent of computational method but not of basis set valence shell type. © 2016 Wiley Periodicals, Inc.     

51V NMR chemical shifts calculated from QM/MM models of vanadium chloroperoxidase

(51)V NMR chemical shifts calculated from QM/MM-optimized (QM=quantum mechanical; MM=molecular mechanical) models of vanadium-dependent chloroperoxidase (VCPO) are presented. An extensive number of protonation states for the vanadium cofactor (active site of the protein) and a number of probable positional isomers for each of the protonation states are considered. The size of the QM region is increased incrementally to observe the convergence behavior of the (51)V NMR chemical shifts. A total of 40 models are assessed by comparison to experimental solid-state (51)V NMR results recently reported in the literature. Isotropic chemical shifts are found to be a poor indicator of the protonation state; however, anisotropic chemical shifts and the nuclear quadrupole tensors appear to be sensitive to changes in the proton environment of the vanadium nuclei. This detailed investigation of the (51)V NMR chemical shifts computed from QM/MM models provides further evidence that the ground state is either a triply protonated (one axial water and one equatorial hydroxyl group) or a doubly protonated vanadate moiety in VCPO. Particular attention is given to the electrostatic and geometric effects of the protein environment. This is the first study to compute anisotropic NMR chemical shifts from QM/MM models of an active metalloprotein for direct comparison with solid-state MAS NMR data. This theoretical approach enhances the potential use of experimental solid-state NMR spectroscopy for the structural determination of metalloproteins.     

Regression formulas for density functional theory calculated 1H and 13C NMR chemical shifts in toluene-d8

This study aimed at investigating the performance of a series of basis sets, density functional theory (DFT) functionals, and the IEF-PCM solvation model in the accurate calculation of (1)H and (13)C NMR chemical shifts in toluene-d(8). We demonstrated that, on a test set of 37 organic species with various functional moieties, linear scaling significantly improved the calculated shifts and was necessary to obtain more accurate results. Inclusion of a solvation model produced larger deviations from the experimental data as compared to the gas-phase calculations. Moreover, we did not find any evidence that very large basis sets were necessary to reproduce the experimental NMR data. Ultimately, we recommend the use of the BMK functional. For the (1)H shifts the use of the 6-311G(d) basis set gave linearly scaled mean unsigned (MU) and root-mean-square (rms) errors of 0.15 ppm and 0.21 ppm, respectively. For the calculation of the (13)C chemical shifts the 6-31G(d) basis set produced MUE of 1.82 ppm and RMSE of 3.29 ppm.     

Fluorooxoperoxo complexes of vanadium(V)

Substances crystallizing under various conditions from the MVO₃(MF, HF)H₂O₂H₂O (M = NH₄, K) systems have been characterized by elemental analysis, infrared and Raman spectra and X-ray powder patterns. Besides the known M₂[VO(O₂)₂F] complexes, complexes of two new types have been obtained: M₃[HV₂O₂(O₂)₃F₄·2H₂O and (NH₄)₃[V₂O₂(O₂)₄F]·nH₂O (n≈2). Vibrational spectra of new complexes are consistent with the presence of dimeric anions containing V(μ₂O₂)V and VFV bridges, respectively.     

A density functional theory study of uranium(VI) nitrate monoamide complexes

Density functional theory calculations were performed on uranyl complexed with nitrate and monoamide ligands (L) [UO(2)(NO(3))(2)·2L]. The obtained results show that the complex stability is mainly governed by two factors: (i) the maximization of the polarizability of the coordinating ligand and (ii) the minimization of the steric hindrance effects. Furthermore, the electrostatic interaction between ligands and uranium(vi) was found to be a crucial parameter for the complex stability. These results pave the way to the definition of (quantitative) property/structure relationships for the in silico screening of monoamide ligands with improved extraction efficiency of uranium(vi) in nitrate acidic solution.     

Synthesis, structure, and spectroscopic properties of chiral oxorhenium(V) complexes incorporating polydentate ligands derived from L-amino acids: a density functional theory/time-dependent density functional theory investigation

The oxorhenium(V) complexes [Re (V)O(L A)Cl 2] bearing the (N-2-pyridylmethyl) of l-valine (HL A (1)), l-leucine (HL A (2)), and l-phenylalanine (HL A (3)) and [Re (V)O(L B)Cl] containing the {(N-2pyridylmethyl)-(N-(5-nitro-2-hydroxybenzyl)} of l-valine (H 2L B (1)), l-leucine (H 2L B (2)), and l-phenylalanine (H 2L B (3)) are presented in this article. The complexes are isolated in enantiomeric pure form examined from X-ray structure determination. The complexes are characterized by spectroscopic and electrochemical methods. The molecular structures observed in the solid state are grossly preserved in solution ( (1)H, (13)C, and circular dichroism spectra). Gas-phase geometry optimization and the electronic structures of [Re (V)O(L A (1))Cl 2], [Re (V)O(L A (2))Cl 2], and [Re (V)O(L B (2))Cl] have been investigated with the framework of density functional theory. The absorption and circular dichroism spectra of the complexes were also calculated applying time-dependent density functional theory (TDDFT) using the conductor-like polarizable continuum solvent model to understand the origin of the electronic excitations. The chemical shift ( (1)H and (13)C) as well as (1)H- (1)H spin-spin coupling constant were also computed by the gauge-independent atomic orbital method, and the computed values are consistent with the experimental data.     

A density functional study of the 13C NMR chemical shifts in functionalized single-walled carbon nanotubes

The 13C NMR chemical shifts for functionalized (7,0), (8,0), (9,0), and (10,0) single-walled carbon nanotubes (SWNTs) have been studied computationally using gauge-including projector-augmented plane-wave (GIPAW) density functional theory (DFT). The functional groups NH, NCH3, NCH2OH, and CH2NHCH2 have been considered, and different sites where covalent addition or substitution may occur have been examined. The shifts of the carbons directly attached to the group are sensitive to the bond which has been functionalized and may, therefore, be used to identify whether the group has reacted with a parallel or a diagonal C-C bond. The addition of NH to a parallel bond renders the functionalized carbons formally sp3-hybridized, yielding shifts of around 44 ppm, independent of the SWNT radius. Reaction with a diagonal bond retains the formal sp2 hybridization of the substituted carbons, and their shifts are slightly lower or higher than those of the unsubstituted carbon atoms. The calculated 1H NMR shifts of protons in the functional groups are also dependent upon the SWNT-group interaction. Upon decreasing the degree of functionalization for the systems where the group is added to a parallel bond, the average chemical shift of the unfunctionalized carbons approaches that of the pristine tube. At the same time, the shifts of the functionalized carbons remain independent upon the degree of functionalization. For the SWNTs where N-R attaches to a parallel bond, the average shift of the sp2 carbons was found to be insensitive to the substituent R. Moreover, the shifts of the functionalized sp3 carbons, as well as of the carbons within the group itself, are independent of the SWNT radius. The results indicate that a wealth of knowledge may be obtained from the 13C NMR of functionalized SWNTs.     

Ethylenediaminetetraacetato-monoperoxo complexes of vanadium(V)

Summary The monoperoxo complexes, M₂[VO(HEDTA)(O₂)]· 4H₂O, where M is K⁺ or NH ₄ ⁺ and H₄EDTA is ethylene-diaminetetraacetic acid, were prepared and characterized by Raman and i.r. spectra in the solid state and in aqueous solution. The single crystal X-ray study revealed a pentagonal bipyramidal anion structure with a tetradentate HEDTA(3—) ligand. The decomposition of complexes in aqueous solution to blue vanadium(IV) complexes as end products proceeds via a nonperoxo complex of vanadium.     

Dinuclear gold (I)-di-N-heterocyclic carbene complexes: syntheses,structure, density functional theory calculation and photoluminescence study

The synthesis and characterizations for a series of dinuclear gold (I)-di-NHC complexes, 1–8 through the trans-metalation method of their respective silver (I)-di-NHC complexes, i–viii are reported (where NHC = N-heterocyclic carbene). The successful complexation of a series of unusual non-symmetrical and symmetrical di-NHC ligands, 3,3'-(ethane-1,2-diyl)-1-alkylbenzimidazolium-1'-butylbenzimidazolium (with alkyl = methyl, ethyl, propyl, butyl, pentyl, hexyl, heptyl, benzyl) with the gold (I) ions are suggested by elemental analysis, Fourier transform-infrared, ¹H- and ¹³C-NMR data. The ¹³C-NMR spectra of 1–8 show a singlet sharp peak in the range of 190.00–192.00 ppm, indicating the presence of a carbene carbon that bonded to the gold (I) ion. From single crystal X-ray diffraction data, the structure of complex 6 with the formula of [di-NHC-Au (I)]₂·2PF₆ is obtained [where NHC = 3,3'-(ethane-1,2-diyl)-1-hexylbenzimidazolium-1'-butylbenzimidazolium]. The photophysical study in solid state of 6 displays an intense photoluminescence with a strong emission maxima, λ_em = 480 nm, upon excitation at 340 nm at room temperature. Interestingly, the emission maximum at 77 K shows a structural character with a strong peak at 410 nm, a medium at 433 nm and a weak at 387 nm, accompanied by a tail band to about 500 nm.     

Synthesis,spectroscopic, theoretical study,and biological activities of vanadium(IV) and vanadium(V) complexes with isonipecotic acid

Russian Journal of General Chemistry - Vanadium(IV) and vanadium(V) complexes 1–7 have been synthesized by the reaction of isonipecotic acid with VOSO4 · 3 H2O, VCl3(THF)3, and NH4VO3 at...