Moving protons with pendant amines: proton mobility in a nickel catalyst for oxidation of hydrogen

Proton transport is ubiquitous in chemical and biological processes, including the reduction of dioxygen to water, the reduction of CO(2) to formate, and the production/oxidation of hydrogen. In this work we describe intramolecular proton transfer between Ni and positioned pendant amines for the hydrogen oxidation electrocatalyst [Ni(P(Cy)(2)N(Bn)(2)H)(2)](2+) (P(Cy)(2)N(Bn)(2) = 1,5-dibenzyl-3,7-dicyclohexyl-1,5-diaza-3,7-diphosphacyclooctane). Rate constants are determined by variable-temperature one-dimensional NMR techniques and two-dimensional EXSY experiments. Computational studies provide insight into the details of the proton movement and energetics of these complexes. Intramolecular proton exchange processes are observed for two of the three experimentally observable isomers of the doubly protonated Ni(0) complex, [Ni(P(Cy)(2)N(Bn)(2)H)(2)](2+), which have N-H bonds but no Ni-H bonds. For these two isomers, with pendant amines positioned endo to the Ni, the rate constants for proton exchange range from 10(4) to 10(5) s(-1) at 25 °C, depending on isomer and solvent. No exchange is observed for protons on pendant amines positioned exo to the Ni. Analysis of the exchange as a function of temperature provides a barrier for proton exchange of ΔG(?) = 11-12 kcal/mol for both isomers, with little dependence on solvent. Density functional theory calculations and molecular dynamics simulations support the experimental observations, suggesting metal-mediated intramolecular proton transfers between nitrogen atoms, with chair-to-boat isomerizations as the rate-limiting steps. Because of the fast rate of proton movement, this catalyst may be considered a metal center surrounded by a cloud of exchanging protons. The high intramolecular proton mobility provides information directly pertinent to the ability of pendant amines to accelerate proton transfers during catalysis of hydrogen oxidation. These results may also have broader implications for proton movement in homogeneous catalysts and enzymes in general, with specific implications for the proton channel in the Ni-Fe hydrogenase enzyme.     

生色团连接的苯骈三氮唑衍生物的激发态分子内质子转移

用从头算和密度泛函理论研究了对硝基二苯乙烯作为生色团连接的2-(2-羟基-苯基)-苯骈三氮唑的衍生物2-羟基-5-[对硝基-二苯乙烯基-氧亚甲基]-苯基-(2H-苯骈三氮唑)(C1)和4′-硝基-3,4-二[2-羟基-(2H-苯骈三氮唑)-苄氧基]-二苯乙烯(C2)发生激发态分子内质子转移(ESIPT)的可能性.系统研究了C1和C2发生ESIPT的互变异构体的基态与激发态的性质变化,包括相关的键长、键角等结构参数,Mulliken电荷和偶极矩,前线轨道以及势能曲线.计算结果表明,对于C1来讲,酮式(keto)的基态(K)不存在稳定结构,因此发生基态分子内质子转移(GSIPT)可能性很小.酮式的激发态(K*)的氢键强度要远强于烯醇式(enol)的激发态(E*)的氢键强度.分子在光致激发后,质子供体所带负电荷减小而质子受体所带负电荷增加.在K*,HOMO→LUMO的电子跃迁导致电子密度从"酚环"向质子化杂环转移.E*→K*跃迁只需要克服较小的能垒(约41 kJ.mol-1).计算结果表明C1发生ESIPT的可能性很大.C2由于具有高能量,其具有基态的单质子转移特征的异构体EK(同时含烯醇E与酮K结构)、具有基态的双质子转移特征的异构体2K(含有双酮结构),以及具有双酮结构特征的激发态2K*均无法获得它们的稳定结构,因此,基态分子内单或双质子转移和激发态分子内双重质子转移发生的可能性极小.然而,由于双烯醇式的激发态(2E*)和EK的激发态(EK*)存在稳定结构,且2E*→EK*跃迁具有低能垒,因此C2有可能发生激发态分子内单重质子转移.本文进一步计算了两个分子的紫外-可见吸收光谱与荧光发射光谱,获得了具有较大斯托克位移的ESIPT的荧光发射峰.     

Dimerization of the keto tautomer of acetohydroxamic acid-infrared matrix isolation and theoretical study

Dimerization of the keto tautomer of acetohydroxamic acid has been studied using FTIR matrix isolation spectroscopy and DFT(B3LYP)/6-31+G(d,p) calculations. Analysis of CH3CONHOH/Ar matrix spectra indicates formation of two dimers in which two intramolecular CO...HON bonds within two interacting acetohydroxamic acid molecules are retained. A chain dimer I is stabilized by the intermolecular CO...HN hydrogen bond, whereas the cyclic dimer II is stabilized by two intermolecular NH...O(H)N bonds. Twelve vibrations were identified for dimer I and six vibrations for dimer II; the observed frequency shifts show a good agreement with the calculated ones for the structures I and II. Both dimers have comparable binding energies (DeltaE(ZPE)(CP)I, II=-7.02, -6.34 kcal mol-1) being less stable than calculated structures III and IV (DeltaE(ZPE)(CP)III, IV=-9.50, -8.87 kcal mol-1) in which one or two intramolecular hydrogen bonds are disrupted. In the most stable 10-membered cyclic dimer III, two intermolecular CO...HON hydrogen bonds are formed at expense of intramolecular hydrogen bonds of the same type. The formation of the less stable (AHA)2 dimers in the studied matrixes indicates that the formation of (AHA)2 is kinetically and not thermodynamically controlled.     

Probing the low-barrier hydrogen bond in hydrogen maleate in the gas phase: a photoelectron spectroscopy and ab initio study

The strength of the low-barrier hydrogen bond in hydrogen maleate in the gas phase was investigated by low-temperature photoelectron spectroscopy and ab initio calculations. Photoelectron spectra of maleic and fumaric acid monoanions (cis-/trans-HO(2)CCH=CHCO(2)(-)) were obtained at low temperatures and at 193 nm photon energy. Vibrational structure was observed for trans-HO(2)CCH=CHCO(2)(-) due to the OCO bending modes; however, cis-HO(2)CCH=CHCO(2)(-) yielded a broad and featureless spectrum. The electron binding energy of cis-HO(2)CCH=CHCO(2)(-) is about 1 eV blue-shifted relative to trans-HO(2)CCH=CHCO(2)(-) due to the formation of intramolecular hydrogen bond in the cis-isomer. Theoretical calculations (CCSD(T)/ aug-cc-pVTZ and B3LYP/aug-cc-pVTZ) were carried out to estimate the strength of the intramolecular hydrogen bond in cis-HO(2)CCH=CHCO(2)(-). Combining experimental and theoretical calculations yields an estimate of 21.5 +/- 2.0 kcal/mol for the intramolecular hydrogen bond strength in hydrogen maleate.     

Potential-energy and free-energy surfaces of glycyl-phenylalanyl-alanine (GFA) tripeptide: experiment and theory

The free-energy surface (FES) of glycyl-phenylalanyl-alanine (GFA) tripeptide was explored by molecular dynamics (MD) simulations in combination with high-level correlated ab initio quantum chemical calculations and metadynamics. Both the MD and metadynamics employed the tight-binding DFT-D method instead of the AMBER force field, which yielded inaccurate results. We classified the minima localised in the FESs as follows: a) the backbone-conformational arrangement; and b) the existence of a COOH...OC intramolecular H-bond (families CO(2)H(free) and CO(2)H(bonded)). Comparison with experimental results showed that the most stable minima in the FES correspond to the experimentally observed structures. Remarkably, however, we did not observe experimentally the CO(2)H(bonded) family (also predicted by metadynamics), although its stability is comparable to that of the CO(2)H(free) structures. This fact was explained by the former's short excited-state lifetime. We also carried out ab initio calculations using DFT-D and the M06-2X functional. The importance of the dispersion energy in stabilising peptide conformers is well reflected by our pioneer analysis using the DFT-SAPT method to explore the nature of the backbone/side-chain interactions.     

Azulene-to-naphthalene rearrangement: the Car-Parrinello metadynamics method explores various reaction mechanisms.

We studied the thermal intramolecular and radical rearrangement of azulene to naphthalene by employing a novel metadynamics method based on Car-Parrinello molecular dynamics. We demonstrate that relatively short simulations can provide us with several possible reaction mechanisms for the rearrangement. We show that different choices of the collective coordinates can steer the reaction along different pathways, thus offering the possibility of choosing the most probable mechanism. We consider herein three intramolecular mechanisms and two radical pathways. We found the norcaradiene pathway to be the preferable intramolecular mechanism, whereas the spiran mechanism is the favored radical route. We obtained high activation energies for all the intramolecular pathways (81.5-98.6 kcal mol(-1)), whereas the radical routes have activation energies of 24-39 kcal mol(-1). The calculations have also resulted in elementary steps and intermediates not yet considered. A few attractive features of the metadynamics method in studying chemical reactions are pointed out.     

Density functional study on the reaction mechanism of proton transfer in 2-pyridone: effect of hydration and self-association

The proton-transfer mechanism in the isolated, mono, dehydrated forms and dimers of 2-pyridone and the effect of hydration or self-assistance on the transition state structures corresponding to proton transfer from the keto form to the enol form have been investigated using B3LYP and BH-LYP hybrid density functional methods at the 6-311++G (2d, 2p) basis set level. The barrier heights for both H2O-assisted and self-assisted reactions are significantly lower than that of the bare tautomerization reaction from 2-pyridone to 2-hydroxypyridine, implying the importance of the superior catalytic effect of H2O and (H2O)2 and the important role of 2-pyridone itself for the intramolecular proton transfer. Long-range solvent effects have also been taken into account by using the continuum model (Onsager model and polarizable continuum model (PCM)) of water. The tautomerization energies and the potential energy barriers are increased both for the water-assisted and for the self-assisted reaction because of the bulk solvent, which imply that the tautomerization of PY becomes less favorable in the polar solvent.     

Cluster phase chemistry: collisions of vibrationally excited cationic dicarboxylic acid clusters with water molecules initiate dissociation of cluster components

A homologous series of cationic gas-phase clusters of dicarboxylic acids (oxalic acid, malonic acid, succinic acid, glutaric acid, and adipic acid) generated via electrospray ionization (ESI) are investigated using collision-induced dissociation (CID). Singly charged cationic clusters with the composition (Na(+))(2n+1)(dicarboxylate(2-))(n), where n = 1-5, are observed as major gas-phase species. Significant abundances of singly charged sodiated hydrogen dicarboxylate clusters with the composition (Na(+))(2n)(dicarboxylate(2-))(n)(H+), where n = 1-6, are observed with oxalic acid, malonic acid, and succinic acid. Isolation of the clusters followed by CID results mainly in sequential loss of disodium dicarboxylate moieties for the clusters of succinic acid, glutaric acid, and adipic acid. However, the dimer of sodiated hydrogen succinate, all malonate clusters, and all oxalate clusters, with the exception of the dimer, exhibit complex chemical reactions initiated by the collision of vibrationally excited clusters with water molecules. Generally, water molecules serve as proton donors for reacting dicarboxylate anions in the cluster, initiating dissociation pathways such as the decomposition of the malonate ion to yield an acetate ion and CO(2). The reactivity of several mixed dicarboxylate clusters is also reported. For example, malonate anion is shown to be more reactive than oxalate anion for decarboxylation when both are present in a cluster. The energetics of several representative cluster phase reactions are evaluated using computational modeling. The present results for cationic clusters are compared and contrasted to earlier studies of anionic sodiated dicarboxylic acid clusters.     

Structure and dynamics of binary and ternary lanthanide(III) and actinide(III) tris[4,4,4-trifluoro-1-(2-thienyl)-1,3-butanedione] (TTA)-tributylphosphate (TBP) complexes. Part 3, the structure, thermodynamics and reaction mechanisms of 8- and 9-coordinated binary and ternary Y-TTA-TBP complexes studied by quantum chemical methods

Possible mechanisms for intermolecular exchange between coordinated and solvent water in the complexes Y(TTA)(3)(OH(2))(2) and Y(TTA)(3)(TBP)(OH(2)) and intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(OH(2))(HTTA) and Y(TTA)(3)(TBP)(HTTA) have been investigated using ab initio quantum chemical methods. The calculations comprise both structures and energies of isomers, intermediates and transition states. Based on these data and experimental NMR data (Part 2) we have suggested intimate reaction mechanisms for water exchange, intramolecular exchange between structure isomers and intermolecular exchange between free HTTA and coordinated TTA. A large number of isomers are possible for the complexes investigated, but only some of them have been investigated, in all of them the most stable geometry is a more or less distorted square anti-prism or bicapped trigonal prism; the energy differences between the various isomers are in general small, less than 10 kJ mol(-1). 9-coordinated intermediates play an important role in all reactions. Y(TTA)(3)(OH(2))(3) has three non-equivalent water ligands that can participate in ligand exchange reactions. The fastest of these exchanging sites has a QM activation energy of 18.1 kJ mol(-1), in good agreement with the experimental activation enthalpy of 19.6 kJ mol(-1). The mechanism for the intramolecular exchange between structure isomers in Y(TTA)(3)(OH(2))(2) involves the opening of a TTA-ring as the rate determining step as suggested by the good agreement between the QM activation energy and the experimental activation enthalpy 47.8 and 58.3 J mol(-1), respectively. The mechanism for the intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(HTTA) and Y(TTA)(3)(TBP)(HTTA) involves the opening of the intramolecular hydrogen bond in coordinated HTTA followed by proton transfer to coordinated TTA. This mechanism is supported by the good agreement between experimental activation enthalpies (within parenthesis) and calculated activation energies 68.7 (71.8) and 35.3 (38.8) kJ mol(-1). The main reason for the difference between the two systems is the much lower energy required to open the intramolecular hydrogen bond in the latter. The accuracy of the QM methods and chemical models used is discussed.     

C-C coupling reactions in the coordination sphere of rhodium(I) and rhodium(III): New routes for the di- and trimerization of terminal alkynes

The alkynyl(vinylidene)rhodium(I) complexes trans-[Rh(C[triple bond, length as m-dash]CR)(=C=CHR)(PiPr3)2] 2, 5, 6 react with CO by migratory insertion to give stereoselectively the butenynyl compounds trans-[Rh{eta1-(Z)-C(=CHR)C[triple bond, length as m-dash]CR}(CO)(PiPr3)2](Z)-7-9, of which (Z)-7 (R=Ph) and (Z)-8 (R=tBu) rearrange upon heating or UV irradiation to the (E) isomers. Similarly, trans-[Rh{eta1-C(=CH2)C[triple bond, length as m-dash]CPh}(CO)(PiPr3)2] 12 and trans-[Rh{eta1-(Z)-C(=CHCO2Me)C[triple bond, length as m-dash]CR}(CO)(PiPr3)2](Z)-15, (Z)-16 have been prepared. At room temperature, the corresponding "non-substituted" derivative trans-[Rh{eta1-C(=CH2)C[triple bond, length as m-dash]CH}(CO)(PiPr3)2] 18 is in equilibrium with the butatrienyl isomer trans-[Rh(eta1-CH=]C=C=CH2)(CO)(PiPr3)2] 19 that rearranges photochemically to the alkynyl complex trans-[Rh(C[triple bond, length as m-dash]CCH=CH2)(CO)(PiPr3)2] 20. Reactions of (Z)-7, (E)-7, (Z)-8 and (E)-8 with carboxylic acids R'CO2H (R'=CH3, CF3) yield either the butenyne (Z)- and/or (E)-RC[triple bond, length as m-dash]CCH=CHR or a mixture of the butenyne and the isomeric butatriene, the ratio of which depends on both R and R'. Treatment of 2 (R=Ph) with HCl at -40 degrees C affords five-coordinate [RhCl(C[triple bond, length as m-dash]CPh){(Z)-CH=CHPh}(PiPr3)2] 23, which at room temperature reacts by C-C coupling to give trans-[RhCl{eta2-(Z)-PhC[triple bond, length as m-dash]CCH=CHPh}(PiPr3)2](Z)-21. The related compound trans-[RhCl(eta2-HC[triple bond, length as m-dash]CCH=CH2)(PiPr3)2] 27, prepared from trans-[Rh(C[triple bond, length as m-dash]CH)(=C=CH2)(PiPr3)2] 17 and HCl, rearranges to the vinylvinylidene isomer trans-[RhCl(=C=CHCH=CH2)(PiPr3)2] 28. While stepwise reaction of 2with CF3CO2H yields, via alkynyl(vinyl)rhodium(III) intermediates (Z)-29 and (E)-29, the alkyne complexes trans-[Rh(kappa1-O2CCF3)(eta2-PhC[triple bond, length as m-dash]CCH=CHPh)(PiPr3)2](Z)-30 and (E)-30, from 2 and CH3CO2H the acetato derivative [Rh(kappa2-O2CCH3)(PiPr3)2] 33 and (Z)-PhC[triple bond, length as m-dash]CCH=]CHPh are obtained. From 6 (R=CO2Me) and HCl or HC[triple bond, length as m-dash]CCO2Me the chelate complexes [RhX(C[triple bond, length as m-dash]CCO2Me){kappa2(C,O)-CH=CHC(OMe)=O}(PiPr3)2] 34 (X=Cl) and 35 (X=C[triple bond, length as m-dash]CCO2Me) have been prepared. In contrast to the reactions of [Rh(kappa2-O2CCH3)(C[triple bond, length as m-dash]CE)(CH=CHE)(PiPr3)2] 37(E=CO2Me) with chloride sources which give, via intramolecular C-C coupling, four-coordinate trans-[RhCl{eta2-(E)-EC[triple bond, length as m-dash]CCH=CHE}(PiPr3)2](E)-36, treatment of 37with HC[triple bond, length as m-dash]CE affords, via insertion of the alkyne into the rhodium-vinyl bond, six-coordinate [Rh(kappa2-O2CCH3)(C[triple bond, length as m-dash]CE){eta1-(E,E)-C(=CHE)CH=CHE}(PiPr3)2] 38. The latter reacts with MgCl2 to yield trans-[RhCl{eta2-(E,E)-EC[triple bond, length as m-dash]CC(=CHE)CH=CHE}(PiPr3)2] 39, which, in the presence of CO, generates the substituted hexadienyne (E,E)-EC[triple bond, length as m-dash]CC(=CHE)CH=CHE 40.     

Proton transfer to nickel-thiolate complexes. 2. Rate-limiting intramolecular proton transfer in the reactions of [Ni(SC(6)H(4)R-4)(PhP[CH(2)CH(2)PPh(2)](2))](+) (R = NO(2), Cl, H, Me, or MeO)

The protonation of [Ni(SC(6)H(4)R-4)(triphos)](+) (triphos = PhP[CH(2)CH(2)PPh(2)](2); R = NO(2), Cl, H, Me, or MeO) by [lutH](+) (lut = 2,6-dimethylpyridine) to form [Ni(S(H)C(6)H(4)R-4)(triphos)](2+) is an equilibrium reaction in MeCN. Kinetic studies, using stopped-flow spectrophotometry, reveal that the reactions occur by a two-step mechanism. Initially, [lutH](+) rapidly binds to the complex (K(2)(R)) in an interaction which probably involves hydrogen-bonding of the acid to the sulfur. Subsequent intramolecular proton transfer from [lutH](+) to sulfur (k(3)(R)) is slow because of both electronic and steric factors. The X-ray crystal structures of [Ni(SC(6)H(4)R-4)(triphos)](+) (R = NO(2), H, Me, or MeO) show that all are best described as square-planar complexes, with the phenyl substituents of the triphos ligand presenting an appreciable barrier to the approach of the sterically demanding [lutH](+) to the sulfur. The kinetic characteristics of the intramolecular proton transfer from [lutH](+) to sulfur have been investigated. The rate of intramolecular proton transfer exhibits a nonlinear dependence on Hammett sigma(+), with both electron-releasing and electron-withdrawing 4-R-substituents on the coordinated thiolate facilitating the rate of proton transfer (NO(2) > Cl > H > Me < MeO). The rate constants for intramolecular proton transfer correlate well with the calculated electron density of the sulfur. The temperature dependence of the rate of the intramolecular proton transfer reactions shows that deltaH() is small but increases as the 4-R-substituent becomes more electron-withdrawing [deltaH = 4.1 (MeO), 6.9 (Me), 11.4 kcal mol(-)(1) (NO(2))], while DeltaS() becomes progressively less negative [deltaS = -50.1 (MeO), -41.2 (Me), -16.4 (NO(2)) cal K(-)(1) mol(-)(1)]. Studies with [lutD](+) show that the rate of intramolecular proton transfer varies with the 4-R-substituent [(k(3)(NO)2)(H)/(k(3)(NO)2)(D) = 0.39; (k(3)(Cl))(H)/(k(3)(Cl))(D) = 0.88; (k(3)(Me))(H)/(k(3)(Me))(D) = 1.3; (k(3)(MeO))(H)/(k(3)(MeO))(D) = 1.2].     

2-(2-巯苯基)苯并噁唑分子内质子转移取代基电子效应的密度泛函理论研究

在B3LYP/6-31G(d,p)和TDB3LYP/6-31++G(d,p)//CIS/6-31G(d,p)水平上研究了2-(2-巯苯基)苯并噁唑及其衍生物基态和激发态分子内质子转移现象,并探讨取代基电子效应对分子内质子转移的影响,研究结果表明,在基态时,硫醇式异构体为优势构象,供电子取代基使基态分子内正向质子转移能垒(烯醇式→酮式)升高;而吸电子取代基则可降低能垒,有利于基态分子内质子转移并有助于硫酮式异构体的稳定.在激发态时,硫酮式结构为优势构象,所研究的2-(2-巯苯基)苯并噁唑化合物及衍生物均可以发生无能垒或低能垒(≤1.5kJ/mol)的激发态分子内质子转移.巯苯基部分是激发态失活的主要活性部分,供电子基团有利于激发态的质子转移,吸电子基团使激发态跃迁困难,不利于激发态的质子转移.     

Coordination sphere flexibility of active-site models for Fe-only hydrogenase: studies in intra- and intermolecular diatomic ligand exchange.

A series of dinuclear complexes, (mu-SRS)Fe(2)(CO)(6) (R = -CH(2)CH(2)-, -CH(2)CH(2)CH(2)-, -CH(2)-C(6)H(4)-CH(2)-; edt, pdt, and o-xyldt, respectively) has been examined for specific characteristics that might relate to structural similarity with the active site of Fe-only hydrogenases. Variable-temperature proton NMR studies display the fluxionality of the iron-dithiocyclohexane unit in (mu-pdt)Fe(2)(CO)(6) while in the (mu-o-xyldt)Fe(2)(CO)(6) compound, the bridge is fixed. Temperature-dependent (13)C NMR spectral studies establish intramolecular CO site exchange localized on discrete Fe(CO)(3) units in all complexes, which is influenced by steric effects of the mu-SRS unit. Kinetic studies of intermolecular CO/CN(-) ligand-exchange reactions establish associative or I(a) mechanisms in sequential steps to form the dicyano dianion, (mu-SRS)[Fe(CO)(2)(CN)](2)(=) with 100% selectivity. Theoretical calculations (DFT) of transition states in the intramolecular site-exchange processes lead to a rationale for the interesting cooperativity in the CN(-)/CO intermolecular ligand-exchange process. The hinge motion of the three light atom S-to-S bridge is related to a possible heterolytic H(2) activation/production process in the enzyme.     

Dinuclear nickel complexes with a Ni(2)O(2) core: a structural and magnetic study

The acetylacetonate complexes [Ni(2)L(1)(acac)(MeOH)] x H(2)O, 1 x H(2)O and [Ni(2)L(3)(acac)(MeOH)] x 1.5H(2)O, 2 x 1.5H(2)O (H(3)L(1) = (2-(2-hydroxyphenyl)-1,3-bis[4-(2-hydroxyphenyl)-3-azabut-3-enyl]-1,3-imidazolidine and H(3)L(3) = (2-(5-bromo-2-hydroxyphenyl)-1,3-bis[4-(5-bromo-2-hydroxyphenyl)-3-azabut-3-enyl]-1,3-imidazolidine) were prepared and fully characterised. Their crystal structures show that they are dinuclear complexes, extended into chains by hydrogen bond interactions. These compounds were used as starting materials for the isolation of the corresponding [Ni(2)HL(x)(o-O(2)CC(6)H(4)CO(2))(H(2)O)] x n MeOH and [Ni(2)HL(x)(O(2)CCH(2)CO(2))(H(2)O)]x nH(2)O dicarboxylate complexes (x = 1, 3; n = 1-3). The crystal structures of [Ni(2)HL(1)(o-O(2)CC(6)H(4)CO(2))(H(2)O)] x MeOH, 3 x MeOH, [Ni(2)HL(3)(o-O(2)CC(6)H(4)CO(2))(H(2)O)] x 3 MeOH, 4 x 3 MeOH and [Ni(2)HL(1)(O(2)CCH(2)CO(2))(H(2)O)] x 2.5H(2)O x 0.25 MeOH x MeCN, 5 x 2.5H(2)O x 0.25 MeOH x MeCN, were solved. Complexes 3-5 show dinuclear [Ni(2)HL(x)(dicarboxylate)(H(2)O)] units, expanded through hydrogen bonds that involve carboxylate and water ligands, as well as solvate molecules. The variable temperature magnetic susceptibilities of all the complexes show an intramolecular ferromagnetic coupling between the Ni(II) ions, which is attempted to be rationalized by comparison with previous results and in the light of molecular orbital treatment. Magnetisation measurements are in accord with a S = 2 ground state in all cases.     

Mechanisms of electron-transfer oxidation of NADH analogues and chemiluminescence. Detection of the keto and enol radical cations

The radical cation of an NADH analogue (BNAH: 1-benzyl-1,4-dihydronicotinamide) has been successfully detected as the transient absorption and ESR spectra in the thermal electron transfer from BNAH to Fe(bpy)(3)(3+) (bpy = 2,2'-bipyridine) and Ru(bpy)(3)(3+). The ESR spectra of the radical cations of BNAH and the dideuterated compound (BNAH-4,4'-d(2)) indicate that the observed radical cation is the keto form rather than the enol form in the tautomerization. The deprotonation rate and the kinetic isotope effects of the keto form of BNAH(*)(+) were determined from the kinetic analysis of the electron-transfer reactions. In the case of electron transfer from BNAH to Ru(bpy)(3)(3+), the chemiluminescence due to Ru(bpy)(3)(2+) was observed in the second electron-transfer step from BNA(*), produced by the deprotonation of the keto form of BNAH(*)(+), to Ru(bpy)(3)(3+). The observation of chemiluminescence due to Ru(bpy)(3)(2+) provides compelling evidence that the Marcus inverted region is observed even for such an intermolecular electron-transfer reaction. When BNAH is replaced by 4-tert-butylated BNAH (4-t-BuBNAH), no chemiluminescence due to Ru(bpy)(3)(2+) has been observed in the electron transfer from 4-t-BuBNAH to Ru(bpy)(3)(3+). This is ascribed to the facile C-C bond cleavage in 4-t-BuBNAH(*)(+). In the laser flash photolysis of a deaerated MeCN solution of BNAH and CHBr(3), the transient absorption spectrum of the enol form of BNAH(*)(+) was detected instead of the keto form of BNAH(*)(+), and the enol form was tautomerized to the keto form. The rate of intramolecular proton transfer in the enol form to produce the keto form of BNAH(*)(+) was determined from the decay of the absorption band due to the enol form and the rise in the absorption band due to the keto form. The kinetic isotope effects were observed for the intramolecular proton-transfer process in the keto form to produce the enol form.     

Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution.

We have studied the keto-enol interconversion of acetone to understand the mechanism of tautomerism relevant to numerous organic and biochemical processes. Applying the ab initio metadynamics method, we simulated the keto-enol isomerism both in the gas phase and in the presence of water. For the gas-phase intramolecular mechanism we show that no other hydrogen-transfer reactions can compete with the simple keto-enol tautomerism. We obtain an intermolecular mechanism and remarkable participation of water when acetone is solvated by neutral water. The simulations reveal that C deprotonation is the kinetic bottleneck of the keto-enol transformation, in agreement with experimental observations. The most interesting finding is the formation of short H-bonded chains of water molecules that provide the route for proton transfer from the carbon to the oxygen atom of acetone. The mechanistic picture that emerged from the present study involves proton migration and emphasizes the importance of active solvent participation in tautomeric interconversion.     

Dissociation of carbonic acid: gas phase energetics and mechanism from ab initio metadynamics simulations

A comprehensive metadynamics study of the energetics, stability, conformational changes, and mechanism of dissociation of gas phase carbonic acid, H2CO3, yields significant new insight into these reactions. The equilibrium geometries, vibrational frequencies, and conformer energies calculated using the density functional theory are in good agreement with the previous theoretical predictions. At 315 K, the cis-cis conformer has a very short life time and transforms easily to the cis-trans conformer through a change in the O=C-O-H dihedral angle. The energy difference between the trans-trans and cis-trans conformers is very small (approximately 1 kcal/mol), but the trans-trans conformer is resistant to dissociation to carbon dioxide and water. The cis-trans conformer has a relatively short path for one of its hydroxyl groups to accept the proton from the other end of the molecule, resulting in a lower activation barrier for dissociation. Comparison of the free and potential energies of dissociation shows that the entropic contribution to the dissociation energy is less than 10%. The potential energy barrier for dissociation of H2CO3 to CO2 and H2O from the metadynamics calculations is 5-6 kcal/mol lower than in previous 0 K studies, possibly due to a combination of a finite temperature and more efficient sampling of the energy landscape in the metadynamics calculations. Gas phase carbonic acid dissociation is triggered by the dehydroxylation of one of the hydroxyl groups, which reorients as it approaches the proton on the other end of the molecule, thus facilitating a favorable H-O-H angle for the formation of a product H2O molecule. The major atomic reorganization of the other part of the molecule is a gradual straightening of the O=C=O bond. The metadynamics results provide a basis for future simulation of the more challenging carbonic acid-water system.     

Theoretical study of intramolecular interaction energies during dynamics simulations of oligopeptides by the fragment molecular orbital-Hamiltonian algorithm method

We analyzed the interaction energies between residues (fragments) in an oligopeptide occurring during dynamic simulations by using the fragment molecular orbital-Hamiltonian algorithm (FMO-HA) method, an ab initio MO-molecular dynamics technique. The FMO method enables not only calculation of large molecules based on ab initio MO but also easy evaluation of interfragment interaction energies. The glycine pentamer [(Gly)(5)] and decamer [(Gly)(10)] were divided into five and ten fragments, respectively. alpha-helix structures of (Gly)(5) and (Gly)(10) were stabilized by attractive interaction energies owing to intramolecular hydrogen bonds between fragments n and n+3 (and n+4), and beta-strand structures were characterized by repulsive interaction energies between fragments n and n+2. We analyzed interfragment interaction energies during dynamics simulations as the peptides' geometries changed from alpha helix to beta strand. Intramolecular hydrogen bonds between fragments 2-4 and 2-5 control the geometrical preference of (Gly)(5) for the beta-strand structure. The pitch of one turn of the alpha helix of (Gly)(10) elongated and thus weakened during dynamics due to a shifting of the intramolecular hydrogen bonds, and enabled the beta-strand structure to form. Changes in interaction energies due to the intramolecular hydrogen bonds controlled the tendency toward alpha-helix or beta-strand structure of (Gly)(5) and (Gly)(10). Evaluation of interfragment interaction energies during dynamics simulations thus enabled detailed analysis of the process of the geometrical changes occurring in oligopeptides.     

Vinyl and carbene ruthenium(II) complexes from hydridoruthenium(II) precursors

The reactions of the hydrido compounds [RuHCl(CO)(L)2][L = PiPr3 (1), PCy3 (2)] with HC(triple bond)CR (R = H, Ph, tBu) afforded by insertion of the alkyne into the Ru-H bond the corresponding vinyl complexes [RuCl(CHCHR)(CO)(L)2], 3-8, which upon protonation with HBF4 gave the cationic five-coordinated ruthenium carbenes [RuCl(CHCH2R)(CO)(L)2]BF4, 9-14. Subsequent reactions of the carbene complexes with PR3(R = Me, iPr) and CH3CN led either to deprotonation and re-generation of the vinyl compounds or to cleavage of the ruthenium-carbene bond and the formation of the six-coordinated complexes [RuCl(CO)(CH3CN)2(PiPr3)2]BF4, 17, and [RuH(CO)(CH3CN)2(PiPr3)2]X, 18a,b. The acetato derivative [RuH(2-O2CCH3)(CO)(PCy3)2], 19, also reacted with acetylene and phenylacetylene by insertion to yield the related vinyl complexes [Ru(CHCHR)(kappa2-O2CCH3)(CO)(PCy3)2], 20, 21, of which that with R = H was protonated with HBF4 to yield the corresponding cationic ruthenium carbene 22. With [RuHCl(H2)(PCy3)2], 25, as the starting material, the five-coordinated chloro(hydrido)ruthenium(II) compounds [RuHCl(PCy3)(dppf)], 26(dppf = [Fe(eta5-C5H4PPh2)2]), [RuHCl[Sb(CH2Ph)3](PCy3)2], 27, and [RuHCl(CH3CN)(PCy3)2], 30, were prepared. The reactions of 27 with HCCR (R = H, Ph) gave the hydrido(vinylidene) complexes [RuHCl(CCHR)(PCy3)2], 28 and 29, whereas treatment of 30 with HC(triple bond)CPh afforded the vinyl compound [RuCl(CHCHPh)(CH3CN)(PCy3)2], 31. The molecular structures of 11(R = tBu, L = PiPr3) and 26 were determined crystallographically.