Vapor pressure osmometry studies of osmolyte-protein interactions: implications for the action of osmoprotectants in vivo and for the interpretation of "osmotic stress" experiments in vitro

Structure and Hydration of l-Proline in Aqueous Solutions

The structure and hydration of L-proline in aqueous solution have been investigated using a combination of neutron diffraction with isotopic substitution, empirical potential structure refinement modeling, and small-angle neutron scattering at three concentrations, 1:10, 1:15, and 1:20 proline/water mole ratios. In each solution the carboxylate oxygen atoms from proline accept less than two hydrogen bonds from the surrounding water solvent and the amine hydrogen atoms donate less than one hydrogen bond to the surrounding water molecules. The solute-solute radial distribution functions indicate relatively weak interactions between proline molecules, and significant clustering or aggregation of proline is absent at all these concentrations. The spatial density distributions for the hydration of the COO- group in proline show a similar shape to that found previously in L-glutamic acid in aqueous solution but with a reduced coordination number.

Local Composition in the Vicinity of a Protein Molecule in an Aqueous Mixed Solvent

This paper is focused on the composition of a cosolvent in the vicinity of a protein surface (local composition) and its dependence on various factors. First, the Kirkwood-Buff theory of solution is used to obtain analytical expressions that connect the excess or deficit number of cosolvent and water molecules in the vicinity of a protein surface with experimentally measurable quantities such as the bulk concentration of the mixed solvent, the preferential binding parameter, and the molar volumes of water and cosolvent. Using these expressions, relations between the preferential binding parameter (at a molal concentration scale) and the above excesses (or deficits) are established. In addition, the obtained expressions are used to examine the effect of the nonideality of the water + cosolvent mixtures and of the molar volume of the cosolvent on the excess (or deficit) number of cosolvent molecules in the vicinity of the protein surface. It is shown that at least for the mixed solvents considered (water + urea and water + glucose) the nonideality of the mixed solvent is not an important factor in the local compositions around a protein molecule and that the main contribution is provided by the nonidealities of the protein-water and protein-cosolvent mixtures. Special attention is paid to urea as cosolvent, because urea is one of only a few compounds with a concentration at the protein surface larger than its concentration in the bulk. The composition dependence of the excess of urea around a protein molecule is calculated for the water + lysozyme + urea mixture at pH = 7.0 and 2.0. At pH = 7.0, the excess of urea becomes almost composition independent at high urea concentrations. Such independence could be explained by assuming that urea totally replaces water in some areas of the protein surface, whereas on the remaining areas of the protein surface both water and urea are present with concentration comparable to those in the bulk. The Schellman exchange model was used to relate the preferential binding parameter in water + lysozyme + urea mixtures to the urea concentration.

Break in the heat capacity change at 303 K for complex binding of netropsin to AATT containing hairpin DNA constructs

Studies performed in our laboratory demonstrated the formation of two thermodynamically distinct complexes on binding of netropsin to a number of hairpin-forming DNA sequences containing AATT-binding regions. These two complexes were proposed to differ only by a bridging water molecule between the drug and the DNA in the lower affinity complex. A temperature-dependent isothermal titration calorimetry (ITC)-binding study was performed using one of these constructs (a 20-mer hairpin of sequence 5'-CGAATTCGTCTCCGAATTCG) and netropsin. This study demonstrated a break in the heat capacity change for the formation of the complex containing the bridging water molecule at approximately 303 K. In the plot of the binding enthalpy change versus temperature, the slope (DeltaCp) was -0.67 kcal mol-1 K-1 steeper after the break at 303 K. Because of the relatively low melting temperature of the 20-mer hairpin (341 K (68 degrees C)), the enthalpy change for complex formation might have included some energy of refolding of the partially denatured hairpin, giving the suggestion of a larger DeltaCp. Studies done on the binding of netropsin to similar constructs, a 24-mer and a 28-mer, with added GC basepairs in the hairpin stem to increase thermal stability, exhibit the same nonlinearity in DeltaCp over the temperature range of from 275 to 333 K. The slopes (DeltaCp) were -0.69 and -0.64 kcal mol-1 K-1 steeper after 303 K for the 24-mer and 28-mer, respectively. This observation strengthens the argument regarding the presence of a bridging water molecule in the lower affinity netropsin/DNA complex. The DeltaCp data seem to infer that because the break in the heat capacity change function for the lower affinity binding occurs at the isoequilibrium temperature for water, water may be included or trapped in the complex. The fact that this break does not occur in the heat capacity change function for formation of the higher affinity complex can similarly be taken as evidence that water is not included in the higher affinity complex.

Effects of Osmolytes on Protein Folding and Aggregation in Cells

Nature has developed many strategies to ensure that the complex and challenging protein folding reaction occurs in vivo with adequate efficiency and fidelity for the success of the organism. Among the strategies widely employed in a huge range of species and cell types is the elaboration of small organic molecules called osmolytes that offset the potentially damaging effects of osmotic stress. While considerable knowledge has been gained in vitro regarding the influence of osmolytes on protein structure and folding, it is of great interest to probe the effects of osmolytes in cells. We have developed an in-cell fluorescent-labeling method that enables the study of protein stability and also protein aggregation in vivo. We utilize a genetically encoded tag called a tetra-Cys motif that binds specifically to a bis-arsenical fluorescein-based dye "FlAsH"; we inserted the tetra-Cys motif into a protein of interest in such a way that the FlAsH signal reported on the state of folding or aggregation of the protein. Then, we designed protocols to assess how various osmolytes influence the stability and propensity to aggregate of our protein of interest. These are described here. Not only are there potential biotechnological applications of osmolytes in the quest to produce greater quantities of well-folded proteins, but also osmolytes may serve as tools and points of departure for therapeutic intervention in protein folding and aggregation diseases. Having in vivo methods to analyze how osmolytes affect folding and aggregation enhances our ability to further these goals greatly.

Application of the Transfer Model to Understand How Naturally Occurring Osmolytes Affect Protein Stability

A primary thermodynamic goal in protein biochemistry is to attain a predictive understanding of the energetic changes responsible for solvent-induced folding and unfolding. This chapter demonstrates the use of Tanford's transfer model to predict solvent-dependent cooperative protein folding/unfolding free energy changes (m values). This approach provides a thermodynamic description of these free energy changes in terms of individual contributions from the peptide backbone and residue side chains. The quantitative success of the transfer model has been hindered for many years because of unresolved issues involving proper measurement of the group transfer-free energies of amino acid side chains and the peptide backbone unit. This chapter demonstrates what is necessary to design experiments properly so that reliable values of group transfer-free energies are obtainable. It then demonstrates how to derive a prediction of the m value for the description of protein folding/unfolding cooperativity and that the calculated values using the transfer model agree quite well with experimentally measured values.

The temperature dependence of salt-protein association is sequence specific

Molecular dynamics (MD) simulations are used to probe the origin of the unexpected temperature dependence of salt accumulation in the C-terminal region of the protein human lymphotactin. As in previous MD simulations, sodium ions accumulate in an enhanced manner near the C-terminal helix at the lower temperature, while the temperature dependence of chloride accumulation is much weaker and slightly positive. In a designed mutant in which all positively charged residues in the C-terminal helix are replaced with neutral polar groups (Ser), the unexpected temperature dependence of the sodium ions is no longer observed. Therefore, these simulations convincingly verified the previous hypothesis that the temperature dependence of ion-protein association is sensitive to the local sequence. This is explained qualitatively in terms of the entropy of association between charged species in solution. These findings have general implications for the interpretation of thermodynamic quantities associated with binding events where ion release is important, such as protein-DNA interactions.

Hydration changes accompanying the binding of minor groove ligands with DNA

4',6-diamidino-2-phenylindole (DAPI), netropsin, and pentamidine are minor groove binders that have terminal -C(NH2)2+ groups. The hydration changes that accompany their binding to the minor groove of the (AATT)2 sequence have been studied using the osmotic stress technique with fluorescence spectroscopy. The affinity of DAPI for the binding site decreases with the increasing osmolality of the solution, resulting in acquisition of 35+/-1 waters upon binding. A competition fluorescence assay was utilized to measure the binding constants and hydration changes of the other two ligands, using the DNA-DAPI complex as the fluorescence reporter. Upon their association to the (AATT)2 binding site, netropsin and pentamidine acquire 26+/-3 and 34+/-2 additional waters of hydration, respectively. The hydration changes are discussed in the context of the terminal functional groups of the ligands and conformational changes in the DNA.

Effect of salts and organic additives on the solubility of proteins in aqueous solutions

The goal of this review is to examine the effect of salts and organic additives on the solubility of proteins in aqueous mixed solvents. The focus is on the correlation between the aqueous protein solubility and the osmotic second virial coefficient or the preferential binding parameter. First, several approaches which connect the solubility and the osmotic second virial coefficient are presented. Most of the experimental and theoretical results correlate the solubility and the osmotic second virial coefficient in the presence of salts. The correlation of the aqueous protein solubility with the osmotic second virial coefficient when the cosolvent is an organic component requires additional research. Second, the aqueous protein solubility is correlated with the preferential binding parameter on the basis of a theory developed by the authors of the present review. This theory can predict (i) the salting-in or -out effect of a cosolvent and (ii) the initial slope of the solubility curve. Good agreement was obtained between theoretical predictions and experimental results.

Differences in hydration coupled to specific and nonspecific competitive binding and to specific DNA Binding of the restriction endonuclease BamHI

Using the osmotic stress technique together with a self-cleavage assay we measure directly differences in sequestered water between specific and nonspecific DNA-BamHI complexes as well as the numbers of water molecules released coupled to specific complex formation. The difference between specific and nonspecific binding free energy of the BamHI scales linearly with solute osmolal concentration for seven neutral solutes used to set water activity. The observed osmotic dependence indicates that the nonspecific DNA-BamHI complex sequesters some 120-150 more water molecules than the specific complex. The weak sensitivity of the difference in number of waters to the solute identity suggests that these waters are sterically inaccessible to solutes. This result is in close agreement with differences in the structures determined by x-ray crystallography. We demonstrate additionally that when the same solutes that were used in competition experiments are used to probe changes accompanying the binding of free BamHI to its specific DNA sequence, the measured number of water molecules released in the binding process is strikingly solute-dependent (with up to 10-fold difference between solutes). This result is expected for reactions resulting in a large change in a surface exposed area.

Inhibition of protein aggregation in vitro and in vivo by a natural osmoprotectant

Small organic molecules termed osmolytes are harnessed by a variety of cell types in a wide range of organisms to counter unfavorable physiological conditions that challenge protein stability and function. Using a well characterized reporter system that we developed to allow in vivo observations, we have explored how the osmolyte proline influences the stability and aggregation of a model aggregation-prone protein, P39A cellular retinoic acid-binding protein. Strikingly, we find that the natural osmolyte proline abrogates aggregation both in vitro and in vivo (in an Escherichia coli expression system). Importantly, proline also prevented aggregation of constructs containing exon 1 of huntingtin with extended polyglutamine tracts. Although compatible osmolytes are known to stabilize the native state, our results point to a destabilizing effect of proline on partially folded states and early aggregates and a solubilizing effect on the native state. Because proline is believed to act through a combination of solvophobic backbone interactions and favorable side-chain interactions that are not specific to a particular sequence or structure, the observed effect is likely to be general. Thus, the osmolyte proline may be protective against biomedically important protein aggregates that are hallmarks of several late-onset neurodegenerative diseases including Huntington's, Alzheimer's, and Parkinson's. In addition, these results should be of practical importance because they may enable protein expression at higher efficiency under conditions where aggregation competes with proper folding.

High concentrations of the compatible solute glycinebetaine destabilize model membranes under stress conditions

Compatible solutes are accumulated by diverse organisms in response to environmental stresses such as drought, salt, or cold. Glycinebetaine (Bet) is such a solute that is accumulated by many plants and microorganisms to high concentrations under stress conditions. It is an osmoprotectant in bacteria and stabilizes both soluble and peripherally membrane-bound proteins in vitro. Here, the effects of Bet on the stability of model lipid membranes are compared to the effects of two other compatible solutes, sucrose and trehalose. Both in the presence of 1M NaCl and during freezing to -20 degrees C, Bet is highly destabilizing to liposomes containing nonbilayer lipids, while the disaccharides are either protective or, in some cases, much less destabilizing. The destabilizing effect of Bet is more pronounced in membranes containing the nonbilayer galactolipid monogalactosyldiacylglycerol from plant chloroplasts than in membranes containing the nonbilayer phospholipid phosphatidylethanolamine. The most dramatic differences between the sugars and Bet were observed in liposomes made from a combination of lipids resembling plant chloroplast thylakoid membranes. Measurements with the dye merocyanine 540 indicate that the water-membrane interface was affected in opposite directions by the presence of high concentrations of sucrose or Bet. The dynamics of the lipids, however, were not differentially affected by the solutes, making direct solute-lipid interactions an unlikely explanation for the different effects on stability. The data offer an explanation, why Bet at high concentrations achieved during exogenous feeding of leaf tissues can be detrimental to cellular stability and survival under stress, while bacterial membranes that contain phosphatidylethanolamine instead of monogalactosyldiacylglycerol, or cyanobacteria that contain highly saturated monogalactosyldiacylglycerol are less susceptible.

Effects of temperature and salt concentration on the structural stability of human lymphotactin: insights from molecular simulations

Extensive molecular dynamics (MD) simulations ( approximately 70 ns total) with explicit solvent molecules and salt ions are carried out to probe the effects of temperature and salt concentration on the structural stability of the human Lymphotactin (hLtn). The distribution of ions near the protein surface and the stability of various structural motifs are observed to exhibit interesting dependence on the local sequence and structure. Whereas chloride association to the protein is overall enhanced as the temperature increases, the sodium distribution in the C-terminal helical region and, to a smaller degree, the chloride distribution in the same region are found higher at the lower temperature. The similar trend is also observed in nonlinear Poisson-Boltzmann calculations with a temperature-dependent water dielectric constant, once conformational averaging over a series of MD snapshots is done. The unexpected temperature dependence in the ion distribution is explained on the basis of the cancellation of association entropy for ion-side chain pairs of opposite-charge and like-charge characters, which have positive and negative contributions, respectively. The C-terminal helix is observed to partially melt whereas a short beta strand forms at the higher temperature with little salt dependence. The N-terminal region, by contrast, develops partial helical structure at a higher salt concentration. These observed behaviors are consistent with solvent and salt screening playing an important role in stabilizing the canonical chemokine fold of hLtn.

The Effect of Salt Stoichiometry on Protein−Salt Interactions Determined by Ternary Diffusion in Aqueous Solutions

We report the four diffusion coefficients for the lysozyme-MgCl2-water ternary system at 25 degrees C and pH 4.5. The comparison with previous results for the lysozyme-NaCl-water ternary system is used to examine the effect of salt stoichiometry on the transport properties of lysozyme-salt aqueous mixtures. We find that the two cross-diffusion coefficients are very sensitive to salt stoichiometry. One of the cross-diffusion coefficients is examined in terms of common-ion, excluded-volume, and protein-preferential hydration effects. We use the four ternary diffusion coefficients to extract chemical-potential cross-derivatives and protein-preferential interaction coefficients. These thermodynamic data characterize the protein-salt thermodynamic interactions. We demonstrate the presence of the common-ion effect (Donnan effect) by analyzing the dependence of the preferential-interaction coefficient not only with respect to salt concentration but also with respect to salt stoichiometry. We conclude that the common-ion effect and the protein-preferential hydration are both important for describing the lysozyme-MgCl2 thermodynamic interaction.

Sequestered water and binding energy are coupled in complexes of lambda Cro repressor with non-consensus binding sequences

We use the osmotic pressure dependence of dissociation rates and relative binding constants to infer differences in sequestered water among complexes of lambda Cro repressor with varied DNA recognition sequences. For over a 1000-fold change in association constant, the number of water molecules sequestered by non-cognate complexes varies linearly with binding free energy. One extra bound water molecule is coupled with the loss of approximately 150 cal/mol complex in binding free energy. Equivalently, every tenfold decrease in binding constant at constant salt and temperature is associated with eight to nine additional water molecules sequestered in the non-cognate complex. The relative insensitivity of the difference in water molecules to the nature of the osmolyte used to probe the reaction suggests that the water is sterically sequestered. If the previously measured changes in heat capacity for lambda Cro binding to different non-cognate sequences are attributed solely to this change in water, then the heat capacity change per incorporated water is almost the same as the difference between ice and water. The associated changes in enthalpies and entropies, however, indicate that the change in complex structure involves more than a simple incorporation of fixed water molecules that act as adaptors between non-complementary surfaces.

Binding of netropsin to several DNA constructs: evidence for at least two different 1:1 complexes formed from an -AATT-containing ds-DNA construct and a single minor groove binding ligand

Isothermal titration calorimetry, ITC, has been used to determine the thermodynamics (DeltaG, DeltaH, and -TDeltaS) for binding netropsin to a number of DNA constructs. The DNA constructs included: six different 20-22mer hairpin forming sequences and an 8-mer DNA forming a duplex dimer. All DNA constructs had a single -AT-rich netropsin binding with one of the following sequences, (A(2)T(2))(2), (ATAT)(2), or (AAAA/TTTT). Binding energetics are less dependent on site sequence than on changes in the neighboring single stranded DNA (hairpin loop size and tail length). All of the 1:1 complexes exhibit an enthalpy change that is dependent on the fractional saturation of the binding site. Later binding ligands interact with a significantly more favorable enthalpy change (partial differential DeltaH(1-2) from 2 to 6 kcal/mol) and a significantly less favorable entropy change (partial differential (-TDeltaS(1-2))) from -4 to -9 kcal/mol). The ITC data could only be fit within expected experimental error by use of a thermodynamic model that includes two independent binding processes with a combined stoichiometry of 1 mol of ligand per 1 mol of oligonucleotide. Based on the biophysical evidence reported here, including theoretical calculations for the energetics of "trapping" or structuring of a single water molecule and molecular docking computations, it is proposed that there are two modes by which flexible ligands can bind in the minor groove of duplex DNA. The higher affinity binding mode is for netropsin to lay along the floor of the minor groove in a bent conformation and exclude all water from the groove. The slightly weaker binding mode is for the netropsin molecule to have a slightly more linear conformation and for the required curvature to be the result of a water molecule that bridges between the floor of the minor groove and two of the amidino nitrogens located at one end of the bound netropsin molecule.

Use of urea and glycine betaine to quantify coupled folding and probe the burial of DNA phosphates in lac repressor-lac operator binding

Thermodynamic analysis of urea-biopolymer interactions and effects of urea on folding of proteins and alpha-helical peptides shows that urea interacts primarily with polar amide surface. Urea is therefore predicted to be a quantitative probe of coupled folding, remodeling, and other large-scale changes in the amount of water-accessible polar amide surface in protein processes. A parallel analysis indicates that glycine betaine [N,N,N-trimethylglycine (GB)] can be used to detect burial or exposure of anionic (carboxylate, phosphate) biopolymer surface. To test these predictions, we have investigated the effects of these solutes (0-3 m) on the formation of 1:1 complexes between lac repressor (LacI) and its symmetric operator site (SymL) at a constant KCl molality. Urea reduces the binding constant K(TO) [initial slope dlnK(TO)/dm(urea) = -1.7 +/- 0.2], and GB increases K(TO) [initial slope dlnK(TO)/dm(GB) = 2.1 +/- 0.2]. For both solutes, this derivative decreases with an increase in solute concentration. Analysis of these initial slopes predicts that (1.5 +/- 0.3) x 10(3) A2 of polar amide surface and (4.5 +/- 1.0) x 10(2) A2 of anionic surface are buried in the association process. Analysis of published structural data, together with modeling of unfolded regions of free LacI as extended chains, indicates that 1.5 x 10(3) A2 of polar amide surface and 6.3 x 10(2) A2 of anionic surface are buried in complexation. Quantitative agreement between structural and thermodynamic results is obtained for amide surface (urea); for anionic surface (GB), the experimental value is approximately 70% of the structural value. For LacI-SymL binding, two-thirds of the structurally predicted change in amide surface (1.0 x 10(3) A2) occurs outside the protein-DNA interface in protein-protein interfaces formed by folding of the hinge helices and interactions of the DNA binding domain (DBD) with the core of the repressor. Since urea interacts principally with amide surface, it is particularly well-suited to detect and quantify the extent of coupled folding and other large-scale remodeling events in the steps of protein-nucleic acid interactions and other protein associations.

Preferential hydration and solubility of proteins in aqueous solutions of polyethylene glycol

This paper is focused on the local composition around a protein molecule in aqueous mixtures containing polyethylene glycol (PEG) and the solubility of proteins in water+PEG mixed solvents. Experimental data from literature regarding the preferential binding parameter were used to calculate the excesses (or deficits) of water and PEG in the vicinity of beta-lactoglobulin, bovine serum albumin, lysozyme, chymotrypsinogen and ribonuclease A. It was concluded that the protein molecule is preferentially hydrated in all cases (for all proteins and PEGs investigated). The excesses of water and deficits of PEG in the vicinity of a protein molecule could be explained by a steric exclusion mechanism, i.e. the large difference in the sizes of water and PEG molecules. The solubility of different proteins in water+PEG mixed solvent was expressed in terms of the preferential binding parameter. The slope of the logarithm of protein (lysozyme, beta-lactoglobulin and bovine serum albumin) solubility versus the PEG concentration could be predicted on the basis of experimental data regarding the preferential binding parameter. For all the cases considered (various proteins, various PEGs molecular weights and various pHs), our theory predicted that PEG acts as a salting-out agent, conclusion in full agreement with experimental observations. The predicted slopes were compared with experimental values and while in some cases good agreement was found, in other cases the agreement was less satisfactory. Because the established equation is a rigorous thermodynamic one, the disagreement might occur because the experimental results used for the solubility and/or the preferential binding parameter do not correspond to thermodynamic equilibrium.

Solutes modify a conformational transition in a membrane transport protein

The bacterial outer-membrane vitamin B(12) transporter, BtuB, undergoes a dramatic order-to-disorder transition in its N-terminal energy-coupling motif (Ton box) upon substrate binding. Here, site-directed spin labeling (SDSL) is used to show that a range of solutes prevents this conformational change when ligand is bound to BtuB, resulting in a more ordered Ton box structure. For each solute examined, the data indicate that solutes effectively block this conformational transition through an osmotic mechanism. The molecular weight dependence of this solute effect has been examined for a series of polyethylene glycols, and a sharp molecular weight cutoff is observed. This cutoff indicates that solutes are preferentially excluded from a cavity within the protein as well as the protein surface. Furthermore, the sensitivity of the conformational change to solution osmolality is consistent with a structural model predicted by SDSL. When the Ton box is unfolded by detergents or mutations (rather than by ligand binding), solutes, such as polyethylene glycols and salts, also induce a more structured compacted conformation. These results suggest that conformational changes in this class of outer membrane transporters, which involve modest energy differences and changes in hydration, may be modulated by a range of solutes, including solutes typically used in protein crystallization.

Hydration changes in the association of Hoechst 33258 with DNA

The role of water in the interaction of Hoechst 33258 with the minor groove binding site of the (AATT)2 sequence was investigated using calorimetric and equilibrium constant measurements. Using isothermal titration calorimetry measurements, the heat capacity change for the reaction is -256 +/- 10 cal/(K mol of Hoechst). Comparison with the heat capacity changes based on area models supports the expulsion of water from the interface of the Hoechst-DNA complex. To further consider the role of water, the osmotic stress method was used to determine if the Hoechst association with DNA was coupled with hydration changes. Using four osmolytes with varying molecular weights and chemical properties, the Hoechst affinity for DNA decreases with increasing osmolyte concentration. From the dependence of the equilibrium constant on the solution osmolality, 60 +/- 13 waters are acquired in the complex relative to the reactants. It is proposed that the osmotic stress technique is measuring weakly bound waters that are not measured via the heat capacity changes.

The Effect of Salt on Protein Chemical Potential Determined by Ternary Diffusion in Aqueous Solutions

We use accurate thermodynamic derivatives extracted from high-precision measurements of the four volume-fixed diffusion coefficients in ternary solutions of lysozyme chloride in aqueous NaCl, NH4Cl, and KCl at pH 4.5 and 25 degrees C to (a) assess the relative contributions of the common-ion and nonideality effects to the protein chemical potential as a function of salt concentration, (b) compare the behavior of the protein chemical potential for the three salts, which we found to be consistent with the Hofmeister series, and (c) discuss our thermodynamic data in relation to the dependence of the protein solubility on salt concentration. The four diffusion coefficients are reported at 0.6 mM lysozyme chloride and 0.25, 0.5, 0.9, 1.2, and 1.5 M KCl and extend into the protein-supersaturated region. The chemical potential cross-derivatives are extracted from diffusion data using the Onsager reciprocal relation and the equality of molal cross-derivatives of solute chemical potentials. They are compared to those calculated previously from diffusion data for lysozyme in aqueous NaCl and NH4Cl. We estimate the effective charge on the diffusing lysozyme cation at the experimental concentrations. Our diffusion measurements on the three salts allowed us to analyze and interpret the four diffusion coefficients for charged proteins in the presence of 1:1 electrolytes. Our results may provide guidance to the understanding of protein crystallization.

Interactions of TrimethylamineN−Oxide and Water withcyclo-Alanylglycine

The osmolyte trimethylamine N-oxide (TMAO) is one of a family of compounds found in living systems that can stabilize biomolecular tertiary structures. As a step in exploring the interactions between this material and polyamino acids, we have determined intermolecular 1H{1H} nuclear Overhauser effects (NOEs) between the protons of cyclo-alanylglycine and protons of solvent components in TMAO-water solutions. Comparison of the results to effects predicted on the basis of the molecular shape of the dipeptide and experimental translational diffusion coefficients suggests that both water and TMAO molecules have properties in the vicinity of the dipeptide that are different from those in the bulk solution. Changes of local concentrations of water and TMAO and changes in the diffusive behavior of these components near the dipeptide are rejected as possible explanations of the discrepancies between observed and calculated Overhauser effects. Rather, it is concluded that TMAO molecules, and the water molecules associated with them, participate to some extent in the formation of long-lived solute-solvent complexes. The aliphatic alcohol tert-butyl alcohol is structurally similar to TMAO. Overhauser effect studies of its interaction with cyclo-alanylglycine in tert-butyl alcohol-water suggest similar kinds of interactions are present in this system but that they are significantly weaker, presumably because of the lower polarity of this alcohol compared to TMAO.

Relationship between preferential interaction of a protein in an aqueous mixed solvent and its solubility

The present paper is devoted to the derivation of a relation between the preferential solvation of a protein in a binary aqueous solution and its solubility. The preferential binding parameter, which is a measure of the preferential solvation (or preferential hydration) is expressed in terms of the derivative of the protein activity coefficient with respect to the water mole fraction, the partial molar volume of protein at infinite dilution and some characteristics of the protein-free mixed solvent. This expression is used as the starting point in the derivation of a relationship between the preferential binding parameter and the solubility of a protein in a binary aqueous solution. The obtained expression is used in two different ways: (1) to produce a simple criterion for the salting-in or salting-out by various cosolvents on the protein solubility in water, (2) to derive equations which predict the solubility of a protein in a binary aqueous solution in terms of the preferential binding parameter. The solubilities of lysozyme in aqueous sodium chloride solutions (pH=4.5 and 7.0), in aqueous sodium acetate (pH=8.3) and in aqueous magnesium chloride (pH=4.1) solutions are predicted in terms of the preferential binding parameter without any adjustable parameter. The results are compared with experiment, and for aqueous sodium chloride mixtures the agreement is excellent, for aqueous sodium acetate and magnesium chloride mixtures the agreement is only satisfactory.

Trehalose and glycerol stabilize and renature yeast inorganic pyrophosphatase inactivated by very high temperatures

A number of naturally occurring small organic molecules, primarily involved in maintaining osmotic pressure in the cell, display chaperone-like activity, stabilizing the native conformation of proteins, and protecting them from various kinds of stress. Most of them are sugars, polyols, amino acids or methylamines. Similar to molecular chaperones, most of these compounds have no substrate specificity, but some specifically stabilize certain proteins. In the present work, the capacity of trehalose and glycerol, two well-known osmolytes, to stabilize and renature inorganic pyrophosphatase is demonstrated. Both trehalose and glycerol significantly protect pyrophosphatase against thermoinactivation achieved by incubating the enzyme at temperatures up to 95 degrees C, and allow the enzyme already inactivated in the presence of these osmolytes to renature upon incubation at low temperatures. To the best of our knowledge, there are no data on the effects of these compounds on renaturation of thermoinactivated proteins. The correlation between the recovery of enzyme activity and structural changes indicated by fluorescence spectroscopy contribute to better understanding of the protein stabilization mechanism.

Evaluation of Hofmeister Effects on the Kinetic Stability of Proteins

Dissolved salts are known to affect properties of proteins in solution including solubility and melting temperature, and the effects of dissolved salts can be ranked qualitatively by the Hofmeister series. We seek a quantitative model to predict the effects of salts in the Hofmeister series on the deactivation kinetics of enzymes. Such a model would allow for a better prediction of useful biocatalyst lifetimes or an improved estimation of protein-based pharmaceutical shelf life. Here we consider a number of salt properties that are proposed indicators of Hofmeister effects in the literature as a means for predicting salt effects on the deactivation of horse liver alcohol dehydrogenase (HL-ADH), alpha-chymotrypsin, and monomeric red fluorescent protein (mRFP). We find that surface tension increments are not accurate predictors of salt effects but find a common trend between observed deactivation constants and B-viscosity coefficients of the Jones-Dole equation, which are indicative of ion hydration. This trend suggests that deactivation constants (log k(d,obs)) vary linearly with chaotropic B-viscosity coefficients but are relatively unchanged in kosmotropic solutions. The invariance with kosmotropic B-viscosity coefficients suggests the existence of a minimum deactivation constant for proteins. Differential scanning calorimetry is used to measure protein melting temperatures and thermodynamic parameters, which are used to calculate the intrinsic irreversible deactivation constant. We find that either the protein unfolding rate or the rate of intrinsic irreversible deactivation can control the observed deactivation rates.

Expanding functionality of RNA: synthesis and properties of RNA containing imidazole modified tandem G-U wobble base pairs

Imidazole modification at C-5 of uridine that is part of tandem G-U wobble base pairs causes slight reduction of thermal stability (DeltaDeltaG(0)(310) < 0.4 kcal mol(-1)) and relatively small change in hydration of short RNA helices.

Trehalose and 6-aminohexanoic acid stabilize and renature glucose-6-phosphate dehydrogenase inactivated by glycation and by guanidinium hydrochloride

A number of naturally occurring small organic molecules, primarily involved in maintaining osmotic pressure in the cell, display chaperone-like activity, stabilizing the native conformation of proteins and protecting them from various kinds of stress. Most of them are sugars, polyols, amino acids or methylamines. In addition to their intrinsic protein-stabilizing activity, these small organic stress molecules regulate the activity of some molecular chaperones, and may stabilize the folded state of proteins involved in unfolding or in misfolding diseases, such as Alzheimer's and Parkinson's diseases, or alpha1-antitrypsin deficiency and cystic fibrosis, respectively. Similar to molecular chaperones, most of these compounds have no substrate specificity, but some specifically stabilize certain proteins, e.g., 6-aminohexanoic acid (AHA) stabilizes apolipoprotein A. In the present work, the capacity of 6-aminohexanoic acid to stabilize non-specifically other proteins is demonstrated. Both trehalose and AHA significantly protect glucose-6-phosphate dehydrogenase (G6PD) against glycation-induced inactivation, and renatured enzyme already inactivated by glycation and by guanidinium hydrochloride (GuHCl). To the best of our knowledge, there are no data on the effect of these compounds on protein glycation. The correlation between the recovery of enzyme activity and structural changes indicated by fluorescence spectroscopy and Western blotting contribute to better understanding of the protein stabilization mechanism.

A contribution to the theory of preferential interaction coefficients

A simple and complete derivation of the relation between concentration-based preferential interaction coefficients and integrals over the relevant pair correlation functions is presented for the first time. Certain omissions from the original treatment of pair correlation functions in multicomponent thermodynamics are also addressed. Connections between these concentration-based quantities and the more common molality-based preferential interaction coefficients are also derived. The pair correlation functions and preferential interaction coefficients of both solvent (water) and cosolvent (osmolyte) in the neighborhood of a macromolecule contain contributions from short-range repulsions and generic long-range attractions originating from the macromolecule, as well as from osmolyte-solvent exchange reactions beyond the macromolecular surface. These contributions are evaluated via a heuristic analysis that leads to simple insightful expressions for the preferential interaction coefficients in terms of the volumes excluded to the centers of the water and osmolyte molecules and a sum over the contributions of exchanging sites in the surrounding solution. The preferential interaction coefficients are predicted to exhibit the experimentally observed dependence on osmolyte concentration. Molality-based preferential interaction coefficients that were reported for seven different osmolytes interacting with bovine serum albumin are analyzed using the this formulation together with geometrical parameters reckoned from the crystal structure of human serum albumin. In all cases, the excluded volume contribution, which is the volume excluded to osmolyte centers minus that excluded to water centers in units of V1, exceeds in magnitude the contribution of the exchange reactions. Under the assumption that the exchange contribution is dominated by sites in the first surface-contiguous layer, the ratio of the average exchange constant to its neutral random value is determined for each osmolyte. These ratios all lie in the range 1.0 +/- 0.15, which indicates rather slight deviations from random occupation near the macromolecular surface. Finally, a mechanism is proposed whereby the chemical identity of an osmolyte might be concealed from partially ordered multilayers of water in clefts, grooves, and pits, and its consequences are noted.

Stability of proteins: Temperature, pressure and the role of the solvent

We focus on the various aspects of the physics related to the stability of proteins. We review the pure thermodynamic aspects of the response of a protein to pressure and temperature variations and discuss the respective stability phase diagram. We relate the experimentally observed shape of this diagram to the low degree of correlation between the fluctuations of enthalpy and volume changes associated with the folding-denaturing transition and draw attention to the fact that one order parameter is not enough to characterize the transition. We discuss in detail microscopic aspects of the various contributions to the free energy gap of proteins and put emphasis on how a cosolvent may either enlarge or diminish this gap. We review briefly the various experimental approaches to measure changes in protein stability induced by cosolvents, denaturants, but also by pressure and temperature. Finally, we discuss in detail our own molecular dynamics simulations on cytochrome c and show what happens under high pressure, how glycerol influences structure and volume fluctuations, and how all this compares with experiments.

Gibbs-Duhem-based relationships among derivatives expressing the concentration dependences of selected chemical potentials for a multicomponent system

For a two-component system, a derivative that specifies the concentration-dependence of one chemical potential can be calculated from the corresponding derivative of the other chemical potential by applying the Gibbs-Duhem Equation. To extend the practical utility of this binary thermodynamic linkage to systems having any number of components, we present a derivation based on a previously unrecognized recursive relationship. Thus, for each independently variable component, kappa, any derivative of its chemical potential, mukappa, with respect to one of the mole ratios {mkappa identical with nkappa/nomega} is related to as a characteristic series of progressively higher order derivatives of muomega for a single "probe" component, omega, with respect to certain of the {mkappa}. For aqueous solutions in which omega is solvent water and one or more of the solutes (kappa) is dilute, under typical conditions each sum of terms expressing a derivative of mukappa consists of at most a few numerically significant contributions, which can be quantified, or at least estimated, by analyzing osmometric data to determine how the single chemical potential muomega depends on the {mkappa} without neglecting any significant contributions from the other components. Expressions derived here also will provide explicit criteria for testing various approximations built into alternative analytic strategies for quantifying derivatives that specify the {mkappa} dependences of mukappa for selected components. Certain quotients of these derivatives are of particular interest in so far as they gauge important thermodynamic effects due to "preferential interactions".

A comparison of the counteracting effects of glycine betaine and TMAO on the activity of RNase A in aqueous urea solution

Trimethylamine-N-oxide (TMAO) and glycine betaine are counteracting osmolytes found in cellular systems under osmotic stress, often in association with high urea concentrations. TMAO is a characteristic component of cartilaginous fish and marine molluscs, while glycine betaine is more widely distributed, occurring in plants, bacteria and the mammalian kidney. As part of a project to explain and understand the action of these methylamines, the RNase A-catalysed degradation of polyuridylic acid in the presence of urea and various osmolytes (0-1.0 M) was studied using (31)P Nuclear Magnetic Resonance spectroscopy. The decrease in reaction rate induced by urea could be fully recovered with 1 molar equivalent of trimethylamine-N-oxide or 1.4 molar equivalents of glycine betaine. These results indicate that the modification of RNase A activity induced by urea is not associated with gross irreversible structural changes and that both glycine betaine and trimethylamine-N-oxide have kinetically detectable counteracting effects.

Solutes probe hydration in specific association of cyclodextrin and adamantane

Using microcalorimetry, we follow changes in the association free energy of beta-cyclodextrin (CD) with the hydrophobic part of adamantane carboxylate (AD) due to added salt or polar (net-neutral) solutes that are excluded from the molecular interacting surfaces. Changes in binding constants with solution osmotic pressure (water activity) translate into changes in the preferential hydration upon complex formation. We find that these changes correspond to a release of 15-25 solute-excluding waters upon CD/AD association. Reflecting the preferential interaction of solute with reactants versus products, we find that changes in hydration depend on the type of solute used. All solutes used here result in a large change in the enthalpy of the CD-AD binding reaction. In one class of solutes, the corresponding entropy change is much smaller, while in the other class, the entropy change almost fully compensates the solute-specific enthalpy. For many of the solutes, the number of waters released correlates well with their effect on air-water surface tensions. We corroborate these results using vapor pressure osmometry to probe individually the hydration of reactants and products of association, and we discuss the possible interactions and forces between cosolute and hydrophobic surfaces responsible for different kinds of solute exclusion.

Hydration Forces Underlie the Exclusion of Salts and of Neutral Polar Solutes from Hydroxypropylcellulose

The distance dependence for the preferential exclusion of several salts and neutral solutes from hydroxypropyl cellulose (HPC) has been measured via the effect of these small molecules on the thermodynamic forces between HPC polymers in ordered arrays. The concentration of salts and neutral solutes decreases exponentially as the spacing between apposing nonpolar HPC surfaces decreases. For all solutes, the spatial decay lengths of this exclusion are remarkably similar to those observed between many macromolecules at close spacings where intermolecular forces have been ascribed to the energetics of water structuring. Exclusion magnitudes depend strongly on the nature and size of the particular salt or solute; for the three potassium salts studied, exclusion follows the anionic Hofmeister series. The change in the number of excess waters associated with HPC polymers is independent of solute concentration suggesting that the dominating interactions are between solutes and the hydrated polymer. These findings further confirm the importance of solvation interactions and reveal an unexpected unity of Hofmeister effects, preferential hydration, and hydration forces.

Rational design of solution additives for the prevention of protein aggregation

We have developed a statistical-mechanical model of the effect of solution additives on protein association reactions. This model incorporates solvent radial distribution functions obtained from all-atom molecular dynamics simulations of particular proteins into simple models of protein interactions. In this way, the effects of additives can be computed along the entire association/dissociation reaction coordinate. We used the model to test our hypothesis that a class of large solution additives, which we term "neutral crowders," can slow protein association and dissociation by being preferentially excluded from protein-protein encounter complexes, in a manner analogous to osmotic stress. The magnitude of this proposed "gap effect" was probed for two simple model systems: the association of two spheres and the association of two planes. Our results suggest that for a protein of 20 A radius, an 8 A additive can increase the free energy barrier for association and dissociation by as much as 3-6 kcal/mol. Because the proposed gap effect is present only for reactions involving multiple molecules, it can be exploited to develop novel additives that affect protein association reactions although having little or no effect on unimolecular reactions such as protein folding. This idea has many potential applications in areas such as the stabilization of proteins against aggregation during folding and in pharmaceutical formulations.

Structural Basis for the Binding of Compatible Solutes by ProX from the Hyperthermophilic Archaeon Archaeoglobus fulgidus

Compatible solutes such as glycine betaine and proline betaine serve as protein stabilizers because of their preferential exclusion from protein surfaces. To use extracellular sources of this class of compounds as osmo-, cryo-, or thermoprotectants, Bacteria and Archaea have developed high affinity uptake systems of the ATP-binding cassette type. These transport systems require periplasmic- or extracellular-binding proteins that are able to bind the transported substance with high affinity. Therefore, binding proteins that bind compatible solutes have to avoid the exclusion of their ligands within the binding pocket. In the present study we addressed the question to how compatible solutes can be effectively bound by a protein at temperatures around 83 degrees C as this is done by the ligand-binding protein ProX from the hyperthermophilic archaeon Archaeoglobus fulgidus. We solved the structures of ProX without ligand and in complex with both of its natural ligands glycine betaine and proline betaine, as well as in complex with the artificial ligand trimethylammonium. Cation-pi interactions and non-classical hydrogen bonds between four tyrosine residues, a main chain carbonyl oxygen, and the ligand have been identified to be the key determinants in binding the quaternary amines of the three investigated ligands. The comparison of the ligand binding sites of ProX from A. fulgidus and the recently solved structure of ProX from Escherichia coli revealed a very similar solution for the problem of compatible solute binding, although both proteins share only a low degree of sequence identity. The residues involved in ligand binding are functionally equivalent but not conserved in the primary sequence.

Bacterial osmosensing: roles of membrane structure and electrostatics in lipid–protein and protein–protein interactions

Bacteria act to maintain their hydration when the osmotic pressure of their environment changes. When the external osmolality decreases (osmotic downshift), mechanosensitive channels are activated to release low molecular weight osmolytes (and hence water) from the cytoplasm. Upon osmotic upshift, osmoregulatory transporters are activated to import osmolytes (and hence water). Osmoregulatory channels and transporters sense and respond to osmotic stress via different mechanisms. Mechanosensitive channel MscL senses the increasing tension in the membrane and appears to gate when the lateral pressure in the acyl chain region of the lipids drops below a threshold value. Transporters OpuA, BetP and ProP are activated when increasing external osmolality causes threshold ionic concentrations in excess of about 0.05 M to be reached in the proteoliposome lumen. The threshold activation concentrations for the OpuA transporter are strongly dependent on the fraction of anionic lipids that surround the cytoplasmic face of the protein. The higher the fraction of anionic lipids, the higher the threshold ionic concentrations. A similar trend is observed for the BetP transporter. The lipid dependence of osmotic activation of OpuA and BetP suggests that osmotic signals are transmitted to the protein via interactions between charged osmosensor domains and the ionic headgroups of the lipids in the membrane. The charged, C-terminal domains of BetP and ProP are important for osmosensing. The C-terminal domain of ProP participates in homodimeric coiled-coil formation and it may interact with the membrane lipids and soluble protein ProQ. The activation of ProP by lumenal, macromolecular solutes at constant ionic strength indicates that its structure and activity may also respond to macromolecular crowding. This excluded volume effect may restrict the range over which the osmosensing domain can electrostatically interact. A simplified version of the dissociative double layer theory is used to explain the activation of the transporters by showing how changes in ion concentration could modulate interactions between charged osmosensor domains and charged lipid or protein surfaces. Importantly, the relatively high ionic concentrations at which osmosensors become activated at different surface charge densities compare well with the predicted dependence of 'critical' ion concentrations on surface charge density. The critical ion concentrations represent transitions in Maxwellian ionic distributions at which the surface potential reaches 25.7 mV for monovalent ions. The osmosensing mechanism is qualitatively described as an "ON/OFF switch" representing thermally relaxed and electrostatically locked protein conformations.

The extended interface: measuring non-local effects in biomolecular interactions

Improvements in the sensitivity and availability of biophysical techniques for the detection of the formation of complexes in solution are revealing that the effects of binding are not restricted to the direct contacts between the biomolecules or even to a localised site. Rather, information about the binding event is transmitted throughout the biomolecules and the surrounding solution through changes in the hydrogen bonding, hydration and electrostatic field as the complex is formed. Calorimetric, volumetric and NMR methods are beginning to provide a quantitative view of the nature and thermodynamic consequences of this extended interface, and the resulting data pose a major challenge for computational models of binding.

Differences between EcoRI nonspecific and "star" sequence complexes revealed by osmotic stress

The binding of the restriction endonuclease EcoRI to DNA is exceptionally specific. Even a single basepair change ("star" sequence) from the recognition sequence, GAATTC, decreases the binding free energy of EcoRI to values nearly indistinguishable from nonspecific binding. The difference in the number of waters sequestered by the protein-DNA complexes of the "star" sequences TAATTC and CAATTC and by the specific sequence complex determined from the dependence of binding free energy on water activity is also practically indistinguishable at low osmotic pressures from the 110 water molecules sequestered by nonspecific sequence complexes. Novel measurements of the dissociation rates of noncognate sequence complexes and competition equilibrium show that sequestered water can be removed from "star" sequence complexes by high osmotic pressure, but not from a nonspecific complex. By 5 Osm, the TAATTC "star" sequence complex has lost almost 90 of the approximately 110 waters initially present. It is more difficult to remove water from the CAATTC "star" sequence complex. The sequence dependence of water loss correlates with the known sequence dependence of "star" cleavage activity.

Preferential interaction of dimethyl sulfoxide and phosphatidyl choline membranes

The interaction free energy of dimethyl sulfoxide (DMSO) and two types phospholipid membranes has been assessed from measurements of vapor pressure. The lipids were phosphatidyl cholines with respectively (14:0/14:0) (DMPC) and (16:0/18:1) (POPC) fatty acid chains. The results were expressed in terms of the iso-osmolal preferential interaction parameter, Gamma(mu1), which remained negative under all experimental conditions investigated here. This shows that water-membrane interactions are more favorable than DMSO-membrane interactions. This condition is known as preferential exclusion of DMSO (or preferential hydration of the membrane), and implies that the local (interfacial) concentration of the solute is reduced compared to the bulk. At room temperature and 1 m DMSO, Gamma(mu1) was -0.3 to -0.4 for both lipids. This corresponds to a sizable reduction in the DMSO concentration in a zone including at least the first two hydration layers of the membrane. Possible origins of the preferential exclusion are discussed. As a direct consequence of the pronounced preferential exclusion, DMSO generates an osmotic stress at the membrane interface. This tends to stabilize lipid phases of low surface areas and to withdraw water from multilamellar stacks of membranes. Based on this, we suggest that the preferential exclusion of DMSO explains both the modulation of phase behavior and the constriction of multilamellar aggregates induced by this solute.

Betaine in human nutrition

Betaine is distributed widely in animals, plants, and microorganisms, and rich dietary sources include seafood, especially marine invertebrates ( approximately 1%); wheat germ or bran ( approximately 1%); and spinach ( approximately 0.7%). The principal physiologic role of betaine is as an osmolyte and methyl donor (transmethylation). As an osmolyte, betaine protects cells, proteins, and enzymes from environmental stress (eg, low water, high salinity, or extreme temperature). As a methyl donor, betaine participates in the methionine cycle-primarily in the human liver and kidneys. Inadequate dietary intake of methyl groups leads to hypomethylation in many important pathways, including 1) disturbed hepatic protein (methionine) metabolism as determined by elevated plasma homocysteine concentrations and decreased S-adenosylmethionine concentrations, and 2) inadequate hepatic fat metabolism, which leads to steatosis (fatty accumulation) and subsequent plasma dyslipidemia. This alteration in liver metabolism may contribute to various diseases, including coronary, cerebral, hepatic, and vascular diseases. Betaine has been shown to protect internal organs, improve vascular risk factors, and enhance performance. Databases of betaine content in food are being developed for correlation with population health studies. The growing body of evidence shows that betaine is an important nutrient for the prevention of chronic disease.

Comparison of the effect of water release on the interaction of the Saccharomyces cerevisiae TATA binding protein (TBP) with "TATA Box" sequences composed of adenosine or inosine

The formation of sequence-specific complexes of TATA binding protein (TBP) with the minor groove of DNA results in the burial of large nonpolar surfaces and the exclusion of water from these interfaces. The release of water is thus expected to provide a significant entropic driving force for formation of the transcription-preinitiated complexes mediated by the binding of TBP to specific sequences. In this article are described equilibrium-binding studies of Saccharomyces cerevisiae TBP to 14 bp oligonucleotides bearing either the tightly bound and efficiently transcribed adenovirus major late promoter (TATAAAAG) or its inosine-substituted derivative (TITIIIIG) as a function of neutral osmolyte concentration. These two DNA sequences present the same pattern of minor groove hydrogen-bond donors and acceptors to the protein. TBP-DNA complex formation was monitored by steady-state fluorescence resonance energy transfer measurements of the oligonucleotides end-labeled with fluorescein (donor) and TAMRA (acceptor). Correct interpretation of the results obtained with the inosine-substituted sequence required careful consideration of the optical properties of the dyes as a function of osmolyte concentration to demonstrate that the relative change in the end-to-end distances for TATAAAAG- and TITIIIIG-bearing oligonucleotides is the same upon TBP binding. Although the affinity of TBP is slightly greater for the adenosine compared with the inosine-substituted TATA sequence in the absence of osmolyte, the end-to-end distances of the bound DNA in complex with TBP, the enthalpic and electrostatic components of binding, are identical within experimental precision. However, approximately 18 additional molecules of water are released upon TBP binding the TATAAAAG as compared with the TITIIIIG sequence resulting in an entropic advantage to the binding of the natural promoter sequence. These results are considered with regard to differences in the flexibility and hydration of the two DNA sequences.

Preferential hydration and the exclusion of cosolvents from protein surfaces

Protein stability is enhanced by the addition of osmolytes, such as sugars and polyols and inert crowders, such as polyethylene glycols. This stability enhancement has been quantified by the preferential hydration parameter which can be determined by experiments. To understand the mechanism of protein stability enhancement, we present a statistical mechanical analysis of the preferential hydration parameter based upon Kirkwood-Buff theory. Previously, the preferential hydration parameter was interpreted in terms of the number of hydration waters, as well as the cosolvent exclusion volume. It was not clear how accurate these interpretations were, nor what the relationship is between the two. By using the Kirkwood-Buff theory and experimental data, we conclude that the contribution from the cosolvent exclusion dominantly determines the preferential hydration parameters for crowders. For osmolytes, although the cosolvent exclusion largely determines the preferential hydration parameters, the contribution from hydration may not be negligible.

Gibbs adsorption isotherm combined with Monte Carlo sampling to see action of cosolutes on protein folding

Driven by conditions set by smaller solutes, proteins fold and unfold. Experimentally, these conditions are stated as intensive variables--pH and other chemical potentials--as though small solutes were infinite resources that come at an externally varied free energy cost. Computationally, the finite spaces of simulation allow only fixed numbers of these solutes. By combining the analytic Gibbs adsorption isotherm with the computational Monte Carlo sampling of polymer configurations, we have been able to overcome an inherent limitation of computer simulation. The idea is to compute analytically the free energy changes wrought by solutes on each particular configuration. Then numerical computation is needed only to sample the set of configurations as efficiently as when no bathing solute is present. For illustration, the procedure is applied to an idealized two-dimensional heteropolymer to yield lessons about the effect of cosolutes on protein stability.

Preferential interactions in aqueous solutions of urea and KCl

A quantitative characterization of the thermodynamic effects due to interactions of salt ions and urea in aqueous solution is needed for rigorous analyses of the effects of changing urea concentration on biopolymer processes in solutions that also contain salt. Therefore, we investigate preferential interactions in aqueous solutions containing KCl and urea by using vapor pressure osmometry (VPO) to measure osmolality as a function of the molality of urea (component 3) over the range 0.09<or=m(3)<or=1.65 m at two fixed molalities of KCl (component 2) (m(2)=0.212 and 0.427 m). With this experimental input and corresponding VPO measurements on solutions that contain only urea or KCl, we evaluate approximately the chemical potential derivative micro(23)=( partial differential micro(KCl)/ partial differential m(urea))(T,P,m(KCl))=( partial differential micro(urea)/ partial differential m(KCl))(T,P,m(urea))= micro(32) and hence the preferential interaction coefficients Gammamicro(3) and Gammamicro(1),micro(3). These results show that for water-KCl-urea solutions neither of these coefficients is determined primarily by contributions from thermodynamic nonideality to micro(23). In aqueous solutions containing a biopolymer and a small solute, the contribution of ideal mixing entropy to micro(23) is negligible in comparison with the experimental uncertainty, whereas in KCl-urea solutions the contribution due to ideal mixing entropy accounts for at least half of the magnitude of micro(23). For comparison, we analyze literature data for NaCl-urea interactions and find again that nonideality makes a smaller contribution to micro(23) than does ideal mixing entropy. In contrast, for aqueous solutions of urea and the protein bovine serum albumin, the experimentally determined contribution of nonideality to micro(23) exceeds the contribution of ideal mixing by a factor of approximately 2 x 10(2).

Osmo-regulation of bacterial transcription via poised RNA polymerase

Adaptation to high-salt environments is critical for the survival of a wide range of cells, especially for pathogenic bacteria that colonize the animal gut and urinary tract. The adaptation strategy involves production of the salt potassium glutamate, which induces a specific gene expression program that produces electro-neutral osmolytes while inhibiting general sigma(70) transcription. These data show that in Escherichia coli potassium glutamate stimulates transcription by disengaging inhibitory polymerase interactions at a sigma(38) promoter. These occur in an upstream region that is marked by an osmotic shock promoter DNA consensus sequence. The disruption activates a poised RNA polymerase to transcribe. This transcription program leads to the production of osmolytes that are shown to have only minor effects on transcription and therefore help to restore normal cell function. An osmotic shock gene expression cycle is discussed.

Stabilization of nucleic acid triplexes by high concentrations of sodium and ammonium salts follows the Hofmeister series

The thermal stability of the triplexes d(C(+)-T)(6):d(A-G)(6);d(C-T)(6) and d(T)(21):d(A)(21);d(T)(21) was studied in the presence of high concentrations of the anions Cl(-), HPO(4)(2-), CH(3)COO(-), SO(4)(2-) and ClO(4)(-). Thermally-induced triplex and duplex transitions were identified by UV- and CD-spectroscopy and T(m) values were determined from melting profiles. A thermodynamic analysis of triplex transitions shows the limitations of commonly used treatments for determining the associated release or uptake of salt, solute or water. Enhancement of the stability of these triplexes follows the rank order of the Hofmeister series for anions of sodium and ammonium salts, whereas water structure-breaking solutes have the opposite effect. The rank order for the Hofmeister series ClO(4)(-)<I(-)<Br(-)<Cl(-)<HPO(4)(2-)<SO(4)(2-)is shown to follow their effective surface charge densities.

Thermoprotection of Bacillus subtilis by exogenously provided glycine betaine and structurally related compatible solutes: involvement of Opu transporters

Bacillus subtilis possesses five osmotically regulated transporters (Opu) for the uptake of various compatible solutes for osmoprotective purposes. We have now found that compatible solutes also function as thermoprotectants for B. subtilis. Low concentrations of glycine betaine enhanced the growth of the B. subtilis wild-type strain JH642 at its maximal growth temperature (52 degrees C) but did not allow an extension of the upper growth limit. A similar enhancement in the growth of B. subtilis was also observed by the addition of several other compatible solutes that are structurally related to glycine betaine or by the addition of proline. Each of these compatible solutes was taken up under heat stress by the cell through the same Opu transporters that are used for their acquisition under osmostress conditions. Northern blot analysis revealed a moderate increase in transcription of the structural genes for each of the Opu transport systems in cells that were propagated at 52 degrees C. In contrast, the uptake level of radiolabeled glycine betaine was very low under high-temperature growth conditions but nevertheless allowed the buildup of an intracellular glycine betaine pool comparable to that found in cells grown at 37 degrees C in the absence of salt stress. Although exogenously added glutamate has only a limited osmoprotective potential for B. subtilis, it was found to be a very effective thermoprotectant. Collectively, our data demonstrate thermoprotection by a variety of compatible solutes in B. subtilis, thus ascribing a new physiological function for this class of compounds in this microorganism and broadening the physiological role of the known osmoprotectant uptake systems (Opu).

Estimating hydration changes upon biomolecular reactions from osmotic stress, high pressure, and preferential hydration experiments

How do we estimate, from thermodynamic measurements, the number of water molecules adsorbed or released from biomolecules as a result of a biochemical process such as binding and allosteric effects? Volumetric and osmotic stress analyses are established methods for estimating water numbers; however, these techniques often yield conflicting results. In contrast, Kirkwood-Buff theory offers a novel way to calculate excess hydration number from volumetric data, provides a quantitative condition to gauge the accuracy of osmotic stress analysis, and clarifies the relationship between osmotic and volumetric analyses. I have applied Kirkwood-Buff theory to calculate water numbers for two processes: (i) the allosteric transition of hemoglobin and (ii) the binding of camphor to cytochrome P450. I show that osmotic stress analysis may overestimate hydration number changes for these processes.

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2'-O-methyl (2'-OMe) oligonucleotides GCGAAUUCGC, (UA)6 and (CG)3 using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Delta(n)W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Delta(n)W was 3-4 for RNA duplexes and 2-3 for DNA and 2'-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2'-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2'-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30? DNA

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition. The variations of l0 and ET with -ln(aw) are similar for acetamide and ethylene glycol at 37 degrees C. Fitting the theory to those data yields the range n = 103-125 for ethylene glycol and n = 71-113 for acetamide, depending on the assumed value of ET for the dehydrated state. The cooperative domain size of the two-state transition is estimated to exceed about 25-30 base pairs (bp). Between 0 and 19.4 w/v % ethylene glycol, the torsion elastic constant, measured by time-resolved fluorescence polarization anisotropy (FPA), increases by 1.37-fold, whereas the measured ET decreases by 1.15-fold over that same range. The implied decrease in bending rigidity over that range is by a factor of about 0.7. The variations of l0 and ET with increasing -ln(aw) due to added ethylene glycol at 37 degrees C are far smaller than the corresponding variations observed previously at 14 and 15 degrees C. However, at 21 degrees C, upon adding either ethylene glycol or acetamide, l0 and ET initially decline steeply with increasing -ln(aw), with slopes possibly comparable to those seen at 14 and 15 degrees C, but then flatten out and follow curves similar to those at 37 degrees C. Possible origins of such mixed behavior are discussed. The effects of betaine at both 37 and 21 degrees C differ qualitatively and quantitatively in various respects from those of ethylene glycol and acetamide. Upon adding sucrose, l0 initially jumps to higher plateaus at both 37 and 21 degrees C, but its effects on ET cannot be reliably assessed, due to the limited range of -ln(aw).