A protein molecule in an aqueous mixed solvent: fluctuation theory outlook

In the present paper a procedure to calculate the properties of proteins in aqueous mixed solvents, particularly the excesses of the constituents of the mixed solvent near the protein molecule and the preferential binding parameters, is suggested. Expressions for the Kirkwood-Buff integrals in ternary mixtures and for the preferential binding parameter were derived and used to calculate various properties of infinitely dilute proteins in aqueous mixed solvents. The derived expressions and experimental information regarding the partial molar volumes and the preferential binding parameters were used to calculate the excesses (deficits) of water and cosolvent (in comparison with the bulk concentrations of protein-free mixed solvent) in the vicinity of ribonuclease A, ribonuclease T1, and lysozyme molecules. The calculations showed that water was in excess in the vicinity of ribonuclease A for water/glycerol and water/trehalose mixtures, and the cosolvent urea was in excess in the vicinity of ribonuclease T1 and lysozyme. The derivative of the activity coefficient of the protein with respect to the mole fraction of water was also calculated. This derivative was negative for the water/glycerol and water/trehalose mixed solvents and positive for the water/urea mixture. The mixture of lysozyme in the water/urea solvent is of particular interest, because the lysozyme at pH 7.0 is in its native state up to 9.3M urea, while at pH 2.0 it is denaturated between 2.5 and 5M and higher concentrations of urea. Our results demonstrated a striking similarity in the hydration of lysozyme at both pHs. It is worthwhile to note that the excesses of urea were only weakly composition dependent on both cases.     

New insights into urea action on proteins: a SANS study of the lysozyme case

We present a study on lysozyme dissolved in mixtures of water and urea, which is ubiquitously used as a protein denaturant. Despite the wide use of urea, the basic molecular mechanisms inducing protein unfolding are not still clarified. Small-angle neutron scattering (SANS) experiments have been performed using little amounts of denaturant in solutions in order to investigate the urea effect on lysozyme preceding the unfolding process. A global fit strategy, applied to analyze SANS experiments, provides an estimation of the average composition of the solvent in the close vicinity of the protein surface and the change of the protein-protein interactions due to the presence of urea. In particular, the thermodynamic equilibrium constant responsible for cosolvent balancing between the bulk and solvation layer has been determined. It turns out that urea is preferentially driven to the protein surface, confirming literature results at infinite dilute conditions. SANS data also reveal a possible variation of the protein net charge as a function of urea concentration, opening new perspectives and questions about the protein surface architecture at the first stages of unfolding processes.     

Preferential hydration of lysozyme in water/glycerol mixtures: a small-angle neutron scattering study

In solution small-angle neutron scattering has been used to study the solvation properties of lysozyme dissolved in water/glycerol mixtures. To detect the characteristics of the protein-solvent interface, 35 different experimental conditions (i.e., protein concentration, water/glycerol fraction in the solvent, content of deuterated compounds) have been considered and a suitable software has been developed to fit simultaneously the whole set of scattering data. The average composition of the solvent in the close vicinity of the protein surface at each experimental condition has been derived. In all the investigated conditions, glycerol resulted especially excluded from the protein surface, confirming that lysozyme is preferentially hydrated. By considering a thermodynamic hydration model based on an equilibrium exchange between water and glycerol from the solvation layer to the bulk, the preferential binding coefficient and the excess solvation number have been estimated. Results were compared with data previously derived for ribonuclease A in the same mixed solvent: even if the investigated solvent compositions were very different, the agreement between data is noticeable, suggesting that a unique mechanism presides over the preferential hydration process. Moreover, the curve describing the excess solvation number as a function of the solvent composition shows the occurrence of a region of maximal hydration, which probably accounts for the changes in protein stability detected in the presence of cosolvents.     

Microcalorimetric study of thermal unfolding of lysozyme in water/glycerol mixtures: an analysis by solvent exchange model

Folded protein stabilization or destabilization induced by cosolvent in mixed aqueous solutions has been studied by differential scanning microcalorimetry and related to difference in preferential solvation of native and denatured states. In particular, the thermal denaturation of a model system formed by lysozyme dissolved in water in the presence of the stabilizing cosolvent glycerol has been considered. Transition temperatures and enthalpies, heat capacity, and standard free energy changes have been determined when applying a two-state denaturation model to microcalorimetric data. Thermodynamic parameters show an unexpected, not linear, trend as a function of solvent composition; in particular, the lysozyme thermodynamic stability shows a maximum centered at water molar fraction of about 0.6. Using a thermodynamic hydration model based on the exchange equilibrium between glycerol and water molecules from the protein solvation layer to the bulk, the contribution of protein-solvent interactions to the unfolding free energy and the changes of this contribution with solvent composition have been derived. The preferential solvation data indicate that lysozyme unfolding involves an increase in the solvation surface, with a small reduction of the protein-preferential hydration. Moreover, the derived changes in the excess solvation numbers at denaturation show that only few solvent molecules are responsible for the variation of lysozyme stability in relation to the solvent composition.     

Various contributions to the osmotic second virial coefficient in protein-water-cosolvent solutions

An analysis of the cosolvent concentration dependence of the osmotic second virial coefficient (OSVC) in water-protein-cosolvent mixtures is developed. The Kirkwood-Buff fluctuation theory for ternary mixtures is used as the main theoretical tool. On its basis, the OSVC is expressed in terms of the thermodynamic properties of infinitely dilute (with respect to the protein) water-protein-cosolvent mixtures. These properties can be divided into two groups: (1) those of infinitely dilute protein solutions (such as the partial molar volume of a protein at infinite dilution and the derivatives of the protein activity coefficient with respect to the protein and water molar fractions) and (2) those of the protein-free water-cosolvent mixture (such as its concentrations, the isothermal compressibility, the partial molar volumes, and the derivative of the water activity coefficient with respect to the water molar fraction). Expressions are derived for the OSVC of ideal mixtures and for a mixture in which only the binary mixed solvent is ideal. The latter expression contains three contributions: (1) one due to the protein-solvent interactions B2(p-s), which is connected to the preferential binding parameter, (2) another one due to protein/protein interactions (B2(p-p)), and (3) a third one representing an ideal mixture contribution (B2(id)). The cosolvent composition dependencies of these three contributions were examined for several water-protein-cosolvent mixtures using experimental data regarding the OSVC and the preferential binding parameter. For the water-lysozyme-arginine mixture, it was found that OSVC exhibits the behavior of an ideal mixture and that B2(id) provides the main contribution to the OSVC. For the other mixtures considered (water-Hm MalDH-NaCl, water-Hm MalDH-(NH4)2SO4, and water-lysozyme-NaCl mixtures), it was found that the contribution of the protein-solvent interactions B2(p-s) is responsible for the composition dependence of the OSVC on the cosolvent concentration, whereas the two remaining contributions (B2(p-p)) and B2(id)) are almost composition independent.     

UV-visible spectroscopic study of solvation in ternary solvent mixtures: ketocyanine dye in methanol + acetone + water and methanol + acetone + benzene

Solvation characteristics of a ketocyanine dye have been studied in completely miscible ternary solvent mixtures, namely, methanol + acetone + water and methanol + acetone + benzene, by monitoring the solvatochromic absorption band of the dye. The maximum energy of absorption (E) of the solute in a ternary solvent mixture differs significantly from the mole fraction average of the E values in the component solvents. Results in the corresponding binary solvent mixtures also show a deviation of the E value from the mole fraction averaged E values. The results have been explained in terms of preferential solvation using a two phase model of solvation. The excess or deficit over the bulk composition of a solvent component in the vicinity of the solute molecule in a ternary solvent mixture has been estimated using the knowledge of solvation in the corresponding binary mixtures.     

Estimation of excess solvation numbers of water and cosolvents from preferential interaction and volumetric experiments

To elucidate, at a molecular level, how cosolvents influence protein stability, it is indispensable to understand the distribution of water and cosolvent molecules around proteins. Calculation of excess solvation numbers of water and cosolvents serves this purpose, and I show that they can be extracted from preferential interaction parameter and volumeric data via the Kirkwood-Buff theory. This scheme was applied to trehalose and glycerol (stabilizers) and urea (denaturant). Important insights from the application include stabilizer-induced enhancement of protein hydration, which, together with the stabilizer's exclusion from protein surfaces, may contribute to protein stabilization at high osmolyte concentrations.     

Influence of the structure of the hydrophobe on the ideality of mixing in micelles in binary mixtures of amphiphilic tricyclic drugs

The composition of the micelles in binary mixtures of the cationic amphiphilic antidepressant drugs nortriptyline, amitriptyline, and doxepin has been determined from an analysis of the variation of the critical micelle concentration from conductivity measurements, as a function of solution composition. Assessment of the nonideality of mixing in terms of the interaction parameter from the regular solution approximation showed small deviations from ideal mixing, with negative interaction parameters for nortriptyline/amitriptyline systems and positive interaction parameters for mixtures of nortriptyline and doxepin. These differences in nonideality have been attributed to differences in the packing of the drugs in the mixed micelles arising from differences in the structure of the hydrophobe.     

Molecular interactions of cationic and anionic surfactants in mixed monolayers and aggregates

The properties of anionic-rich and cationic-rich mixtures of CTAB (cetyltrimethylammonium bromide) and SDS (sodium dodecyl sulfate) were investigated with conductometry and surface tension measurements and by determining the surfactant NMR self-diffusion coefficients. The critical aggregate concentration (CAC), surface tension reduction effectiveness(gamma(CAC)), surface excess(Gamma(max)), and mean molecular surface area (A(min)) were determined from plots of the surface tension (gamma) as a function of the total surfactant concentration. The compositions of the adsorbed films (Z) and aggregates (chi) were estimated by using regular solution theory, and then the interaction parameters in the aggregates (beta) and the adsorbed film phases (beta(sigma)) were calculated. The results showed that the synergism between the surfactants enhances the formation of mixed aggregates and reduces the surface tension. Further, the nature and strength of the interaction between the surfactants in the mixtures were obtained by calculating the values of the following parameters: the interaction parameter, beta, the size parameter, rho, and the nonrandom mixing parameter, P*. These results indicate that in ionic surfactant mixtures the optimized packing parameter has the highest value and that the size parameter can be used to account for deviations from the predictions of regular solution theory. It was concluded that, for planar air/aqueous interfaces and aggregation systems, this nonideality increases as the temperature increases. This trend is attributed to the increased dehydration of the surfactant head groups that results from increases in temperature. Further, our conductometry measurements show that the counterion binding number of mixed micelles formed in mixtures with a high CTAB content is different to those with a high SDS content. This difference is due to either their different aggregation sizes or the different interactions between the head groups and the counterions.     

Solubility of phenobarbital in aqueous cosolvent mixtures revisited: IKBI preferential solvation analysis

The preferential solvation parameters of phenobarbital in aqueous binary mixtures of 1,4-dioxane, t-butanol, n-propanol, ethanol, propylene glycol and glycerol were derived from solution thermodynamic properties by using the IKBI method. This drug is sensitive to preferential solvation effects in all these mixtures. The preferential solvation parameter by the cosolvent (δx_1,3) is negative in almost all the water-rich mixtures but positive in mixtures with similar proportions of solvents and cosolvent-rich mixtures, except in 1-propanol + water mixtures, where negative values are also found in mixtures with x₁ ≥ 0.70. Hydrophobic hydration around the non-polar ethyl and phenyl groups of this drug in water-rich mixtures could play a relevant role in drug solvation. Otherwise, in mixtures of similar solvent compositions and in cosolvent-rich mixtures the preferential solvation by cosolvent could be due to the acidic behaviour of the drug.     

Preferential solvation of some antiepileptic drugs in {cosolvent (1) + water (2)} mixtures at 298.15 K

The solubility of lamotrigine (LTG), clonazepam (CZP) and diazepam (DZP) in some {cosolvent (1) + water (2)} mixtures expressed in mole fraction at 298.15 K was calculated from reported solubility values expressed in molarity by using the densities of the saturated solutions. Aqueous binary mixtures of ethanol, propylene glycol and N-methyl-2-pyrrolidone were considered. From mole fraction solubilities and some thermodynamic properties of the solvent mixtures, the preferential solvation of these drugs by both solvents in the mixtures was analysed by using the inverse Kirkwood–Buff integrals. It is observed that LTG, CZP and DZP are preferentially solvated by water in water-rich mixtures in all the three binary systems analysed. In {ethanol (1) + water (2)} mixtures, preferential solvation by water is also observed in ethanol-rich mixtures. Nevertheless, in {propylene glycol (1) + water (2)} and {N-methyl-2-pyrrolidone (1) + water (2)} mixtures preferential solvation by the cosolvent was observed in cosolvent-rich mixtures.     

Interfacial and micellar properties of binary mixtures of surfactant and phospholipid in an aqueous medium

Mixed micelles formed by zwitterionic surfactant dimethyldodecylammniopropane sulfonate and short-chain phospholipid 1,2-diheptanoyl-sn-glycero-3-phosphocholine in different proportions in an aqueous medium have been studied physicochemically at an air/water interface and in the bulk by using interfacial tension and pyrene fluorescence intensity measurements, respectively. The critical micellar concentration and free energies of micellization and of interfacial adsorption have been determined. The interfacial study reveals that a mixed monolayer is formed at the air/water interface by the adsorption of surfactant and phospholipid monomers. This has been confirmed by evaluating the interfacial parameters; the maximum surface excess, the minimum area per molecule of a surface-active compound, and the Gibbs surface excess related to surface pressure. The nonideality of mixing, expressed in the terms of the regular solution interaction parameter, #, has negative values over the whole mole fraction range. The negative # values indicate the mutual synergism between the surfactant and phospholipid monomers. The equilibrium distribution of components between micelle and monomer phases was evaluated using a theoretical treatment based on excess thermodynamics quantities evaluated by Motomura's formulation.     

Molecular solvation in water-methanol and water-sorbitol mixtures: the roles of preferential hydration, hydrophobicity, and the equation of state

Molecular dynamics simulations of aqueous mixtures of methanol and sorbitol were performed over a wide range of binary composition, density (pressure), and temperature to study the equation of state and solvation of small apolar solutes. Experimentally, methanol is a canonical solubilizing agent for apolar solutes and a protein denaturant in mixed-aqueous solvents; sorbitol represents a canonical "salting-out" or protein-stabilizing cosolvent. The results reported here show increasing sorbitol concentration under isothermal, isobaric conditions results in monotonic increases in apolar solute excess chemical potential (mu2ex) over the range of experimentally relevant temperatures. For methanol at elevated temperatures, increasing cosolvent composition results in monotonically decreasing mu2ex. However, at lower temperatures mu2ex exhibits a maximum versus cosolvent concentration, as seen experimentally for Ar in ethanol-water solutions. Both density anomalies and hydrophobic effects--characterized by temperatures of density maxima and apolar solute solubility minima, respectively--are suppressed upon addition of either sorbitol or methanol at all temperatures and compositions simulated here. Thus, the contrasting effects of sorbitol and methanol on solute chemical potential cannot be explained by qualitative differences in their ability to enhance or suppress hydrophobic effects. Rather, we find mu2ex values across a broad range of temperatures and cosolvent composition can be quantitatively explained in terms of isobaric changes in solvent density--i.e., the equation of state--along with the corresponding packing fraction of the solvent. Analysis in terms of truncated preferential interaction parameters highlights that care must be taken in interpreting cosolvent effects on solvation in terms of local preferential hydration.     

Water rotational relaxation and diffusion in hydrated lysozyme

This paper is concerned with the dynamics of water around a small globular protein. Dipolar second-rank relaxation time and diffusion properties of surface water were computed by extensive molecular dynamics simulations of lysozyme in water which lasted a total of 28 ns. Our results indicate that the rotational relaxation of water in the vicinity of lysozyme is 3-7 times slower than that in the bulk depending on how the hydration shell is defined in the calculation. We have also verified that the dynamics of water translational diffusion in the vicinity of lysozyme have retardations similar to rotational relaxation. This is a common assumption in nuclear magnetic relaxation dispersion (NMRD) studies to derive residence times. In contrast to bulk water dynamics, surface water is in a dispersive diffusion regime or subdiffusion. Very good agreement of dipolar second-rank relaxation time with NMRD estimates is obtained by using appropriate dimensions of the hydration shell. Although our computed second-rank dipolar retardations are independent of the water model, SPC/E describes more realistically the time scale of the water dynamics around lysozyme than does TIP3P.     

New description of heterogeneity effects in adsorption from multicomponent solutions and liquid adsorption chromatography

General equations describing adsorption from solutions on solids and liquid adsorption chromatography with mixed mobile phases are formulated in terms of the bulk and surface activity coefficeints. Definition of the surface activity coefficients is extended; they describe nonideality of surface solution due to difference in molecular interactions of the components as well as nonideality of this solution generated by the adsorbent heterogeneity. It is shown that the above general equations predict simpler expressions known already in the literature.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Thermochemistry of electrolyte solutions

The results of calorimetric investigations of electrolyte solutions in the mixtures of water, methanol, N,N-dimethylformamide, and acetonitrile with numerous organic cosolvents are discussed with regard to the intermolecular interactions that occur in the solution. Particular attention is given to answer the questions how and to what extent the properties of the systems examined are modified by the cosolvent added and how much the properties of the cosolvent are revealed in the mixtures with the solvents mentioned above. To this goal, the analysis of the electrolyte dissolution enthalpies, single ionic transfer enthalpies, and enthalpic pair interaction coefficients as well as the preferential solvation (PS) model are applied. The analysis performed shows that in the case of the dissolution enthalpies of simple inorganic electrolytes in water–organic solvent mixtures, the shape of the dependence of the standard dissolution enthalpy on the mixed solvent composition reflects to a large extent the hydrophobic properties of the organic cosolvent. In the mixtures of methanol with organic cosolvents, the ions are preferentially solvated either by methanol molecules or by molecules of the cosolvent, depending on the properties of the mixed solvent components. The behavior of inorganic salts in the mixtures containing N,N-dimethylformamide is mostly influenced by the DMF which is a relatively strongly ion solvating solvent, whereas in acetonitrile mixtures, the thermochemical behavior of electrolyte solutions is influenced to a large extent by the properties of the cosolvent particularly due to the PS of cation by the cosolvent molecules.     

Revealing the influence of the strength of coulomb interactions on the viscosity and interfacial tension of ionic liquid cosolvent mixtures

Within this work, viscosity and interfacial tension of selected ionic liquid cosolvent mixtures, [EMIM][EtOSO3] (1-ethyl-3-methyl-immidazolium ethyl sulfate) with water and ethanol, were studied as a function of composition by surface light scattering (SLS) and the pendant drop method in a consistent manner, allowing a close insight into the nature of interactions. Here, we show that the viscosity behavior clearly reflects the bulk structure of the ionic liquid cosolvent mixtures and correlates to the fluid structure at the phase boundary. In contrast to former work, we found the viscosity of ionic liquid [EMIM][EtOSO3] to be decreasing the stronger by small amounts of the cosolvents and the lower their dielectric constant. Furthermore, two distinct trends for the dependence of the viscosity on the cosolvent concentration were resolved. These were assigned to ion-dipole interactions dominating in the salt-rich region and to dipole-dipole interactions in the diluted one. A crossover between both regions is reflected by the interfacial tension data, where it seems that up to a "critical" concentration almost no cosolvent is present at the phase boundary.     

The Kirkwood-Buff theory of solutions and the local composition of liquid mixtures

The present paper is devoted to the local composition of liquid mixtures calculated in the framework of the Kirkwood-Buff theory of solutions. A new method is suggested to calculate the excess (or deficit) number of various molecules around a selected (central) molecule in binary and multicomponent liquid mixtures in terms of measurable macroscopic thermodynamic quantities, such as the derivatives of the chemical potentials with respect to concentrations, the isothermal compressibility, and the partial molar volumes. This method accounts for an inaccessible volume due to the presence of a central molecule and is applied to binary and ternary mixtures. For the ideal binary mixture it is shown that because of the difference in the volumes of the pure components there is an excess (or deficit) number of different molecules around a central molecule. The excess (or deficit) becomes zero when the components of the ideal binary mixture have the same volume. The new method is also applied to methanol + water and 2-propanol + water mixtures. In the case of the 2-propanol + water mixture, the new method, in contrast to the other ones, indicates that clusters dominated by 2-propanol disappear at high alcohol mole fractions, in agreement with experimental observations. Finally, it is shown that the application of the new procedure to the ternary mixture water/protein/cosolvent at infinite dilution of the protein led to almost the same results as the methods involving a reference state.