Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Thiothionyl bromide (S=SBr2): The elusive isomer of dibromodisulfane (BrSSBr)

CCSD(T) with a series of correlation consistent basis up to quadruple-zeta is used to investigate the structures, vibrational spectra, relative stability, heats of formation, and barrier to isomerization of S=SBr2 and BrSSBr. It represents the most accurate and detailed characterization of these molecules to date. We show that the frequency mode at 302 cm(-1), detected in various studies and assigned to impurities by some authors, and to the anti-symmetric SBr stretch in BrSSBr by others, thus in fact corresponds to the anti-symmetric SBr stretch in the elusive S=SBr2 species; it thus corroborates and complements an earlier partial IR spectra study attributable to S=SBr2. Including corrections for relativistic and core-valence correlation effects, we also predict 26.33 (12.74) kcal/mol for Delta H-f (298.15 K) of S=SBr2 (BrSSBr). For the S=SBr2 -> BrSSBr reaction, our best estimates for the Gibbs free energy and the enthalpy of the reaction at 298.15 K are -13.71 and -13.44 kcal/mol, respectively. For a value of Delta G(#) equal to 23.52 kcal/mol, we estimate a TST rate constant, at 298.15 K, of 3.57 x 10(-5) s(-1). (c) 2007 Elsevier B.V. All rights reserved.

Ab Initio Analysis of Monomers and Dimers of Trialkylphosphine Oxides: Structural and Thermodynamic Stability

Structural and thermodynamic stabilities of monomers and dimers of trialkylphosphine oxides (TRPO) were Studied using quantum chemistry calculations. Density functional theory calculations were carried Out and the structures Of four TRPO have been determined: TMPO (methyl; R = CH(3)), TEPO (ethyl; R = CH(3)CH(2)), TBPO (n-butyl; R = CH(3)(CH(2))(3)), and TOPO (n-octyl; R = CH(3)(CH(2))(7)). TRPO homodimers were investigated considering two isomeric possibilities for each dimer. Relative binding energies and the enthalpic and entropic contributions to the Gibbs free energy were Calculated for all dimers. The formation of dimers from the individual monomeric TRPO species as a function of temperature was also analyzed. (C) 2008 Wiley Periodicals, Inc. Int J Quantum Chem 109: 250-258, 2009

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

Pozzolanic behavior of bamboo leaf ash: Characterization and determination of the kinetic parameters

The paper presents a characterization and study of the pozzolanic behavior between calcium hydroxide (CH) and bamboo leaf ash (BLAsh), which was obtained by calcining bamboo leaves at 600 degrees C for 2 h in a laboratory electric furnace. To evaluate the pozzolanic behavior the conductometric method was used, which is based on the measurement of the electrical conductivity in a BLAsh/CH solution with the reaction time. Later, the kinetic parameters are quantified by applying a kinetic-diffusive model. The kinetic parameters that characterize the process (in particular, the reaction rate constant and free energy of activation) were determined with relative accuracy in the fitting process of the model. The pozzolanic activity is quantitatively evaluated according to the values obtained of the kinetic parameters. Other experimental techniques, such as X-ray diffraction (XRD) and scanning electron microscopy (SEM), were also employed. The results show that this kind of ash is formed by silica with a completely amorphous nature and a high pozzolanic activity. The correlation between the values of free energy of activation (Delta G(#)) and the reaction rate constants (K) are in correspondence with the theoretical studies about the rate processes reported in the literature. (C) 2010 Elsevier Ltd. All rights reserved.

Selection of operational parameters for the production of instant soy protein isolate by pulsed fluid bed agglomeration

The objective of this study was to select the optimal operational conditions for the production of instant soy protein isolate (SPI) by pulsed fluid bed agglomeration. The spray-dried SPI was characterized as being a cohesive powder, presenting cracks and channeling formation during its fluidization (Geldart type A). The process was carried out in a pulsed fluid bed, and aqueous maltodextrin solution was used as liquid binder. Air pulsation, at a frequency of 600 rpm, was used to fluidize the cohesive SPI particles and to allow agglomeration to occur. Seventeen tests were performed according to a central composite design. Independent variables were (i) feed flow rate (0.5-3.5 g/min), (ii) atomizing air pressure (0.5-1.5 bar) and (iii) binder concentration (10-50%). Mean particle diameter, process yield and product moisture were analyzed as responses. Surface response analysis led to the selection of optimal operational parameters, following which larger granules with low moisture content and high process yield were produced. Product transformations were also evaluated by the analysis of size distribution, flowability, cohesiveness and wettability. When compared to raw material, agglomerated particles were more porous and had a more irregular shape, presenting a wetting time decrease, free-flow improvement and cohesiveness reduction. (C) 2010 Elsevier B.V. All rights reserved.

Study of the pozzolanic reaction kinetics in sugar cane bagasse-clay ash/calcium hydroxide system: kinetic parameters and pozzolanic activity

This paper presents a study of the pozzolanic reaction kinetics between calcium hydroxide and a mixture of sugar cane bagasse with 20 and 30% of clay, burned at 800 and 1000 degrees C (SCBCA) by electrical conductivity measurements. A kinetic-diffusive model produced in previous studies by some of the authors was used. The model was fitted to the experimental data, which allowed the computation of the kinetic parameters of the pozzolanic reaction (reaction rate constant and free energy of activation) that rigorously characterised the pozzolanic activity of the materials. The results show that SCBCA demonstrated reactivity and good pozzolanic qualities in the range 800-1000 degrees C.

Thermodynamic and statistical analysis of soybean oil extraction process using renewable solvent

The main goal of this work was to evaluate thermodynamic parameters of the soybean oil extraction process using ethanol as solvent. The experimental treatments were as follows: aqueous solvents with water contents varying from 0 to 13% (mass basis) and extraction temperature varying from 50 to 100 degrees C. The distribution coefficients of oil at equilibrium have been used to calculate enthalpy, entropy and free energy changes. The results indicate that oil extraction process with ethanol is feasible and spontaneous, mainly under higher temperature. Also, the influence of water level in the solvent and temperature were analysed using the response surface methodology (RSM). It can be noted that the extraction yield was highly affected by both independent variables. A joint analysis of thermodynamic and RSM indicates the optimal level of solvent hydration and temperature to perform the extraction process.

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Electronic Structure and Spectro-Structural Correlations of Fe(III)Zn(II) Biomimetics for Purple Acid Phosphatases: Relevance to DNA Cleavage and Cytotoxic Activity

Purple acid phosphatases (PAPs) are a group of metallohydrolases that contain a dinuclear Fe(II)M(II) center (M(II) = Fe, Mn, Zn) in the active site and are able to catalyze the hydrolysis of a variety of phosphoric acid esters. The dinuclear complex [(H(2)O)Fe(III)(mu-OH)Zn(II)(L-H)](CIO(4))(2) (2) with the ligand 2-[N-bis(2-pyridylmethyl)aminomethyl]-4-methyl-6-[N-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl]phenol (H(2)L-H) has recently been prepared and is found to closely mimic the coordination environment of the Fe(III)Zn(II) active site found in red kidney bean PAP (Neves et al. J. Am. Chem. Soc. 2007, 129, 7486). The biomimetic shows significant catalytic activity in hydrolytic reactions. By using a variety of structural, spectroscopic, and computational techniques the electronic structure of the Fe(III) center of this biomimetic complex was determined. In the solid state the electronic ground state reflects the rhombically distorted Fe(III)N(2)O(4) octahedron with a dominant tetragonal compression align ad along the mu-OH-Fe-O(phenolate) direction. To probe the role of the Fe-O(phenolate) bond, the phenolate moiety was modified to contain electron-donating or -withdrawing groups (-CH(3), -H, -Br, -NO(2)) in the 5-position. Tie effects of the substituents on the electronic properties of the biomimetic complexes were studied with a range of experimental and computational techniques. This study establishes benchmarks against accurate crystallographic struck ral information using spectroscopic techniques that are not restricted to single crystals. Kinetic studies on the hydrolysis reaction revealed that the phosphodiesterase activity increases in the order -NO(2)<- Br <- H <- CH(3) when 2,4-bis(dinitrophenyl)phosphate (2,4-bdnpp) was used as substrate, and a linear free energy relationship is found when log(k(cat)/k(0)) is plotted against the Hammett parameter a. However, nuclease activity measurements in the cleavage of double stranded DNA showed that the complexes containing the electron-withdrawing -NO(2) and electron-donating CH3 groups are the most active while the cytotoxic activity of the biomimetics on leukemia and lung tumoral cells is highest for complexes with electron-donating groups.

Structural modeling of high-affinity thyroid receptor-ligand complexes

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors` ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

The number of open paths in an oriented rho-percolation model

We study the asymptotic properties of the number of open paths of length n in an oriented rho-percolation model. We show that this number is e(n alpha(rho)(1+o(1))) as n ->infinity. The exponent alpha is deterministic, it can be expressed in terms of the free energy of a polymer model, and it can be explicitly computed in some range of the parameters. Moreover, in a restricted range of the parameters, we even show that the number of such paths is n(-1/2)We (n alpha(rho))(1+o(1)) for some nondegenerate random variable W. We build on connections with the model of directed polymers in random environment, and we use techniques and results developed in this context.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved