Role of solvation in the reduction of the uranyl dication by water: a density functional study

We have studied the solvation of uranyl, UO(2)(2+), and the reduced species UO(OH)(2+) and U(OH)(2)(2+) systematically using three levels of approximation: direct application of a continuum model (M1); explicit quantum-chemical treatment of the first hydration sphere (M2); a combined quantum-chemical/continuum model approach (M3). We have optimized complexes with varying numbers of aquo ligands (n = 4-6) and compared their free energies of solvation. Models M1 and M2 have been found to recover the solvation energy only partially, underestimating it by approximately 100 kcal/mol or more. With our best model M3, the calculated hydration free energy Delta(h)G degrees of UO(2)(2+) is about -420 kcal/mol, which shifts to about -370 kcal/mol when corrected for the expected error of the model. This value agrees well with the experimentally determined interval, -437 kcal/mol < Delta(h)G degrees < -318 kcal/mol. Complexes with 5 and 6 aquo ligands have been found to be about equally favored with models M2 and M3. The same solvation models have been applied to a two-step reduction of UO(2)(2+) by water, previously theoretically studied in the gas phase. Our results show that the solvation contribution to the reaction free energy, about 60 kcal/mol, dominates the endoergicity of the reduction.     

Study of Hg22+ and complexes of NpO2+ and UO22+ in solution. examples of cation-cation interactions

Density functional theory (BPW91/TZ2P) is used to explore the nature of cation-cation interactions (CCIs) that exist between two actinyl cations in solution. Solvation, which is modeled using COSMO, favors the complexes (ONpO-ONpO)2+ and (ONpO-OUO)3+ over separated NpO2+(aq) and UO2(2+)(aq) cations because of the quadratic dependence of solvation on charge. For (OUO-OUO)4+, solvation effects, even though very large, are unable to overcome intrinsic electrostatic repulsion between the units. The actinyl-actinyl complexes are T-shaped, with the oxygen of one unit coordinated to the actinide metal of the other unit. The association free energies of (ONpO-ONpO)2+ and (ONpO-OUO)3+ are calculated as -42.1 and -29.2 kcal/mol. Explicit consideration of the first solvation shell at the B3LYP/LANL2DZ level suggests that the free energies of binding may be overestimated. The Hg2(2+) dication, though not considered a "traditional" CCI, is very similar to the actinyl-actinyl interaction. The binding free energy of Hg2(2+) in solution is calculated as -16.0 kcal/mol.     

Using ligand bite angles to control the hydricity of palladium diphosphine complexes

A series of [Pd(diphosphine)(2)](BF(4))(2) and Pd(diphosphine)(2) complexes have been prepared for which the natural bite angle of the diphosphine ligand varies from 78 degrees to 111 degrees. Structural studies have been completed for 7 of the 10 new complexes described. These structural studies indicate that the dihedral angle between the two planes formed by the two phosphorus atoms of the diphosphine ligands and palladium increases by over 50 degrees as the natural bite angle increases for the [Pd(diphosphine)(2)](BF(4))(2) complexes. The dihedral angle for the Pd(diphosphine)(2) complexes varies less than 10 degrees for the same range of natural bite angles. Equilibrium reactions of the Pd(diphosphine)(2) complexes with protonated bases to form the corresponding [HPd(diphosphine)(2)](+) complexes were used to determine the pK(a) values of the corresponding hydrides. Cyclic voltammetry studies of the [Pd(diphosphine)(2)](BF(4))(2) complexes were used to determine the half-wave potentials of the Pd(II/I) and Pd(I/0) couples. Thermochemical cycles, half-wave potentials, and measured pK(a) values were used to determine both the homolytic ([HPd(diphosphine)(2)](+) --> [Pd(diphosphine)(2)](+) + H*) and the heterolytic ([HPd(diphosphine)(2)](+) --> [Pd(diphosphine)(2)](2+) + H(-)) bond-dissociation free energies, Delta G(H*)* and Delta G(H-)*, respectively. Linear free-energy relationships are observed between pK(a) and the Pd(I/0) couple and between Delta G(H-)* and the Pd(II/I) couple. The measured values for Delta G(H*)* were all 57 kcal/mol, whereas the values of Delta G(H-)* ranged from 43 kcal/mol for [HPd(depe)(2)](+) (where depe is bis(diethylphosphino)ethane) to 70 kcal/mol for [HPd(EtXantphos)(2)](+) (where EtXantphos is 9,9-dimethyl-4,5-bis(diethylphosphino)xanthene). It is estimated that the natural bite angle of the ligand contributes approximately 20 kcal/mol to the observed difference of 27 kcal/mol for Delta G(H-)*.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Dependence of pKa on solute cavity for diprotic and triprotic acids

A systematic study of ΔG(aq)/pK(a) for monoprotic, diprotic, and triprotic acids has been carried out based on DFT/aug-cc-pVTZ combined with CPCM and SMD solvation modeling. All DFT/cavity set combinations considered showed similar accuracy for ΔG(aq)(1)/pK(a1) (70% within ±2.5 kcal mol(-1) of experiment) while only the M05-2X/Pauling cavity combination gave reasonable results for ΔG(aq)(2)/pK(a2) when both pK(a) values are separated by more than three units (70% within ±5.0 kcal mol(-1) of experiment). The choice of experimental data is critical to the interpretation of the calculated accuracy especially for several inorganic acids. For the calculation of ΔG(aq)(3)/pK(a3), the larger experimental uncertainty and an unrealistic orbital population of diffuse function for trianions in the gas phase hinders an evaluation of the predictive performance. We find the M05-2X functional with the Pauling cavity set is the best choice for ΔG(aq)(2)/pK(a2) prediction in aqueous media while all DFT/cavity sets considered were competitive for ΔG(aq)(1)/pK(a1).     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Thermodynamics of

Formation of cobalt(II)-thiocyanato complexes in nonionic surfactant solutions of poly(ethylene oxide) type with varying poly(ethylene oxide) chain lengths of 7.5 (Triton X-114), 30 (Triton X-305), and 40 (Triton X-405) has been studied by titration spectrophotometry and calorimetry at 298 K. Data were analyzed by assuming formation of a series of ternary complexes Co(NCS)(n)Y(m)((2-n)+) (Y=surfactant) with an overall formation constant beta(nm). In all the surfactant systems examined, data obtained can be explained well in terms of formation of Co(NCS)(+) and Co(NCS)(2) in an aqueous phase (aq), and Co(NCS)(4)Y(2-) in micelles, and their formation constants, enthalpies, and entropies have been determined. The beta(41)/beta(20) ratio increases and the corresponding enthalpy becomes significantly less negative with an increasing number of ethylene oxide groups. This suggests that micelles of these nonionic surfactants have a heterogeneous inner structure consisting of ethylene oxide and octylphenyl moieties. Indeed, on the basis of molar volumes of ethylene oxide and octylphenyl groups, intrinsic thermodynamic parameters have been extracted for the reaction Co(NCS)(2)(aq)+2NCS(-)(aq)=Co(NCS)(4)Y(2-) (Delta(r)G degrees, Delta(r)H degrees, and Delta(r)S degrees ) at each moiety. The Delta(r)G degrees, Delta(r)H degrees, and Delta(r)S degrees values are -16 kJ mol(-1), -15 kJ mol(-1), and 3 J K(-1) mol(-1), respectively, for the ethylene oxide moiety, and -15 kJ mol(-1), -70 kJ mol(-1), and -183 J K(-1) mol(-1) for octylphenyl. Significantly less negative Delta(r)H degrees and Delta(r)S degrees values for ethylene oxide imply that the hydrogen-bonded network structure of water is extensively formed at the ethylene oxide moiety, and the structure is thus broken around the Co(NCS)(4)(2-) complex with weak hydrogen-bonding ability. Copyright 2001 Academic Press.     

Thermochemistry of Lewis adducts of BH3 and nucleophilic substitution of triethylamine on NH3BH3 in tetrahydrofuran

The thermochemistry of the formation of Lewis base adducts of BH(3) in tetrahydrofuran (THF) solution and the gas phase and the kinetics of substitution on ammonia borane by triethylamine are reported. The dative bond energy of Lewis adducts were predicted using density functional theory at the B3LYP/DZVP2 and B3LYP/6-311+G** levels and correlated ab initio molecular orbital theories, including MP2, G3(MP2), and G3(MP2)B3LYP, and compared with available experimental data and accurate CCSD(T)/CBS theory results. The analysis showed that the G3 methods using either the MP2 or the B3LYP geometries reproduce the benchmark results usually to within ~1 kcal/mol. Energies calculated at the MP2/aug-cc-pVTZ level for geometries optimized at the B3LYP/DZVP2 or B3LYP/6-311+G** levels give dative bond energies 2-4 kcal/mol larger than benchmark values. The enthalpies for forming adducts in THF were determined by calorimetry and compared with the calculated energies for the gas phase reaction: THFBH(3) + L → LBH(3) + THF. The formation of NH(3)BH(3) in THF was observed to yield significantly more heat than gas phase dative bond energies predict, consistent with strong solvation of NH(3)BH(3). Substitution of NEt(3) on NH(3)BH(3) is an equilibrium process in THF solution (K ≈ 0.2 at 25 °C). The reaction obeys a reversible bimolecular kinetic rate law with the Arrhenius parameters: log A = 14.7 ± 1.1 and E(a) = 28.1 ± 1.5 kcal/mol. Simulation of the mechanism using the SM8 continuum solvation model shows the reaction most likely proceeds primarily by a classical S(N)2 mechanism.     

Molecule-induced alkane homolysis with dioxiranes.

The mechanisms of C-H and C-C bond activations with dimethyldioxirane (DMD) were studied experimentally and computationally at the B3LYP/6-311+G**//B3LYP/6-31G* density functional theory level for the propellanes 3,6-dehydrohomoadamantane (2) and 1,3-dehydroadamantane (3). The sigma(C-C) activation of 3 with DMD (Delta G(*) = 23.9 kcal mol(-1) and Delta G(r) = -5.4 kcal mol(-1)) is the first example of a molecule-induced homolytic C-C bond cleavage. The C-H bond hydroxylation observed for 2 is highly exergonic (Delta G(r) = -74.4 kcal mol(-1)) and follows a concerted pathway (Delta G(*) = 34.8 kcal mol(-1)), in contrast to its endergonic molecule-induced homolysis (Delta G(*) = 28.8 kcal mol(-1) and Delta G(r) = +9.2 kcal mol(-1)). The reactivities of 2 and 3 with CrO(2)Cl(2), which follow a molecule-induced homolytic activation mechanism, parallel the DMD results only for highly reactive 3, but differ considerably for more stable propellanes such as 4-phenyl-3,6-dehydrohomoadamantane (1) and 2.     

Quantum chemical simulation of the reaction of methane with gold(I) acetylacetonate and aquaacetylacetonate complexes

The interaction between methane and gold(I) acetylacetonate via electrophilic substitution (reaction (I)) and oxidative addition (reaction (II)) is simulated. In both cases, the formation of the products is thermodynamically favorable: the decrease in energy is 31 kcal/mol for reaction (I) and 26 kcal/mol for reaction (II). The product of reaction (II) is additionally stabilized by Au-H interaction. Both reactions have a low activation barrier and proceed via the formation of structurally different methane complexes reducing the energy of the system by 9.3 kcal/mol for reaction (I) and by 10.9 kcal/mol for reaction (II). The complex [Au(H₂O)(acac)] is also capable of forming methane complexes. These complexes result from a thermally neutral reaction and turn into products after overcoming a low energy barrier. The structure of the complex activating methane in the gold-rutin system is deduced from the data obtained.     

Mechanism of olefin metathesis with catalysis by ruthenium carbene complexes: density functional studies on model systems

Gradient-corrected (BP86) density functional calculations were used to study alternative mechanisms of the metathesis reactions between ethene and model catalysts [(PH(3))(L)Cl(2)Ru[double bond]CH(2)] with L=PH3 (I) and L=C(3)N(2)H(4)=imidazol-2-ylidene (II). On the associative pathway, the initial addition of ethene is calculated to be rate-determining for both catalysts (Delta G(22-25)*[double bond] kcal mol(-1)). The dissociative pathway starts with the dissociation of phosphane, which is rather facile (Delta G(298)* is approximately equal to 5-10 kcal mol(-1)). The resulting active species (L)Cl(2)Ru[double bond]CH(2) can coordinate ethene cis or trans to L. The cis addition is unfavorable and mechanistically irrelevant (Delta G(298)* is approximately equal to 21-25 kcal mol(-1)). The trans coordination is barrierless, and the rate-determining step in the subsequent catalytic cycle is either ring closure of the complex to yield the ruthenacyclobutane (catalyst I, Delta G(298)*=12 kcal mol(-1)), or the reverse reaction (catalyst II, ring opening, Delta G(298)*=10 kcal mol(-1)), that is, II is slightly more active than I. For both catalysts, the dissociative mechanism with trans olefin coordination is favored. The relative energies of the species on this pathway can be tuned by ligand variation, as seen in (PMe(3))(2)Cl(2)Ru[double bond]CH(2) (III), in which phosphane dissociation is impeded and olefin insertion is facilitated relative to I. The differences in calculated relative energies for the model catalysts I-III can be rationalized in terms of electronic effects. Comparisons with experiment indicate that steric effects must also be considered for real catalysts containing bulky substituents.     

Interaction of adenine adducts with thymine: a computational study

The existence of DNA adducts bring the danger of carcinogenesis because of mispairing with normal DNA bases. 1,N6-ethenoadenine adducts (epsilonA) and 1,N6-ethanoadenine adducts (EA) have been considered as DNA adducts to study the interaction with thymine, as DNA base. Several different stable conformers for each type of adenine adduct with thymine, [epsilonA(1)-T(I), epsilonA(2)-T(I), epsilonA(3)-T(I) and EA(1)-T(I), EA(2)-T(I), EA(3)-T(I)] and [epsilonA(1)-T(II), epsilonA(2)-T(II), epsilonA(3)-T(II) and EA(1)-T(II), EA(2)-T(II), EA(3)-T(II)], have been considered with regard to their interactions. The differences in their geometrical structures, energetic properties, and hydrogen-bonding strengths have also been compared with Watson-Crick adenine-thymine base pair (A-T). Single-point energy calculations at MP2/6-311++G** levels on B3LYP/6-31+G* optimized geometries have also been carried out to better estimate the hydrogen-bonding strengths. The basis set superposition error corrected hydrogen-bonding strength sequence at MP2/6-311++G**//B3LYP/6-31+G* for the most stable complexes is found to be EA(2)-T(I) (15.30 kcal/mol) > EA(1)-T(II) (14.98 kcal/mol) > EA(3)-T(II) (14.68 kcal/mol) > epsilonA(2)-T(I) (14.54 kcal/mol) > epsilonA(3)-T(II) (14.22 kcal/mol) > epsilonA(3)-T(II) (13.64 kcal/mol) > A-T (13.62 kcal/mol). The calculated reaction enthalpy value for epsilonA(2)-T(I) is 10.05 kcal/mol, which is the highest among the etheno adduct-thymine complexes and about 1.55 kcal/mol more than those obtained for Watson-Crick A-T base pair and the reaction enthalpy value for EA(1)- T(II) is 10.22 kcal/mol, which is highest among the ethano addcut-thymine complexes and about 1.72 kcal/mol more than those obtained for Watson-Crick A-T base pair. The aim of this research is to provide fundamental understanding of adenine adduct and thymine interaction at the molecular level and to aid in future experimental studies toward finding the possible cause of DNA damage.     

Computational studies of Cu(II)/Met and Cu(I)/Met binding motifs relevant for the chemistry of Alzheimer's disease

A systematic study of the binding motifs of Cu(II) and Cu(I) to a methionine model peptide, namely, N-formylmethioninamide 1, has been carried out by quantum chemical computations. Geometries of the coordination modes obtained at the B3LYP/6-31G(d) level of theory are discussed in the context of copper coordination by the peptide backbone and the S atom of a methionine residue in peptides with special emphasis on Met35 of the amyloid-beta peptide (Abeta) of Alzheimer's disease. The relative binding free energies in the gas phase, DeltaG(g), are calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level of theory, and the solvation affects are included by means of the COSMO model to obtain the relative binding energies in solution, DeltaG(aq). A free energy of binding, DeltaG(aq) = -19.4 kJ mol(-1), relative to aqueous Cu(II) and the free peptide is found for the most stable Cu(II)/Met complex, 12. The most stable Cu(I)/Met complex, 23, is bound by -15.6 kJ mol(-1) relative to the separated species. The reduction potential relative to the standard hydrogen electrode is estimated to be E degrees (12/23) = 0.41 V. On the basis of these results, the participation of Met35 as a low affinity binding site of Cu(II) in Abeta, and its role in the redox chemistry underlying Alzheimer's disease is discussed.     

Binding affinities for models of biologically available potential Cu(II) ligands relevant to Alzheimer's disease: an ab initio study

A systematic study of the binding affinities of the model biological ligands X: = (CH3)2S, CH3S-, CH3NH2, 4-CH3-imidazole (MeImid), C6H5O-, and CH3CO2- to (NH3)i(H2O)3-iCu(II)-H2O (i = 3, 2, 1, 0) complexes has been carried out using quantum chemical calculations. Geometries have been obtained at the B3LYP/ 6-31G(d) level of theory, and binding energies, Delta, relative to H2O as a ligand, have been calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level. Solvation effects have been included using the COSMO model, and the relative binding free energies in aqueous solution (Delta) have been determined at pH 7 for processes that are pH dependent. CH3S- (Delta = -16.0 to -53.5 kJ mol(-1)) and MeImid (Delta = -18.5 to -35.2 kJ mol(-1)) give the largest binding affinities for Cu(II). PhO- and (CH3)2S are poor ligands for Cu(II), Delta = 20.6 to -9.7 and 19.8 to -3.7 kJ mol(-1), respectively. The binding affinities for CH3NH2 range from -0.8 to -15.0 kJ mol(-1). CH3CO2- has Cu(II) binding affinities in the ranges Delta = -13.5 to -32.4 kJ mol(-1) if an adjacent OH bond is available for hydrogen bonding and Delta = 10.1 to -4.6 kJ mol(-1) if this interaction is not present. In the context of copper coordination by the Abeta peptide of Alzheimer's disease, the binding affinities suggest preferential binding of Cu(II) to the three histidine residues plus a lysine or the N-terminus. For a 3N1O Cu(II) ligand arrangement, it is more probable that the oxygen ligand comes from an aspartate/glutamate residue side chain than from the tyrosine at position 10. Methionine appears unlikely to be a Cu(II) ligand in Abeta.     

Determination of the gas-phase acidities of cysteine-polyalanine peptides using the extended kinetic method

We determined the gas-phase acidities of two cysteine-polyalanine peptides, HSCA3 and HSCA4, using a triple-quadrupole mass spectrometer through application of the extended kinetic method with full entropy analysis. Five halogenated carboxylic acids were used as the reference acids. The negatively charged proton-bound dimers of the deprotonated peptides with the conjugate bases of the reference acids were generated by electrospray ionization. Collision-induced dissociation (CID) experiments were carried out at three collision energies. The enthalpies of deprotonation (Delta(acid)H) of the peptides were derived according to the linear relationship between the logarithms of the CID product ion branching ratios and the differences of the gas-phase acidities. The values were determined to be Delta(acid)H(HSCA3) = 317.3 +/- 2.4 kcal/mol and Delta(acid)H (HSCA4) = 316.2 +/- 3.9 kcal/mol. Large entropy effects (Delta(DeltaS) = 13-16 cal/mol K) were observed for these systems. Combining the enthalpies of deprotonation with the entropy term yielded the apparent gas-phase acidities (Delta(acid)G(app)) of 322.1 +/- 2.4 kcal/mol (HSCA3) and 320.1 +/- 3.9 kcal/mol (HSCA4), in agreement with the results obtained from the CID-bracketing experiments. Compared with that in the isolated cysteine residue, the thiol group in HSCA3,4 has a stronger gas-phase acidity by about 20 kcal/mol. This increased acidity is likely due to the stabilization of the negatively charged thiolate group through internal solvation.     

Exploring a reaction mechanism for acetato ligand replacement in paddlewheel tetrakisacetatodirhodium (II,II) complex by ammonia: computational density functional theory study

This study focuses on the first step of interaction between DNA and the paddle-wheel dirhodium complex. The ammonia molecule was used to model the oligonucleotide sequence. The reaction was considered in neutral and acidic conditions, in gas phase, and in solvent, using the COSMO model. Molecular structures of the complexes were optimized in both models at the B3PW91/6-31G(d) level. The B3LYP functional and aug-cc-pvdz basis set were employed for single-point energy determination and electron distribution analyses. It was shown that in neutral solution the replacement of axial aqua ligand is mildly exoergic. The reaction is characterized by a relatively low activation barrier (10-12 kcal/mol), and, according to Eyring transition state theory, it proceeds very quickly. The breaking of the Rh-O(ac) bond in neutral solution is mildly endoergic (less than 1 kcal/mol) with an activation barrier of about 21 kcal/mol. However, this process can occur much more spontaneously (ΔG of -14 kcal/mol) when the dirhodium complex is protonated at the acetyl oxygen in remote position.     

Coordination mode of nitrate in uranyl(VI) complexes: a first-principles molecular dynamics study

According to Car-Parrinello molecular dynamics simulations for [UO(2)(NO(3))(3)](-), [UO(2)(NO(3))(4)](2-), and [UO(2)(OH(2))(4-)(NO(3))](+) complexes in the gas phase and in aqueous solution, the nitrate coordination mode to uranyl depends on the interplay between ligand-metal attractions, interligand repulsions, and solvation. In the trinitrate, the eta(2)-coordination is clearly favored in water and in the gas phase, leading to a coordination number (CN) of 6. According to pointwise thermodynamic integration involving constrained molecular dynamics simulations, a change in free energy of +6 kcal/mol is predicted for eta(2)- to eta(1)-transition of one of the three nitrate ligands in the gas phase. In the gas phase, the mononitrate-hydrate complex also prefers a eta(2)-binding mode but with a CN of 5, one H(2)O molecule being in the second shell. This contrasts with the aqueous solution where the nitrate binds in a eta(1)-fashion and uranyl coordinates to four H2O ligands. A driving force of ca. -3 kcal/mol is predicted for the eta(2)- to eta(1)- transition in water. This structural preference is interpreted in terms of steric arguments and differential solvation of terminal vs uranyl-coordinated O atoms of the nitrate ligands. The [UO(2)(NO(3))(4)](2-) complex with two eta(2)- and two eta(1)- coordinated nitrates, observed in the solid state, is stable for 1-2 ps in the gas phase and in solution. In the studied series, the modulation of uranyl-ligand distances upon immersion of the complex in water is found to depend on the nature of the ligand and the composition of the complex.     

Supramolecular binding thermodynamics by dispersion-corrected density functional theory

The equilibrium association free enthalpies ΔG(a) for typical supramolecular complexes in solution are calculated by ab initio quantum chemical methods. Ten neutral and three positively charged complexes with experimental ΔG(a) values in the range 0 to -21?kcal?mol(-1) (on average -6?kcal?mol(-1) ) are investigated. The theoretical approach employs a (nondynamic) single-structure model, but computes the various energy terms accurately without any special empirical adjustments. Dispersion corrected density functional theory (DFT-D3) with extended basis sets (triple-ζ and quadruple-ζ quality) is used to determine structures and gas-phase interaction energies (ΔE), the COSMO-RS continuum solvation model (based on DFT data) provides solvation free enthalpies and the remaining ro-vibrational enthalpic/entropic contributions are obtained from harmonic frequency calculations. Low-lying vibrational modes are treated by a free-rotor approximation. The accurate account of London dispersion interactions is mandatory with contributions in the range -5 to -60?kcal?mol(-1) (up to 200?% of ΔE). Inclusion of three-body dispersion effects improves the results considerably. A semilocal (TPSS) and a hybrid density functional (PW6B95) have been tested. Although the ΔG(a) values result as a sum of individually large terms with opposite sign (ΔE vs. solvation and entropy change), the approach provides unprecedented accuracy for ΔG(a) values with errors of only 2?kcal?mol(-1) on average. Relative affinities for different guests inside the same host are always obtained correctly. The procedure is suggested as a predictive tool in supramolecular chemistry and can be applied routinely to semirigid systems with 300-400 atoms. The various contributions to binding and enthalpy-entropy compensations are discussed.     

Activation barriers and rate constants for hydration of platinum and palladium square-planar complexes: an ab initio study

In the present work, an ab initio study on hydration (a metal-ligand replacement by water molecule or OH- group) of cis- and transplatin and their palladium analogs was performed within a neutral pseudomolecule approach (e.g., metal-complex+water as reactant complex). Subsequent replacement of the second ligand was considered. Optimizations were performed at the MP2/6-31+G(d) level with single-point energy evaluation using the CCSD(T)/6-31++G(d,p) approach. For the obtained structures of reactants, transition states (TS's), and products, both thermodynamic (reaction energies and Gibbs energies) and kinetic (rate constants) characteristics were estimated. It was found that all the hydration processes are mildly endothermic reactions-in the first step they require 8.7 and 10.2 kcal/mol for ammonium and chloride replacement in cisplatin and 13.8 and 17.8 kcal/mol in the transplatin case, respectively. Corresponding energies for cispalladium amount to 5.2 and 9.8 kcal/mol, and 11.0 and 17.7 kcal/mol for transpalladium. Based on vibrational analyses at MP2/6-31+G(d) level, transition state theory rate constants were computed for all the hydration reactions. A qualitative agreement between the predicted and known experimental data was achieved. It was also found that the close similarities in reaction thermodynamics of both Pd(II) and Pt(II) complexes (average difference for all the hydration reactions are approximately 1.8 kcal/mol) do not correspond to the TS characteristics. The TS energies for examined Pd(II) complexes are about 9.7 kcal/mol lower in comparison with the Pt analogs. This leads to 10(6) times faster reaction course in the Pd cases. This is by 1 or 2 orders of magnitude more than the results based on experimental measurements.     

Solvation of Au+ versus Au0 in aqueous solution: electronic structure governs solvation shell patterns

The solvation behavior of Au(+) and Au(0) in liquid water under ambient conditions has been studied using ab initio molecular dynamics. The Au(+) aqua ion forms a rigid and well-defined quasi-linear structure in the sense of ligand field theory, where two water molecules are tightly bound to the gold cation through oxygen atoms ("cationic solvation"). Yet, transient charge accumulation in the direction perpendicular to the O-Au(+)-O linear core structure leads occasionally to the formation of a short Au(+)-H contact within the distance range of the first solvation shell, which is typical of "anionic solvation". Upon adding an electron to Au(+), the resulting solvation pattern of Au(0)(aq) has nothing in common with that of Au(+)(aq). Quite surprisingly we discover that the first solvation shell of Au(0)(aq) consists of a single water molecule and features both "anionic" and "cationic" solvation patterns depending on fluctuation and polarization effects. Thus, charging/decharging of metals dissolved in water, M(0)? M(+) + e(-), as occurring e.g. during elementary electrochemical steps, is expected to change dramatically their solvation behavior in the sense of re-solvation processes.