Equilibrium isotope effects as a probe of nonbonding attractions

The equilibrium isotope effects (EIEs) in the capsule complex (12 + p-xylene + CCl4) were computed with two simple models that mimic two types of contacts, -CH3.Ar and -CH3.CCl4. By using the MP2/6-311++G(d,p) method, the calculated equilibrium constants are consistent with the reported experimental data. Harmonic frequency analyses indicate that vibrational frequency changes caused by nonbonding attractions lead to the observed EIE.     

Equilibrium isotope effects on noncovalent interactions in a supramolecular host-guest system

The self-assembled supramolecular complex [Ga(4)L(6)](12-) (1; L = 1,5-bis[2,3-dihydroxybenzamido]naphthalene) can act as a molecular host in aqueous solution and bind cationic guest molecules to its highly charged exterior surface or within its hydrophobic interior cavity. The distinct internal cavity of host 1 modifies the physical properties and reactivity of bound guest molecules and can be used to catalyze a variety of chemical transformations. Noncovalent host-guest interactions in large part control guest binding, molecular recognition and the chemical reactivity of bound guests. Herein we examine equilibrium isotope effects (EIEs) on both exterior and interior guest binding to host 1 and use these effects to probe the details of noncovalent host-guest interactions. For both interior and exterior binding of a benzylphosphonium guest in aqueous solution, protiated guests are found to bind more strongly to host 1 (K(H)/K(D) > 1) and the preferred association of protiated guests is driven by enthalpy and opposed by entropy. Deuteration of guest methyl and benzyl C-H bonds results in a larger EIE than deuteration of guest aromatic C-H bonds. The observed EIEs can be well explained by considering changes in guest vibrational force constants and zero-point energies. DFT calculations further confirm the origins of these EIEs and suggest that changes in low-frequency guest C-H/D vibrational motions (bends, wags, etc.) are primarily responsible for the observed EIEs.     

Comparisons of phosphorothioate with phosphate transfer reactions for a monoester,diester, and triester: isotope effect studies

Phosphorothioate esters are sometimes used as surrogates for phosphate ester substrates in studies of enzymatic phosphoryl transfer reactions. To gain better understanding of the comparative inherent chemistry of the two types of esters, we have measured equilibrium and kinetic isotope effects for several phosphorothioate esters of p-nitrophenol (pNPPT) and compared the results with data from phosphate esters. The primary (18)O isotope effect at the phenolic group ((18)k(bridge)), the secondary nitrogen-15 isotope effect ((15)k) in the nitro group, and (for the monoester and diester) the secondary oxygen-18 isotope effect ((18)k(nonbridge)) in the phosphoryl oxygens were measured. The equilibrium isotope effect (EIE) (18)k(nonbridge) for the deprotonation of the monoanion of pNPPT is 1.015 +/- 0.002, very similar to values previously reported for phosphate monoesters. The EIEs for complexation of Zn(2+) and Cd(2+) with the dianion pNPPT(2-) were both unity. The mechanism of the aqueous hydrolysis of the monoanion and dianion of pNPPT, the diester ethyl pNPPT, and the triester dimethyl pNPPT was probed using heavy atom kinetic isotope effects. The results were compared with the data reported for analogous phosphate monoester, diester, and triester reactions. The results suggest that leaving group bond fission in the transition state of reactions of the monoester pNPPT is more advanced than for its phosphate counterpart pNPP, while alkaline hydrolysis of the phosphorothioate diester and triester exhibits somewhat less advanced bond fission than that of their phosphate ester counterparts.     

Experimental evidence for a temperature dependent transition between normal and inverse equilibrium isotope effects for oxidative addition of H2 to Ir(PMe2Ph)2(CO)Cl

The equilibrium isotope effect (EIE) for oxidative addition of H(2) and D(2) to Ir(PMe(2)Ph)(2)(CO)Cl has been measured over a large temperature range, thereby demonstrating that the inverse (<1) EIE previously observed at ambient temperature becomes normal (>1) at high temperature (>90 degrees C). The temperature dependence of the EIE for oxidative addition of H(2) and D(2) to Ir(PH(3))(2)(CO)Cl has been calculated using the geometry and vibrational frequencies obtained from DFT (B3LYP) calculations on Ir(PH(3))(2)(CO)ClH(2) and Ir(PH(3))(2)(CO)ClD(2), and is in accord with the experimentally observed transition from an inverse to normal EIE for oxidative addition of H(2) and D(2) to Ir(PMe(2)Ph)(2)(CO)Cl: the EIE is calculated to be inverse between 0 and 510 K, reach a maximum value of 1.15 at 867 K and then slowly decrease to unity as the temperature approaches infinity. This deviation from simple van't Hoff behavior, and the occurrence of a maximum in the EIE, is the result of the entropy term being temperature dependent. At low temperature, the enthalpy term dominates and the EIE is inverse, whereas at high temperatures the entropy term dominates and the EIE is normal. The observation of both normal and inverse EIEs for the same system indicate that inferences pertaining to the magnitude of an isotope effect at a single temperature may require more detailed consideration than previously realized.     

Temperature-dependent transitions between normal and inverse equilibrium isotope effects for coordination and oxidative addition of C-H and H-H bonds to a transition metal center

The temperature dependence of the equilibrium isotope effects (EIEs) for coordination and oxidative addition of C-H and H-H bonds to the tungstenocene species {[H2Si(C5H4)2]W} has been determined with the aid of DFT (B3LYP) calculations. The EIE for coordination of CH4 and CD4 does not exhibit typical van't Hoff type behavior in which there is a monotonic variation of EIE with temperature; rather, the temperature dependence of the EIE exhibits a maximum, with inverse values (<1) at low temperature and normal values (>1) at high temperatures. The temperature dependence of the EIE for oxidative addition of CH4 and CD4 differs significantly from that for coordination, with the EIE being normal at all temperatures and approaching infinity at 0 K. In contrast to oxidative addition of methane which is normal at all temperatures, the EIE for oxidative addition of H2 and D2 exhibits a transition from inverse to normal upon raising the temperature. The existence of inverse EIEs in these systems at low temperatures is a result of the zero point energy changes for the products upon isotopic substitution being greater than those for the reactants (H2 or CH4).     

A steric deuterium isotope effect in 1,1,3,3-tetramethylcyclohexane

The equilibrium isotope effect (EIE) for the interconversion of the two chair isotopomers of 1-trideutero-1,3,3-trimethylcyclohexane was predicted using geometry and vibrational force constants derived from electronic structure theory at HF, B3LYP, and MP2 levels as input for the program THERMISTP. Agreement between theory and previously reported NMR results is very good (experimental K(eq) = 1.042 +/- 0.001 vs K(eq) = 1.0409 at MP2/6-311G* level, K(eq) = 1.0503 at HF/6-311G*, and K(eq) = 1.0417 at B3LYP/6-311G* level, all at 17 degrees C). In order to investigate the origin of this isotope effect, the calculated EIEs for the monodeuterated isotopomers were analyzed. It has been shown that the hydrogen atom on an axial methyl group which is unusually close to its counterpart on the other axial methyl is responsible for the large (steric) isotope effect in the compound studied.     

Investigation of hydroxyl-mediated spectral interference from easily ionized elements in laser-enhanced ionization spectrometry in the 281–285 nm range

The spectral background from 281 to 285 nm in the laser-enhanced ionization (LEI) spectrum of aqueous samples containing easily ionized elements (EIE) at concentrations similar to those found in blood was investigated. A complex, structured spectral background was observed, which appears in the presence of Na or K, but does not match the spectral signature of either element. The same behavior was also observed for Li. It was established that this background originates from an energy transfer between laser-excited hydroxyl (OH) molecules and ground-state EIEs. The intensity of this spectral feature was found to increase with EIE concentration and applied electrode voltage. This unexpected source of spectral interference may complicate the determination of trace metals by LEI in the presence of EIEs, since it can not be prevented by simply avoiding interference from atomic lines.     

Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase

It has been suggested that the magnitudes of secondary kinetic isotope effects (2 degrees KIEs) of enzyme-catalyzed reactions are an indicator of the extent of reaction-center rehybridization at the transition state. A 2 degrees KIE value close to the corresponding secondary equilibrium isotope effects (2 degrees EIE) is conventionally interpreted as indicating a late transition state that resembles the final product. The reliability of using this criterion to infer the structure of the transition state is examined by carrying out a theoretical investigation of the hybridization states of the hydride donor and acceptor in the Escherichia coli dihydrofolate reductase (ecDHFR)-catalyzed reaction for which a 2 degrees KIE close to the 2 degrees EIE was reported. Our results show that the donor carbon at the hydride transfer transition state resembles the reactant state more than the product state, whereas the acceptor carbon is more productlike, which is a symptom of transition state imbalance. The conclusion that the isotopically substituted carbon is reactant-like disagrees with the conclusion that would have been derived from the criterion of 2 degrees KIEs and 2 degrees EIEs, but the breakdown of the correlation with the equilibrium isotope effect can be explained by considering the effect of tunneling.     

18O kinetic isotope effects in non-heme iron enzymes: probing the nature of Fe/O2 intermediates

Contrasted here are the competitive 18O/16O kinetic isotope effects (18O KIEs) on kcat/Km(O2) for three non-heme iron enzymes that activate O2 at an iron center coordinated by a 2-His-1-carboxylate facial triad: taurine dioxygenase (TauD), (S)-(2)-hydroxypropylphosphonic acid epoxidase (HppE), and 1-aminocyclopropyl-1-carboxylic acid oxidase (ACCO). Measured 18O KIEs of 1.0102 +/- 0.0002 (TauD), 1.0120 +/- 0.0002 (HppE), and 1.0215 +/- 0.0005 (ACCO) suggest the formation in the rate-limiting step of O2 activation of an FeIII-peroxohemiketal, FeIII-OOH, and FeIV O species, respectively. The comparison of the measured 18O KIEs with calculated or experimental 18O equilibrium isotope effects (18O EIEs) provides new insights into the O2 activation through an inner-sphere mechanism at a non-heme iron center.     

Cavity Catalysis by Cooperative Vibrational Strong Coupling of Reactant and Solvent Molecules

Here, we report the catalytic effect of vibrational strong coupling (VSC) on the solvolysis of para‐nitrophenyl acetate (PNPA), which increases the reaction rate by an order of magnitude. This is observed when the microfluidic Fabry–Perot cavity in which the VSC is generated is tuned to the C=O vibrational stretching mode of both the reactant and solvent molecules. Thermodynamic experiments confirm the catalytic nature of VSC in the system. The change in the reaction rate follows an exponential relation with respect to the coupling strength of the solvent, indicating a cooperative effect between the solvent molecules and the reactant. Furthermore, the study of the solvent kinetic isotope effect clearly shows that the vibrational overlap of the C=O vibrational bands of the reactant and the strongly coupled solvent molecules is critical for the catalysis in this reaction. The combination of cooperative effects and cavity catalysis confirms the potential of VSC as a new frontier in chemistry.     

Vibrational analysis of free and hydrogen bonded complexes of nicotinamide and picolinamide

Harmonic and anharmonic vibrations of free nicotinamide (NIA) and picolinamide (PIA) molecules together with their hydrogen bonded complexes H₂O–NIA and H₂O–PIA have been studied by means of density functional method. The calculation results of the vibrational spectra of free molecules have been investigated and are compared to the available experimental spectra. The vibrational wavenumbers of both molecules have also been calculated by polarizable continuum model (PCM) that represents the solvent as a polarizable continuum and places the solute in a cavity within the solvent (water is chosen as the solvent in this study). The results of PCM calculations and the H₂O–NIA, H₂O–PIA complexes, are used to investigate the H-bonding interactions of both molecules with the water molecule. The harmonic wavenumbers have been scaled by proper factors obtained by comparing the observed versus calculated wavenumbers and it is shown that anharmonic corrections on the vibrational spectra provided a better agreement between the observed and calculated wavenumbers compared to the results obtained by scaling factor method.     

Molecular analysis of water clusters: Calculation of the cluster structures and vibrational spectrum using density functional theory

Density functional theory (DFT) is applied to obtain absorption spectra at THz frequencies for molecular clusters of H₂O. The vibrational modes of the clusters are calculated. Coupling among molecular vibrational modes explains their spectral features associated with THz excitation. THz excitation is associated with vibrational frequencies which are here calculated within the DFT approximation of electronic states. This is done for both isolated molecules and collections of molecules in a cluster. The principal result of the paper is that a crystal-like cluster of 38 water molecules together with a continuum solvent background is sufficient to replicate well the experimental vibrational frequencies.     

Hydration of inorganic phosphates in crystal lattices and in aqueous solution. An experimental and theoretical study

B3LYP/6-31G* and 6-311++G** calculations have been carried out in order to study the hydration of phosphates in aqueous media. Optimized geometries and relative stabilities for PO4(-3), HPO(4)-2, H2PO4(-1) have been calculated considering the interaction with one, two, three, four and five discrete water molecules and taking into account the solvent effect by using the self-consistent reaction field theory (Onsager and PCM methods). The role of specific and bulk contributions of solvent effect on the observable properties of phosphate compounds is analysed. Good agreement between theoretical and available experimental results of harmonic vibration frequencies is found. Significant effects on the geometrical and vibrational frequencies are found for those studied phosphate anions. The results presented here provide a first step toward the understanding of the phosphate group as a hydration sensor in lipid bilayers.     

Electronic structure studies on deprotonation of dithiophosphinic acids in water clusters

We report herein a computational study of proton transfer reactions between dithiophosphinic acids (HAs) and water clusters using B3LYP and MP2 methods. The ground-state and transition-state structures of HA-(H(2)O)(n) (n = 1, 2, 3) cluster complexes have been calculated. The influence of water molecules on energy barrier heights of proton transfer reactions has been examined in the gas phase and solution for bis[o-(trifluoromethyl)phenyl]- and bis(2,4,4-trimethylpentyl)dithiophosphinic acids (HA1 and HA2, respectively). Gas-phase calculations indicate that electron-withdrawing substituents and trifluoromethyl groups in the ortho position favor deprotonation of HA1 when three water molecules are included in the cluster. This suggests that at least three water molecules are necessary to solvate the abstracted proton in the presence of the anion. In the case of HA2, the electron-donating groups favor the reverse proton transfer reaction, namely, protonation of dithiophosphinate anion. Bulk solvent effects have been modeled for aqueous and organic media with the CPCM model. The calculated results show that polar solvents can lower the activation energy for less energetically stable transition states that have more localized charges.     

Hydride reduction of NAD+ analogues by isopropyl alcohol: kinetics, deuterium isotope effects and mechanism

Observed pseudo-first-order rate constants (k(obs)) of the hydride-transfer reactions from isopropyl alcohol (i-PrOH) to two NAD(+) analogues, 9-phenylxanthylium ion (PhXn(+)) and 10-methylacridinium ion (MA(+)), were determined at temperatures ranging from 49 to 82 degrees C in i-PrOH containing various amounts of AN or water. Formations of the alcohol-cation ether adducts (ROPr-i) were observed as side equilibria. The equilibrium constants for the conversion of PhXn(+) to PhXnOPr-i in i-PrOH/AN (v/v = 1) were determined, and the equilibrium isotope effect (EIE = K(i-PrOH)/K(i-PrOD)) at 62 degrees C was calculated to be 2.67. The k(H) of the hydride-transfer step for both reactions were calculated on the basis of the k(obs) and K. The corresponding deuterium kinetic isotope effects (e.g., KIE(OD)(H) = k(H)(i-PrOH)/k(H)(i-PrOD) and KIE(beta-D6)(H) = k(obs)(i-PrOH)/k(obs)((CD3)2CHOH)), as well as the activation parameters, were derived. For the reaction of PhXn(+) (62 degrees C) and MA(+) (67 degrees C), primary KIE(alpha-D)(H) (4.4 and 2.1, respectively) as well as secondary KIE(OD)(H) (1.07 and 1.18) and KIE(beta-D6)(H) (1.1 and 1.5) were observed. The observed EIE and KIE(OD)(H) were explained in terms of the fractionation factors for deuterium between OH and OH(+)(OH(delta+)) sites. The observed inverse kinetic solvent isotope effect for the reaction of PhXn(+) (k(obs)(i-PrOH)/k(obs)(i-PrOD) = 0.39) is consistent with the intermolecular hydride-transfer mechanism. The dramatic reduction of the reaction rate for MA(+), when the water or i-PrOH cosolvent was replaced by AN, suggests that the hydride-transfer T.S. is stabilized by H-bonding between O of the solvent OH and the substrate alcohol OH(delta+). This result suggests an H-bonding stabilization effect on the T.S. of the alcohol dehydrogenase reactions.     

High-precision determination of 18O/16O ratios of silver phosphate by EA-pyrolysis-IRMS continuous flow technique

A high-precision, and rapid on-line method for oxygen isotope analysis of silver phosphate is presented. The technique uses high-temperature elemental analyzer (EA)-pyrolysis interfaced in continuous flow (CF) mode to an isotopic ratio mass spectrometer (IRMS). Calibration curves were generated by synthesizing silver phosphate with a 13 per thousand spread in delta(18)O values. Calibration materials were obtained by reacting dissolved potassium dihydrogen phosphate (KH(2)PO(4)) with water samples of various oxygen isotope compositions at 373 K. Validity of the method was tested by comparing the on-line results with those obtained by classical off-line sample preparation and dual inlet isotope measurement. In addition, silver phosphate precipitates were prepared from a collection of biogenic apatites with known delta(18)O values ranging from 12.8 to 29.9 per thousand (V-SMOW). Reproducibility of +/- 0.2 per thousand was obtained by the EA-Py-CF-IRMS method for sample sizes in the range 400-500 microg. Both natural and synthetic samples are remarkably well correlated with conventional (18)O/(16)O determinations. Silver phosphate is a very stable material and easy to degas and, thus, could be considered as a good candidate to become a reference material for the determination of (18)O/(16)O ratios of phosphate by high-temperature pyrolysis.     

Theoretical studies of UO2(H2O)n(2+), NpO2(H2O)n(+), and PuO2(H2O)n(2+) complexes (n=4-6) in aqueous solution and gas phase

Extensive ab initio calculations both in gas phase and solution have been carried out to study the equilibrium structure, vibrational frequencies, and bonding characteristics of various actinyl (UO2(2+), NpO2(+), and PuO2(2+)) and their hydrated forms, AnO2(H2O)n(z+) (n=4, 5, and 6). Bulk solvent effects were studied using a continuum method. The geometries were fully optimized at the coupled-cluster singles + doubles (CCSD), density-functional theory (DFT), and M?ller-Plesset (MP2) level of theories. In addition vibrational frequencies have been obtained at the CCSD as well as MP2/DFT levels. The results show that both the short-range and long-range solvent effects are important. The combined discrete-continuum model, in which the ionic solute and the solvent molecules in the first and second solvation shells are treated quantum mechanically while the solvent is simulated by a continuum model, can predict accurately the bonding characteristics. Moreover, our values of solvation free energies suggest that five- and six-coordinations are equally preferred for UO2(2+), and five-coordinated species are preferred for NpO2(+) and PuO2(2+). On the basis of combined quantum-chemical and continuum treatments of the hydrated complexes, we are able to determine the optimal cavity radii for the solvation models. The coupled-cluster computations with large basis sets were employed for the vibrational spectra and equilibrium geometries both of which compare quite favorably with experiment. Our most accurate computations reveal that both five- and six-coordination complexes are important for these species.     

A theoretical investigation of the plausibility of reactions between ammonia and carbonyl species (formaldehyde, acetaldehyde, and acetone) in interstellar ice analogs at ultracold temperatures

We have reexamined the reaction between formaldehyde and ammonia, which was previously studied by us and other workers in modestly sized cluster calculations. Larger model systems with up to 12H(2)O were employed, and reactions of two more carbonyl species, acetaldehyde and acetone, were also carried out. Calculations were performed at the B3LYP/6-31+G** level with bulk solvent effects treated with a polarizable continuum model; limited MP2/6-31+G** calculations were also performed. We found that while the barrier for the concerted proton relay mechanism described in previous work remains modest, it is still prohibitively high for the reaction to occur under the ultracold conditions that prevail in dense interstellar clouds. However, a new pathway emerged in more realistic clusters that involves at least one barrierless step for two of the carbonyl species considered here: ammonia reacts with formaldehyde and acetaldehyde to form a partial charge transfer species in small clusters (4H(2)O) and a protonated hydroxyamino intermediate species in large clusters (9H(2)O, 12H(2)O); modest barriers that decrease sharply with cluster size are found for the analogous processes for the acetone-NH(3) reaction. Furthermore, if a second ammonia replaces one of the water molecules in calculations in the 9H(2)O clusters, deprotonation can occur to yield the same neutral hydroxyamino species that is formed via the original concerted proton relay mechanism. In at least one position, deprotonation is barrierless when zero-point energy is included. In addition to describing the structures and energetics of the reactions between formaldehyde, acetaldehyde, and acetone with ammonia, we report spectroscopic predictions of the observable vibrational features that are expected to be present in ice mixtures of different composition.     

Protonated benzofuran,anthracene, naphthalene,benzene, ethene,and ethyne: measurements and estimates of pK(a) and pK(R)

Aqueous solvolyses of acyl derivatives of hydrates (water adducts) of anthracene and benzofuran yield carbocations which undergo competitive deprotonation to form the aromatic molecules and nucleophilic reaction with water to give the aromatic hydrates. Trapping experiments with azide ions yield rate constants k(p) for the deprotonation and k(H2O) for the nucleophilic reaction based on the "azide clock". Combining these with rate constants for (a) the H(+)-catalyzed reaction of the hydrate to form the carbocation and (b) hydrogen isotope exchange of the aromatic molecule (from the literature) yields pK(R) = -6.0 and -9.4 and pK(a) = -13.5 and -16.3 for the protonated anthracene and protonated benzofuran, respectively. These pK values may be compared with pK(R) = -6.7 for naphthalene hydrate (1-hydroxy-1,2-dihydronaphthalene), extrapolated to water from measurements by Pirinccioglu and Thibblin for acetonitrile-water mixtures, and pK(a) = -20.4 for the 2-protonated naphthalene from combining k(p) with an exchange rate constant. The differences between pK(R) and pK(a) correspond to pK(H2O), the equilibrium constant for hydration of the aromatic molecule (pK(H2O) = pK(R) - pK(a)). For naphthalene and anthracene values of pK(H2O) = +13.7 and +7.5 compare with independent estimates of +14.2 and +7.4. For benzene, pK(a) = -24.3 is derived from an exchange rate constant and an assigned value for the reverse rate constant close to the limit for solvent relaxation. Combining this pK(a) with calculated values of pK(H2O) gives pK(R) = -2.4 and -2.1 for protonated benzenes forming 1,2- and 1,4-hydrates, respectively. Coincidentally, the rate constant for protonation of benzene is similar to those for protonation of ethylene and acetylene (Lucchini, V.; Modena, G. J. Am. Chem Soc. 1990, 112, 6291). Values of pK(a) for the ethyl and vinyl cations (-24.8) may thus be derived in the same way as that for the benzenonium ion. Combining these with appropriate values of pK(H2O) then yields pK(R) = -39.8 and -29.6 for the vinyl and ethyl cations, respectively.