DFT study on the mechanism of water-assisted dihydrogen elimination in group 6 octahedral metal hydride complexes

The reaction of water with octahedral bis-, tris- and tetrakis-(phosphine)tungsten, (phosphine)molybdenum and (phosphine)chromium complexes has been studied using B3LYP/def2-TZVPP level of DFT to elucidate dissociative, associative and hydride migratory insertion mechanisms for hydrogen elimination. In the dissociative mechanism, phosphine dissociation requires 19.3-28.5 kcal mol(-1) of energy. The phosphine-water ligand exchange is endergonic due to poor coordination ability of water to group 6 metals (binding energy 8.8-15.5 kcal mol(-1)). The ligand exchange leads to intermolecular M-HH(2)O dihydrogen interaction and facilitates dihydrogen elimination (G(act) = 6.8-15.5 kcal mol(-1)). In the associative mechanism, a water molecule in the first solvation shell interacts with the M-H bond through a dihydrogen bond (interaction energy 2.7-4.0 kcal mol(-1)) and leads to the elimination of H(2) by forming a hydroxide complex. Compared to the dissociative mechanism, G(act) of associative mechanisms are ～22 kcal mol(-1) higher. In the hydride migratory insertion mechanism, the hydride ligand shifts to the CO ligand (G(act) = 25.4-30.4 kcal mol(-1)) to afford a formyl complex and subsequently the H-H bond coupling occurs between formyl and water ligand (G(act) = 2.8-4.4 kcal mol(-1)). In many cases, the migratory insertion mechanism can simultaneously operate with the dissociative mechanism as a minor pathway, whereas owing to high G(act) value, the associative mechanism can be described as inactive in the reaction. The general argument that dihydrogen elimination is preceded by the formation of a dihydrogen intermediate is not applicable for the systems studied herein as the most favoured dissociative mechanism does not pass through such an intermediate. On the other hand, irrespective of the mechanisms, dihydrogen elimination invariably occurs with the formation of a dihydrogen bonded transition state. Our results also suggest that group 6 octahedral metal hydride complexes are attractive targets for the design of water splitting reactions.     

The first example of a nitrile hydratase model complex that reversibly binds nitriles

Nitrile hydratase (NHase) is an iron-containing metalloenzyme that converts nitriles to amides. The mechanism by which this biochemical reaction occurs is unknown. One mechanism that has been proposed involves nucleophilic attack of an Fe-bound nitrile by water (or hydroxide). Reported herein is a five-coordinate model compound ([Fe(III)(S(2)(Me2)N(3)(Et,Pr))](+)) containing Fe(III) in an environment resembling that of NHase, which reversibly binds a variety of nitriles, alcohols, amines, and thiocyanate. XAS shows that five-coordinate [Fe(III)(S(2)(Me2)N(3)(Et,Pr))](+) reacts with both methanol and acetonitrile to afford a six-coordinate solvent-bound complex. Competitive binding studies demonstrate that MeCN preferentially binds over ROH, suggesting that nitriles would be capable of displacing the H(2)O coordinated to the iron site of NHase. Thermodynamic parameters were determined for acetonitrile (DeltaH = -6.2(+/-0.2) kcal/mol, DeltaS = -29.4(+/-0.8) eu), benzonitrile (-4.2(+/-0.6) kcal/mol, DeltaS = -18(+/-3) eu), and pyridine (DeltaH = -8(+/-1) kcal/mol, DeltaS = -41(+/-6) eu) binding to [Fe(III)(S(2)(Me2)N(3)(Et,Pr))](+) using variable-temperature electronic absorption spectroscopy. Ligand exchange kinetics were examined for acetonitrile, iso-propylnitrile, benzonitrile, and 4-tert-butylpyridine using (13)C NMR line-broadening analysis, at a variety of temperatures. Activation parameters for ligand exchange were determined to be DeltaH(+ +) = 7.1(+/-0.8) kcal/mol, DeltaS(+ +) = -10(+/-1) eu (acetonitrile), DeltaH(+ +) = 5.4(+/-0.6) kcal/mol, DeltaS(+ +) = -17(+/-2) eu (iso-propionitrile), DeltaH(+ +) = 4.9(+/-0.8) kcal/mol, DeltaS(+ +) = -20(+/-3) eu (benzonitrile), and DeltaH(+ +) = 4.7(+/-1.4) kcal/mol DeltaS(+ +) = -18(+/-2) eu (4-tert-butylpyridine). The thermodynamic parameters for pyridine binding to a related complex, [Fe(III)(S(2)(Me2)N(3)(Pr,Pr))](+) (DeltaH = -5.9(+/-0.8) kcal/mol, DeltaS = -24(+/-3) eu), are also reported, as well as kinetic parameters for 4-tert-butylpyridine exchange (DeltaH(+ +) = 3.1(+/-0.8) kcal/mol, DeltaS(+ +) = -25(+/-3) eu). These data show for the first time that, when it is contained in a ligand environment similar to that of NHase, Fe(III) is capable of forming a stable complex with nitriles. Also, the rates of ligand exchange demonstrate that low-spin Fe(III) in this ligand environment is more labile than expected. Furthermore, comparison of [Fe(III)(S(2)(Me2)N(3)(Et,Pr))](+) and [Fe(III)(S(2)(Me2)N(3)(Pr,Pr))](+) demonstrates how minor distortions induced by ligand constraints can dramatically alter the reactivity of a metal complex.     

Gas-phase acidity studies of multiple sites of adenine and adenine derivatives

The acidities of multiple sites in the purine nucleobase adenine (1) and adenine alkyl derivatives 9-ethyladenine (2), 3-methyladenine (3), 1-methyladenine (4), and N,N-dimethyladenine (5) have been investigated for the first time, using computational and experimental methods to provide an understanding of adenine reactivity. We have previously measured two acidic sites on adenine, with the N9 site being 19 kcal mol(-)(1) more acidic than the N10 site (333 +/- 2 versus 352 +/- 4 kcal mol(-)(1), respectively). In this work, we have established that 9-ethyladenine has two sites more acidic than water: the N10 (352 +/- 4 kcal mol(-)(1)) and the C8 (374 +/- 2 kcal mol(-)(1)). We have likewise measured two acidities for 3-methyladenine, the N10 (347 +/- 4 kcal mol(-)(1)) and the C2 (370 +/- 3 kcal mol(-)(1)). For 1-methyladenine and N,N-dimethyladenine, we measure the N9H acidity to be 331 +/- 2 and 333 +/- 2 kcal mol(-)(1), respectively. We believe that the bracketing of only one site for the latter species is a kinetic effect, which we discuss further in the paper. Computationally, we have found the interesting result that some of the vinylic C-H sites in these purine bases are predicted to be much more acidic than water (DeltaH(acid) = 390.7 kcal mol(-)(1)) in the gas phase, on the order of 373 kcal mol(-)(1). The acidic vinylic C-H sites are always adjacent to an N-R group, and this pattern is maintained regardless of whether the site is on the five- or six-membered ring of the purine. Vinylic C-H sites elsewhere on the purine have calculated acidities of about 400 kcal mol(-)(1). The differing acidities are interpreted through electrostatic potential calculations. We also relate our results to the intriguing biochemical decarboxylation of orotate ribose monophosphate, which involves a vinylic anion adjacent to an N-R group; this decarboxylation is the last step in the de novo biosynthesis of pyrimidine nucleotides, and the enzyme that catalyzes the reaction, orotate ribose monophosphate decarboxylase, has been the subject of intense study recently, as its mechanism remains elusive.     

A DFT study of SiH(4) activation by Cp(2)LnH

A theoretical study of SiH(4) activation by Cp(2)LnH complexes for the entire series of lanthanides has been carried out at the DFT-B3PW91 level of theory. The reaction paths corresponding to H/H exchange and silylation, formation of Cp(2)Ln(SiH(3)), have been computed. They both occur via a single-step sigma-bond metathesis mechanism. For the athermal H/H exchange reaction, the calculated activation barrier averages 1.8 kcal.mol(-)(1) relative to the precursor adduct Cp(2)LnH(eta(2)-SiH(4)) for all lanthanide elements. The silylation path is slightly exogenic (DeltaE approximately -6.5 kcal.mol(-1)) with an activation barrier averaging 5.2 kcal.mol(-1) relative to the precursor adduct where SiH(4) is bonded by two Si-H bonds. Both pathways are therefore thermally accessible. The H/H exchange path is calculated to be kinetically more favorable whereas the silylation reaction is thermodynamically preferred. The reactivity of this familly of lanthanide complexes with SiH(4) contrasts strongly with that obtained previously with CH(4). The considerably lower activation barrier for silylation relative to methylation is attributed to the ability of Si to become hypervalent.     

Theoretical study of the oxygen exchange in uranyl hydroxide. An old riddle solved?

A multistep mechanism for the experimentally observed oxygen exchange [Inorg. Chem. 1999, 38, 1456] of UO2(2+) cations in highly alkaline solutions is suggested and probed computationally. It involves an equilibrium between [UO2(OH)4](2-) and [UO2(OH)5](3-), followed by formation of the stable [UO3(OH)3 x H2O](3-) intermediate that forms from [UO2(OH)5](3-) through intramolecular water elimination. The [UO3(OH)3 x H2O](3-) intermediate facilitates oxygen exchange through proton shuttling, retaining trans-uranyl structures throughout, without formation of the cis-uranyl intermediates proposed earliar. Alternative cis-uranyl pathways have been explored but were found to have activation energies that are too high. Relativistic density functional theory (DFT) has been applied to obtain geometries and vibrational frequencies of the different species (reactants, intermediates, transition states, products) and to calculate reaction paths. Two different relativistic methods were used: a scalar four-component all-electron relativistic method and the zeroeth-order regular approximation. Calculations were conducted for both gas phase and condensed phase, the latter treated using the COSMO continuum model. An activation energy of 12.5 kcal/mol is found in solution for the rate-determining step, the reaction of changing the four-coordinated uranyl hydroxide to the five-coordinated one. This compares favorably to the experimental value of 9.8 +/- 0.7 kcal/mol. Activation energies of 7.8 and 5.1 kcal/mol are found for the hydrogen transfer between equatorial and axial oxygens through a water molecule in [UO3(OH)3 x H2O](3-) in the gas phase and condensed phase, respectively. Contrary to previously proposed mechanisms that resulted in high activation barriers, we find energies that are low enough to facilitate the reaction at room temperature. For the activation energies, two approximate DFT methods, B3LYP and PBE, are compared. The differences in activation energies are only about 1-2 kcal/mol for these methods.     

Effective monkey saddle points and berry and lever mechanisms in the topomerization of SF(4) and related tetracoordinated AX(4) species

The topomerization mechanisms of the SF(4) and SCl(2)F(2) sulfuranes, as well as their higher (SeF(4), TeF(4)) and isoelectronic analogues PF(4)(-), AsF(4)(-), SbF(4)(-), SbCl(4)(-), ClF(4)(+), BrF(4)(+), BrCl(2)F(2)(+), and IF(4)(+)), have been computed at B3LYP/6-31+G and at B3LYP/6-311+G. All species have trigonal bipyramidal (TBP) C(2)(v)() ground states. In such four-coordinated molecules, Berry rotation exchanges both axial with two equatorial ligands simultaneously while the alternative "lever" mechanism exchanges only one axial ligand with one equatorial ligand. While the barrier for the lever exchange in SF(4) (18.8 kcal mol(-1)) is much higher than that for the Berry process (8.1 kcal mol(-1)), both mechanisms are needed for complete ligand exchange. The F(ax)F(ax) and F(eq)F(eq) isomers of SF(2)Cl(2) have nearly the same energy and readily interconvert by BPR with a barrier of 7.6 kcal mol(-1). The enantiomerization of the F(ax)F(eq) chiral isomer can occur by either the Berry process (transition state barrier 8.3 kcal mol(-1)) or the "lever" mechanism via either of two C(s)() transition states, based on the TBP geometry: Cl(ax) <--> Cl(eq) or F(ax) <--> F(eq) exchanges with barriers of 6.3 and 15.7 kcal mol(-1), respectively. Full scrambling of all ligand sites is possible only by inclusion of the lever mechanism. Planar, "tetrahedral", and triplet forms are much higher in energy. The TBP C(3)(v) structures of AX(4) either have two imaginary frequencies (NIMAG = 2) for the X = F, Cl species or are minima (NIMAG = 0) for the X = Br, I compounds. These "effective monkey saddle points" have degenerate modes with two small frequencies, imaginary or real. Although a strictly defined "monkey saddle" (with degenerate frequencies exactly zero) is not allowed, the flat C(3)(v) symmetry region serves as a "transition state" for trifurcation of the pathways. The BPR mechanism also is preferred over the alternative lever process in the topomerization of the selenurane SeF(4) (barriers 5.9 vs. 12.1 kcal mol(-1)), the tellurane TeF(4) (2.1 vs. 6.4), and the interhalogen cations ClF(4)(+) (2.5 vs 14.8), BrF(4)(+) (4.7 vs. 11.3), BrF(2)Cl(2)(+) (14.6 vs. 17.4), and IF(4)(+) (1.4 vs. 6.0), as well as for the series PF(4)(-) (7.0 vs. 9.0), AsF(4)(-) (9.3 vs. 17.2), and SbF(4)(-) (3.8 vs. 5.3 kcal mol(-1)), all computed at B3LYP/6-311+G with the inclusion of quasirelativistic pseudopotentials for Te, I, and Sb. The heavier halogens increasingly favor the lever process, where the barrier (2.6 kcal mol(-1)) pertaining to the effective monkey saddle point (C(3)(v) minimum for SbCl(4)(-)) is less than that for the Berry process (8.2 kcal mol(-1)).     

New model for a theoretical density functional theory investigation of the mechanism of the carbonic anhydrase: how does the internal bicarbonate rearrangement occur?

A theoretical density functional theory (DFT, B3LYP) investigation has been carried out on the catalytic cycle of the carbonic anhydrase. A model system including the Glu106 and Thr199 residues and the "deep" water molecule has been used. It has been found that the nucleophilic attack of the zinc-bound OH on the CO(2) molecule has a negligible barrier (only 1.2 kcal mol(-1)). This small value is due to a hydrogen-bond network involving Glu106, Thr199, and the deep water molecule. The two usually proposed mechanisms for the internal bicarbonate rearrangement have been carefully examined. In the presence of the two Glu106 and Thr199 residues, the direct proton transfer (Lipscomb mechanism) is a two-step process, which proceeds via a proton relay network characterized by two activation barriers of 4.4 and 9.0 kcal mol(-1). This pathway can effectively compete with a rotational mechanism (Lindskog mechanism), which has a barrier of 13.2 kcal mol(-1). The fast proton transfer found here is basically due to the effect of the Glu106 residue, which stabilizes an intermediate situation where the Glu106 fragment is protonated. In the absence of Glu106, the barrier for the proton transfer is much larger (32.3 kcal mol(-1)) and the Lindskog mechanism becomes favored.     

Predicting the energy of the water exchange reaction and free energy of solvation for the uranyl ion in aqueous solution

The structures and vibrational frequencies of UO2(H2O)4(2+) and UO2(H2O)5(2+) have been calculated using density functional theory and are in reasonable agreement with experiment. The energies of various reactions were calculated at the density functional theory (DFT) and MP2 levels; the latter provides the best results. Self-consistent reaction field calculations in the PCM and SCIPCM approximations predicted the free energy of the water exchange reaction, UO2(H2O)4(2+) + H2O <--> UO2(H2O)5(2+). The calculated free energies of reaction are very sensitive to the choice of radii (O and H) and isodensity values in the PCM and SCIPCM models, respectively. Results consistent with the experimental HEXS value of -1.19 +/- 0.42 kcal/mol (within 1-3 kcal/mol) are obtained with small cavities. The structures and vibrational frequencies of the clusters with second solvation shell waters: UO2(H2O)4(H2O)8(2+), UO2(H2O)4(H2O)10(2+), UO2(H2O)4(H2O)11(2+), UO2(H2O)5(H2O)7(2+), and UO2(H2O)5(H2O)10(2+), were calculated and are in better agreement with experiment as compared to reactions involving only UO2(H2O)4(2+) and UO2(H2O)5(2+). The MP2 reaction energies for water exchange gave gas-phase results that agreed with experiment in the range -5.5 to +3.3 kcal/mol. The results were improved by inclusion of a standard PCM model with differences of -1.2 to +2.7 kcal/mol. Rearrangement reactions based on an intramolecular isomerization leading to a redistribution of water in the two shells provide good values in comparison to experiment with values of Delta G(exchange) from -2.2 to -0.5 kcal/mol so the inclusion of a second hydration sphere accounts for most solvation effects. Calculation of the free energy of solvation of the uranyl cation yielded an upper bound to the solvation energy of -410 +/- 5 kcal/mol, consistent with the best experimental value of -421 +/- 15 kcal/mol.     

DFT study of the glutathione peroxidase-like activity of phenylselenol incorporating solvent-assisted proton exchange

Modeling of the glutathione peroxidase-like activity of phenylselenol has been accomplished using density-functional theory and solvent-assisted proton exchange (SAPE). SAPE is a modeling technique intended to mimic solvent participation in proton transfer associated with chemical reaction. Within this method, explicit water molecules incorporated into the gas-phase model allow relay of a proton through the water molecules from the site of protonation in the reactant to that in the product. The activation barriers obtained by SAPE for the three steps of the GPx-like mechanism of PhSeH fall within the limits expected for a catalytic system at physiological temperatures (DeltaG(1)++ = 19.1 kcal/mol; DeltaG(2)++= 6.6 kcal/mol; G(3)++ = 21.7 kcal/mol) and are significantly lower than studies which require direct proton transfer. The size of the SAPE network is also considered for the model of the reduction of the selenenic acid, step 2 of the GPx-like cycle. Use of a four-water network better accommodates the reaction pathway and reduces the activation barrier by 5 kcal/mol over the two-water model.     

Intermolecular interaction in water hexamer

The origin of the intermolecular interaction, especially the many-body interaction, in eight low-lying water hexamer structures (prism, cage, book-1, book-2, cyclic-chair, bag, cyclic-boat-1, and cyclic-boat-2) is unraveled using the localized molecular orbital energy decomposition analysis (LMO-EDA) method at the second-order M?ller-Plesset perturbation (MP2) level of theory with a large basis set. It is found that the relative stabilities of these hexamer structures are determined by delicate balances between different types of interaction. According to LMO-EDA, electrostatic and exchange interactions are strictly pairwise additive. Dispersion interaction in these water hexamer structures is almost pairwise additive, with many-body effects varying from -0.13 to +0.05 kcal/mol. Repulsion interaction is roughly pairwise additive, with many-body effects varying from -0.84 to -0.62 kcal/mol. Polarization interaction is not pairwise additive, with many-body effects varying from -13.10 to -8.85 kcal/mol.     

Searching for stable, five-coordinate aquated Al(III) species. Water exchange mechanism and effect of pH

Density functional theory calculations have been performed for the water exchange mechanism of aquated Al(III). The effect of pH was considered by studying the exchange processes for [Al(H2O)6]3+ and its conjugated base, [Al(H2O)5OH]2+. Both complexes were found to exchange water in a dissociative way with activation energies (EA) of 15.9 and 10.2 kcal/mol, respectively. The influence of solvent molecules on the gas-phase cluster model was considered by the addition of up to four water molecules to the model system. The stabilizing effect of the solvent on the transition state decreases EA to 8.6 (hexa-aqua complex) and 7.6 (monohydroxo complex) kcal/mol, whereas EA for all hydroxo species is consistently significantly lower than those for the related aqua systems, which indicates a much faster water exchange rate. For the hydroxo complex, all calculated five-coordinate intermediates, nH2O.[Al(H2O)4(OH)]2+ (n = 1, 2, 3, 4, 5), are more stable than the corresponding six-coordinate reactants. Our results therefore suggest the presence of a stable five-coordinate species of aquated Al(III), namely, the [Al(H2O)4(OH)]2+ complex.     

Mechanism of epoxide hydrolysis in microsolvated nucleotide bases adenine, guanine and cytosine: a DFT study

Six water molecules have been used for microsolvation to outline a hydrogen bonded network around complexes of ethylene epoxide with nucleotide bases adenine (EAw), guanine (EGw) and cytosine (ECw). These models have been developed with the MPWB1K-PCM/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) level of DFT method and calculated S(N)2 type ring opening of the epoxide due to amino group of the nucleotide bases, viz. the N6 position of adenine, N2 position of guanine and N4 position of cytosine. Activation energy (E(act)) for the ring opening was found to be 28.06, 28.64, and 28.37 kcal mol(-1) respectively for EAw, EGw and ECw. If water molecules were not used, the reactions occurred at considerably high value of E(act), viz. 53.51 kcal mol(-1) for EA, 55.76 kcal mol(-1) for EG and 56.93 kcal mol(-1) for EC. The ring opening led to accumulation of negative charge on the developing alkoxide moiety and the water molecules around the charge localized regions showed strong hydrogen bond interactions to provide stability to the intermediate systems EAw-1, EGw-1 and ECw-1. This led to an easy migration of a proton from an activated water molecule to the alkoxide moiety to generate a hydroxide. Almost simultaneously, a proton transfer chain reaction occurred through the hydrogen bonded network of water molecules and resulted in the rupture of one of the N-H bonds of the quaternized amino group. The highest value of E(act) for the proton transfer step of the reaction was 2.17 kcal mol(-1) for EAw, 2.93 kcal mol(-1) for EGw and 0.02 kcal mol(-1) for ECw. Further, the overall reaction was exothermic by 17.99, 22.49 and 13.18 kcal mol(-1) for EAw, EGw and ECw, respectively, suggesting that the reaction is irreversible. Based on geometric features of the epoxide-nucleotide base complexes and the energetics, the highest reactivity is assigned for adenine followed by cytosine and guanine. Epoxide-mediated damage of DNA is reported in the literature and the present results suggest that hydrated DNA bases become highly S(N)2 active on epoxide systems and the occurrence of such reactions can inflict permanent damage to the DNA.     

Theoretical study on the reaction of VS+(3∑-, 1Γ) with CO

Two possible reaction mechanisms of VS (3Ε-, 1Γ) with CO in the gas phase have been studied by using B3LYP/TZVP and CCSD(T)/6-311 G (3df, 3pd) methods: the O/S exchange reaction (VS CO → VO CS) and the S-transfer reaction (VS CO → V COS). The two reactions proceed via two-step and one-step mechanism, respectively. The barriers of the triplet and singlet PESs are 30.6 and 50.9 kcal/mol, respectively, for O/S exchange reaction and 7.3 and 50.2 kcal/mol, respectively, for the S-transfer reaction. The results indicate that the triplet ground state reaction is more favorable, and the S-transfer reaction is more favorable than the O/S exchange reaction, which is in good agreement with the experimental observation.     

Beyond the Roger Brown rearrangement: long-range atom topomerization in conjugated polyynes

Long-range carbon atom topomerization in a 1,3-diyne has been demonstrated for the first time. 1-Phenyl-4-p-tolyl-1,3-butadiyne, (13)C-enriched at C-1, was synthesized and subjected to flash vacuum pyrolysis. At 800 degrees C and 0.01 Torr, this resulted in nearly complete (13)C label equilibration between C-1 and C-2, as seen by NMR analysis. Pyrolysis at 900 degrees C further led to ca. 35% of the label migrating about equally to C-3 and C-4. These results demonstrate that both intrabond and interbond atom exchange processes are operative, with the former having a lower activation barrier. DFT and Moller-Plesset calculations support a mechanism that passes through Brown rearrangement (1,2-shift), closure to trialene (bicyclo[1.1.0]-1,3-butadiene), bond-shift isomerization to exchange C-2 and C-3, and ring opening. The resulting vinylidene can rearrange to a butadiyne with the isotopic label at C-3 or C-4. Consistent with earlier calculations, trialene is predicted to have alternating peripheral bonds, with a weak central sigma bond and significant diradical character. Trialene is predicted [(B3LYP/6-311+G(2d,p)] to lie 64.6 kcal/mol above butadiyne, with barriers of 2.2 and 4.4 kcal/mol, respectively, for ring opening or bond-shift isomerization. Other potential rearrangement mechanisms which pass through tetrahedrene (E(rel) = 167.2 kcal/mol) or 1,2,3-cyclobutatriene (E(rel) = 161.1 kcal/mol) lie at much higher energies.     

Synthesis,structure, and dynamic behavior of cyclopentadienyl-lithium, -sodium,and -potassium annelated with bicyclo[2.2.2]octene units: a systematic study on site exchange of alkali metals on a cyclopentadienyl ring in tetrahydrofuran

Novel cyclopentadienyl (Cp)-alkali metal complexes 1-M and 2-M (M = Li, Na, K), in which the Cp ring is annelated with two bicyclo[2.2.2]octene units and substituted with a phenyl group for 1 and a tert-butyl group for 2, were synthesized, and their structures and dynamic behaviors were investigated by means of X-ray crystallography, dynamic (13)C NMR, and DFT calculations. The X-ray crystallography results indicated that 1-Li, 1-Na, and 2-Na form monomeric contact ion pairs (CIP) with three THF molecules coordinated to the metal atom. Also, in THF-d(8), all of the 1-M and 2-M form monomeric CIP in the ground state. However, variable-temperature (13)C NMR measurements of 1-M and 2-M in THF-d(8) demonstrated dynamic behavior in which the metal ion exchanges positions between the upper and lower faces of the Cp ring. From a study of the concentration dependence of the dynamic behavior, the exchange was found to proceed principally as an intramolecular process at concentration ranges lower than 0.2 M. The experimentally observed deltaG values for the intramolecular exchange process for all the 1-M and 2-M (except for 2-Li, whose intramolecular process was too slow to observe) were found to be quite similar in THF-d(8) solution and to fall within the range of 12-14 kcal mol(-)(1). Within this range, a tendency was observed for the deltaG values to increase as the size of the metal decreased. Theoretical calculations (B3LYP/6-31G(d)) afforded considerably large values as the gas-phase dissociation energy for 1-M (162.7 kcal mol(-)(1) for M = Li; 131.6 kcal mol(-)(1) for M = Na; 110.9 kcal mol(-)(1) for M = K) and for 2-M (170.0 kcal mol(-)(1) for M = Li; 137.5 kcal mol(-)(1) for M = Na; 115.4 kcal mol(-)(1) for M = K). These values should be compensated for by a decrease in the solvation energies for the metal ions with increasing size, as exemplified by the calculated solvation energy for M(+)(Me(2)O)(4), which serves as a model for metal ions solvated with four molecules of THF (-122.9 kcal mol(-)(1) for M = Li; -94.7 kcal mol(-)(1) for M = Na; -67.7 kcal mol(-)(1) for M = K). This compensation results in a small difference in the overall energy for dissociation of 1-M or 2-M in ethereal solutions, thus supporting the similar deltaG values observed for the intramolecular metal exchange.     

Improved density functionals for water

The accuracy of existing density functional methods for describing the noncovalent interaction energies in small water clusters is investigated by testing 25 density functionals against a data set of 28 water dimers and 8 water trimers whose structures are taken from the literature and from simulations. The most accurate functionals are found to be PW6B95 with a mean unsigned error of 0.13 kcal/mol and MPWB1K and B98 with mean unsigned errors of 0.15 kcal/mol; the best functional with no Hartree-Fock exchange is mPWLYP, which is a GGA with a mean unsigned error of 0.28 kcal/mol. In comparison, the most popular GGA functionals, PBE and BLYP, have mean unsigned errors of 0.52 and 1.03 kcal/mol, respectively. Since GGAs are very cost efficient for both condensed-phase simulations and electronic structure calculations on large systems, we optimized four new GGAs for water. The best of these, PBE1W and MPWLYP1W, have mean unsigned errors of 0.12 and 0.17 kcal/mol, respectively. These new functionals are well suited for use in condensed-phase simulations of water and ice.     

On the hydrolysis mechanism of the second-generation anticancer drug carboplatin

The hydrolysis reaction mechanisms of carboplatin, a second-generation anticancer drug, have been explored by combining density functional theory (DFT) with the conductor-like dielectric continuum model (CPCM) approach. The decomposition of carboplatin in water is expected to take place through a biphasic mechanism with a ring-opening process followed by the loss of the malonato ligand. We have investigated this reaction in water and acid conditions and established that the number of protons present in the malonato ligand has a direct effect on the energetics of this system. Close observation of the optimised structures revealed a necessary systematic water molecule in the vicinity of the amino groups of carboplatin. For this reason we have also investigated this reaction with an explicit water molecule. From the computed potential-energy surfaces it is established that the water hydrolysis takes place with an activation barrier of 30 kcal mol(-1), confirming the very slow reaction observed experimentally. The decomposition of carboplatin upon acidification was also investigated and we have computed a 21 kcal mol(-1) barrier to be overcome (experimental value 23 kcal mol(-1)). We have also established that the rate-limiting process is the first hydration, and ascertained the importance of a water molecule close to the two amine groups in lowering the activation barriers for the ring-opening reaction.     

Equilibrium thermodynamic studies in water: reactions of dihydrogen with rhodium(III) porphyrins relevant to Rh-Rh, Rh-H, and Rh-OH bond energetics

Aqueous solutions of rhodium(III) tetra p-sulfonatophenyl porphyrin ((TSPP)Rh(III)) complexes react with dihydrogen to produce equilibrium distributions between six rhodium species including rhodium hydride, rhodium(I), and rhodium(II) dimer complexes. Equilibrium thermodynamic studies (298 K) for this system establish the quantitative relationships that define the distribution of species in aqueous solution as a function of the dihydrogen and hydrogen ion concentrations through direct measurement of five equilibrium constants along with dissociation energies of D(2)O and dihydrogen in water. The hydride complex ([(TSPP)Rh-D(D(2)O)](-4)) is a weak acid (K(a)(298 K) = (8.0 +/- 0.5) x 10(-8)). Equilibrium constants and free energy changes for a series of reactions that could not be directly determined including homolysis reactions of the Rh(II)-Rh(II) dimer with water (D(2)O) and dihydrogen (D(2)) are derived from the directly measured equilibria. The rhodium hydride (Rh-D)(aq) and rhodium hydroxide (Rh-OD)(aq) bond dissociation free energies for [(TSPP)Rh-D(D(2)O)](-4) and [(TSPP)Rh-OD(D(2)O)](-4) in water are nearly equal (Rh-D = 60 +/- 3 kcal mol(-1), Rh-OD = 62 +/- 3 kcal mol(-1)). Free energy changes in aqueous media are reported for reactions that substitute hydroxide (OD(-)) (-11.9 +/- 0.1 kcal mol(-1)), hydride (D(-)) (-54.9 kcal mol(-1)), and (TSPP)Rh(I): (-7.3 +/- 0.1 kcal mol(-1)) for a water in [(TSPP)Rh(III)(D(2)O)(2)](-3) and for the rhodium hydride [(TSPP)Rh-D(D(2)O)](-4) to dissociate to produce a proton (9.7 +/- 0.1 kcal mol(-1)), a hydrogen atom (approximately 60 +/- 3 kcal mol(-1)), and a hydride (D(-)) (54.9 kcal mol(-1)) in water.     

NH tautomerism in the natural chlorin derivatives

The NH tautomerism of five Mg-free chlorophyll a and b derivatives 2-6 was studied utilizing NMR spectroscopy and molecular modeling. The results from the dynamic NMR measurements of the chlorins revealed that substituent effects contribute crucially to the free energy of activation (DeltaG(double dagger)) in the NH tautomeric processes. An intermediate tautomer for the total tautomeric NH exchange in a chlorin was observed for the first time, when the (1)H NMR spectra of chlorin e(6) TME (3) and rhodin g(7) TME (4) (TME = trimethyl ester) were measured at lower temperatures. The lower energy barriers (DeltaG(1)(double dagger)) obtained for the formation of the intermediate tautomers of 3 and 4, assigned to the N(22)-H, N(24)-H trans-tautomer, were 10.8 and 10.6 kcal/mol, respectively. The energy barrier (DeltaG(2)(double dagger) value) for the total tautomeric NH exchange in the five chlorins was found to vary from 13.6 kcal/mol to values higher than 18 kcal/mol. The lowest DeltaG(2)(double dagger) value (13.6 kcal/mol) was obtained for rhodochlorin XV dimethyl ester (2), which was the only chlorophyll derivative lacking the C(15) substituent. In the case of chlorins 4 and 5, the steric crowding around the methoxycarbonylmethyl group at C(15) raised the DeltaG(2)(double dagger) activation free-energy to 17.1 kcal/mol. However, the highest energy barrier with DeltaG(2)(double dagger) > 18 kcal/mol was observed for the NH exchange of pyropheophorbide a methyl ester (6), possessing the macrocycle rigidifying isocyclic ring E. Our results demonstrate that the steric strain, arising either from the steric crowding around the bulky substituent at C(15) or the macrocycle rigidifying isocyclic ring E, slows down the NH tautomeric process. We suggest that deformations in the chlorin skeleton are closely connected to the NH tautomeric exchange and that the exchange occurs by a stepwise proton-transfer mechanism via a hydrogen bridge.     

o-Quinone methide as alkylating agent of nitrogen, oxygen, and sulfur nucleophiles. The role of H-bonding and solvent effects on the reactivity through a DFT computational study

The reactivity of the alkylating agent o-quinone methide (o-QM) toward NH(3), H(2)O, and H(2)S, prototypes of nitrogen-, oxygen-, and sulfur-centered nucleophiles, has been studied by quantum chemical methods in the frame of DF theory (B3LYP) in reactions modeling its reactivity in water with biological nucleophiles. The computational analysis explores the reaction of NH(3), H(2)O, and H(2)S with o-QM, both free and H-bonded to a discrete water molecule, with the aim to rationalize the specific and general effect of the solvent on o-QM reactivity. Optimizations of stationary points were done at the B3LYP level using several basis sets [6-31G(d), 6-311+G(d,p), adding d and f functions to the S atom, 6-311+G(d,p),S(2df), and AUG-cc-pVTZ]. The activation energies calculated for the addition reactions were found to be reduced by the assistance of a water molecule, which makes easier the proton-transfer process in these alkylation reactions by at least 12.9, 10.5, and 6.0 kcal mol(-1) [at the B3LYP/AUG-cc-pVTZ//B3LYP/6-311+G(d,p) level], for ammonia, water, and hydrogen sulfide, respectively. A proper comparison of an uncatalyzed with a water-catalyzed reaction mechanism has been made on the basis of activation Gibbs free energies. In gas-phase alkylation of ammonia and water by o-QM, reactions assisted by an additional water molecule H-bonded to o-QM (water-catalyzed mechanism) are favored over their uncatalyzed counterparts by 5.6 and 4.0 kcal mol(-1) [at the B3LYP/6-311+G(d,p) level], respectively. In contrast, the hydrogen sulfide alkylation reaction in the gas phase shows a slight preference for a direct alkylation without water assistance, even though the free energy difference (DeltaDeltaG(#)) between the two reaction mechanisms is very small (by 1.0 kcal mol(-1) at the B3LYP/6-311+G(d,p),S(2df) level of theory). The bulk solvent effect, evaluated by the C-PCM model, significantly modifies the relative importance of the uncatalyzed and water-assisted alkylation mechanism by o-QM in comparison to the case in the gas phase. Unexpectedly, the uncatalyzed mechanism becomes highly favored over the catalyzed one in the alkylation reaction of ammonia (by 7.0 kcal mol(-1)) and hydrogen sulfide (by 4.0 kcal mol(-1)). In contrast, activation induced by water complexation still plays an important role in the o-QM hydration reaction in water as solvent.