Ferric and ferrous iron in nitroso-myoglobin: computer simulations of stable and metastable states and their infrared spectra

The binding of NO to iron is involved in the biological function of many heme proteins. Contrary to ligands like CO and O-2, which only bind to ferrous (Fe-II) iron, NO binds to both ferrous and ferric (Fe-II) iron. In a particular protein, the natural oxidation state can therefore be expected to be tailored to the required function. Herein, we present an ob initio potential-energy surface for ferric iron interacting with NO. This potential-energy surface exhibits three minima corresponding to eta'-NO coordination (the global minimum), eta(1)-ON coordination and eta(2) coordination. This contrasts with the potential-energy surface for Fe-II-NO, which ex- hibits only two minima (the eta(2) coordination mode for Fe-II is a transition state, not a minimum). In addition, the binding energies of NO are substantially larger for Fe-III than for Fe-II. We have performed molecular dynamics simulations for NO bound to ferric myoglobin (Mb(III)) and compare these with results obtained for Mb(II). Over the duration of our simulations (1.5 ns), all three binding modes are found to be stable at 200 K and transiently stable at 300 K, with eventual transformation to the eta(1)-NO global-minimum conformation. We discuss the implication of these results related to studies of rebinding processes in myoglobin.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

Spectro-electrochemical and DFT studies of a planar Cu(II)–phenolate complex active in the aerobic oxidation of primary alcohols

A square-planar compound [Cu(pyrimol)Cl] (pyrimol = 4-methyl-2-N-(2-pyridylmethylene)aminophenolate) abbreviated as CuL–Cl) is described as a biomimetic model of the enzyme galactose oxidase (GOase). This copper(II) compound is capable of stoichiometric aerobic oxidation of activated primary alcohols in acetonitrile/water to the corresponding aldehydes. It can be obtained either from Hpyrimol (HL) or its reduced/hydrogenated form Hpyramol (4-methyl-2-N-(2-pyridylmethyl)aminophenol; H2L) readily converting to pyrimol (L-) on coordination to the copper(II) ion. Crystalline CuL–Cl and its bromide derivative exhibit a perfect square-planar geometry with Cu–O(phenolate) bond lengths of 1.944(2) and 1.938(2) Å. The cyclic voltammogram of CuL–Cl exhibits an irreversible anodic wave at +0.50 and +0.57 V versus ferrocene/ferrocenium (Fc/Fc+) in dry dichloromethane and acetonitrile, respectively, corresponding to oxidation of the phenolate ligand to the corresponding phenoxyl radical. In the strongly donating acetonitrile the oxidation path involves reversible solvent coordination at the Cu(II) centre. The presence of the dominant CuII–L. chromophore in the electrochemically and chemically oxidised species is evident from a new fairly intense electronic absorption at 400–480 nm ascribed to a several electronic transitions having a mixed pi-pi(L.) intraligand and Cu–Cl -> L. charge transfer character. The EPR signal of CuL–Cl disappears on oxidation due to strong intramolecular antiferromagnetic exchange coupling between the phenoxyl radical ligand (L.) and the copper(II) centre, giving rise to a singlet ground state (S = 0). The key step in the mechanism of the primary alcohol oxidation by CuL–Cl is probably the alpha-hydrogen abstraction from the equatorially bound alcoholate by the phenoxyl moiety in the oxidised pyrimol ligand, Cu–L., through a five-membered cyclic transition state.

Time-resolved gas-phase kinetic, quantum chemical, and RRKM studies of reactions of silylene with alcohols

Time-resolved kinetic studies of silylene, SiH2, generated by laser flash photolysis of 1-silacyclopent-3-ene and phenylsilane, have been carried out to obtain rate constants for its bimolecular reactions with methanol, ethanol, 1-propanol, 1-butanol and 2-methyl-1-butanol. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at room temperature. In the study with methanol several buffer gases were used. All five reactions showed pressure dependences characteristic of third body assisted association reactions. The rate constant pressure dependences were modelled using RRKM theory, based on Eo values of the association complexes obtained by ab initio calculation (G3 level). Transition state models were adjusted to fit experimental fall-off curves and extrapolated to obtain k∞ values in the range 1.9 to 4.5 × 10-10 cm3 molecule-1 s-1. These numbers, corresponding to the true bimolecular rate constants, indicate efficiencies of between 16 and 67% of the collision rates for these reactions. In the reaction of SiH2 + MeOH there is a small kinetic component to the rate which is second order in MeOH (at low total pressures). This suggests an additional catalysed reaction pathway, which is supported by the ab initio calculations. These calculations have been used to define specific MeOH-for-H2O substitution effects on this catalytic pathway. Where possible our experimental and theoretical results are compared with those of previous studies.

Mode-specific reactivity of CH4 on Pt(110)-(1Ã2): The concerted role of stretch and bend excitation

The state-resolved reaction probability of CH4 on Pt�110-�1�2 was measured as a function of CH4 translational energy for four vibrational eigenstates comprising different amounts of C-H stretch and bend excitation. Mode-specific reactivity is observed both between states from different polyads and between isoenergetic states belonging to the same polyad of CH4. For the stretch/bend combination states, the vibrational efficacy of reaction activation is observed to be higher than for either pure C-H stretching or pure bending states, demonstrating a concerted role of stretch and bend excitation in C-H bond scission. This concerted role, reflected by the nonadditivity of the vibrational efficacies, is consistent with transition state structures found by ab initio calculations and indicates that current dynamical models of CH4 chemisorption neglect an important degree of freedom by including only C-H stretching motion.

Warranty cost analysis for products with a dormant state

Most building services products are installed while a building is constructed, but they are not operated until the building is commissioned. The warranty of the products may cover the time starting from their installation to the end of the warranty period. Prior to the commissioning of the building, the products are at a dormant mode (i.e., not operated) but protected by the warranty. For such products, both the usage intensity and the failure patterns are different from those with continuous usage intensity and failure patterns. This paper develops warranty cost models for repairable products with a dormant mode from both the manufacturer's and buyer's perspectives. Relationships between the failure patterns at the dormant mode and at the operational mode are also discussed. Numerical examples and sensitivity analysis are used to demonstrate the applicability of the methodology derived in the paper.

Revealing ensemble state transition patterns in multi-electrode neuronal recordings using hidden Markov models

In order to harness the computational capacity of dissociated cultured neuronal networks, it is necessary to understand neuronal dynamics and connectivity on a mesoscopic scale. To this end, this paper uncovers dynamic spatiotemporal patterns emerging from electrically stimulated neuronal cultures using hidden Markov models (HMMs) to characterize multi-channel spike trains as a progression of patterns of underlying states of neuronal activity. However, experimentation aimed at optimal choice of parameters for such models is essential and results are reported in detail. Results derived from ensemble neuronal data revealed highly repeatable patterns of state transitions in the order of milliseconds in response to probing stimuli.

Reversible state transition in nanoconfined aqueous solutions

Using molecular dynamics simulations, we find a reversible transition between the dispersion and aggregation states of solute molecules in aqueous solutions confined in nanoscale geometry, which is not observed in macroscopic systems. The nanoscale confinement also leads to a significant increase of the critical aggregation concentration (CAC). A theoretical model based on Gibbs free energy calculation is developed to describe the simulation results. It indicates that the reversible state transition is attributed to the low free energy barrier (of order kBT) in between two energy minima corresponding to the dispersion and aggregation states, and the enhancement of the CAC results from the fact that at lower concentrations the number of solute molecules is not large enough to allow the formation of a stable cluster in the confined systems.

Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach

The analysis step of the (ensemble) Kalman filter is optimal when (1) the distribution of the background is Gaussian, (2) state variables and observations are related via a linear operator, and (3) the observational error is of additive nature and has Gaussian distribution. When these conditions are largely violated, a pre-processing step known as Gaussian anamorphosis (GA) can be applied. The objective of this procedure is to obtain state variables and observations that better fulfil the Gaussianity conditions in some sense. In this work we analyse GA from a joint perspective, paying attention to the effects of transformations in the joint state variable/observation space. First, we study transformations for state variables and observations that are independent from each other. Then, we introduce a targeted joint transformation with the objective to obtain joint Gaussianity in the transformed space. We focus primarily in the univariate case, and briefly comment on the multivariate one. A key point of this paper is that, when (1)-(3) are violated, using the analysis step of the EnKF will not recover the exact posterior density in spite of any transformations one may perform. These transformations, however, provide approximations of different quality to the Bayesian solution of the problem. Using an example in which the Bayesian posterior can be analytically computed, we assess the quality of the analysis distributions generated after applying the EnKF analysis step in conjunction with different GA options. The value of the targeted joint transformation is particularly clear for the case when the prior is Gaussian, the marginal density for the observations is close to Gaussian, and the likelihood is a Gaussian mixture.

Hidden Markov analysis describes dive patterns in semiaquatic animals

Existing methods of dive analysis, developed for fully aquatic animals, tend to focus on frequency of behaviors rather than transitions between them. They, therefore, do not account for the variability of behavior of semiaquatic animals, and the switching between terrestrial and aquatic environments. This is the first study to use hidden Markov models (HMM) to divide dives of a semiaquatic animal into clusters and thus identify the environmental predictors of transition between behavioral modes. We used 18 existing data sets of the dives of 14 American mink (Neovison vison) fitted with time-depth recorders in lowland England. Using HMM, we identified 3 behavioral states (1, temporal cluster of dives; 2, more loosely aggregated diving within aquatic activity; and 3, terminal dive of a cluster or a single, isolated dive). Based on the higher than expected proportion of dives in State 1, we conclude that mink tend to dive in clusters. We found no relationship between temperature and the proportion of dives in each state or between temperature and the rate of transition between states, meaning that in our study area, mink are apparently not adopting different diving strategies at different temperatures. Transition analysis between states has shown that there is no correlation between ambient temperature and the likelihood of mink switching from one state to another, that is, changing foraging modes. The variables provided good discrimination and grouped into consistent states well, indicating promise for further application of HMM and other state transition analyses in studies of semiaquatic animals.