Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Dopant-vacancy binding effects in Li-doped magnesium hydride

We use a combination of ab initio calculations and statistical mechanics to investigate the substitution of Li+ for Mg2+ in magnesium hydride (MgH2) accompanied by the formation of hydrogen vacancies with positive charge (with respect to the original ion at the site). We show that the binding energy between dopants and vacancy defects leads to a significant fraction of trapped vacancies and therefore a dramatic reduction in the number of free vacancies available for diffusion. The concentration of free vacancies initially increases with dopant concentration but reaches a maximum at around 1 mol % Li doping and slowly decreases with further doping. At the optimal level of doping, the corresponding concentration of free vacancies is much higher than the equilibrium concentrations of charged and neutral vacancies in pure MgH2 at typical hydrogen storage conditions. We also show that Li-doped MgH2 is thermodynamically metastable with respect to phase separation into pure magnesium and lithium hydrides at any significant Li concentration, even after considering the stabilization provided by dopant-vacancy interactions and configurational entropic effects. Our results suggest that lithium doping may enhance hydrogen diffusion hydride but only to a limited extent determined by an optimal dopant concentration and conditioned to the stability of the doped phase.

Phosphorus sorption and availability in soils amended with animal manures and sewage sludge

Soils that receive large applications of animal wastes and sewage sludge are vulnerable to releasing environmentally significant concentrations of dissolved P available to subsurface flow owing to the gradual saturation of the soil's P sorption capacity. This study evaluated P sorption (calculated from Langmuir isotherms) and availability of P (as CaCl2-P and resin P) in soils incubated for 20 d with poultry litter, poultry manure, cattle slurry, municipal sewage sludge, or KH2PO4, added on a P-equivalent basis (100 mg P kg(-1)). All the P sources had a marked negative effect on P sorption and a positive effect on P availability in all soils. In the cattle slurry- and KH2PO4- treated soils, the decreases in P sorption maximum (19-66%) and binding energy (25-89%) were consistently larger than the corresponding decreases (7-41% and 11-30%) in poultry litter-, poultry manure-, and sewage sludge-treated soils. The effects of cattle slurry and KH2PO4 on P availability were, in most cases, larger than those of the other P sources. In the poultry litter, poultry manure, and sewage sludge treatments, the increase in soil solution P was inversely related (R-2 = 0.75) to the input of Ca from these relatively high Ca (13.5-42 g kg(-1)) sources. Correlation analyses implied that the magnitude of the changes in P sorption and availability was not related to the water-extractable P content of the P sources. Future research on the sustainable application of organic wastes to agricultural soils needs to consider the non-P- as well as P-containing components of the waste.

Dissociation of water on oxygen-covered Rh{111}

The adsorption of water and coadsorption with oxygen on Rh{111} under ultrahigh vacuum conditions was studied using synchrotron-based photoemission and photoabsorption spectroscopy. Water adsorbs intact on the clean surface at temperatures below 154 K. Irradiation with x-rays, however, induces fast dissociation and the formation of a mixed OH+H(2)O layer indicating that the partially dissociated layer is thermodynamically more stable. Coadsorption of water and oxygen at a coverage below 0.3 monolayers has a similar effect, leading to the formation of a hydrogen-bonded network of water and hydroxyl molecules at a ratio of 3:2. The partially dissociated layers are more stable than chemisorbed intact water with the maximum desorption temperatures up to 30 K higher. For higher oxygen coverage, up to 0.5 monolayers, water does not dissociate and an intact water species is observed above 160 K, which is characterized by an O 1s binding energy 0.6 eV higher than that of chemisorbed water and a high desorption temperature similar to the partially dissociated layer. The extra stabilization is most likely due to hydrogen bonds with atomic oxygen.

Water adsorption on O-covered Ru{0001}: coverage-dependent change from dissociation to molecular adsorption

When water is coadsorbed with oxygen at coverages above 0.25ML an intact water species is observed in high resolution X-ray photoelectron spectroscopy up to 220 K, which is significantly more stable than intact water on the clean surface. The presence of this species causes a shift in the O 1s binding energy of the pre-adsorbed oxygen, which indicates the formation of hydrogen bonds between the two adsorbates. Low coverages of oxygen induce partial dissociation and recombinative desorption in the same temperature range, which illustrates that desorption temperatures alone cannot be used to determine whether water is molecularly intact or not.

Non-stoichiometric oxide and metal interfaces and reactions

We have employed a combination of experimental surface science techniques and density functional calculations to study the reduction of TiO2(110) surfaces through the doping with submonolayer transition metals. We concentrate on the role of Ti adatoms in self doping of rutile and contrast the behaviour to that of Cr. DFT+U calculations enable identification of probable adsorption structures and their spectroscopic characteristics. Adsorption of both metals leads to a broken symmetry and an asymmetric charge transfer localised around the defect site of a mixed localised/delocalised character. Charge transfer creates defect states with Ti 3d character in the band gap at similar to 1-eV binding energy. Cr adsorption, however, leads to a very large shift in the valence-band edge to higher binding energy and the creation of Cr 3d states at 2.8-eV binding energy. Low-temperature oxidation lifts the Ti-derived band-gap states and modifies the intensity of the Cr features, indicative of a change of oxidation state from Cr3+ to Cr4+. Higher temperature processing leads to a loss of Cr from the surface region, indicative of its substitution into the bulk.

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

Time-resolved kinetic studies of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reactions with oxirane, oxetane, and tetrahydrofuran (THF). The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at four or five temperatures in the range 294-605 K. All three reactions showed pressure dependences characteristic of third-body-assisted association reactions with, surprisingly, SiH2 + oxirane showing the least and SiH2 + THF showing the most pressure dependence. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equations where the error limits are single standard deviations: log(k(oxirane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.03 +/- 0.07) + (5.70 +/- 0.51) kJ mol(-1)/RT In 10 log(k(oxetane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.17 +/- 0.11) + (9.04 +/- 0.78) kJ mol(-1)/RT In 10 log(k(THF)(infinity)/cm(3) molecule(-1) s(-1)) = (-10.59 +/- 0.10) + (5.76 +/- 0.65) kJ mol(-1)/RT In 10 Binding-energy values of 77, 97, and 92 kJ mol(-1) have been obtained for the donor-acceptor complexes of SiH2 with oxirane, oxetane, and THF, respectively, by means of quantum chemical (ab initio) calculations carried Out at the G3 level. The use of these values to model the pressure dependences of these reactions, via RRKM theory, provided a good fit only in the case of SiH2 + THF. The lack of fit in the other two cases is attributed to further reaction pathways for the association complexes of SiH2 with oxirane and oxetane. The finding of ethene as a product of the SiH2 + oxirane reaction supports a pathway leading to H2Si=O + C2H4 predicted by the theoretical calculations of Apeloig and Sklenak.

Highly proton-ordered water structures on oxygen precovered Ru{0001}

We present helium scattering measurements of a water ad-layer grown on a O(2 1)/Ru(0001) surface. The adsorbed water layer results in a well ordered helium diffraction pattern with systematic extinctions of diffraction spots due to glide line symmetries. The data reflects a well-defined surface structure that maintains proton order even at surprisingly high temperatures of 140 K. The diffraction data we measure is consistent with a structure recently derived from STM measurements performed at 6 K. Comparison with recent DFT calculation is in partial agreement, suggesting that these calculations might be underestimating the contribution of relative water molecule orientations to the binding energy.

Effect of an ancillary ligand on single helix-double helix interconversion in copper complexes. Copper(I)-water bond

Reaction of Cu(ClO(4))(2)center dot 6H(2)O with the 1:2 condensate of benzildihydrazone and 2-acetylpyridine, in methanol in equimolar ratio yields a green compound which upon recrystallisation from 1:1 CH(2)Cl(2)-C(6)H(6) mixture affords [CuL(H(2)O)](ClO(4))(2)center dot 1/2C(6)H(6). The complex crystallises in the space group P-1 with a = 8.028(11) angstrom, b = 12.316(17) angstrom, c = 18.14(3) angstrom, alpha = 97.191(10)degrees, beta = 94.657(10)degrees and gamma = 108.039(10)degrees. It is single helical with the metal having a distorted trigonal bipyramidal N(4)O coordination sphere. The acid dissociation constant of the Cu(I) complex in CH(3)CN is 3.34 +/- 0.19. The X band EPR spectrum of the compound is rhombic with g(1) = 2.43, g(2) = 2.10 g(3) = 2.02 and A(1) = 79.3 x 10(-4) cm(-1). The Cu(II/I) potential of the complex in CH(2)Cl(2) at a glassy carbon electrode is 0.43 V vs SCE. It is argued that the copper-water bond persists in the corresponding copper(I) species. Its implications on the single helix-double helix interconversion in copper helicates are discussed. DFT calculations at the B3LYP/6-311G** level shows that the binding energy of water in the single helicol live-coordinate copper(I) species [CuL(H(2)O)](+) is similar to 40 kJ mol(-1).

Emulation of chiral hydrogenation catalysts by adsorbtion of small molecules onto Ni surfaces

Adsorption of small molecules on the Ni{111} and NiO{111} surfaces is investigated under UHV and elevated pressures (~10-1 mbar) of hydrogen and water. The molecules considered are chosen for their relevance to understanding the mechanism of enantioselective hydrogenation on Raney Nickel modified by chiral molecules. Adsorption of water onto, and its subsequent reaction with, oxygen-covered Ni{111} is dependent on the initial atomic oxygen coverage. An OH species (O1s binding energy 531.5eV), oriented perpendicular to the surface, forms at atomic oxygen coverages <0.25ML. The reaction does not consume all the adsorbed oxygen for coverages ≥0.12ML. The p(2×2) atomic oxygen uperstructure is unreactive, while an OH species is formed on the p(√3×√3) superstructure at binding energy 530.9eV. L-alanine is adsorbed on Ni{111} as a model chiral modifier molecule. At low coverages, alanine forms a presumed tridentate alaninate species for coverages ≥0.11ML at 250K. A minority, bidentate zwitterionic species forms at coverages >0.11ML, but was not observed at 300K. Saturation occurs at 0.25ML. At high alanine coverages (≥0.19ML) and H2 pressure (≥1×10-2 mbar), the tridentate L-alaninate converts to bidentate zwitterionic L-alanine at 300K. Thermal evolution of L-alanine on Ni{111} under varying hydrogen pressures is examined. Adsorption of L-alanine onto hydroxylated NiO{111} at 300K in UHV, mimicking a catalyst surface under aqueous conditions, yields the tridentate alaninate which is immune to the effects of elevated hydrogen pressure. Exposing the L-alanine/Ni{111} adsorption system to water (≤10-1 mbar) oxidises the surface and recreates the L-alanine/hydroxylated NiO{111} system. Pyruvic acid on Ni{111} is examined as a model for hydrogenation substrate adsorption. Behaviour is coverage dependent and several conformations are possible at low coverages (≤0.1ML). Annealing at coverages <0.2ML causes a condensation reaction, releasing water onto the surface. High coverages do not condense and a saturation coverage of ~0.35ML is found.

Accumulated metal speciation in earthworm populations with multi-generational exposure to metalliferous soils: cell fractionation and high energy synchrotron analyses

Predicting metal bioaccumulation and toxicity in soil organisms is complicated by site-specific biotic and abiotic parameters. In this study we exploited tissue fractionation and digestion techniques, combined with X-ray absorption spectroscopy (XAS), to investigate the whole-body and subcellular distributions, ligand affinities, and coordination chemistry of accumulated Pb and Zn in field populations of the epigeic earthworm Lumbricus rubellus inhabiting three contrasting metalliferous and two unpolluted soils. Our main findings were (i) earthworms were resident in soils with concentrations of Pb and Zn ranging from 1200 to 27 000 mg kg(-1) and 200 to 34 000 mg kg(-1), respectively; (ii) Pb and Zn primarily accumulated in the posterior alimentary canal in nonsoluble subcellular fractions of earthworms; (iii) site-specific differences in the tissue and subcellular partitioning profiles of populations were observed, with earthworms from a calcareous site partitioning proportionally more Pb to their anterior body segments and Zn to the chloragosome-rich subcellular fraction than their acidic-soil inhabiting counterparts; (iv) XAS indicated that the interpopulation differences in metal partitioning between organs were not accompanied by qualitative differences in ligand-binding speciation, because crystalline phosphate-containing pyromorphite was a predominant chemical species in the whole-worm tissues of all mine soil residents. Differences in metal (Pb, Zn) partitioning at both organ and cellular levels displayed by field populations with protracted histories of metal exposures may reflect their innate ecophysiological responses to essential edaphic variables, such as Ca2+ status. These observations are highly significant in the challenging exercise of interpreting holistic biomarker data delivered by "omic" technologies.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background: MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results: A large dataset comprising MHC-peptide structural complexes was created by remodelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion: The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

Dual role of the beta(2)-adrenergic receptor C terminus for the binding of beta-arrestin and receptor internalization

Homologous desensitization of beta(2)-adrenergic and other G-protein-coupled receptors is a two-step process. After phosphorylation of agonist-occupied receptors by G-protein-coupled receptor kinases, they bind beta-arrestins, which triggers desensitization and internalization of the receptors. Because it is not known which regions of the receptor are recognized by beta-arrestins, we have investigated beta-arrestin interaction and internalization of a set of mutants of the human beta(2)-adrenergic receptor. Mutation of the four serine/threonine residues between residues 355 and 364 led to the loss of agonist-induced receptor-beta-arrestin2 interaction as revealed by fluorescence resonance energy transfer (FRET), translocation of beta-arrestin2 to the plasma membrane, and receptor internalization. Mutation of all seven serine/threonine residues distal to residue 381 did not affect agonist-induced receptor internalization and beta-arrestin2 translocation. A beta(2)-adrenergic receptor truncated distal to residue 381 interacted normally with beta-arrestin2, whereas its ability to internalize in an agonist-dependent manner was compromised. A similar impairment of internalization was observed when only the last eight residues of the C terminus were deleted. Our experiments show that the C terminus distal to residue 381 does not affect the initial interaction between receptor and beta-arrestin, but its last eight amino acids facilitate receptor internalization in concert with beta-arrestin2.