The Gas-Phase Acidities of Substituted Hydroxamic and Silahydroxamic Acids: A Comparative ab Initio Study

Conformational and structural studies of N-methylacetohydroxamic acid and of its mono- and bis-chelated uranium(VI) complexes

The thermodynamics and kinetics of the cis/trans isomerism of N-methylacetohydroxamic acid (NMAH) and its conjugated base (NMA(-)) have been reinvestigated in aqueous media by (1)H NMR spectroscopy. Hindered rotation around the central C-N bond due to electronic delocalization becomes slow enough on the NMR time scale to observe both rotamers in equilibrium in D2O at room temperature. By properly assigning the methyl group resonances, evidence for the prevalence of the E over the Z form is unambiguously provided [K300=[E]/[Z]=2.86(2) and 9.63(5) for NMAH and NMA(-), respectively], closing thereby a long-lasting dispute about the most stable conformer. To that end, calculations of the chemical shifts by density functional theory (DFT), which accurately reproduced the experimental data, turned out to be a much more reliable method than the direct computation of the relative energy for each conformer. The Z ⇌ E interconversion dynamics was probed at 300 K in D2O by 2D exchange-correlated spectroscopy (EXSY), affording the associated rate constants [kZE=9.0(2) s(-1) and kEZ=3.14(5) s(-1) for NMAH, kZE=0.96(3) s(-1) and kEZ=0.10(2) s(-1) for NMA(-)] and activation barriers at 300 K [ΔG(≠)ZE=68.0 kJ mol(-1) and ΔG(≠)EZ=70.6 kJ mol(-1) for NMAH, ΔG(≠)ZE=73.6 kJ mol(-1) and ΔG(≠)EZ=79.2 kJ mol(-1) for NMA(-)]. For the first time, mono- and bis-chelated uranium(VI) complexes of NMA(-) have been isolated. Crystals of [UO2(NMA)(NO3)(H2O)2] and [UO2(NMA)2(H2O)] have been characterized by X-ray diffractometry, infrared and Raman spectroscopies.

Investigation of simple and water assisted tautomerism in a derivative of 1,3,4-oxadiazole: a DFT study

Investigation of tautomerism and transition states in a derivative of 1,3,4-oxadiazole (A, B, C and D) in the gas phase and in solution and in a micro hydrated environment with 1-3 water molecules was performed by calculations at the DFT-B3LYP/6-311++G(d,p) level of theory. The solvent effect is taken into account via the self-consistent reaction field (SCRF) method. The geometries of four possible tautomers of 5-amino-1,3,4-oxadiazole-2(3H)-one were optimized in the gas phase and solution with polarized continuum model (PCM). It was found that in the gas phase and different solvents, A and C tautomers are the most stable and unstable forms, respectively. The results show that the tautomeric interconversion C to D has the lowest Gibbs free energy changes and so the highest equilibrium constant in the gas phase and solution. The equilibrium and rate constants of intermolecular tautomerism in the absence and presence of 1-3 molecules of water were also calculated. The calculated results show that the presence of water molecules considerably reduces the barrier energy of the various reactions. Therefore, this water-assisted tautomerism can be performed fast, especially, with the assistance of two molecules of water.

Toward the accurate calculation of pKa values in water and acetonitrile

We present a simple approach for the calculation of accurate pKa values in water and acetonitrile based on the straightforward calculation of the gas-phase absolute free energies of the acid and conjugate base with use of only a continuum solvation model to obtain the corresponding solution-phase free energies. Most of the error in such an approach arises from inaccurate differential solvation free energies of the acid and conjugate base which is removed in our approach using a correction based on the realization that the gas-phase acidities have only a small systematic error relative to the dominant systematic error in the differential solvation. The methodology is outlined in the context of the calculation of a set of neutral acids with water as the solvent for a reasonably accurate electronic structure level of theory (DFT), basis set, and implicit solvation model. It is then applied to the comparison of results for three different hybrid density functionals to illustrate the insensitivity to the functional. Finally, the approach is applied to the comparison of results for sets of neutral acids and protonated amine cationic acids in both aqueous (water) and nonaqueous (acetonitrile) solvents. The methodology is shown to generally predict the pKa values for all the cases investigated to within 1 pH unit so long as the differential solvation error is larger than the systematic error in the gas-phase acidity calculations. Such an approach is rather general and does not have additional complications that would arise in a cluster-continuum method, thus giving it strength as a simple high-throughput means to calculate absolute pKa values. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

An ONIOM investigation on anion recognition of alkali-metal complexes with diurea calix[4]arene receptor

The ONIOM(B3LYP/6-31G(d):AM1) optimized structures of complexes of diurea calix[4]arene receptor (L) with alkali metals Li(+), Na(+) and K(+) and their complexes with halide ions F(-), Cl(-), Br(-), oxygen-containing anions HCO(3)(-), HSO(4)(-) and CH(3)COO(-) ions were obtained. Binding energies and thermodynamic properties of complex receptors LiL(+), NaL(+) and KL(+) with these anions were determined. The binding stabilities according to binding energies of LiL(+), NaL(+) and KL(+) associated with anions computed either at the ZPVE-corrected ONIOM(B3LYP/6-31G(d):AM1) or BSSE-corrected B3LYP/6-31 + G(d,p)//ONIOM(B3LYP/6-31G(d):AM1) are in the same order: F(-) >> CH(3)COO(-) ≈ HCO(3)(-) > Br(-) ≈ HSO(4)(-) ≈ Cl(-). All the receptors LiL(+), NaL(+) and KL(+) were found to be selective toward fluoride ion.

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.

Molecular modeling of nitrosamines adsorbed on H-ZSM-5 zeolite: an ONIOM study

Binding energies of nitrosamine compounds, N-nitrosamine (NA), N-methyl-N-nitrososamine (NMA), N-ethyl-N-nitrososamine (NEA), N,N-dimethyl-N-nitrosoamine (NDMA), N-ethyl-N-methyl-N-nitrosoamine (NEMA) and N,N-diethyl-N-nitrosoamine (NDEA) on the H-ZSM-5 zeolite were obtained using the ONIOM(B3LYP/6-31G(d):AM1) approach. Based on amino and imino isomers of nitrosamines, there are two adsorption configurations on the H-ZSM-5 for NA (as NA_a and NA_i), NMA (as NMA_a and NMA_i) and NEA (as NEA_a and NEA_i). The relative binding energies of nitrosamines are in order: NA_a > NMA_a approximately NEA_a > NA_i > NMA_i approximately NEA_i > NEMA approximately NDEA > NDMA. The order of adsorption selectivity for nitrosamines of the H-ZSM-5 is NA_a approximately NA_i >> NMA_a approximately NEA_a > NDMA approximately NMA_i approximately NEMA > NDEA approximately NEA_i. The selective recognition of the NA by the H-ZSM-5 was obviously found.

Computation of pKaValues of Substituted Aniline Radical Cations in Dimethylsulfoxide Solution

A newly developed computation strategy was used to calculate the absolute pKa values of 18 substituted aniline radical cations in dimethylsulfoxide (DMSO) solution with the error origin elucidated and deviation minimized. The B3LYP/6-311++G(2df,2p) method was applied and was found to be capable of reproducing the gas-phase proton-transfer free energies of substituted anilines with a precision of 0.83 kcal/mol. The IEF-PCM solvation model with gas-phase optimized structures was adopted in calculating the pKa values of the substituted neutral anilines in DMSO, regenerating the experimental results within a standard deviation of 0.4 pKa unit. When the IEF-PCM solvation model was applied to calculate the standard redox potentials of anilide anions, it showed that the computed values agreed well with experiment, but the redox potentials of substituted anilines were systematically overestimated by 0.304 eV. The cause of this deviation was found to be related to the inaccuracy of the calculated solvation free energies of aniline radical cations. By adjusting the size of the cavity in the IEF-PCM method, we derived a reliable procedure that can reproduce the experimental pKa values of aniline radical cations within 1.2 pKa units to those from experiment.

Ab initio molecular dynamics study of the hydration of the formohydroxamate anion

We apply ab initio molecular dynamics (AIMD) to study the hydration structures and electronic properties of the formohydroxamate anion in liquid water. We consider the cis- nitrogen-deprotonated, cis- oxygen-deprotonated, and trans- oxygen-deprotonated formohydroxamate tautomers. They form an average of 6.3, 6.9, and 6.0 hydrogen bonds with water molecules, respectively. The predicted pair correlation functions and time dependence of the hydration numbers suggest that water is highly structured around the nominally negatively charged oxime oxygen in O-deprotonated tautomers but significantly less so around the nitrogen atom in the N-deprotonated species. Wannier function analysis suggests that, in the O-deprotonated anions, the negative charge is concentrated on the oxime oxygen, while in the N-deprotonated case, it is partially delocalized between the nitrogen and the adjoining oxime oxygen atom.

Molecular model for host–guest interaction of tetraamino-tert-butylthiacalix[4]arene and tetraamino-tert-butylcalix[4]arene receptors with carboxylate and dicarboxylate guests: an ONIOM study

Geometry optimizations of tetraamino-tert-butylthiacalix[4]arene (tatbtc4a) and tetraamino-tert-butylcalix[4]arene (tatbc4a) complexes with acetate, oxalate,malonate, succinate, glutarate, adipate, and pimelate were carried out using the integrated MO:MO method. Thermodynamic quantities, preorganization energies and complexation energies of these complexes were obtained at the ONIOM(B3LYP/6-31G(d):AM1) level of theory. The relative stabilities of the tatbtc4a and tatbc4a complexes with carboxylate guests are reported. The complexes tatbtc4a/malonate and tatbc4a/oxalate were found to be the most stable species. The selectivity of the tatbtc4a receptor toward to malonate with respect to oxalate, in terms of selectivity coefficient, K(oxalatemalonate) is 9.90×10(2).

First-Principle Predictions of Absolute pKa's of Organic Acids in Dimethyl Sulfoxide Solution

MP2/6-311++G(d,p) and B3LYP/6-311++G(2df,p) methods were found to be able to predict the gas-phase acidities of various organic acids with a precision of 2.2 and 2.3 kcal/mol. A PCM cluster-continuum solvation method was developed that could predict the solvation free energies of various neutral, cationic, and anionic organic species in DMSO with a precision of about 2.0 kcal/mol. Using these carefully tested methods, we successfully predicted the pKa's of 105 organic acids in DMSO with a precision of 1.7-1.8 pKa units. We also predicted the pKa's of a variety of organosilanes in DMSO for the first time using the newly developed methods. This study was one of the first that employed first-principle methods for calculating pKa's of unrelated compounds in organic solutions.

Deprotonation sites of acetohydroxamic acid isomers. A theoretical and experimental study

Theoretical (ab initio calculations) and experimental (NMR, spectrophotometric, and potentiometric measurements) investigations of the isomers of acetohydroxamic acid (AHA) and their deprotonation processes have been performed. Calculations with the Gaussian 98 package, refined at the MP2(FC)/AUG-cc-pVDZ level considering the molecule isolated, indicate that the Z(cis) amide is the most stable form of the neutral molecule. This species and the less stable (Z)-imide form undergo deprotonation, giving rise to two stable anions. Upon deprotonation, the E(trans) forms give three stable anions. The ab initio calculations were performed in solution as well, regarding water as a continuous dielectric; on the basis of the relative energies of the most stable anion and neutral forms, calculated with MP2/PCM/AUG-cc-pVDZ, N-deprotonation of the amide (Z or E) structure appeared to be the most likely process in solution. NMR measurements provided evidence for the existence of (Z)- and (E)-isomers of both the neutral and anion forms in solution. Comparisons of the dynamic NMR and NOESY (one-dimensional) results obtained for the neutral species and their anions were consistent with N-deprotonation, which occurred preferentially to O-deprotonation. The (microscopic) acid dissociation constants of the two isomers determined at 25 degrees C from the pH dependence of the relevant chemical shifts, pK(E) = 9.01 and pK(Z) = 9.35, were consistent with the spectrophotometric and potentiometric evaluations (pK(HA) = 9.31).