51V NMR Chemical Shifts Calculated from QM/MM Models of Vanadium Chloroperoxidase

Predicting 51V nuclear magnetic resonance observables in molecular crystals

Solid-state nuclear magnetic resonance (NMR) spectroscopy and quantum chemical density functional theory (DFT) calculations are widely used to characterize vanadium centers in biological and pharmaceutically relevant compounds. Several techniques have been recently developed to improve the accuracy of predicted NMR parameters obtained from DFT. Fragment-based and planewave-corrected methods employing hybrid density functionals are particularly effective tools for solid-state applications. A recent benchmark study involving molecular crystal compounds found that fragment-based NMR calculations using hybrid density functionals improve the accuracy of predicted 51V chemical shieldings by 20% relative to traditional planewave methods. This work extends the previous study, including a careful analysis of 51V chemical shift anisotropy, electric field gradient calculations, and a more extensive test set. The accuracy of planewave-corrected techniques and recently developed fragment-based methods using electrostatic embedding based on the polarized continuum model (PCM) are found to be highly competitive with previous methods. Planewave-corrected methods achieve a 34% improvement in the errors of predicted 51V chemical shieldings relative to planewave. Additionally, planewave-corrected and fragment-based calculations were performed using PCM embedding, improving the accuracy of predicted 51V chemical shielding (CS) tensor principal values by 30% andC q values by 15% relative to traditional planewave methods. The performance of these methods is further examined using a redox-active oxovandium complex and a common 51V NMR reference compound.

Accurate fragment-based 51-V chemical shift predictions in molecular crystals

Nuclear magnetic resonance (NMR) spectroscopy plays a crucial role in determining molecular structure for complex biological and pharmaceutical compounds. NMR investigations are increasingly reliant on computation for mapping spectral features to chemical structures. Here we benchmark the accuracy of fragment-based ⁵¹V chemical shielding tensor calculations using a training set comprised of 10 biologically and pharmaceutically relevant oxovanadium complexes. Using our self-consistent reproduction of the Madelung potential (SCRMP) electrostatic embedding model, we demonstrate comparable performance between fragment methods and computationally demanding cluster-based techniques. Specifically, fragment methods employing hybrid density functionals are capable of reproducing the experimental ⁵¹V isotropic chemical shifts with a training set rms error of ~9 ​ppm, representing a 20% improvement over traditional plane wave techniques. We provide training set-derived linear regression models for mapping the absolute shieldings obtained from computation to the experimentally determined chemical shifts using four common density functionals; PBE0, B3LYP, PBE, and BLYP. Finally, we establish the utility of fragment methods and the reported regression parameters examining four oxovanadium structures excluded from the training set including the tetracoordinate oxovanadium silicate [Formula: see text] , VO¹⁵NGlySalbz which contains redox-active ligands, and the solid-state form of the common ⁵¹V NMR reference compound VOCl₃.

Calculation of higher protonation states and of a new resting state for vanadium chloroperoxidase using QM/MM, with an Atom-in-Molecules analysis

Earlier QM/MM studies of the resting state of vanadium chloroperoxidase (VCPO) focused on the diprotonated states of the vanadate cofactor. Herein, we report a new extensive QM/MM study that includes the tri- and quadprotonated states of VCPO at neutral pH. We identify certain di- and triprotonated states as being candidates for the resting state based on a comparison of relative energies. The quadprotonated states as well as some of the triprotonated states are ruled out as the resting state. An Atoms-in-Molecules (AIM) analysis of the complex hydrogen bonding around the vanadate cofactor helps to explain the relative energies of the protonation states considered herein, and it also indicates new hydrogen bonding which has not been recognized previously. A Natural Bond Orbital (NBO) study is presented to give a better understanding of the electronic structure of the vanadate co-factor.

Metal-Metal (MM) Bond Distances and Bond Orders in Binuclear Metal Complexes of the First Row Transition Metals Titanium Through Zinc

This survey of metal-metal (MM) bond distances in binuclear complexes of the first row 3d-block elements reviews experimental and computational research on a wide range of such systems. The metals surveyed are titanium, vanadium, chromium, manganese, iron, cobalt, nickel, copper, and zinc, representing the only comprehensive presentation of such results to date. Factors impacting MM bond lengths that are discussed here include (a) the formal MM bond order, (b) size of the metal ion present in the bimetallic core (M₂) ⁿ⁺, (c) the metal oxidation state, (d) effects of ligand basicity, coordination mode and number, and (e) steric effects of bulky ligands. Correlations between experimental and computational findings are examined wherever possible, often yielding good agreement for MM bond lengths. The formal bond order provides a key basis for assessing experimental and computationally derived MM bond lengths. The effects of change in the metal upon MM bond length ranges in binuclear complexes suggest trends for single, double, triple, and quadruple MM bonds which are related to the available information on metal atomic radii. It emerges that while specific factors for a limited range of complexes are found to have their expected impact in many cases, the assessment of the net effect of these factors is challenging. The combination of experimental and computational results leads us to propose for the first time the ranges and "best" estimates for MM bond distances of all types (Ti-Ti through Zn-Zn, single through quintuple).

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

Unsaturated trinuclear iron fluoroborylene complexes

The unsaturated trinuclear iron fluoroborylene complexes Fe₃(BF)₃(CO) _n (n = 7, 6) have been studied using density functional theory (DFT). Relatively complicated potential energy surfaces are found with nine and eight structures within 15 kcal mol^-1 of the lowest energy structures for the Fe₃(BF)₃(CO)₇ and Fe₃(BF)₃(CO)₆ systems, respectively. In each of these low-energy structures all three BF groups are either edge-bridging or face-bridging but never terminal groups. Some, but not all, of the low-energy structures also have edge-bridging and/or face-bridging CO groups leading to some structures with as many as five bridging groups. The relatively narrow range of Fe-Fe distances in the central Fe₃ triangles of the Fe₃(BF)₃(CO) _n (n = 7, 6) structures, mainly between 2.37 and 2.55 Å, suggests considerable delocalization in these unsaturated systems. Graphical Abstract The lowest energy Fe₃(BF)₃(CO)₇ and Fe₃(BF)₃(CO)₆ structures have a face-bridging μ₃-BF group with the two remaining BF groups bridging edges. The lowest energy Fe₃(BF)₃(CO)₆ structure also has one four-electron donor bridging η²-μ-CO group.

51V magic angle spinning NMR spectroscopy and quantum chemical calculations in vanadium bio-inorganic systems: current perspective

In recent years, 51V magic angle spinning (MAS) NMR spectroscopy has been widely used to characterize vanadium centers in biology, biomimetic complexes, and inorganic compounds of medicinal and industrial relevance. It has been demonstrated that 51V NMR parameters are sensitive probes of the coordination geometry and chemical environment of the metal center, beyond the first coordination sphere. To establish the relationships between NMR parameters and structure of the vanadium centers, over the past decade a large series of coordination complexes have been analyzed by MAS NMR spectroscopy. It has been demonstrated that the interpretation of the NMR parameters requires the use of theoretical methods, such as density functional (DFT) theory, whereby the experimental NMR observables are linked to the electronic and structural properties of a molecule. DFT calculations have been successfully employed to not only predict NMR parameters but to also yield valuable information regarding the structure and function of various vanadium compounds. In this report, we review the current state of the field, and present a survey of bioinorganic vanadium complexes as well as vanadium-dependent haloperoxidases analyzed using 51V MAS NMR spectroscopy and DFT calculations, to illustrate the rich information content available from such a combined approach.

Enzymatic Halogenases and Haloperoxidases: Computational Studies on Mechanism and Function

Despite the fact that halogenated compounds are rare in biology, a number of organisms have developed processes to utilize halogens and in recent years, a string of enzymes have been identified that selectively insert halogen atoms into, for instance, a CH aliphatic bond. Thus, a number of natural products, including antibiotics, contain halogenated functional groups. This unusual process has great relevance to the chemical industry for stereoselective and regiospecific synthesis of haloalkanes. Currently, however, industry utilizes few applications of biological haloperoxidases and halogenases, but efforts are being worked on to understand their catalytic mechanism, so that their catalytic function can be upscaled. In this review, we summarize experimental and computational studies on the catalytic mechanism of a range of haloperoxidases and halogenases with structurally very different catalytic features and cofactors. This chapter gives an overview of heme-dependent haloperoxidases, nonheme vanadium-dependent haloperoxidases, and flavin adenine dinucleotide-dependent haloperoxidases. In addition, we discuss the S-adenosyl-l-methionine fluoridase and nonheme iron/α-ketoglutarate-dependent halogenases. In particular, computational efforts have been applied extensively for several of these haloperoxidases and halogenases and have given insight into the essential structural features that enable these enzymes to perform the unusual halogen atom transfer to substrates.

Insights into enzymatic halogenation from computational studies

The halogenases are a group of enzymes that have only come to the fore over the last 10 years thanks to the discovery and characterization of several novel representatives. They have revealed the fascinating variety of distinct chemical mechanisms that nature utilizes to activate halogens and introduce them into organic substrates. Computational studies using a range of approaches have already elucidated many details of the mechanisms of these enzymes, often in synergistic combination with experiment. This Review summarizes the main insights gained from these studies. It also seeks to identify open questions that are amenable to computational investigations. The studies discussed herein serve to illustrate some of the limitations of the current computational approaches and the challenges encountered in computational mechanistic enzymology.

Isomerization mechanism of the HcRed fluorescent protein chromophore

To understand how the protein achieves fluorescence, the isomerization mechanism of the HcRed chromophore is studied both under vacuum and in the solvated red fluorescent protein. Quantum mechanical (QM) and quantum mechanical/molecular mechanical (QM/MM) methods are applied both for the ground and the first excited state. The photoinduced processes in the chromophore mainly involve torsions around the imidazolinone-bridge bond (τ) and the phenoxy-bridge bond (φ). Under vacuum, the isomerization of the cis-trans chromophore essentially proceeds by τ twisting, while the radiationless decay requires φ torsion. By contrast, the isomerization of the cis-trans chromophore in HcRed occurs via simultaneous τ and φ twisting. The protein environment significantly reduces the barrier of this hula twist motion compared with vacuum. The excited-state isomerization barrier via the φ rotation of the cis-coplanar conformer in HcRed is computed to be significantly higher than that of the trans-non-coplanar conformer. This is consistent with the experimental observation that the cis-coplanar-conformation of the chromophore is related to the fluorescent properties of HcRed, while the trans-non-planar conformation is weakly fluorescent or non-fluorescent. Our study shows how the protein modifies the isomerization mechanism, notably by interactions involving the nearby residue Ile197, which keeps the chromophore coplanar and blocks the twisting motion that leads to photoinduced radiationless decay.

Platinated DNA affects zinc finger conformation. Interaction of a platinated single-stranded oligonucleotide and the C-terminal zinc finger of nucleocapsid protein HIVNCp7

This paper describes for the first time the intimate molecular details of the association between a platinated oligonucleotide and a zinc finger peptide. Site-specific platination of the guanine in a single-stranded hexanucleotide gave {[Pt(dien)d(5'-TACGCC-3')], Pt(dien)(6-mer)} (II) characterized by mass spectrometry and (1)H nuclear magnetic resonance (NMR) spectroscopy. The work extends the study of platinum-nucleobase complex-zinc finger interactions using small molecules such as [Pt(dien)(9-EtGua)](2+) (I). The structure of the (34-52) C-terminal finger of HIV nucleocapsid protein HIVNCp7 (ZF1) was characterized by (1)H NMR spectroscopy and compared with that of the N-terminal single finger and the two-finger "intact" NCp7. Interaction of II with ZF1 results in significant changes in comparison to the "free" uncomplexed hexanucleotide; the major changes occurring for Trp37 resonances that are broadened and moved upfield, and other major shifts are for Gln45 (Hε21, Hγ3, Qβ), Met46 (NH, Hγ2), Lys47 (NH, Qγ), and Glu50 (Hγ2, Hγ3). The Zn-Cys/His chemical shifts show only marginal deviations. The solution structures of ZF1 and the 6-mer-ZF1 and II-ZF1 adducts were calculated from the nuclear Overhauser effect spectroscopy-derived distance constraints. The DNA position in the II-ZF1 adduct is completely different than in the absence of platinum. Major differences are the appearance of new Met46-Cyt6 H5 and Trp37-Cyt5 H5 contacts but severe weakening of the Trp37-Gua4 contact, attributed to the steric effects caused by Gua4 platination, accompanied by a change in the position of the aromatic ring. The results demonstrate the feasibility of targeting specific ZF motifs with DNA-tethered coordination compounds, such as Pt compounds and Co macrocycles, with implications for drug targetting and indeed the intimate mechanisms of DNA repair of platinated DNA.

Thermal and solvent effects on NMR indirect spin-spin coupling constants of a prototypical Chagas disease drug

NMR J-couplings across hydrogen bonds reflect the static and dynamic character of hydrogen bonding. They are affected by thermal and solvent effects and can therefore be used to probe such effects. We have applied density functional theory (DFT) to compute the NMR (n)J(N,H) scalar couplings of a prototypical Chagas disease drug (metronidazole). The calculations were done for the molecule in vacuo, in microsolvated cluster models with one or few water molecules, in snapshots obtained from molecular dynamics simulations with explicit water solvent, and in a polarizable dielectric continuum. Hyperconjugative and electrostatic effects on spin-spin coupling constants were assessed through DFT calculations using natural bond orbital (NBO) analysis and atoms in molecules (AIM) theory. In the calculations with explicit solvent molecules, special attention was given to the nature of the hydrogen bonds formed with the solvent molecules. The results highlight the importance of properly incorporating thermal and solvent effects into NMR calculations in the condensed phase.

A simple method for analyzing 51V solid-state NMR spectra of complex systems

Five vanadium complexes as models for biological systems were investigated using (51)V-MAS-NMR spectroscopy. All spectra show an uncommon line shape, which can be attributed to a shorter relaxation time of the satellite transition in contrast to the central one. A method for the reliable analysis of such kind of spectra is presented for the first time and the most important NMR parameters of the investigated complexes (quadrupolar coupling constant C(Q), asymmetry of the EFG tensor η(Q), isotropic chemical shift δ(iso), chemical shift anisotropy δ(σ) and asymmetry of the CSA tensor η(σ)) are presented. These results are of particular importance with respect to the analysis of the (51)V-MAS-NMR spectra of vanadium moieties in biological matrices such as vanadium chloroperoxidase, which show hitherto unexplained low intensity of the satellite sideband pattern.

Diverse roles of hydrogen in rhenium carbonyl chemistry: hydrides, dihydrogen complexes, and a formyl derivative

Rhenium carbonyl hydride chemistry dates back to the 1959 synthesis of HRe(CO)₅ by Hieber and Braun. The binuclear H₂Re₂(CO)₈ was subsequently synthesized as a stable compound with a central Re₂(μ-H)₂ unit analogous to the B₂(μ-H)₂ unit in diborane. The complete series of HRe(CO)(n) (n = 5, 4, 3) and H₂Re₂(CO)(n) (n = 9, 8, 7, 6) derivatives have now been investigated by density functional theory. In contrast to the corresponding manganese derivatives, all of the triplet rhenium structures are found to lie at relatively high energies compared with the corresponding singlet structures consistent with the higher ligand field splitting of rhenium relative to manganese. The lowest energy HRe(CO)₅ structure is the expected octahedral structure. Low-energy structures for HRe(CO)(n) (n = 4, 3) are singlet structures derived from the octahedral HRe(CO)₅ structure by removal of one or two carbonyl groups. For H₂Re₂(CO)₉ a structure HRe₂(CO)₉(μ-H), with one terminal and one bridging hydrogen atom, lies within 3 kcal/mol of the structure Re₂(CO)₉(η²-H₂), similar to that of Re₂(CO)₁₀. For H₂Re₂(CO)(n) (n = 8, 7, 6) the only low-energy structures are doubly bridged singlet Re₂(μ-H)₂(CO)(n) structures. Higher energy dihydrogen complex structures are also found.

Metal-metal bonding in binuclear metal carbonyls with three bridging methylaminobis(difluorophosphine) ligands: iron, cobalt, and nickel derivatives

The [CH3N(PF2)2]3M2(CO)n (M = Fe, Co, Ni; n = 2, 1, 0) derivatives, in which a metal-metal bond is bridged by three CH3N(PF2)2 groups similar to the known very stable [CH3N(PF2)2]3Co2(CO)2, have been investigated by density functional theory. The lowest energy structures for the dicarbonyls [CH3N(PF2)2]3M2(CO)2 (M = Ni, Co, Fe) are predicted to be closed-shell singlets with metal-metal distances of approximately 3.94 A, approximately 2.80 A, and approximately 2.63 A, respectively, corresponding to the metal-metal bonds of orders 0, 1, and 2, required to give both metal atoms the favored 18-electron configurations. For the monocarbonyls [CH3N(PF2)2]3M2(CO) (M = Ni, Co, Fe), unsymmetrical structures with a carbonyl group terminally bonded to one metal are greatly preferred energetically over structures in which the carbonyl group bridges the pair of metal atoms. The lowest energy structures for [CH3N(PF2)2]3M2(CO) (M = Ni, Co, Fe) are singlet, triplet, and singlet states, respectively, with metal-metal distances of approximately 2.77 A, approximately 2.38 A, and approximately 2.13 A, respectively. These distances correspond to the metal-metal dative bonds of orders 1, 2, and 3, required to give both metal atoms the favored 18-electron configuration. Among the carbonyl-free species the iron compound [CH3N(PF2)2]3Fe2 is of particular interest since its very short Fe[double double bonds]Fe distance of approximately 2.02 A suggests the formal quadruple bond required to give both iron atoms the favored 18-electron configuration.