Christoph Held

Compatible solutes: Thermodynamic properties and biological impact of ectoines and prolines

Compatible solutes like ectoine and its derivatives are deployed by halophile organisms as osmolytes to sustain the high salt concentration in the environment. This work investigates the relation of the thermodynamic properties of compatible solutes and their impact as osmolytes. The ectoines considered in this work are ectoine, hydroxyectoine, and homoectoine. Besides solution densities (15-45°C) and solubilities in water (3-80°C), component activity coefficients in the aqueous solutions were determined in the temperature range between 0 and 50°C. The latter is important for adjusting a certain water activity and therewith a respective osmotic pressure within a cell. The characteristic effect of ectoines is compared to that of prolines, as well as to that of incompatible solutes as salts and urea. The experimental results show that the influence on the activity (coefficient) of water is quite different for compatible and incompatible solutes: whereas compatible solutes cause decreasing water activity coefficients, incompatible solutes lead to an increase in water activity coefficients. Based on this quantity, the paper discusses the impact of various osmolytes on biological systems and contributes to the explanation why some osmolytes are more often and at other temperatures used than others. Moreover, it was found that the anti-stress effect of an osmolyte is weakened in the presence of a salt. Finally, it is shown that the thermodynamic properties of compatible solutes can be modeled and even predicted using the thermodynamic model PC-SAFT (Perturbed-Chain Statistical Associating Fluid Theory).

Solubility of Sugars and Sugar Alcohols in Ionic Liquids: Measurement and PC-SAFT Modeling

Biorefining processes using ionic liquids (ILs) require proper solubility data of biomass-based compounds in ILs, as well as an appropriate thermodynamic approach for the modeling of such data. Carbohydrates and their derivatives such as sugar alcohols represent a class of compounds that could play an important role in biorefining. Thus, in this work, the pure IL density and solubility of xylitol and sorbitol in five different ILs were measured between 288 and 339 K. The ILs under consideration were 1-ethyl-3-methylimidazolium dicyanamide, 1-butyl-3-methylimidazolium dicyanamide ([bmim][DCA]), Aliquat dicyanamide, trihexyltetradecylphosphonium dicyanamide, and 1-ethyl-3-methylimidazolium trifluoroacetate. Comparison with the literature data was performed, showing good agreement. With the exception of [bmim][DCA], the solubility of these sugar alcohols in the other ILs is presented for the first time. The measured data as well as previously published solubility data of glucose and fructose in these ILs were modeled by means of PC-SAFT using a molecular-based associative approach for ILs. PC-SAFT was used in this work as it has shown to be applicable to model the solubility of xylitol and sorbitol in ILs (Paduszyński; et al. J. Phys. Chem. B 2013, 117, 7034-7046). For this purpose, three pure IL parameters were fitted to pure IL densities, activity coefficients of 1-propanol at infinite dilution in ILs, and/or xylitol solubility in ILs. This approach allows accurate modeling of the pure IL data and the mixture data with only one binary interaction parameter k(ij) between sugar and the IL or sugar alcohol and the IL. In cases where only the pure IL density and activity coefficients of 1-propanol at infinite dilution in ILs were used for the IL parameter estimation, the solubility of the sugars and sugar alcohols in the ILs could be predicted (k(ij) = 0 between sugar and the IL or sugar alcohol and the IL) with reasonable accuracy.

Molecular interactions in 1-butanol + IL solutions by measuring and modeling activity coefficients.

Molecular interactions in 1-butanol + ionic liquid (IL) solutions have been investigated by measuring and modeling activity-coefficient data. The activity coefficients in binary solutions containing 1-butanol and an IL were determined experimentally: the ILs studied were 1-decyl-3-methyl-imidazolium tetracyanoborate ([Im10.1](+)[tcb](-)), 4-decyl-4-methyl-morpholinium tetracyanoborate ([Mo10.1](+)[tcb](-)), 1-decyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ([Im10.1](+)[ntf2](-)), and 4-decyl-4-methyl-morpholinium bis(trifluoromethylsulfonyl)imide ([Mo10.1](+)[ntf2](-)). The methods used to determine the activity coefficients included vapor-pressure osmometry, headspace-gas chromatography, and gas-liquid chromatography. The results from all of these techniques were combined to obtain activity-coefficient data over the entire IL concentration range, and the ion-specific interactions of the ILs investigated were identified with 1-butanol. The highest (1-butanol)-IL interactions of the ILs considered in this work were found for [Im10.1](+)[tcb](-); thus, [Im10.1](+)[tcb](-) showed the highest affinity for 1-butanol in a binary mixture. The experimental data were modeled with the Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT). PC-SAFT was able to accurately describe the pure IL and (1-butanol)-IL data. Moreover, the model was shown to be predictive and extrapolative with respect to concentration and temperature.

The role of activity coefficients in bioreaction equilibria: Thermodynamics of methyl ferulate hydrolysis

The Gibbs energy of reaction (Δ(R)g) is the key quantity in the thermodynamic characterization of biological reactions. Its calculation requires precise standard Gibbs energy of reaction (Δ(R)g(+)) values. The value of Δ(R)g(+) is usually determined by measuring the apparent (concentration-dependent) equilibrium constants K, e.g., the molality-based Km. However, the thermodynamically consistent determination of Δ(R)g(+) requires the thermodynamic (activity-based) equilibrium constant Ka. These values (Km and Ka) are equal only if the ratio of the activity coefficients of the reactants to the activity coefficients of the products (Kγ) is equal to unity. In this work, the impact of Kγ on the estimation of Ka for biological reactions was investigated using methyl ferulate (MF) hydrolysis as a model reaction. The value of Kγ was experimentally determined from Km values that were measured at different reactant concentrations. Moreover, Kγ was independently predicted using the thermodynamic model ePC-SAFT. Both the experimentally determined and the predicted Kγ values indicate that this value cannot be assumed to be unity in the considered reaction. In fact, in the reaction conditions considered in this work, Kγ was shown to be in the range of 3<Kγ<6 for different reactant molalities (2<mmol MF kg(-1)<10). The inclusion of Kγ and thus the use of the thermodynamically correct Ka value instead of Km lead to remarkable differences (almost 40%) in the determination of Δ(R)g(+). Moreover, the new value for Δ(R)g(+) increases the concentration window at which the reaction can thermodynamically occur. The influence of additives was also investigated both experimentally and theoretically. Both procedures consistently indicated that the addition of NaCl (0 to 1molkg(-1) water) moderately decreased the value of Kγ, which means that the values of Km increase and that a higher amount of products is obtained as a result of the addition of salt. Additionally, Km was found to strongly depend on pH. A ten-fold increase in the Km values was observed in the pH range of 6 to 7; this increase corresponds to a change of more than 100% in the value of Δ(R)g(+).

PC-SAFT Modeling of CO2 Solubilities in Deep Eutectic Solvents

Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT), a physically based model that accounts for different molecular interactions explicitly, was applied to describe for the first time the phase behavior of deep eutectic solvents (DESs) with CO2 at temperatures from 298.15 to 318.15 K and pressures up to 2 MPa. DESs are mixtures of two solid compounds, a hydrogen bond donor (HBD) and a hydrogen bond acceptor (HBA), which form liquids upon mixing with melting points far below that of the individual compounds. In this work, the HBD is lactic acid and the HBAs are tetramethylammonium chloride, tetraethylammonium chloride, and tetrabutylammonium chloride. Two different modeling strategies were considered for the PC-SAFT modeling. In the first strategy, the so-called pseudo-pure component approach, a DES was considered as a pseudo-pure compound, and its pure-component parameters were obtained by fitting to pure DES density data. In the second strategy, the so-called individual-component approach, a DES was considered to consist of two individual components (HBA and HBD), and the pure-component parameters of the HBA and HBD were obtained by fitting to the density of aqueous solutions containing only the individual compounds of the DES. In order to model vapor-liquid equilibria (VLE) of DES + CO2 systems, binary interaction parameters were adjusted to experimental data from the literature and to new data measured in this work. It was concluded that the individual-component strategy allows quantitative prediction of the phase behavior of DES + CO2 systems containing those HBD:HBA molar ratios that were not used for k(ij) fitting. In contrast, applying the pseudo-pure component strategy required DES-composition specific k(ij) parameters.

A thermodynamic investigation of the glucose-6-phosphate isomerization

In this work, Δ(R)g(+) values for the enzymatic G6P isomerization were determined as a function of the G6P equilibrium molality between 25 °C and 37 °C. The reaction mixtures were buffered at pH=8.5. In contrast to standard literature work, Δ(R)g(+) values were determined from activity-based equilibrium constants instead of molality-based apparent values. This yielded a Δ(R)g(+) value of 2.55±0.05 kJ mol(-1) at 37 °C, independent of the solution pH between 7.5 and 8.5. Furthermore, Δ(R)h(+) was measured at pH=8.5 and 25 °C yielding 12.05±0.2 kJ mol(-1). Accounting for activity coefficients turned out to influence Δ(R)g(+) up to 30% upon increasing the G6P molality. This result was confirmed by predictions using the thermodynamic model ePC-SAFT. Finally, the influence of the buffer and of potassium glutamate as an additive on the reaction equilibrium was measured and predicted with ePC-SAFT in good agreement.

Benchmark Thermochemistry for Biologically Relevant Adenine and Cytosine. A Combined Experimental and Theoretical Study.

The thermochemical properties available in the literature for adenine and cytosine are in disarray. A new condensed phase standard (p° = 0.1 MPa) molar enthalpy of formation at T = 298.15 K was measured by using combustion calorimetry. New molar enthalpies of sublimation were derived from the temperature dependence of vapor pressure measured by transpiration and by the quarz-crystal microbalance technique. The heat capacities of crystalline adenine and cytosine were measured by temperature-modulated DSC. Thermodynamic data on adenine and cytosine available in the literature were collected, evaluated, and combined with our experimental results. Thus, the evaluated collection of data together with the new experimental results reported here has helped to resolve contradictions in the available enthalpies of formation. A set of reliable thermochemical data is recommended for adenine and cytosine for further thermochemical calculations. Quantum-chemical calculations of the gas phase molar enthalpies of formation of adenine and cytosine have been performed by using the G4 method and results were in excellent agreement with the recommended experimental data. The standard molar entropies of formation and the standard molar Gibbs functions of formation in crystal and gas state have been calculated. Experimental vapor-pressure data measured in this work were used to estimate pure-component PC-SAFT parameters. This allowed modeling solubility of adenine and cytosine in water over the temperature interval 278-310 K.

Thermodynamics of Bioreactions

Thermodynamic principles have been applied to enzyme-catalyzed reactions since the beginning of the 1930s in an attempt to understand metabolic pathways. Currently, thermodynamics is also applied to the design and analysis of biotechnological processes. The key thermodynamic quantity is the Gibbs energy of reaction, which must be negative for a reaction to occur spontaneously. However, the application of thermodynamic feasibility studies sometimes yields positive Gibbs energies of reaction even for reactions that are known to occur spontaneously, such as glycolysis. This article reviews the application of thermodynamics in enzyme-catalyzed reactions. It summarizes the basic thermodynamic relationships used for describing the Gibbs energy of reaction and also refers to the nonuniform application of these relationships in the literature. The review summarizes state-of-the-art approaches that describe the influence of temperature, pH, electrolytes, solvents, and concentrations of reacting agents on the Gibbs energy of reaction and, therefore, on the feasibility and yield of biological reactions.

Crowders and Cosolvents - Major Contributors to the Cellular Milieu and Efficient Means to Counteract Environmental Stresses

The free energy and conformational landscape of biomolecular systems as well as biochemical reactions depend not only on temperature and pressure, but also on the particular solution conditions. Such conditions include the effects of cosolvents (for example osmolytes) and macromolecular crowding, which are crucial components to understand the energetics and kinetics of biological processes in living system. Such conditions are also important for the understanding of many debilitating diseases, such as those where misfolding and amyloid formation of proteins are involved. Moreover, understanding their effects on biomolecular processes is prerequisite for designing industrially relevant enzymatic reactions, which seldom take place under neat conditions. Here, we review and discuss experimental and theoretical studies on the characterization of cosolvent and crowding induced effects in biologically relevant systems, approaching even the complexity of living organisms. In particular, we focus on cosolvent and crowding effects on the conformational equilibrium and folding kinetics of proteins and nucleic acids as well as on enzymatic reactions, including their effects on the temperature and pressure dependence of these processes. By presenting a few representative examples, we show how such effects are unveiled and described in thermodynamic and kinetic terms.

Benzoic Acid and Chlorobenzoic Acids: Thermodynamic Study of the Pure Compounds and Binary Mixtures With Water

Benzoic acid is a model compound for drug substances in pharmaceutical research. Process design requires information about thermodynamic phase behavior of benzoic acid and its mixtures with water and organic solvents. This work addresses phase equilibria that determine stability and solubility. In this work, Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) was used to model the phase behavior of aqueous and organic solutions containing benzoic acid and chlorobenzoic acids. Absolute vapor pressures of benzoic acid and 2-, 3-, and 4-chlorobenzoic acid from literature and from our own measurements were used to determine pure-component PC-SAFT parameters. Two binary interaction parameters between water and/or benzoic acid were used to model vapor-liquid and liquid-liquid equilibria of water and/or benzoic acid between 280 and 413 K. The PC-SAFT parameters and 1 binary interaction parameter were used to model aqueous solubility of the chlorobenzoic acids. Additionally, solubility of benzoic acid in organic solvents was predicted without using binary parameters. All results showed that pure-component parameters for benzoic acid and for the chlorobenzoic acids allowed for satisfying modeling phase equilibria. The modeling approach established in this work is a further step to screen solubility and to predict the whole phase region of mixtures containing pharmaceuticals.