Outer-Sphere Contributions to the Electronic Structure of Type Zero Copper Proteins

Bioinorganic canon states that active-site thiolate coordination promotes rapid electron transfer (ET) to and from type 1 copper proteins. In recent work, we have found that copper ET sites in proteins also can be constructed without thiolate ligation (called “type zero” sites). Here we report multifrequency electron paramagnetic resonance (EPR), magnetic circular dichroism (MCD), and nuclear magnetic resonance (NMR) spectroscopic data together with density functional theory (DFT) and spectroscopy-oriented configuration interaction (SORCI) calculations for type zero Pseudomonas aeruginosa azurin variants. Wild-type (type 1) and type zero copper centers experience virtually identical ligand fields. Moreover, O-donor covalency is enhanced in type zero centers relative that in the C112D (type 2) protein. At the same time, N-donor covalency is reduced in a similar fashion to type 1 centers. QM/MM and SORCI calculations show that the electronic structures of type zero and type 2 are intimately linked to the orientation and coordination mode of the carboxylate ligand, which in turn is influenced by outer-sphere hydrogen bonding

Magnetic Properties of Monomer and Dimer Tetrahedral VOx Entities Dispersed on Amorphous Silica-based Materials: Prediction of EPR Parameters from Relativistic DFT Calculations and Broken Symmetry Approach to Exchange Couplings

Molecular structures of the isolated tetrahedral oxovanadium(IV) and bridged μ-oxo-divanadium(IV) complexes hosted by the clusters mimicking surfaces of amorphous silica-based materials were investigated using density functional theory (DFT) calculations. Principal values of the g and A tensors for the monomer vanadyl species were obtained using the coupled-perturbed DFT level of theory and the spin–orbit mean-field approximation (SOMF). Magnetic exchange interaction for the μ-oxo bridged vanadium(IV) dimer was investigated within the broken symmetry approach. An antiferromagnetic coupling of the individual magnetic moments of the vanadium(IV) centers in the [VO–O–VO]2+ bridges was revealed and discussed in detail. The coupling explains pronounced decrease of the electron paramagnetic resonance signal (EPR) intensity, observed for the reduced VOx/SiO2 samples with the increasing coverage of vanadia, in terms of transformation of the paramagnetic monomer species into the dimers with S = 0 ground state

Absence of Insect Juvenile Hormones in the American Dog Tick, Dermacentor veriabilis (Say) (Acari: Ixodidae), and in Ornithodoros parkeri Cooley (Acari: Argasidae)

Synganglia, salivary gland, midgut, ovary, fat body and muscle alone and in combination from the ixodid tick, Dermacentor variabilis (Say), or the argasid tick, Ornithodoros parkeri Cooley, were incubated in vitro in separate experiments with L-[methyl-3H]methionine and farnesoic acid or with [1-14C]acetate. Life stages examined in D. variabilis were 3 and 72 h old (after ecdysis) unfed nymphs, partially fed nymphs (18 and 72 h after attachment to the host), fully engorged nymphs (2 d after detachment from host), 3 and 72 h old (after eclosion) unfed females, partially fed unmated females (12–168 h after attachment to host) and mated replete females (2 d after detachment from the host). Those from O. parkeri were third and fourth stadium nymphs and female O. parkeri, 1–2 d after detachment. Corpora allata from Diploptera punctata, Periplaneta americana and Gromphadorina portentosa were used as positive controls in these experiments. No farnesol, methyl farnesoate, JH I, JH II, JH III, or JHIII bisepoxide was detected by radio HPLC from any tick analysis while JH III, methyl farnesoate, and farnesol were detected in the positive controls. To examine further for the presence of a tick, insect-juvenilizing agent, Galleria pupal–cuticle bioassays were conducted on lipid extracts from 10 and 15 d old eggs, unfed larvae (1–5 d after ecdysis), unfed nymphs (1–7 d after ecdysis), and partially fed, unmated female adults (completed slow feeding phase) of D. variabilis. Whole body extracts of fourth stadium D. punctata and JH III standard were used as positive controls. No juvenilizing activity in any of the tick extracts could be detected. Electron impact, gas chromatography–mass spectrometry of hemolymph extracts from fed, virgin (forcibly detached 7 d after attachment) and mated, replete (allowed to drop naturally) D. variabilis and fully engorged (1–2 d after detachment) O. parkeri females also failed to identify the common insect juvenile hormones. The same procedures were successful in the identification of JH III in hemolymph of fourth stadium D. punctata. Last stadium nymphal (female) O. parkeri implanted with synganglia from second nymphal instars underwent normal eclosion to the adult. The above studies in toto suggest that D. variabilis and O. parkeri do not have the ability to make the common insect juvenile hormones, and these juvenile hormones do not regulate tick metamorphosis or reproduction as hypothesized in the literature

Density functional theory

Density functional theory (DFT) finds increasing use in applications related to biological systems. Advancements in methodology and implementations have reached a point where predicted properties of reasonable to high quality can be obtained. Thus, DFT studies can complement experimental investigations, or even venture with some confidence into experimentally unexplored territory. In the present contribution, we provide an overview of the properties that can be calculated with DFT, such as geometries, energies, reaction mechanisms, and spectroscopic properties. A wide range of spectroscopic parameters is nowadays accessible with DFT, including quantities related to infrared and optical spectra, X-ray absorption and Mössbauer, as well as all of the magnetic properties connected with electron paramagnetic resonance spectroscopy except relaxation times. We highlight each of these fields of application with selected examples from the recent literature and comment on the capabilities and limitations of current methods

On the feasibility of N2 fixation via a single-site FeI/FeIV cycle: Spectroscopic studies of FeI(N2)FeI, FeIV=N, and related species

The electronic properties of an unusually redox-rich iron system, [PhBPR 3]FeNx (where [PhBPR 3] is [PhB(CH2PR2)3]−), are explored by Mössbauer, EPR, magnetization, and density-functional methods to gain a detailed picture regarding their oxidation states and electronic structures. The complexes of primary interest in this article are the two terminal iron(IV) nitride species, [PhBPiPr 3]FeN (3a) and [PhBPCH2Cy 3]FeN (3b), and the formally diiron(I) bridged-Fe(μ-N2)Fe species, {[PhBPiPr 3]Fe}2(μ-N2) (4). Complex 4 is chemically related to 3a via a spontaneous nitride coupling reaction. The diamagnetic iron(IV) nitrides 3a and 3b exhibit unique electronic environments that are reflected in their unusual Mössbauer parameters, including quadrupole-splitting values of 6.01(1) mm/s and isomer shift values of −0.34(1) mm/s. The data for 4 suggest that this complex can be described by a weak ferromagnetic interaction (J/D < 1) between two iron(I) centers. For comparison, four other relevant complexes also are characterized: a diamagnetic iron(IV) trihydride [PhBPiPr 3]Fe(H)3(PMe3) (5), an S = 3/2 iron(I) phosphine adduct [PhBPiPr 3]FePMe3 (6), and the S = 2 iron(II) precursors to 3a, [PhBPiPr 3]FeCl and [PhBPiPr 3]Fe-2,3:5,6-dibenzo-7-aza bicyclo[2.2.1]hepta-2,5-diene (dbabh). The electronic properties of these respective complexes also have been explored by density-functional methods to help corroborate our spectral assignments and to probe their electronic structures further

Search for LBV Candidates in the M33 Galaxy

A total of 185 luminous blue variable (LBV) candidates with V < 18.5 and B-V < 0.35 are selected based on the photometrical Survey of Local Group Galaxies made by P. Massey et al. 2006. The candidates were selected using aperture photometry of H-alpha images. The primary selection criterion is that the prospective candidate should be a blue star with H-aplha emission. In order not to miss appreciably reddened LBV candidates, we compose an additional list of 25 presumably reddened (0.35 < B-V < 1.2, V < 18.5) emission star candidates. A comparison with the list of known variables in the M33 galaxy showed 29% of our selected candidates to be photometrically variable. We also find our list to agree well with the lists of emission-line objects obtained in earlier papers using different methods.Comment: 6 figure

Photochemistry of Furyl- and Thienyldiazomethanes: Spectroscopic Characterization of Triplet 3-Thienylcarbene

Photolysis (λ \u3e 543 nm) of 3-thienyldiazomethane (1), matrix isolated in Ar or N2 at 10 K, yields triplet 3-thienylcarbene (13) and α-thial-methylenecyclopropene (9). Carbene 13 was characterized by IR, UV/vis, and EPR spectroscopy. The conformational isomers of 3-thienylcarbene (s-E and s-Z) exhibit an unusually large difference in zero-field splitting parameters in the triplet EPR spectrum (|D/hc| = 0.508 cm–1, |E/hc| = 0.0554 cm–1; |D/hc| = 0.579 cm–1, |E/hc| = 0.0315 cm–1). Natural Bond Orbital (NBO) calculations reveal substantially differing spin densities in the 3-thienyl ring at the positions adjacent to the carbene center, which is one factor contributing to the large difference in D values. NBO calculations also reveal a stabilizing interaction between the sp orbital of the carbene carbon in the s-Z rotamer of 13 and the antibonding σ orbital between sulfur and the neighboring carbon—an interaction that is not observed in the s-E rotamer of 13. In contrast to the EPR spectra, the electronic absorption spectra of the rotamers of triplet 3-thienylcarbene (13) are indistinguishable under our experimental conditions. The carbene exhibits a weak electronic absorption in the visible spectrum (λmax = 467 nm) that is characteristic of triplet arylcarbenes. Although studies of 2-thienyldiazomethane (2), 3-furyldiazomethane (3), or 2-furyldiazomethane (4) provided further insight into the photochemical interconversions among C5H4S or C5H4O isomers, these studies did not lead to the spectroscopic detection of the corresponding triplet carbenes (2-thienylcarbene (11), 3-furylcarbene (23), or 2-furylcarbene (22), respectively)

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.Financial support by the “Ministerio de Economía y Competitividad” (MINECO) of Spain and European FEDER funds through projects CTQ2011-27253 and CTQ2012-31914 is acknowledged. The support of the Generalitat Valenciana (Prometeo/2012/053) is also acknowledged. J.A. thanks the EU for the FP7-PEOPLE-2012-IEF-329513 grant. J.C. acknowledges the “Ministerio de Educación, Cultura y Deporte” (MECD) of Spain for a predoctoral FPU grant

Accurate ab initio spin densities

We present an approach for the calculation of spin density distributions for molecules that require very large active spaces for a qualitatively correct description of their electronic structure. Our approach is based on the density-matrix renormalization group (DMRG) algorithm to calculate the spin density matrix elements as basic quantity for the spatially resolved spin density distribution. The spin density matrix elements are directly determined from the second-quantized elementary operators optimized by the DMRG algorithm. As an analytic convergence criterion for the spin density distribution, we employ our recently developed sampling-reconstruction scheme [J. Chem. Phys. 2011, 134, 224101] to build an accurate complete-active-space configuration-interaction (CASCI) wave function from the optimized matrix product states. The spin density matrix elements can then also be determined as an expectation value employing the reconstructed wave function expansion. Furthermore, the explicit reconstruction of a CASCI-type wave function provides insights into chemically interesting features of the molecule under study such as the distribution of $\alpha$- and $\beta$-electrons in terms of Slater determinants, CI coefficients, and natural orbitals. The methodology is applied to an iron nitrosyl complex which we have identified as a challenging system for standard approaches [J. Chem. Theory Comput. 2011, 7, 2740].Comment: 37 pages, 13 figure

Theoretical investigation of the electronic structure of Fe(II) complexes at spin-state transitions

The electronic structure relevant to low spin (LS)high spin (HS) transitions in Fe(II) coordination compounds with a FeN6 core are studied. The selected [Fe(tz)6]2+(1) (tz=1H-tetrazole), [Fe(bipy)3]2+(2) (bipy=2,2’-bipyridine) and [Fe(terpy)2]2+ (3) (terpy=2,2’:6’,2’’-terpyridine) complexes have been actively studied experimentally, and with their respective mono-, bi-, and tridentate ligands, they constitute a comprehensive set for theoretical case studies. The methods in this work include density functional theory (DFT), time-dependent DFT (TD-DFT) and multiconfigurational second order perturbation theory (CASPT2). We determine the structural parameters as well as the energy splitting of the LS-HS states (ΔEHL) applying the above methods, and comparing their performance. We also determine the potential energy curves representing the ground and low-energy excited singlet, triplet, and quintet d6 states along the mode(s) that connect the LS and HS states. The results indicate that while DFT is well suited for the prediction of structural parameters, an accurate multiconfigurational approach is essential for the quantitative determination of ΔEHL. In addition, a good qualitative agreement is found between the TD-DFT and CASPT2 potential energy curves. Although the TD-DFT results might differ in some respect (in our case, we found a discrepancy at the triplet states), our results suggest that this approach, with due care, is very promising as an alternative for the very expensive CASPT2 method. Finally, the two dimensional (2D) potential energy surfaces above the plane spanned by the two relevant configuration coordinates in [Fe(terpy)2]2+ were computed both at the DFT and CASPT2 levels. These 2D surfaces indicate that the singlet-triplet and triplet-quintet states are separated along different coordinates, i.e. different vibration modes. Our results confirm that in contrast to the case of complexes with mono- and bidentate ligands, the singlet-quintet transitions in [Fe(terpy)2]2+ cannot be described using a single configuration coordinate