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.

Dinuclear copper and zinc complexes of a hexaazamacrocycle containing p-xylyl spacers and bridging anions: theoretical and spectroscopic studies

The hexaazamacrocycle 7,22-dimethyl-3,7,11,18,22,26-hexaazatricyclo[26.2.2.2(13,16)] tetratriaconta-1(30), 13,15,28,31,33- hexaene (Me-2[30] pbz(2)N(6)) was synthesized and characterised by single crystal X-ray diffraction. The macrocycle adopts a conformation with the two aromatic rings almost parallel at a distance of ca. 4.24 Angstrom, but displaced relative to each other by ca. 1.51 Angstrom. The protonation constants of this compound and the stability constants of its complexes with Cu2+ and Zn2+, were determined in water - methanol (9 : 1 v/v) at 25 degreesC with ionic strength 0.10 mol dm(-3) in KCl. The potentiometric and spectroscopic studies (NMR of zinc, cadmium and lead complexes, and EPR of the copper complexes) indicate the formation of only dinuclear complexes. The association constants of the dinuclear copper complex with anions ( thiocyanate, terephthalate and glyphosate) and neutral molecules (1,4-benzenedimethanol, p-xylylenediamine and terephthalic acid) were determined at 20 degreesC in methanol. The structural preferences of this ligand and of its dinuclear copper(II) complex with a variety of bridging ligands were evaluated theoretically by molecular mechanics calculations (MM) and molecular dynamics (MD) using quenching techniques.

Geometrical isomers of [TEAH][Co(L-Se)(2)]center dot xH(2)O: synthesis, structural, spectroscopic and computational studies

Trans-1, [HNEt3][Co-III(L-Se)(2)]center dot H2O and cis-1, [HNEt3][Co-III(L-Se)(2)]center dot 3H(2)O have been synthesized and characterized by single-crystal X-ray studies. The counter ion Et3NH+ plays a crucial role in the crystal packing leading to the formation of two distinctly different supramolecular assemblies in the two complexes. In trans-1, Co-bisphenolate units and triethylamine molecules are arranged in a linear fashion leading to a supramolecular columnar assembly along the crystallographic a-axis. In this assembly, triethylammonium ions are sandwiched between successive Co-bisphenolate units and act as gluing agents joining Co-bisphenolate units on either side through C-H center dot center dot center dot pi interactions. In sharp contrast to trans-1, Co-bisphenolate units and triethylammonium ions in cis-1 are arranged in a helical supramolecular assembly through similar C-H center dot center dot center dot pi interactions along the crystallographic b-axis. The Se center dot center dot center dot Se van der Waals interactions may be responsible for the predominant occurrence of the cis-isomer. The cyclic voltammetric studies showed quasi-reversible waves for the cobalt(III) -> cobalt(II) reductions with E-1/2 = 0.635 and 0.628 V vs. Ag/AgCl for cis-1 (at similar to 5 degrees C) and trans-1 (at similar to 25 degrees C), respectively. DFT calculations show that the trans-form is the thermodynamic product with higher stability than the cis-one, which is consistent with the variable temperature H-1 NMR studies

C-H...O hydrogen bonding in 4-phenyl-benzaldehyde: a comprehensive crystallographic, spectroscopic and computational study

The crystal structure of 4-phenyl-benzaldehyde reveals the presence of a dimer linked by the C=O and C( 9)-H groups of adjacent molecules. In the liquid phase, the presence of C-(HO)-O-... bonded forms is revealed by both vibrational and NMR spectroscopy. A Delta H value of - 8.2 +/- 0.5 kJ mol(-1) for the dimerisation equilibrium is established from the temperature-dependent intensities of the bands assigned to the carbonyl-stretching modes. The NMR data suggest the preferential engagement of the C(2,6)-H and C(10/12)/C(11)-H groups as hydrogen bond donors, instead of the C(9)-H group. While ab initio calculations for the isolated dimers are unable to corroborate these NMR results, the radial distribution functions obtained from molecular dynamics simulations show a preference for C(2,6)-H and C(10/12)/C(11)-(HO)-O-... contacts relative to the C(9)-(HO)-O-... ones.

Solvent-dependent modulation of metal–metal electronic interactions in a dinuclear cyanoruthenate complex: a detailed electrochemical, spectroscopic and computational study

The dinuclear complex [{Ru(CN)4}2(μ-bppz)]4− shows a strongly solvent-dependent metal–metal electronic interaction which allows the mixed-valence state to be switched from class 2 to class 3 by changing solvent from water to CH2Cl2. In CH2Cl2 the separation between the successive Ru(II)/Ru(III) redox couples is 350 mVand the IVCT band (from the UV/Vis/NIR spectroelectrochemistry) is characteristic of a borderline class II/III or class III mixed valence state. In water, the redox separation is only 110 mVand the much broader IVCT transition is characteristic of a class II mixed-valence state. This is consistent with the observation that raising and lowering the energy of the d(π) orbitals in CH2Cl2 or water, respectively, will decrease or increase the energy gap to the LUMO of the bppz bridging ligand, which provides the delocalisation pathway via electron-transfer. IR spectroelectrochemistry could only be carried out successfully in CH2Cl2 and revealed class III mixed-valence behaviour on the fast IR timescale. In contrast to this, time-resolved IR spectroscopy showed that the MLCTexcited state, which is formulated as RuIII(bppz˙−)RuII and can therefore be considered as a mixed-valence Ru(II)/Ru(III) complex with an intermediate bridging radical anion ligand, is localised on the IR timescale with spectroscopically distinct Ru(II) and Ru(III) termini. This is because the necessary electron-transfer via the bppz ligand is more difficult because of the additional electron on bppz˙− which raises the orbital through which electron exchange occurs in energy. DFT calculations reproduce the electronic spectra of the complex in all three Ru(II)/Ru(II), Ru(II)/Ru(III) and Ru(III)/Ru(III) calculations in both water and CH2Cl2 well as long as an explicit allowance is made for the presence of water molecules hydrogen-bonded to the cyanides in the model used. They also reproduce the excited-state IR spectra of both [Ru(CN)4(μ-bppz)]2– and [{Ru(CN)4}2(μ-bppz)]4− very well in both solvents. The reorganization of the water solvent shell indicates a possible dynamical reason for the longer life time of the triplet state in water compared to CH2Cl2.

Bridge-localized HOMO-binding character of divinylanthracene-bridged dinuclear ruthenium carbonyl complexes: spectroscopic, spectroelectrochemical, and computational studies

The electronic properties of four divinylanthracene-bridged diruthenium carbonyl complexes [{RuCl(CO)(PMe3)3}2(μ[BOND]CH[DOUBLE BOND]CHArCH[DOUBLE BOND]CH)] (Ar=9,10-anthracene (1), 1,5-anthracene (2), 2,6-anthracene (3), 1,8-anthracene (4)) obtained by molecular spectroscopic methods (IR, UV/Vis/near-IR, and EPR spectroscopy) and DFT calculations are reported. IR spectroelectrochemical studies have revealed that these complexes are first oxidized at the noninnocent bridging ligand, which is in line with the very small ν(C[TRIPLE BOND]O) wavenumber shift that accompanies this process and also supported by DFT calculations. Because of poor conjugation in complex 1, except oxidized 1+, the electronic absorption spectra of complexes 2+, 3+, and 4+ all display the characteristic near-IR band envelopes that have been deconvoluted into three Gaussian sub-bands. Two of the sub-bands belong mainly to metal-to-ligand charge-transfer (MLCT) transitions according to results from time-dependent DFT calculations. EPR spectroscopy of chemically generated 1+–4+ proves largely ligand-centered spin density, again in accordance with IR spectra and DFT calculations results.

Synthesis, structure and solution chemistry of mixed-ligand oxovanadium(IV) and oxovanadium(V) complexes incorporating tridentate ONO donor hydrazone ligands

In the title family, the ONO donor ligands are the acetylhydrazones of salicylaidehyde (H2L1) and 2-hydroxyacetophenone (H2L2) (general abbreviation, H2L). The reaction of bis(acetylacetonato)oxovanadium(IV) with a mixture of tridentate H2L and a bidentate NN donor [e.g., 2,2'-bipyridine(bpy) or 1,10-phenanthroline(phen), hereafter B] ligands in equimolar ratio afforded the tetravalent complexes of the type [(VO)-O-IV(L)(B)]; complexes (1)-(4) whereas, if B is replaced by 8-hydroxyquinoline(Hhq) (which is a bidentate ON donor ligand), the above reaction mixture yielded the pentavalent complexes of the type [(VO)-O-V(L)(hq)]; complexes (5) and (6). Aerial oxygen is most likely the oxidant (for the oxidation of V-IV -> V-V) in the synthesis of pentavalent complexes (5) and (6). [(VO)-O-IV(L)(B)] complexes are one electron paramagnetic and display axial EPR spectra, while the [(VO)-O-V(L)(hq)] complexes are diamagnetic. The X-ray structure of [(VO)-O-V(L-2)(hq)] (6) indicates that H2L2 ligand is bonded with the vanadium meridionally in a tridentate dinegative fashion through its phenolic-O, enolic-O and imine-N atoms. The general bond length order is: oxo < phenolato < enolato. The V-O (enolato) bond is longer than V-O (phenolato) bond by similar to 0.07 angstrom and is identical with V-O (carboxylate) bond. H-1 NMR spectrum of (6) in CDCl3 solution indicates that the binding nature in the solid state is also retained in solution. Complexes (1)(4) display two ligand-field transitions in the visible region near 820 and 480 nm in DMF solution and exhibit irreversible oxidation peak near +0.60 V versus SCE in DMSO solution, while complexes (5) and (6) exhibit only LMCT band near 535 nm and display quasi-reversible one electron reduction peak near -0.10 V versus SCE in CH2Cl2 solution. The VO3+-VO2+ E-1/2 values shift considerably to more negative values when neutral NN donor is replaced by anionic ON donor species and it also provides better VO3+ binding via phenolato oxygen. For a given bidentate ligand, E-1/2 increases in the order: (L-2)(2-) < (L-1)(2-). (c) 2004 Elsevier B.V. All rights reserved.

Design, synthesis and computational modelling of aromatic tweezer-molecules as models for chain-folding polymer blends

Novel 'tweezer-type' complexes that exploit the interactions between pi-electron-rich pyrenyl groups and pi-electron deficient diimide units have been designed and synthesised. The component molecules leading to complex formation were accessed readily from commercially available starting materials through short and efficient syntheses. Analysis of the resulting complexes, using the visible charge-transfer band, revealed association constants that increased sequentially from 130 to 11,000 M-1 as increasing numbers of pi-pi-stacking interactions were introduced into the systems. Computational modelling was used to analyse the structures of these complexes, revealing low-energy chain-folded conformations for both components, which readily allow close, multiple pi-pi-stacking and hydrogen bonding to be achieved. In this paper, we give details of our initial studies of these complexes and outline how their behaviour could provide a basis for designing self-healing polymer blends for use in adaptive coating systems. (C) 2008 Elsevier Ltd. All rights reserved.

Synthesis, structure and solution chemistry of quaternary oxovanadium(V) complexes incorporating hydrazone ligands

[VIVO(acac)(2)] reacts with an equimolar amount of benzoyl hydrazone of 2-hydroxyacetophenone (H2L1) or 5-chloro-2-hydroxyacetophenone (H2L2) in the presence of excess pyridine (py) in methanol to produce the quaternary [(VO)-O-V(L-1)(OCH3)(py)] (1) and [(VO)-O-V(L-2)(OCH3)(py)] (2) complexes, respectively, while under similar condition, the benzoyl hydrazones of 2-hydroxy-5-methylacetophenone (H2L3) and 2-hydroxy-5-methoxyacetophenone (H2L4) afforded only the methoxy bridged dimeric [(VO)-O-V(L-3/L-4)(OCH3)](2) complexes. The X-ray structural analysis of 1 and 2 indicates that the geometry around the metal is distorted octahedral where the three equatorial positions are occupied by the phenolate-O, enolate-O and the imine-N of the fully deprotonated hydrazone ligand in its enolic form and the fourth one by a methoxide-O atom. An oxo-O and a pyridine-N atom occupy two axial positions. Quaternary complexes exhibit one quasi-reversible one-electron reduction peak near 0.25 V versus SCE in CH2Cl2 and they decompose appreciably to the corresponding methoxy bridged dimeric complex in CDCl3 solution as indicated by their H-1 NMR spectra. These quaternary VO3+ complexes are converted to the corresponding V2O34+-complexes simply on refluxing them in acetone and to the VO2+-complexes on reaction with KOH in methanol. An equimolar amount of 8-hydroxyquinoline (Hhq) converts these quaternary complexes to the ternary [(VO)-O-V(L)(hq)] complexes in CHCl3. (C) 2009 Elsevier B. V. All rights reserved.

Synthesis, structure and solution chemistry of a family of dinuclear hydrazonato-vanadium(V) complexes with [OV(mu-O)VO](4+) core

Dinuclear trioxidic [{VOL}(2)mu-O] (1-4) complexes were synthesized from the reaction of [(VO)-O-IV(acac)(2)] with an equimolar amount of H2L [H2L is the general abbreviation of hydrazone ligands (H2L1-4) derived from the condensation of benzoyl hydrazine with either 2-hydroxyacetophenone or its para substituted derivatives] in acetone or CH2Cl2 or acetonitrile. These V2O3L2 complexes were also obtained from the reaction of VOSO4 with H2L in the presence of two equivalents sodium acetate in aqueous-methanolic (50% V/V) medium and also from the decomposition of [(VO)-O-IV(L)(bipy/phen)] complexes in CH2Cl2 Solution. Black monoclinic crystals of 2 and 4 with C2/c space group were obtained from the reaction of [(VO)-O-IV(acac)(2)], respectively, with H2L2 and H2L4 in acetone in which the respective ligands are bonded meridionally to vanadium in their fully deprotonated enol forms. The V-O bond lengths follow the order: V-O(oxo) < V-O(oxo-bridged) < V-O(phenolate) < V-O(enolate). Complexes (1-4) are diamagnetic exhibiting LMCT transition band near 380 nm in CH2Cl2 solution and they are electroactive displaying a quasi-reversible reduction peak in the 0.14-0.30 V versus SCE region. The and the reduction peak potential (E-p(c)) values show linear relationships with the Hammett constant (sigma) of the substituents in the hydrazone ligands. These dinuclear complexes are converted to the corresponding mononuclear cis dioxo complexes K(H2O)(+)[(VO2)-O-V(L)](-) (5-8) and mixed-ligand [(VO)-O-V(L)(hq)] complexes on reaction, respectively, with two equivalents KOH in methanol and two equivalents 8-hydroxyquinoline (Hhq) in CHCl3. Ascorbic acid reduces the dioxovanadium(V) complexes reversibly under aerobic condition. (C) 2008 Elsevier Ltd. All rights reserved.

Theoretical and experimental investigation of wall confluent jets ventilation and comparison with wall displacement ventilation

This paper reports CFD and experimental results of the characteristics of wall confluent jets in a room. The results presented show the behaviour of wall confluent jets in the form of velocity profiles, the spreading rate of jets on the surface, jets decay, etc. The empirical equations derived are compared with other types of air jets. In addition, the flow in wall confluent jets is compared with the flow in displacement ventilation supply, with regards to the vertical and horizontal spreading on the floor. It is concluded that the jet momentum of wall confluent jets can be more conserved than other jets. Thus, wall confluent jets have a greater spread over the floor than displacement flow.

Theoretical and experimental approaches towards the determination of solute effective diffusivities in foods

The aim of this review paper is to present experimental methodologies and the mathematical approaches used to determine effective diffusivities of solutes in food materials. The paper commences by describing the diffusion phenomena related to solute mass transfer in foods and effective diffusivities. It then focuses on the mathematical formulation for the calculation of effective diffusivities considering different diffusion models based on Fick's second law of diffusion. Finally, experimental considerations for effective diffusivity determination are elucidated primarily based on the acquirement of a series of solute content versus time curves appropriate to the equation model chosen. Different factors contributing to the determination of the effective diffusivities such as the structure of food material, temperature, diffusion solvent, agitation, sampling, concentration and different techniques used are considered. (c) 2005 Elsevier Inc. All rights reserved.

Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.