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

Can domain-based local pair natural orbitals approaches accurately predict phosphorescence energies?

Since the discovery of the peculiar conducting and optical properties of aromatics, many efforts have been made to characterize and predict their phosphorescence. This physical process is exploited in modern Organic Emitting Light Diodes (OLEDs), and it is also one of the processes decreasing the efficiency of Dye-sensitized solar cells (DSSCs). Herein, we propose a computational strategy for the accurate calculation of singlet–triplet gaps of aromatic compounds, which provides results that are in excellent agreement with available experimental data. Our approach relies on the domain-based local pair natural orbital (DLPNO) variant of the “gold standard” CCSD(T) method. The convergence of our results with respect to the key technical parameters of the calculation, such as the basis set used, the approximations employed in the perturbative triples correction, and the dimension of the PNOs space, was thoroughly discussed

Laboratory Characterization and Astrophysical Detection of Vibrationally Excited States of Vinyl Cyanide in Orion-KL

New laboratory data of CH$_2$CHCN (vinyl cyanide) in its ground and vibrationally excited states at the microwave to THz domain allow searching for these excited state transitions in the Orion-KL line survey. Frequency-modulated spectrometers combined into a single broadband 50-1900 GHz spectrum provided measurements of CH$_2$CHCN covering a spectral range of 18-1893 GHz, whose assignments was confirmed by Stark modulation spectra in the 18-40 GHz region and by ab-initio anharmonic force field calculations. For analyzing the emission lines of CH$_2$CHCN species detected in Orion-KL we used the excitation and radiative transfer code (MADEX) at LTE conditions. The rotational transitions of the ground state of this molecule emerge from four cloud components of hot core nature which trace the physical and chemical conditions of high mass star forming regions in the Orion-KL Nebula. The total column density of CH$_2$CHCN in the ground state is (3.0$\pm$0.9)x10$^{15}$ cm$^{-2}$. We report on the first interstellar detection of transitions in the v10=1/(v11=1,v15=1) dyad in space, and in the v11=2 and v11=3 states in Orion-KL. The lowest energy vibrationally excited states of vinyl cyanide such as v11=1 (at 328.5 K), v15=1 (at 478.6 K), v11=2 (at 657.8 K), the v10=1/(v11=1,v15=1) dyad (at 806.4/809.9 K), and v11=3 (at 987.9 K) are populated under warm and dense conditions, so they probe the hottest parts of the Orion-KL source. Column density and rotational and vibrational temperatures for CH$_2$CHCN in their ground and excited states, as well as for the isotopologues, have been constrained by means of a sample of more than 1000 lines in this survey. Moreover, we present the detection of methyl isocyanide (CH$_3$NC) for the first time in Orion-KL and a tentative detection of vinyl isocyanide (CH$_2$CHNC) and give column density ratios between the cyanide and isocyanide isomers.Comment: 46 pages, 22 figures, 14 tables, 9 online table

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of “vertically desolvated” spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run

Explicit kinetic heterogeneity: mechanistic models for interpretation of labeling data of heterogeneous cell populations

Estimation of division and death rates of lymphocytes in different conditions is vital for quantitative understanding of the immune system. Deuterium, in the form of deuterated glucose or heavy water, can be used to measure rates of proliferation and death of lymphocytes in vivo. Inferring these rates from labeling and delabeling curves has been subject to considerable debate with different groups suggesting different mathematical models for that purpose. We show that the three models that are most commonly used are in fact mathematically identical and differ only in their interpretation of the estimated parameters. By extending these previous models, we here propose a more mechanistic approach for the analysis of data from deuterium labeling experiments. We construct a model of "kinetic heterogeneity" in which the total cell population consists of many sub-populations with different rates of cell turnover. In this model, for a given distribution of the rates of turnover, the predicted fraction of labeled DNA accumulated and lost can be calculated. Our model reproduces several previously made experimental observations, such as a negative correlation between the length of the labeling period and the rate at which labeled DNA is lost after label cessation. We demonstrate the reliability of the new explicit kinetic heterogeneity model by applying it to artificially generated datasets, and illustrate its usefulness by fitting experimental data. In contrast to previous models, the explicit kinetic heterogeneity model 1) provides a mechanistic way of interpreting labeling data; 2) allows for a non-exponential loss of labeled cells during delabeling, and 3) can be used to describe data with variable labeling length

Reduction of CO\u3csub\u3e2\u3c/sub\u3e By a Masked Two-Coordinate Cobalt(I) Complex and Characterization of a Proposed Oxodicobalt(II) Intermediate

Fixation and chemical reduction of CO2 are important for utilization of this abundant resource, and understanding the detailed mechanism of C-O cleavage is needed for rational development of CO2 reduction methods. Here, we describe a detailed analysis of the mechanism of the reaction of a masked two-coordinate cobalt(i) complex, LtBuCo (where LtBu = 2,2,6,6-tetramethyl-3,5-bis[(2,6-diisopropylphenyl)imino]hept-4-yl), with CO2, which yields two products of C-O cleavage, the cobalt(i) monocarbonyl complex LtBuCo(CO) and the dicobalt(ii) carbonate complex (LtBuCo)2(μ-CO3). Kinetic studies and computations show that the κN,η6-arene isomer of LtBuCo rearranges to the κ2N,N′ binding mode prior to binding of CO2, which contrasts with the mechanism of binding of other substrates to LtBuCo. Density functional theory (DFT) studies show that the only low-energy pathways for cleavage of CO2 proceed through bimetallic mechanisms, and DFT and highly correlated domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) calculations reveal the cooperative effects of the two metal centers during facile C-O bond rupture. A plausible intermediate in the reaction of CO2 with LtBuCo is the oxodicobalt(II) complex LtBuCoOCoLtBu, which has been independently synthesized through the reaction of LtBuCo with N2O. The rapid reaction of LtBuCoOCoLtBu with CO2 to form the carbonate product indicates that the oxo species is kinetically competent to be an intermediate during CO2 cleavage by LtBuCo. LtBuCoOCoLtBu is a novel example of a thoroughly characterized molecular cobalt-oxo complex where the cobalt ions are clearly in the +2 oxidation state. Its nucleophilic reactivity is a consequence of high charge localization on the μ-oxo ligand between two antiferromagnetically coupled high-spin cobalt(ii) centers, as characterized by DFT and multireference complete active space self-consistent field (CASSCF) calculations

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

Using Cepheids to determine the galactic abundance gradient I. The solar neighbourhood

A number of studies of abundance gradients in the galactic disk have been performed in recent years. The results obtained are rather disparate: from no detectable gradient to a rather significant slope of about -0.1 dex kpc -1. The present study concerns the abundance gradient based on the spectroscopic analysis of a sample of classical Cepheids. These stars enable one to obtain reliable abundances of a variety of chemical elements. Additionally, they have well determined distances which allow an accurate determination of abundance distributions in the galactic disc. Using 236 high resolution spectra of 77 galactic Cepheids, the radial elemental distribution in the galactic disc between galactocentric distances in the range 6-11 kpc has been investigated. Gradients for 25 chemical elements (from carbon to gadolinium) are derived...Comment: 28 pages, 14 postscript figures, LaTeX, uses Astronomy and Astrophysics macro aa.cls, graphicx package, to be published in Astronomy and Astrophysics (2002) also available at http://www.iagusp.usp.br/~maciel/index.htm