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.

Solvation in Pure Liquids: What Can Be Learned from the Use of Pairs of Indicators?

The solvation of six solvatochromic probes in a large number of solvents (33-68) was examined at 25 degrees C. The probes employed were the following: 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-triphenyl pyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), respectively. Of these, MePMBr is a novel compound. They can be grouped in three pairs, each with similar pK(a) in water but with different molecular properties, for example, lipophilicity and dipole moment. These pairs are formed by RB and MePM; QB and MePMBr; WB and MePMBr(2), respectively. Theoretical calculations were carried out in order to calculate their physicochemical properties including bond lengths, dihedral angles, dipole moments, and wavelength of absorption of the intramolecular charge-transfer band in four solvents, water, methanol, acetone, and DMSO, respectively. The data calculated were in excellent agreement with available experimental data, for example, bond length and dihedral angles. This gives credence to the use of the calculated properties in explaining the solvatochromic behaviors observed. The dependence of an empirical solvent polarity scale E(T)(probe) in kcal/mol on the physicochemical properties of the solvent (acidity, basicity, and dipolarity/polarizability) and those of the probes (pK(a), and dipole moment) was analyzed by using known multiparameter solvation equations. For each pair of probes, values of E(T)(probe) (for example, E(T)(MePM) versus E(T)(RB)) were found to be linearly correlated with correlation coefficients, r, between 0.9548 and 0.9860. For the mercyanine series, the values of E(T)(probe) also correlated linearly, with (r) of 0.9772 (MePMBr versus MePM) and 0.9919 (MePMBr(2) versus MePM). The response of each pair of probes (of similar pK(a)) to solvent acidity is the same, provided that solute-solvent hydrogen-bonding is not seriously affected by steric crowding (as in case of RB). We show, for the first time, that the response to solvent dipolarity/polarizability is linearly correlated to the dipole moment of the probes. The successive introduction of bromine atoms in MePM (to give MePMBr, then MePMBr(2)) leads to the following linear decrease: pK(a) in water, length of the phenolate oxygen-carbon bond, length of the central ethylenic bond, susceptibility to solvent acidity, and susceptibility to solvent dipolarity/polarizability. Thus studying the solvation of probes whose molecular structures are varied systematically produces a wealth of information on the effect of solute structure on its solvation. The results of solvation of the present probes were employed in order to test the goodness of fit of two independent sets of solvent solvatochromic parameters.

First study on the thermo-solvatochromism in aqueous 1-(1-butyl)-3-methylimidazolium tetrafluoroborate: A comparison between the solvation by an ionic liquid and by aqueous alcohols

The thermo-solvatochrornic behaviors of 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate, RB; 2,6-dichloro-4-(2,4,6-triphenyloyridinium-1-yl) phenolate, WB; 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2); 2,6-dibromo-4-[(E)-2-(1-n-octylpyridinium-4-yl)ethenyl] phenolate, OcPMBr(2), have been investigated in binary mixtures of the ionic liquid, IL, 1-(1-butyl)-3-methylimidazolium tetrafluorborate, [BuMeIm][BF(4)], and water (W), in the temperature range from 10 to 60 degrees C. Plots of the empirical solvent polarities, ET (probe) in kcal mol(-1), versus the mole fraction of water in the binary mixture, chi(w) showed nonlinear, i.e., nonideal behavior. Solvation by these IL-W mixtures shows the following similarities to that by aqueous aliphatic alcohols: The same solvation model can be conveniently employed to treat the data obtained; it is based on the presence in the system-bulk medium and probe solvation shell of IL, W, and the ""complex"" solvent 1:1 IL-W. The origin of the nonideal solvation behavior appears to be the same, preferential solvation of the probe, in particular by the complex solvent. The strength of association of the IL-W complex, and the polarity of the IL are situated between the corresponding values of aqueous methanol and aqueous ethanol. Temperature increase causes a gradual desolvation of all probes employed. A difference between solvation by IL-W and aqueous alcohols is that probe-solvent hydrophobic interactions appear to play a minor role in case of the former mixture, probably because solvation is dominated by hydrogen-bonding and Coulombic interactions between the ions of the IL and the zwitterionic probes.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Use of genetic algorithms and solvation potential to study peptides structure

In this work, genetic algorithms concepts along with a rotamer library for proteins side chains and implicit solvation potential are used to optimize the tertiary structure of peptides. We starting from the known PDB structure of its backbone which is kept fixed while the side chains allowed adopting the conformations present in the rotamer library. It was used rotamer library independent of backbone and a implicit solvation potential. The structure of Mastoporan-X was predicted using several force fields with a growing complexity; we started it with a field where the only present interaction was Lennard-Jones. We added the Coulombian term and we considered the solvation effects through a term proportional to the solvent accessible area. This paper present good and interesting results obtained using the potential with solvation term and rotamer library. Hence, the algorithm (called YODA) presented here can be a good tool to the prediction problem. (c) 2007 Elsevier B.V. All rights reserved.

Predicting peptides structure with solvation potential and rotamer library dependent of the backbone

In this work, genetic algorithms concepts along with a rotamer library dependent of backbone and implicit solvation potential are used to study the tertiary structure of peptides. We starting from known primary sequence and optimize the structure of the backbone while the side chains allowed adopting the conformations present in a rotamer library. The GA, implemented with two force fields with a growing complexity, was used predict the structure of a polyalanine and a poly-isolueucine. This paper presents good and interesting results about the study of peptides structures and about the development of computational tools to study peptides structures. (C) 2007 Elsevier B.V. All rights reserved.

Correlation between solvation of peptide-resins and solvent properties

The solvation properties of model resin and peptide-resins measured in ca. 30 solvent systems correlated better with the sum of solvent electron acceptor (AN) and electron donor (DN) numbers, in 1:1 proportion, than with other solvent polarity parameters. The high sensitivity of the (AN+DN) term to detect differentiated solvation behaviors of peptide-resins, taken as model of heterogeneous and complex solutes, seems to be in agreement with the previously proposed two-parameter model, where the sum of the Lewis acidity and Lewis basicity characters of solvent are proposed for scaling solvent effect. Besides these physicochemical aspects regarding solute-solvent interactions, important implications of this study for the solid phase peptide synthesis were also observed. Each class of peptide-resin displayed a specific salvation profile that was dependent on the amount and the nature of the resin-bound peptide sequence. Plots of resin swelling versus solvent (AN+DN) values allowed the visualization of a maximum salvation region characteristic for each class of resin. This strategy facilitates the selection of solvent systems for optimal solvation conditions of peptide chains in every step of the entire synthesis cycle. Moreover, only the AN and DN concepts allow the understanding of rules for solvation/shrinking of peptide-resins when in homogeneous or in heterogeneous mixed solvents.

Correlation between the mobility of spin-labeled peptide chains and resin solvation: An approach to optimize the synthesis of aggregating sequences

Resin solvation properties affect the efficiency of the coupling reactions in solid-phase peptide synthesis. Here we report a novel approach to evaluate resin solvation properties, making use of spin label electron paramagnetic resonance (EPR) spectroscopy. The aggregating VVLGAAIV and ING sequences were assembled in benzhydrylamine-resin with different amino group contents (up to 2.6 mmol/g) to examine the extent of chain association within the beads. These model peptidyl-resins were first labeled at their N-terminus with the amino acid spin label 2,2,6,6-tetramethylpiperidine-N-oxyl-4-amino-4-carboxylic acid (Toac). Their solvation properties in different solvents were estimated, either by bead swelling measurement or by assessing the dynamics of their polymeric matrixes through the analysis of Toac EPR spectra, and were correlated with the yield of the acylation reaction. In most cases the coupling rate was found to depend on bead swelling. Comparatively, the EPR approach was more effective. Line shape analysis allowed the detection of more than one peptide chain population, which influenced the reaction. The results demonstrated the unique potential of EPR spectroscopy not only for improving the yield of peptide synthesis, even in challenging conditions, but also for other relevant polymer-supported methodologies in chemistry and biology.

EPR investigation of the influence of side chain protecting groups on peptide-resin solvation of the Asx and Glx model containing peptides

In spite of all progressive efforts aiming to optimize SPPS, serious problems mainly affecting the assembly of aggregating sequences have persisted. Following the study intended to unravel the complex solvation phenomenon of peptide-resin beads, the XING and XAAAA model aggregating segments were labeled with a paramagnetic probe and studied via EPR spectroscopy. Low and high substituted resins were also comparatively used, with the X residue being Asx or Glx containing the main protecting groups used in the SPPS. Notably, the cyclo-hexyl group used for Asp and Glu residues in Boc-chemistry induced greater chain immobilization than its tert-butyl partner-protecting group of the Fmoc strategy. Otherwise, the most impressive peptide chain immobilization occurred when the large trytil group was used for Asn and Gln protection in Fmoc-chemistry. These surprising results thus seem to stress the possibility of the relevant influence of the amino-acid side chain protecting groups in the overall peptide synthesis yield. (C) 2007 Elsevier Ltd. All rights reserved.

A hybrid neutron diffraction and computer simulation study on the solvation of N-methylformamide in dimethylsulfoxide

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Solvation in aqueous binary mixtures: consequences of the hydrophobic character of the ionic liquids and the solvatochromic probes

Information on the solvation in mixtures of water, W, and the ionic liquids, ILs, 1-allyl-3-R-imidazolium chlorides; R = methyl, 1-butyl, and 1-hexyl, has been obtained from the responses of the following solvatochromic probes: 2,6-dibromo-4-[(E)-2-(1-R-pyridinium-4-yl)ethenyl] phenolate, R = methyl, MePMBr2; 1-octyl, OcPMBr(2), and the corresponding quinolinium derivative, MeQMBr(2). A model developed for solvation in binary mixtures of W and molecular solvents has been extended to the present mixtures. Our objective is to assess the relevance to solvation of hydrogen-bonding and the hydrophobic character of the IL and the solvatochromic probe. Plots of the medium empirical polarity, E-T(probe) versus its composition revealed non-ideal behavior, attributed to preferential solvation by the IL and, more efficiently, by the IL-W hydrogen-bonded complex. The deviation from linearity increases as a function of increasing number of carbon atoms in the alkyl group of the IL, and is larger than that observed for solvation by W plus molecular solvents (1-propanol and 2-(1-butoxy)ethanol) that are more hydrophobic than the ILs investigated. This enhanced deviation is attributed to the more organized structure of the ILs proper, which persists in their aqueous solutions. MeQMBr(2) is more susceptible to solvent lipophilicity than OcPMBr(2), although the former probe is less lipophilic. This enhanced susceptibility agrees with the important effect of annelation on the contributions of the quinonoid and zwitterionic limiting structures to the ground and excited states of the probe, hence on its response to both medium composition and lipophilicity of the IL.

Chiroptical properties of bioactive molecules: sensitivity to conformation and solvation

Chiroptical spectroscopies play a fundamental role in pharmaceutical analysis for the stereochemical characterisation of bioactive molecules, due to the close relationship between chirality and optical activity and the increasing evidence of stereoselectivity in the pharmacological and toxicological profiles of chiral drugs. The correlation between chiroptical properties and absolute stereochemistry, however, requires the development of accurate and reliable theoretical models. The present thesis will report the application of theoretical chiroptical spectroscopies in the field of drug analysis, with particular emphasis on the huge influence of conformational flexibility and solvation on chiroptical properties and on the main computational strategies available to describe their effects by means of electronic circular dichroism (ECD) spectroscopy and time-dependent density functional theory (TD-DFT) calculations. The combination of experimental chiroptical spectroscopies with state-of-the-art computational methods proved to be very efficient at predicting the absolute configuration of a wide range of bioactive molecules (fluorinated 2-arylpropionic acids, β-lactam derivatives, difenoconazole, fenoterol, mycoleptones, austdiol). The results obtained for the investigated systems showed that great care must be taken in describing the molecular system in the most accurate fashion, since chiroptical properties are very sensitive to small electronic and conformational perturbations. In the future, the improvement of theoretical models and methods, such as ab initio molecular dynamics, will benefit pharmaceutical analysis in the investigation of non-trivial effects on the chiroptical properties of solvated systems and in the characterisation of the stereochemistry of complex chiral drugs.

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using Complete Basis Set and Gaussian-n Models combined with CPCM Continuum Solvation Methods

Complete basis set and Gaussian-n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pKa values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol, to calculate pKa values with cycle 1. The complete basis set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001

Accurate pKa Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pKa values for six carboxylic acids. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol to calculate pKa values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit.

Accurate relative pKa calculations for carboxylic acids using complete basis set and Gaussian-n models combined with continuum solvation methods

The complete basis set methods CBS-4, CBS-QB3, and CBS-APNO, and the Gaussian methods G2 and G3 were used to calculate the gas phase energy differences between six different carboxylic acids and their respective anions. Two different continuum methods, SM5.42R and CPCM, were used to calculate the free energy differences of solvation for the acids and their anions. Relative pKa values were calculated for each acid using one of the acids as a reference point. The CBS-QB3 and CBS-APNO gas phase calculations, combined with the CPCM/HF/6-31+G(d)//HF/6-31G(d) or CPCM/HF/6-31+G(d)//HF/6-31+G(d) continuum solvation calculations on the lowest energy gas phase conformer, and with the conformationally averaged values, give results accurate to ½ pKa unit. © 2001 American Institute of Physics.