Accurate gas phase acidities of carboxylic acids estimated by scaling the vibrational contribution of ab initio gibbs free energies

The gas phase Gibbs free energies deltaG(T) of dissociation reaction of 14 carboxylic acids were calculated on the SCF, as well as G3 and CBS-Q levels. Corresponding accuracies were critically compared with experimental data. Since all of the results suffer from systematic errors, the procedure of scaling of thermal contribution to Gibbs free energy was applied for minimizing differences between theoretical and experimental values of deltaG(T). Two parameters were adjusted, namely the scaling of thermal contribution to Gibbs free energy of neutral and anionic forms. The presented results suggest the great effectiveness of such a procedure since for all applied basis sets within the SCF framework the achieved accuracy was below the experimental error. Besides, the proposed low-cost approximation method leads to precision comparable to or even exceeding the quality offered by more sophisticated composite quantum chemistry methods. The extension of the set of training molecules up to 82 has an insignificant impact on the overall quality of deltaG(T) estimation, which suggests that a wisely chosen set of reference data may be used for the characteristics of the whole class of compounds. There is a straightforward way for the analysis of acidities/basicities of other classes of chemicals such as DNA bases, alcohols, phenols, amines, amino acids, etc.     

Thermodynamic analysis of trinitrotoluene biodegradation and mineralization pathways

Biodegradation of 2,4,6-trinitrotoluene (TNT) proceeds through several different metabolic pathways. However, the reaction steps which are considered rate-controlling have not been fully determined. Glycolysis and other biological pathways contain biochemical reactions which are acutely rate-limiting due to enzyme control. These rate-limiting steps also have large negative Gibbs free energy changes. Because xenobiotic compounds such as TNT can be used by biological systems as nitrogen, carbon, and energy sources, it is likely that their degradation pathways also contain acutely rate-limiting steps. Identification of these rate-controlling reactions will enhance and better direct genetic engineering techniques to increase specific enzyme levels.This article identifies likely rate-controlling steps (or sets of steps) in reported TNT biodegradation pathways by estimating the Gibbs free energy change for each step and for the overall pathways. The biological standard Gibbs free energy change of reaction was calculated for each pathway step using a group contribution method specifically tailored for biomolecules. The method was also applied to hypothetical "pathways" constructed to mineralize TNT using several different microorganisms. Pathways steps that have large negative Gibbs free energy changes are postulated to be potentially rate-controlling. The microorganisms which utilize degradation pathways with the largest overall (from TNT to citrate) negatiave Gibbs free energy changes were also determined. Such microorganisms can extract more energy from the starting substrate and are thus assumed to have a competitive advantage over other microorganisms. Results from this modeling-based research are consistent with much experimental work available in the literature. (c) 1996 John Wiley & Sons, Inc.     

DFT study of the V(IV)/V(V) oxidation mechanism in the presence of N-hydroxyacetamide

The oxidation mechanism of V(IV)/V(V) in the presence of N-hydroxyacetamide (acetohydroxamic acid, HL) in aqueous solution has been investigated using density functional theory (DFT) calculations aiming to contribute to the understanding of this process at a molecular level. The mechanism proposed involves formation of the *OH, *OOH, H2O2 radicals and complexes formed from the interaction of these species with VOL2 complex. The Gibbs free energy of each step of the mechanism has been evaluated. The solvation energy has been estimated by means of united atoms Hartree-Fock/polarizable continuum method (UAHF/PCM). The Gibbs free energy of the global reaction of the V(IV)/V(V) oxidation has been estimated and compared with the available experimental equilibrium constant. The difference between the calculated and experimental estimates for the reaction energy of the global equation is about 1.5 kcal mol(-1). The thermodynamic profile of the reaction mechanism has been provided and discussed in terms of the possible intermediates. The influence of the ligand and the reaction rate in terms of the steady-state approximation has been briefly discussed.     

On the calculation of Gibbs energy corresponding to enantioselective interactions at a direct HRGC separation of enantiomers

A novel procedure is proposed for the calculation of Gibbs energy corresponding to enantiospecific interactions of 2-(2, 4-dinitrophenoxy)-, 2-phenoxy-, and 2-halogen-n-pentane enantiomers with a beta-cyclodextrin (ChirasilDex) stationary phase under gas chromatographic conditions. This energy is calculated from retention data as a difference between the Gibbs energy of an enantiomer and its corresponding achiral congener. The procedure for the determination of 2-(2,4-dinitrophenoxy)-, 2-phenoxy- and 2-halogen- n-pentane achiral congener retention data is discussed in detail.     

Computation of the bond dissociation enthalpies and free energies of hydroxylic antioxidants using the ab initio Hartree–Fock method

Abstract

Introduction

A new method for calculating theoretical bond dissociation enthalpy (BDE) and bond dissociation free energy (BDFE) of hydroxylic antioxidants is forwarded. BDE and BDFE may be understood as activation energies accompanying the formation of transition states, which may undergo downhill homolytic dissociation. The new method does not involve the complete fission of O–H bonds.

Method

Theoretical gas phase BDE values were calculated with the ab initio unrestricted Hartree–Fock (UHF) method, as changes in enthalpy between ground singlet states (GS) and triplet dissociative states (DS). Similarly, gas phase BDFEs were estimated from the corresponding changes in Gibbs free energy. The results were then compared with reliable experimental reports.

Results

The proposed theoretical approach of BDE and BDFE determination was tested using 10 simple phenols, 5 flavonoids, and l-ascorbic acid derivatives. The agreement between our calculated gas phase results and the adopted experimental values were generally within 0.5 kcal mol^?1, with a very few exceptions.

Discussion

Generally, steric interactions as well as intramolecular hydrogen bonding involving the dissociating OH group should be minimized in the GS. The DS are both electronically and vibrationally exited transition states. They have one unpaired electron on the carbon atom, which bears the homolytically dissociating OH group and are second order saddle points with a fixed <C–O–H bond angel of 180°.

Conclusion

It was concluded that ab initio UHF was well suited for the estimation of gas phase BDE and BDFE. The method presented has a good potential for application across a range of hydroxylic antioxidants. Currently, work is underway to extend its application in other class of antioxidants.     

Redox Potential as a Parameter To Predict the Reductive Dechlorination Pathway of Chloroanilines in Anaerobic Environments

Abstract The hypothesis that the microbially catalyzed pathway proceeds with a step that would yield the highest energy was examined for the reductive dechlorination of chloroanilines (CAs) under anaerobic conditions. The Gibbs free energy of formation was estimated with Benson's method, then the redox potentials were determined using an H⁺/H₂ couple as the reference system. The observed pathways were compared using the redox potential of the reaction, and the results showed that the redox potential correctly predicts the pathway that yields the highest energy for the dechlorination step. Received: 5 April 1996; Accepted: 9 August 1996     

DFT and MP2 study of low barrier proton transfer in hydrazide schiff base tautomers via water bridges and in the gas

MP2 and DFT studies were performed on the tautomers of N′-ethylideneacetohydrazide in different environments including gas phase, continuum solvent and microhydrated environment. The ground electronic state structures of the tautomers were optimized at the MP2 and B3LYP levels of theory using 6-311++G(d,p), separately. The optimized geometries of the transition states of different tautomerism processes, which occur through the proton transfer (PT) reaction, were determined using the QST3 approach at the same levels of theory. The same stability order as T1Z〉 T1E〉 T2ZE〉 T2ZZ〉 T2EE〉 T2EZ〉 T6〉 T4E was found for the tautomers in the gas phase, continuum solvent and microhydrated environment for both B3LYP and MP2 levels of theory. In addition, the variations of the Gibbs free energies of tautomeric processes, the activation Gibbs free energies of the forward and reverse tautomeric processes with the polarity of the solvent (in continuum solvent model) and the number of water molecules (in microhydrated environment) were investigated. It was found that the reverse tautomeric process is more favorable in all considered environments thermodynamically and kinetically. In addition, it was shown that the rate constants of the reverse and forward considered tautomeric processes decrease with the solvent polarity in the continuum solvent model and the process becomes more difficult than the gas phase. The opposite trend is seen in the microhydrated environment.     

Oxidative addition of methyl iodide to [Rh(CH3COCHCOCH3)(CO)(P(OCH2)3CCH3)]

The reaction rate of the oxidative addition and the following CO insertion step of methyl iodide with [Rh(acac)(CO)(P(OCH₂)₃CCH₃)] is determined. The key finding is that while [Rh(acac)(CO)(P(OCH₂)₃CCH₃)] oxidatively adds methyl iodide ca 300 times faster than the Monsanto catalyst, the CO insertion step is much slower. However, the rate-determining step of the oxidative addition reaction of the phosphorus-containing acetylacetonato-rhodium(I) complex, the carbonyl insertion step, is still in the same order or faster than the rate-determining oxidative addition step of iodomethane to [Rh(CO)₂I₂]⁻.     

Biomass as an energy source: thermodynamic constraints on the performance of the conversion process

In the present work an equilibrium model (gas-solid), based on the minimization of the Gibbs energy, has been used in order to estimate the theoretical yield and the equilibrium composition of the reaction products (syngas and char) of biomass thermochemical conversion processes (pyrolysis and gasification). The data obtained from this model have also been used to calculate the heating value of the fuel gas, in order to evaluate the overall energy efficiency of the thermal conversion stage. The proposed model has been applied both to partial oxidation and steam gasification processes with varying air to biomass (ER) and steam to carbon (SC) ratio values and using different feedstocks; the obtained results have been compared with experimental data and with other model predictions obtaining a satisfactory agreement.     

Quantitative assignment of reaction directionality in a multicompartmental human metabolic reconstruction

Reaction directionality is a key constraint in the modeling of genome-scale metabolic networks. We thermodynamically constrained reaction directionality in a multicompartmental genome-scale model of human metabolism, Recon 1, by calculating, in vivo, standard transformed reaction Gibbs energy as a function of compartment-specific pH, electrical potential, and ionic strength. We show that compartmental pH is an important determinant of thermodynamically determined reaction directionality. The effects of pH on transport reaction thermodynamics are only seen to their full extent when metabolites are represented as pseudoisomer groups of multiple protonated species. We accurately predict the irreversibility of 387 reactions, with detailed propagation of uncertainty in input data, and manually curate the literature to resolve conflicting directionality assignments. In at least half of all cases, a prediction of a reversible reaction directionality is due to the paucity of compartment-specific quantitative metabolomic data, with remaining cases due to uncertainty in estimation of standard reaction Gibbs energy. This study points to the pressing need for 1), quantitative metabolomic data, and 2), experimental measurement of thermochemical properties for human metabolites.     

Determination of the thermodynamic parameters of the complex formation between malvidin-3-O-glucoside and polyphenols. Copigmentation effect in red wines

The thermodynamics of the molecular association process between the malvidin-3-O-glucoside and a series of polyphenol derivatives (called 'copigmentation' in food chemistry) were studied in aqueous media. The Gibbs free energy, enthalpy and entropy values were determined by the fluorometric method. A combination of the Job's method with the van't Hoff theory was applied for data evaluation. The results show the exothermic character of the copigmentation process. The change of the enthalpy seems to be the same in every complexation step. However, the decreasing of the entropy term is higher at higher stoichiometries. As a result, the Gibbs free energy changes and, thus, the complex stability decreases quickly with increasing stoichiometry. Quantum-chemical investigation reveals the complexity of molecular interactions between malvidin and polyphenols, which is preferably based on pi-pi and OH-pi interaction moderated by repulsive Coulomb-type interactions.     

Reductive unfolding of serum albumins uncovered by Raman spectroscopy

The reductive unfolding of bovine serum albumin (BSA) and human serum albumin (HSA) induced by dithiothreitol (DTT) is investigated using Raman spectroscopy. The resolution of the S-S Raman band into both protein and oxidized DTT contributions provides a reliable basis for directly monitoring the S-S bridge exchange reaction. The related changes in the protein secondary structure are identified by analyzing the protein amide I Raman band. For the reduction of one S-S bridge of BSA, a mean Gibbs free energy of -7 kJ mol(-1) is derived by studying the reaction equilibrium. The corresponding value for the HSA S-S bridge reduction is -2 kJ mol(-1). The reaction kinetics observed via the S-S or amide I Raman bands are identical giving a reaction rate constant of (1.02 +/- 0.11) M(-1) s(-1) for BSA. The contribution of the conformational Gibbs free energy to the overall Gibbs free energy of reaction is further estimated by combining experimental data with ab initio calculations.     

Optimization of the allylic oxidation in the synthesis of 7-keto-delta5-steroidal substrates

A variety of delta5-steroids were converted into alpha, beta-unsaturated 7-ketones using a modification of the already known method of t-butyl hydroperoxide in the presence of copper iodide in acetonitrile. The same alteration was applied to another oxidative procedure, which had never been used before on steroidal substrates. The same oxidative agent was used in the presence of copper iodide, and tetra-n-butylammonium bromide was used as a phase-transfer catalyst in a two-phase system of water/methylene chloride. It was found that the allylic oxidation proceeded more efficiently when t-butyl hydroperoxide was added to the reaction mixture in portions. The initial addition of the total amount of oxidant or its dropwise addition afforded low yields. This observation contributes to the investigation of the reaction mechanism, and high-yield conversions of steroidal 5,6-enes into the corresponding conjugated 7-ones in short reaction times are reported.     

Analysis of ligand binding curves on basis of mean intrinsic thermodynamic quantities

The mean intrinsic thermodynamic quantity can be defined by considering the relative population of complex species in the solution and the value of intrinsic thermodynamic quantity corresponds to each step of ligation. In the present study a new method is introduced for analysis of experimental ligand binding data on basis of mean intrinsic thermodynamic quantities. In this regard, a deviation parameter was defined by comparing the non-interacting system with the cooperative interactive one. This parameter can be calculated just by estimation of the first binding constant. A set of relations between this deviation parameter and other binding characteristics, such as mean intrinsic Gibbs free energy of binding and mean Gibbs free energy of site-site interaction, have been developed. This model presents binding mechanism in a unified way that is simple, yet stringent, more straightforward, more reliable and informative. This analyzing method has been successfully applied for evaluation of various systems such as oxygen binding to hemoglobin, laurate and warfarin binding to human serum albumin, and reveals some new biological features of these binding systems.     

Gas-liquid chromatographic determination of antibiotic X-537A, lasalocid

A method is described for the chemical analysis of lasalocid, formerly known as antibiotic X-537A. The method is based on quantitative thermolysis of the antibiotic to a retroaldol ketone which is assayed by gas chromatography. No catalyst or special pyrolysis attachment is needed for the assay, as cleavage of the antibiotic is complete in the injection port of the gas chromatograph. This technique has the advantage that no derivatization step is necessary prior to analysis.     

5-Enolpyruvylshikimate-3-phosphate synthase: determination of the protonation state of active site residues by the semiempirical method

EPSP synthase (EPSPS) catalyzes the addition of shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP) to form a tetrahedral intermediate (TI) that is converted to 5-enolpyruvylshikimate-3-phosphate (EPSP) and inorganic phosphate. A semiempirical molecular modeling study of the EPSPS active site containing the TI was implemented for the assignment of the protonation states of four basic residues, Lys22, Lys340, His385, and Lys411, based on the evaluation of 16 different protonation states and comparison of the resulting energy minimized heavy atoms coordinates with available X-ray crystallographic data of the D313A mutant of EPSPS. The results, employing both gas phase and continuum solvent models, are indicative that after the TI formation the histidine residue is most probably in neutral form (N^ε-protonated) and the lysine residues are in protonated form, which suggests that none of the presently proposed assignments of aminoacid residues involved in the reaction mechanism could be completely correct. The protonated state of Lys22 in the presence of the TI supports the proposal that this residue is a general acid catalyst for TI breakdown. Modeling of the native enzyme active site suggests that Asp313 residue has only minor effects on the definition of the TI position inside the active site. Hydrogen-bonds distances suggest that, in order to act as a base, Asp313 needs the intermediacy of a hydroxyl group of the TI for effecting the attack on the TI methyl group in the elimination step leading to EPSP, as suggested previously in the literature.     

Water plays a different role on activation thermodynamic parameters of alcoholysis reaction catalyzed by lipase in gaseous and organic media

The effect of water on the alcoholysis of methyl propionate and n-propanol catalyzed by immobilized Candida antarctica lipase B (CALB) has been compared in a continuous solid-gas reactor and in an organic liquid medium. The enthalpic and entropic contributions of water to the Gibbs free energy of activation in the gas phase were different from the ones in the organic phase, the inverse trends being observed for the variation of both DeltaH* and DeltaS* with water activity.Different phenomena were identified for their influence on the thermodynamic parameters. When increasing a(w), the enhanced flexibility of the enzyme was predominant in the gas phase whereas substrate-solvent interactions due to an increased polarity of the solvent affected mainly the thermodynamic parameters in the organic phase. The observed variations of DeltaG* with water activity were in accordance with kinetics results previously obtained in both reaction media.     

Equilibrium calculations on systems of biochemical reactions at specified pH and pMg.

The use of G' in discussing the thermodynamics of biochemical reactions at a specified pH and pMg is justified by use of a Legendre transform of the Gibbs energy G. When several enzymatic reactions occur simultaneously in a system, the standard transformed Gibbs energies of reaction delta rG'0 can be used in a computer program to calculate the equilibrium composition that minimizes the transformed Gibbs energy at the specified pH and pMg. The calculation of standard transformed Gibbs energies of formation of reactants at pH 7 and pMg 3 is described. In addition a method for calculating the equilibrium concentrations of reactants is illustrated for a system with steady state concentrations of some reactants like ATP and NAD.     

Thermodynamic analysis of biodegradation pathways

Microorganisms provide a wealth of biodegradative potential in the reduction and elimination of xenobiotic compounds in the environment. One useful metric to evaluate potential biodegradation pathways is thermodynamic feasibility. However, experimental data for the thermodynamic properties of xenobiotics is scarce. The present work uses a group contribution method to study the thermodynamic properties of the University of Minnesota Biocatalysis/Biodegradation Database. The Gibbs free energies of formation and reaction are estimated for 914 compounds (81%) and 902 reactions (75%), respectively, in the database. The reactions are classified based on the minimum and maximum Gibbs free energy values, which accounts for uncertainty in the free energy estimates and a feasible concentration range relevant to biodegradation. Using the free energy estimates, the cumulative free energy change of 89 biodegradation pathways (51%) in the database could be estimated. A comparison of the likelihood of the biotransformation rules in the Pathway Prediction System and their thermodynamic feasibility was then carried out. This analysis revealed that when evaluating the feasibility of biodegradation pathways, it is important to consider the thermodynamic topology of the reactions in the context of the complete pathway. Group contribution is shown to be a viable tool for estimating, a priori, the thermodynamic feasibility and the relative likelihood of alternative biodegradation reactions. This work offers a useful tool to a broad range of researchers interested in estimating the feasibility of the reactions in existing or novel biodegradation pathways. Biotechnol. Bioeng. 2009;103: 532–541. © 2009 Wiley Periodicals, Inc.