Record Mt

Quantifying why urea is a protein denaturant, whereas glycine betaine is a protein stabilizer

To explain the large, opposite effects of urea and glycine betaine (GB) on stability of folded proteins and protein complexes, we quantify and interpret preferential interactions of urea with 45 model compounds displaying protein functional groups and compare with a previous analysis of GB. This information is needed to use urea as a probe of coupled folding in protein processes and to tune molecular dynamics force fields. Preferential interactions between urea and model compounds relative to their interactions with water are determined by osmometry or solubility and dissected using a unique coarse-grained analysis to obtain interaction potentials quantifying the interaction of urea with each significant type of protein surface (aliphatic, aromatic hydrocarbon (C); polar and charged N and O). Microscopic local-bulk partition coefficients Kp for the accumulation or exclusion of urea in the water of hydration of these surfaces relative to bulk water are obtained. Kp values reveal that urea accumulates moderately at amide O and weakly at aliphatic C, whereas GB is excluded from both. These results provide both thermodynamic and molecular explanations for the opposite effects of urea and glycine betaine on protein stability, as well as deductions about strengths of amide NH—amide O and amide NH—amide N hydrogen bonds relative to hydrogen bonds to water. Interestingly, urea, like GB, is moderately accumulated at aromatic C surface. Urea m-values for protein folding and other protein processes are quantitatively interpreted and predicted using these urea interaction potentials or Kp values.

Interpretation of preferential interaction coefficients of nonelectrolytes and of electrolyte ions in terms of a two-domain model.

For a three-component system consisting of solvent (1), polymer or polyelectrolyte (2J), and a nonelectrolyte or electrolyte solute (3), a two-domain description is developed to describe thermodynamic effects of interactions between solute components (2J) and (3). Equilibrium dialysis, which for an electrolyte solute produces the Donnan distribution of ions across a semipermeable membrane, provides a fundamental basis for this two-domain description whose applicability is not restricted, however, to systems where dialysis equilibrium is established. Explicit expressions are obtained for the solute-polymer preferential interaction coefficient gamma 3,2J (nonelectrolyte case) and for gamma +,2J and gamma -,2J, which are corresponding coefficients defined for single (univalent) cations and anions, respectively: gamma +,2J = magnitude of ZJ + gamma -,2J = 0.5(magnitude of ZJ + B-,2J + B+,2J) - B1,2Jm3/m1 Here B+,2J, B-,2J, and B1,2J are defined per mole of species J, respectively, as the number of moles of cation, anion, and water included within the local domains that surround isolated molecules of J; ZJ is the charge on J; m3 is the molal concentration of uniunivalent electrolyte, and m1 = 55.5 mol/kg for water. Incorporating this result into a general thermodynamic description(derived by us elsewhere) of the effects of the activity a+ of excess uniunivalent salt on an equilibrium involving two or more charged species J (each of which is dilute in comparison with the salt) yields:SaKobs bS/d a+ A(r+2J r 2j) A(B+2J B-2 2B12Jm3/m1)where KObS is an equilibrium quotient defined in terms of the molar concentrations of the participants, J, and A denotes astoichio metrically weighted combination of terms pertaining to the reactant(s) and product(s). The derivation presented here does not depend on any particular molecular model for salt-polyelectrolyte (or solute-polymer) interactions; it therefore generalizes our earlier (1978) derivation.

Conformational transitions of duplex and triplex nucleic acid helices: thermodynamic analysis of effects of salt concentration on stability using preferential interaction coefficients.

For order-disorder transitions of double- and triple-stranded nucleic acid helices, the midpoint temperatures Tm depend strongly on a +/-, the mean ionic activity of uniunivalent salt. Experimental determinations of dTm/d ln a +/- and of the enthalpy change (delta H(o)) accompanying the transition in excess salt permit evaluation of delta gamma, the stoichiometrically weighted combination of preferential interaction coefficients, each of which reflects thermodynamic effects of interactions of salt ions with a reactant or product of the conformational transition (formula; see text) Here delta H(o) is defined per mole of nucleotide by analogy to delta gamma. Application of Eq. 1 to experimental values of delta H(o) and Tm yields values of delta gamma for the denaturation of B-DNA over the range of NaCl concentrations 0.01-0.20 M (Privalov et al. (1969), Biopolymers 8,559) and for each of four order-disorder transitions of poly rA.(poly rU)n, n = 1, 2 over the range of NaCl concentrations 0.01-1.0 M (Krakauer and Sturtevant (1968), Biopolymers 6, 491). For denaturation of duplexes and triplexes, delta gamma is negative and not significantly dependent on a +/-, but delta gamma is positive and dependent on a +/- for the disproportionation transition of poly rA.poly rU duplexes. Quantitative interpretations of these trends and magnitudes of delta gamma in terms of coulombic and excluded volume effects are obtained by fitting separately each of the two sets of thermodynamic data using Eq. 1 with delta gamma PB evaluated from the cylindrically symmetric Poisson-Boltzmann (PB) equation for a standard model of salt-polyelectrolyte solutions. The only structural parameters required by this model are: b, the mean axial distance between the projections of adjacent polyion charges onto the cylindrical axis; and a, the mean distance of closest approach between a salt ion center and the cylindrical axis. Fixing bMS and aMS for the multi-stranded (ordered) conformations, we determined the corresponding best fitted values of bSS and aSS for single-stranded RNA and DNA. The resulting best fitted values of aSS are systematically less than aDS by 2-4 A. Uncertainty in the best-fitted values of bSS is significantly lower than in the aSS, because bMS is known with relatively high precision and because the larger uncertainty in aMS has a relatively small effect on the best-fitted values of bSS:bSS = 3.2 +/- 0.6 A for single-stranded poly rA and poly rU; and bSS = 3.4 +/- 0.2 A for single-stranded DNA. These values are approximately one-halt of those expected for a fully extended single-stranded conformation. With the best fitted values of ass and bss, our calculations of delta gamma PB are in close quantitative agreement with experimental observations on each of five nucleic acid order-disorder transitions.

Grand canonical Monte Carlo molecular and thermodynamic predictions of ion effects on binding of an oligocation (L8+) to the center of DNA oligomers

Grand canonical Monte Carlo (GCMC) simulations are reported for aqueous solutions containing excess univalent salt (activities a +/- = 1.76-12.3 mM) and one of the following species: an octacationic rod-like ligand, L8+; a B-DNA oligomer with N phosphate charges (8 < or = N < or = 100); or a complex resulting from the binding of L8+ at the center of an N-mer (24 < or = N < or = 250). Simplified models of these multiply charged species are used in the GCMC simulations to predict the fundamental coulombic contributions to the following experimentally relevant properties: 1) the axial distance over which ligand binding affects local counterion concentrations at the surface of the N-mer; 2) the dependence on N of GCMC preferential interaction coefficients, gamma 32MC identical to delta C3/delta C2l a +/-, T, where C3 and C2 are, respectively, the molar concentrations of salt and the multiply charged species (ligand, N-mer or complex); and 3) the dependence on N of SaKobs identical to d in Kobs/d in a +/- = delta (magnitude of ZJ + 2 gamma 32J), where Kobs is the equilibrium concentration quotient for the binding of L8+ to the center of an N-mer and delta denotes the stoichiometric combination of terms, each of which pertains to a reactant or product J having magnitude of ZJ charges. The participation of electrolyte ions in the ligand binding interaction is quantified by the magnitude of SaKobs, which reflects the net (stoichiometrically weighted) difference in the extent of thermodynamic binding of salt ions to the products and reactants. Results obtained here from GCMC simulations yield a picture of the salient molecular consequences of binding a cationic ligand, as well as thermodynamic predictions whose applicability can be tested experimentally. Formation of the central complex is predicted to cause a dramatic reduction in the surface counterion (e.g., Na+) concentration over a region including but extending well beyond the location of the ligand binding site. For binding a cationic ligand, SaKobs is predicted to be negative, indicating net electrolyte ion release in the binding process. At small enough N, -SaKobs is predicted to decrease strongly toward zero with decreasing N. At intermediate N, -SaKobs appears to exceed its limiting value as N-->infinity.

Preferential interactions between small solutes and the protein backbone: A computational analysis

To improve our understanding of the effects of small solutes on protein stability, we conducted atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), and a triglycine peptide, which is a good model for a protein backbone. Multiple solute concentrations were analyzed, and each solute-peptide-water ternary system was studied with approximately 200-300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Gamma(23)) and experimentally measured values suggests that semiquantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein, and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Gamma(23)) results can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from an existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Gamma(23) into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results with GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.

Analyses of thermodynamic data for concentrated hemoglobin solutions using scaled particle theory: implications for a simple two-state model of water in thermodynamic analyses of crowding in vitro and in vivo.

Quantitative description of the thermodynamic consequences of macromolecular crowding (excluded volume nonideality) is an important component of analyses of the thermodynamics and kinetics of noncovalent interactions of biopolymers in vivo and in concentrated polymer solutions in vitro. By analyzing previously published thermodynamic data, we have investigated extensively the comparative applicability of two forms of scaled particle theory (SPT). In both forms, macromolecules are treated as hard spheres, but MSPT, introduced by Ross and Minton, treats the solvent as a structureless continuum, whereas bulk water molecules are included explicitly as hard spheres in BSPT, an approach developed by Berg. Here we use both MSPT and BSPT to calculate the excluded volume component of the macromolecular activity coefficient of hemoglobin (Hb) at concentrations up to 509 mg/ml by fitting osmotic pressure data for Hb and sedimentation equilibrium data for Hb and sickle-cell Hb (HbS). Both forms of SPT also are used here to analyze the effects of other globular proteins (BSA and Hb) on the solubility of HbS. In applying MSPT and BSPT to analyze macromolecular crowding, the extent of hydration delta Hb (in gH2O/gprotein) is introduced as an adjustable parameter to specify the effective (hard sphere) radius of hydrated Hb. In our nonlinear least-squares fittings based on BSPT, the hard sphere radius of bulk water molecules is either fixed at 1.375 A or floated. Although both forms of SPT yield good fittings (with different values of delta Hb) at Hb concentrations up to 350 mg/ml, only BSPT gives good fittings of all available Hb osmotic pressure data as well as of the sedimentation equilibrium and solubility data. Only BSPT predicts values for delta Hb (approximately 0.5-0.6 g/g) in the range obtained for Hb from hydrodynamic measurements (approximately 0.36-0.78 g/g). These findings indicate the applicability, at least in the context of BSPT, of a simple two-state classification of water (bulk water and water of macromolecular hydration) as a basis for interpreting excluded volume nonideality in concentrated solutions of globular proteins.

Prediction of Hofmeister ion effects on biopolymer processes

At moderate to high concentrations, salt ions exert a wide range of effects on protein folding and other protein processes, from extremely destabilizing (GuH+, SCN-) to very stabilizing (SO42-). The Hofmeister series is a qualitative ranking of these effects, originally based on the effectiveness of salts as protein precipitants and subsequently observed for other processes including creating air-water surface and dissolving hydrocarbons in water. Recently surface spectroscopy, molecular dynamics simulations and molecular thermodynamic analysis of surface tension and model compound solubility data have all provided evidence that local accumulation or exclusion of individual salt ions, relative to bulk concentrations, is responsible for their Hofmeister effects. In particular, application of a novel two-domain salt ion partitioning model (SPM) to analyze effects of Hofmeister salts on the surface tension of water and on hydrocarbon and peptide solubility (Pegram & Record, J. Phys. Chem. B 112, 9428 (2008); 111, 5411 (2007)) provides a quantitative molecular thermodynamic description of the individual partitioning of salt cations and anions at uncharged interfaces, with predictive capability. This analysis shows that the Hofmeister rank order of ions arises from their interactions with nonpolar surface. Surface-bulk partition coefficients of ions obtained from hydrocarbon and amide model compound solubility data, together with a coarse-grained description of functional groups that make up molecular surfaces, allow the quantitative prediction of Hofmeister (noncoulombic) salt effects on micelle formation, protein folding, protein crystallization and DNA helix formation.This work is supported by NIH GM47022.