Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

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.

Infrared absorption spectra, radiative efficiencies, and global warming potentials of perfluorocarbons: Comparison between experiment and theory

Experimentally and theoretically determined infrared spectra are reported for a series of straight-chain perfluorocarbons: C2F6, C3F8, C4F10, C5F12, C6F14, and C8F18. Theoretical spectra were determined using both density functional (DFT) and ab initio methods. Radiative efficiencies (REs) were determined using the method of Pinnock et al. (1995) and combined with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). Theoretically determined absorption cross sections were within 10% of experimentally determined values. Despite being much less computationally expensive, DFT calculations were generally found to perform better than ab initio methods. There is a strong wavenumber dependence of radiative forcing in the region of the fundamental C-F vibration, and small differences in wavelength between band positions determined by theory and experiment have a significant impact on the REs. We apply an empirical correction to the theoretical spectra and then test this correction on a number of branched chain and cyclic perfluoroalkanes. We then compute absorption cross sections, REs, and GWPs for an additional set of perfluoroalkenes.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, 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(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Self-assembled molecular complexes and coordination polymers of CdII, hexamine, and monocarboxylates: structural analysis and theoretical studies of supramolecular interactions

Four new cadmium(II) complexes [Cd-2(bz)(4)(H2O)(4)(mu 2-hmt)]center dot Hbz center dot H2O (1), [Cd-3(bz)(6)(H2O)(6)(mu 2-hmt)(2)]center dot 6H(2)O (2), [Cd(pa)(2)(H2O)(mu(2)-hmt)](n) (3), and {[Cd-3(ac)(6)(H2O)(3)(mu(3)-hmt)(2)]center dot 6H(2)O}(n) (4) with hexamine (hmt) and monocarboxylate ions, benzoate (bz), phenylacetate (pa), or acetate (ac) have been synthesized and characterized structurally. Structure determinations reveal that 1 is dinuclear, 2 is trinuclear, 3 is a one-dimensional (1D) infinite chain, and 4 is a two-dimensional (2D) polymer with fused hexagonal rings consisting of Cd-II and hmt. All the Cd-II atoms in the four complexes (except one CdII in 2) possess seven-coordinate pentagonal bipyramidal geometry with the various chelating bidentate carboxylate groups in equatorial sites. One of the CdII ions in 2, a complex that contains two monodentate carboxylates is in a distorted octahedral environment. The bridging mode of hmt is mu 2- in complexes 1-3 but is mu 3- in complex 4. In all complexes, there are significant numbers of H-bonds, C-H/pi, and pi-pi interactions which play crucial roles in forming the supramolecular networks. The importance of the noncovalent interactions in terms of energies and geometries has been analyzed using high level ab initio calculations. The effect of the cadmium coordinated to hmt on the energetic features of the C-H/pi interaction is analyzed. Finally, the interplay between C-H/pi and pi-pi interactions observed in the crystal structure of 3 is also studied.

Proteogenomics and in silico structural and functional annotation of the barley powdery mildew Blumeria graminis f. sp. hordei

Blumeria graminis is an economically important obligate plant-pathogenic fungus, whose entire genome was recently sequenced and manually annotated using ab initio in silico predictions [7]. Employing large scale proteogenomic analysis we are now able to verify independently the existence of proteins predicted by 24% of open reading frame models. We compared the haustoria and sporulating hyphae proteomes and identified 71 proteins exclusively in haustoria, the feeding and effector-delivery organs of the pathogen. These proteins are ‘significantly smaller than the rest of the protein pool and predicted to be secreted. Most do not share any similarities with Swiss–Prot or Trembl entries nor possess any identifiable Pfam domains. We used a novel automated prediction pipeline to model the 3D structures of the proteins, identify putative ligand binding sites and predict regions of intrinsic disorder. This revealed that the protein set found exclusively in haustoria is significantly less disordered than the rest of the identified Blumeria proteins or random (and representative) protein sets generated from the yeast proteome. For most of the haustorial proteins with unknown functions no good templates could be found, from which to generate high quality models. Thus, these unknown proteins present potentially new protein folds that can be specific to the interaction of the pathogen with its host.

Assessment of the consistency of H2O line intensities over the near-infrared using sun-pointing ground-based Fourier transform spectroscopy

We report on the consistency of water vapour line intensities in selected spectral regions between 800–12,000 cm−1 under atmospheric conditions using sun-pointing Fourier transform infrared spectroscopy. Measurements were made across a number of days at both a low and high altitude field site, sampling a relatively moist and relatively dry atmosphere. Our data suggests that across most of the 800–12,000 cm−1 spectral region water vapour line intensities in recent spectral line databases are generally consistent with what was observed. However, we find that HITRAN-2008 water vapour line intensities are systematically lower by up to 20% in the 8000–9200 cm−1 spectral interval relative to other spectral regions. This discrepancy is essentially removed when two new linelists (UCL08, a compilation of linelists and ab-initio calculations, and one based on recent laboratory measurements by Oudot et al. (2010) [10] in the 8000–9200 cm−1 spectral region) are used. This strongly suggests that the H2O line strengths in the HITRAN-2008 database are indeed underestimated in this spectral region and in need of revision. The calculated global-mean clear-sky absorption of solar radiation is increased by about 0.3 W m−2 when using either the UCL08 or Oudot line parameters in the 8000–9200 cm−1 region, instead of HITRAN-2008. We also found that the effect of isotopic fractionation of HDO is evident in the 2500–2900 cm−1 region in the observations.

The molecular structure of propylene sulphide (methylthiirane) by gas-phase electron diffraction and theoretical calculations: a molecule used in MOCVD

Gas-phase electron-diffraction (GED) data together with results from ab initio molecular orbital calculations have been used to determine the structure of propylene sulphide. Values found for the main structural parameters for the molecule are consistent with those obtained from microwave studies and are compared here with those found for similar sulphur containing rings of general formula S(CH2)n (n = 2–5). A high ring strain enthalpy was calculated for propylene sulphide which is consistent with the small C–S–C angle (48.2(6)degrees) and the relatively long C–S bond lengths (ra = 1.831(2) Å). This is thought to account for the ease of ring opening in propylene sulphide observed in MOCVD reactions and the ready polymerisation of the molecule.

Steric effects on uranyl complexation: synthetic, structural, and theoretical studies of carbamoyl pyrazole compounds of the uranyl(VI) ion

New bifunctional pyrazole based ligands of the type [C3HR2N2CONR'] (where R = H or CH3; R' = CH3, C2H5, or (C3H7)-C-i) were prepared and characterized. The coordination chemistry of these ligands with uranyl nitrate and uranyl bis(dibenzoyl methanate) was studied with infrared (IR), H-1 NMR, electrospray-mass spectrometry (ES-MS), elemental analysis, and single crystal X-ray diffraction methods. The structure of compound [UO2(NO3)(2)(C3H3N2CON{C2H5}(2))] (2) shows that the uranium(VI) ion is surrounded by one nitrogen atom and seven oxygen atoms in a hexagonal bipyramidal geometry with the ligand acting as a bidentate chelating ligand and bonds through both the carbamoyl oxygen and pyrazolyl nitrogen atoms. In the structure of [UO2(NO3)(2)(H2O)(2)(C5H7N2CON {C2H5}(2))(2)], (5) the pyrazole figand acts as a second sphere ligand and hydrogen bonds to the water molecules through carbamoyl oxygen and pyrazolyl nitrogen atoms. The structure of [UO2(DBM)(2)C3H3N2CON{C2H5}(2)] (8) (where DBM = C6H5COCHCOC6H5) shows that the pyrazole ligand acts as a monodentate ligand and bonds through the carbamoyl oxygen to the uranyl group. The ES-MS spectra of 2 and 8 show that the ligand is similarly bonded to the metal ion in solution. Ab initio quantum chemical studies show that the steric effect plays the key role in complexation behavior.

The interplay of secondary Hg⋯S, Hg⋯N and Hg⋯π bonding interactions in supramolecular structures of phenylmercury(ii) dithiocarbamates

Three new phenylmercury(II) and one mercury(II) dithiocarbamate complexes viz. PhHg S2CN(PyCH2) Bz (1), PhHg S2CN(PyCH2)CH3 (2), PhHg S2CN(Bz)CH3 (3), and [Hg (NCS2(PyCH2)Bz)(2)] (4) (Py = pyridine; Bz = benzyl) have been synthesized and characterized by elemental analyses, IR, electronic absorption, H-1 and C-13 NMR spectroscopy. The crystal structures of 1, 2 and 3 showed a linear S-Hg-C core at the centre of the molecule, in which the metal atom is bound to the sulfur atom of the dithiocarbamate ligand and a carbon atom of the aromatic ring. In contrast the crystal structure of 4 showed a linear S-Hg-S core at the Hg(II) centre of the molecule. Weak intermolecular Hg center dot center dot center dot N (Py) interactions link molecules into a linear chain in the case of 1, whereas chains of dimers are formed in 2 through intermolecular Hg center dot center dot center dot N (Py) and Hg center dot center dot center dot S interactions. 3 forms a conventional face-to-edge dimeric structure through intermolecular Hg center dot center dot center dot S secondary bonding and 4 forms a linear chain of dimers through face-to-face Hg center dot center dot center dot S secondary bonding. In order to elucidate the nature of these secondary bonding interactions and the electronic absorption spectra of the complexes, ab initio quantum chemical calculations at the MP2 level and density functional theory calculations were carried out for 1-3. Complexes 1 and 2 exhibited photoluminescent properties in the solid state as well as in the solution phase. Studies indicate that Hg center dot center dot center dot S interactions decrease and Hg center dot center dot center dot N interactions increase the chances of photoluminescence in the solid phase

Global warming potentials and radiative efficiencies of halocarbons and related compounds: a comprehensive review

In the mid-1970s it was recognized that, as well as being substances that deplete stratospheric ozone, chlorofluorocarbons (CFCs) were strong greenhouse gases that could have substantial impacts on radiative forcing of climate change. Around a decade later, this group of radiatively active compounds was expanded to include a large number of replacements for ozone-depleting substances such as chlorocarbons, hydrochlorocarbons, hydrochlorofluorocarbons (HCFCs), hydrofluorocarbons (HFCs), perfluorocarbons (PFCs), bromofluorocarbons, and bromochlorofluorocarbons. This paper systematically reviews the published literature concerning the radiative efficiencies (REs) of CFCs, bromofluorocarbons and bromochlorofluorocarbons (halons), HCFCs, HFCs, PFCs, SF6, NF3, and related halogen containing compounds. In addition we provide a comprehensive and self-consistent set of new calculations of REs and global warming potentials (GWPs) for these compounds, mostly employing atmospheric lifetimes taken from the available literature. We also present Global Temperature change Potentials (GTPs) for selected gases. Infrared absorption spectra used in the RE calculations were taken from databases and individual studies, and from experimental and ab initio computational studies. Evaluations of REs and GWPs are presented for more than 200 compounds. Our calculations yield REs significantly (> 5%) different from those in the Intergovernmental Panel on Climate Change Fourth Assessment Report (AR4) for 49 compounds. We present new RE values for more than 100 gases which were not included in AR4. A widely-used simple method to calculate REs and GWPs from absorption spectra and atmospheric lifetimes is assessed and updated. This is the most comprehensive review of the radiative efficiencies and global warming potentials of halogenated compounds performed to date.

The reaction between silylene and ammonia: some gas-phase kinetic and quantum chemical studies

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by 193 nm laser flash photolysis of silacyclopent-3-ene, have been carried out in the presence of ammonia, NH3. Second order kinetics were observed. The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas at each of the three temperatures, 299, 340 and 400 K. The second order rate constants (laser pulse energy of 60 mJ/pulse) fitted the Arrhenius equation: log(k/cm3 molecule-1 s-1) = (-10.37 ± 0.17) + (0.36 ± 1.12 kJ mol-1)/RTln10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range 10-100 Torr, but showed small decreases in value at 3 and 1 Torr. There was also a weak intensity dependence, with rate constants decreasing at laser pulse energies of 30 mJ/pulse. Ab initio calculations at the G3 level of theory, show that SiH2 + NH3 should form an initial adduct (donor-acceptor complex), but that energy barriers are too great for further reaction of the adduct. This implies that SiH2 + NH3 should be a pressure dependent association reaction. The experimental data are inconsistent with this and we conclude that SiH2 decays are better explained by reaction of SiH2 with the amino radical, NH2, formed by photodissociation of NH3 at 193 nm. The mechanism of this previously unstudied reaction is discussed.

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by Heyd-Scuseria-Ernzerhof screened hybrid functional calculations.

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.