Comparison of Methods for Characterizing Nonideal Solute Self-Association by Sedimentation Equilibrium

We have examined in detail analytical solutions of expressions for sedimentation equilibrium in the analytical ultracentrifuge to describe self-association under nonideal conditions. We find that those containing the radial dependence of total solute concentration that incorporate the Adams-Fujita assumption for composition-dependence of activity coefficients reveal potential shortcomings for characterizing such systems. Similar deficiencies are shown in the use of the NONLIN software incorporating the same assumption about the interrelationship between activity coefficients for monomer and polymer species. These difficulties can be overcome by iterative analyses incorporating expressions for the composition-dependence of activity coefficients predicted by excluded volume considerations. A recommendation is therefore made for the replacement of current software packages by programs that incorporate rigorous statistical-mechanical allowance for thermodynamic nonideality in sedimentation equilibrium distributions reflecting solute self-association.     

Allowance for effects of thermodynamic nonideality in sedimentation equilibrium distributions reflecting protein dimerization

This reexamination of a high-speed sedimentation equilibrium distribution for α-chymotrypsin under slightly acidic conditions (pH 4.1, I(M) 0.05) has provided experimental support for the adequacy of nearest-neighbor considerations in the allowance for effects of thermodynamic nonideality in the characterization of protein self-association over a moderate concentration range (up to 8 mg/mL). A widely held but previously untested notion about allowance for thermodynamic nonideality effects is thereby verified experimentally. However, it has also been shown that a greater obstacle to better characterization of protein self-association is likely to be the lack of a reliable estimate of monomer net charge, a parameter that has a far more profound effect on the magnitude of the measured equilibrium constant than any deficiency in current procedures for incorporating the effects of thermodynamic nonideality into the analysis of sedimentation equilibrium distributions reflecting reversible protein self-association.     

Ultracentrifugal studies of the effect of molecular crowding by trimethylamine N-oxide on the self-association of muscle glycogen phosphorylase b.

The suitability of sedimentation equilibrium for characterizing the self-association of muscle glycogen phosphorylase b has been reappraised. Whereas sedimentation equilibrium distributions for phosphorylase b in 40 mM Hepes buffer (pH 6.8) supplemented with 1 mM AMP signify a lack of chemical equilibrium attainment, those in buffer supplemented additionally with potassium sulfate conform with the requirements of a dimerizing system in chemical as well as sedimentation equilibrium. Because the rate of attainment of chemical equilibrium under the former conditions is sufficiently slow to allow resolution of the dimeric and tetrameric enzyme species by sedimentation velocity, this procedure has been used to examine the effects of thermodynamic nonideality arising from molecular crowding by trimethylamine N-oxide on the self-association behaviour of phosphorylase b. In those terms the marginally enhanced extent of phosphorylase b self-association observed in the presence of high concentrations of the cosolute is taken to imply that the effects of thermodynamic nonideality on the dimer-tetramer equilibrium are being countered by those displacing the T<==>R isomerization equilibrium for dimer towards the smaller, nonassociating T state. Because the R state is the enzymically active form, an inhibitory effect is the predicted consequence of molecular crowding by high concentrations of unrelated solutes. Thermodynamic nonideality thus provides an alternative explanation for the inhibitory effects of high concentrations of glycerol, sucrose and ethylene glycol on phosphorylase b activity, phenomena that have been attributed to extremely weak interaction of these cryoprotectants with the T state of the enzyme.     

Elimination of self-association as the source of the thermodynamic nonideality in aqueous proline solutions

The effect of high concentrations of proline on the diffusion coefficient of water has been examined to assess the extent to which the resulting thermodynamic nonideality could be explained on the statistical-mechanical basis of excluded volume. In fact, such a space-filling role not only accounts for the proline concentration-dependence of the diffusion coefficient of water but it also accounts for the nonideality of proline in freezing point depression and isopiestic measurements. These findings refute the conclusion (Schobert, B. and Tschesche, H. (1978) Biochim. Biophys. Acta 541, 270–277) that the stabilization of enzyme structure by high concentrations of proline stems from self-association of the imino acid via intermolecular hydrogen bonding; and thereby support the concept that the protective effect of proline on enzyme stability must reside mainly in its action as an inert, space-filling solute.     

Analysis of sedimentation equilibrium distributions reflecting nonideal macromolecular associations

A rigorous statistical-mechanical approach is adopted to derive general quantitative expressions that allow for the effects of thermodynamic nonideality in equilibrium measurements reflecting interaction between dissimilar macromolecular reactants. An analytical procedure based on these expressions is then formulated for obtaining global estimates of equilibrium constants and the corresponding reference thermodynamic activities of the free reactants in each of several sedimentation equilibrium experiments. The method is demonstrated by application to results from an ultracentrifugal study of an electrostatic interaction between ovalbumin and cytochrome c (Winzor, D. J., M. P. Jacobsen, and P. R. Wills. 1998. Biochemistry. 37:2226-2233). It is demonstrated that reliable estimates of relevant thermodynamic parameters are extracted from the data through statistical analysis by means of a simple nonlinear fitting procedure.     

Effect of sucrose on the dimerization of alpha-chymotrypsin. Allowance for thermodynamic nonideality arising from the presence of a small inert solute

The space-filling effects of sucrose on the dimerization of alpha-chymotrypsin have been investigated by sedimentation equilibrium studies on the enzyme in acetate-chloride buffer, pH 3.9, I 0.2. From the extent of enhancement of the apparent dimerization constant in the presence of 0.05-0.16 M sucrose, it is concluded that this effect of thermodynamic nonideality finds quantitative explanation in terms of excluded volume. However, the suggested approximation that the radius of an inert small solute would be sufficiently small to be neglected in the calculation of covolumes (D.J. Winzor and P.R. Wills, Biophys. Chem. 25 (1986) 243) has not withstood the more stringent test afforded by the present study of alpha-chymotrypsin dimerization. A value of 0.34 nm for the effective thermodynamic radius of sucrose was inferred from the covolume for self-interaction obtained by frontal gel chromatography on Sephadex G-10 under the conditions of the ultracentrifugal studies. Finally, results of sedimentation equilibrium experiments on alpha-chymotrypsin in the presence of 0.1 M glycerol were also shown to be consistent with interpretation in terms of the model of space-filling effects entailing complete exclusion of small solute from the hydrated protein domain.     

Negative second virial coefficients as predictors of protein crystal growth: evidence from sedimentation equilibrium studies that refutes the designation of those light scattering parameters as osmotic virial coefficients

The effects of ammonium sulphate concentration on the osmotic second virial coefficient (BAA/MA) for equine serum albumin (pH 5.6, 20 degrees C) have been examined by sedimentation equilibrium. After an initial steep decrease with increasing ammonium sulphate concentration, BAA/MA assumes an essentially concentration-independent magnitude of 8-9 ml/g. Such behaviour conforms with the statistical-mechanical prediction that a sufficient increase in ionic strength should effectively eliminate the contributions of charge interactions to BAA/MA but have no effect on the covolume contribution (8.4 ml/g for serum albumin). A similar situation is shown to apply to published sedimentation equilibrium data for lysozyme (pH 4.5). Although termed osmotic second virial coefficients and designated as such (B22), the negative values obtained in published light scattering studies of both systems have been described incorrectly because of the concomitant inclusion of the protein-salt contribution to thermodynamic nonideality of the protein. Those negative values are still valid predictors of conditions conducive to crystal growth inasmuch as they do reflect situations in which there is net attraction between protein molecules. However, the source of attraction responsible for the negative virial coefficient stems from the protein-salt rather than the protein-protein contribution, which is necessarily positive.     

Nonequivalence of second virial coefficients from sedimentation equilibrium and static light scattering studies of protein solutions

Experimental data for ovalbumin and lysozyme are presented to highlight the nonequivalence of second virial coefficients obtained for proteins by sedimentation equilibrium and light scattering. Theoretical considerations confirm that the quantity deduced from sedimentation equilibrium distributions is B(22), the osmotic second virial coefficient describing thermodynamic nonideality arising solely from protein self-interaction. On the other hand, the virial coefficient determined by light scattering is shown to reflect the combined contributions of protein-protein and protein-buffer interactions to thermodynamic nonideality of the protein solution. Misidentification of the light scattering parameter as B(22) accounts for published reports of negative osmotic second virial coefficients as indicators of conditions conducive to protein crystal growth. Finally, textbook assertions about the equivalence of second virial coefficients obtained by sedimentation equilibrium and light scattering reflect the restriction of consideration to single-solute systems. Although sedimentation equilibrium distributions for buffered protein solutions are, indeed, amenable to interpretation in such terms, the same situation does not apply to light scattering measurements because buffer constituents cannot be regarded as part of the solvent: instead they must be treated as non-scattering cosolutes.     

Effects of thermodynamic nonideality on protein interactions. Equivalence of interpretations based on excluded volume and preferential solvation

Published results on the stabilization of proteins by sucrose (J.C. Lee and S.N. Timasheff, J. Biol. Chem. 256 (1981) 7193) have been reexamined and interpreted in terms of thermodynamic nonideality. The composition dependence of activity coefficients may be accounted for on a statistical-mechanical basis using the concept of excluded volume. An expression is derived in which the effect of sucrose on determination of the partial specific volume of a protein, previously interpreted in terms of preferential protein solvation, is also seen to be attributable to excluded volume. Gel chromatographic studies of the reversible unfolding of alpha-chymotrypsin are presented which demonstrate temperature- and sucrose-mediated changes in the effective volume of the enzyme. These measurements support the quantitative interpretation of the stabilization in terms of thermodynamic nonideality arising from the difference between covolumes for sucrose and the two isomeric states of alpha-chymotrypsin. By establishing the equivalence of the two approaches that have been used to account for the effects of inert solutes on protein transitions, the present investigation eliminates the need for any distinction between such solutes on the basis of molecular size; and also enhances greatly the potential sensitivity of thermodynamic nonideality as a means of probing protein isomerizations, since greater displacement of the equilibrium position may be effected by small rather than by macromolecular solutes present at the same weight concentrations.     

Chromatographic evidence of the self-association of oxyhemoglobin in concentrated solutions: its biological implications.

Expressions that take into account the effects of thermodynamic non-ideality, described in terms of a high-order virial expansion, are derived for the concentration-dependence of the weight-average partition coefficient in exclusion chromatography of a single solute and of a solute undergoing reversible self-association. Comparison of the concentration-dependences predicted by those expressions with results obtained for bovine and human oxyhemoglobins on CPG-10-120 porous glass beads in 0.156 I phosphate-chloride buffer, pH 7.3, shows that neither oxyhemoglobin conforms with the concept of it being a single alpha 2 beta 2 entity with Stokes radius of 3.13 nm, the experimental value. Previously published osmotic pressure and sedimentation equilibrium results are also shown to be inconsistent with this concept. On the other hand, both sets of exclusion chromatography results are consistent with the joint operation of thermodynamic non-ideality and reversible association of the alpha 2 beta 2 species. From the magnitude of the equilibrium constant, derived for either of two possible modes of association, it is calculated that only half of the oxyhemoglobin would be in the alpha 2 beta 2 states under conditions of oxygen saturation and a concentration of 320 g/liter, that pertaining in the red blood cell. The consequences of this association phenomenon are discussed in relation to the oxygen binding curves obtained by others in the presence and absence of 2,3-diphosphoglycerate (DPG). An explanation is provided of the observed dependence on hemoglobin concentration of oxygen-binding in the presence of DPG, and of the absence of such an effect in DPG-free solutions. It is concluded that the control of oxygen binding to hemoglobin in the physiological situation involves the joint operation of self-association and allosteric effects.     

Assessment of the validity of the Adams and Fujita approximation for the higher oligomers of human spectrin

Analysing the self-association behaviour of human erythrocyte spectrin is complicated by a large degree of nonideality. Adams and Fujita [1] proposed that, as a first order approximation, the logarithm of the activity coefficient of the protomer of a self-associating system can be considered to be linearly dependent on the total concentration of the protein, and that the same second virial coefficient could be considered to apply to all species. As a consequence of the Adams and Fujita approximation, the apparent equilibrium constant is equal to the thermodynamic equilibrium constant. The equilibrium concentrations at 30°C of each oligomer spectrin species up to the 14-mer were determined after electrophoresis at low temperature. An apparent equilibrium constant for forming tetramer (K_2,4) of (1.2 ± 0.1) × 10⁶ l/mol was obtained, a value of (9.4 ± 0.7) × 10⁴ l/mol was obtained for K_4,6 and for all reactions forming oligomers higher than the hexamer an average approximate value of (2.7 ± 0.4) × 10⁵ l/mol was obtained. The apparent equilibrium constants for the formation of all oligomer species of spectrin up to the tetrakaidecamer (14-mer) remain relatively independent of total spectrin concentration, and indicate that within the precision of the measurements a single virial coefficient is sufficient to account for the nonideality of spectrin self-association over the range 2–42 g/l, thus further justifying the use of the Adams and Fujita approximation for this protein over this concentration range.     

Evaluation of nonideality from gel chromatographic partition coefficients. A technique with greater versatility than equilibrium dialysis

Frontal gel chromatography has been used to measure partition coefficients which enable a quantitative evaluation of the thermodynamic nonideality of small solutes generated by the presence of high concentrations of macromolecular solutes. Equivalence of results obtained by the present method and by equilibrium dialysis is demonstrated in a comparison of results for dextran sulfate-NaCl and dextran-sorbitol systems. Interaction coefficients obtained for dextran-sorbitol and protein-polyethylene glycol 4000 systems yields results which are in reasonable agreement with those predicted on the statistical-mechanical basis of excluded volume. Because of its greater versatility in regard to the range of systems that may be studied, the frontal gel chromatographic procedure is likely to be of particular value for the quantitative characterization of thermodynamic nonideality arising from excluded volume effects in concentrated mixtures of macromolecular solutes.     

Interactions of lens proteins. Self-association and mixed-association studies of bovine alpha-crystallin and gamma-crystallin

Concentrated solutions of calf alpha-crystallin (up to 45 g/l) and gamma-crystallin (up to 67 g/l) were subjected to frontal exclusion chromatography at pH 7.3, ionic strength 0.17 and 20 degrees C. The experimental concentration dependence of the weight-average partition coefficient was compared with theoretical expressions, which include considerations of thermodynamic non-ideality effects, for the concentration dependence of a single solute and of a solute undergoing reversible self-association. Two types of association pattern were examined, discrete dimerization and indefinite self-association. The partition chromatography results are consistent with an indefinite self-association of gamma-crystallin, governed by an isodesmic association constant of 6.7 X 10(-3) l/g. alpha-Crystallin appears to self-associate either very weakly, with a maximal association constant of 0.9 X 10(-3) l/g, or not at all; the distinction depends on the assessment of the non-ideality coefficients. The consequences of excluded volume effects on these self-association equilibria at high total protein concentration are discussed. Mixtures of alpha-crystallin and gamma-crystallin were analyzed by frontal exclusion chromatography (up to 14 g/l) and sedimentation velocity (up to 115 g/l): no interaction was observed.     

The indefinite self-association of lysozyme: consideration of composition-dependent activity coefficients

An improved iterative method for computing association constants from sedimentation equilibrium results obtained with self-interacting protein systems is presented which accounts for the composition-dependence of the activity coefficients of all oligomeric species. The method is based on the calculation of virial coefficients from covolume and charge considerations, the statistical mechanical basis of which is discussed in relation to the DLVO theory. The method is applied to results obtained with lysozyme in diethylbarbiturate buffer of pH 8.0 and ionic strength 0.15 at 15°C. It is shown that these results, encompassing a range of total solute concentration up to 19.7 g/liter are consistent with self-association patterns comprising either a monomer-dimer-trimer system or an isodesmic indefinite self-association of the monomer, the latter being favored. A firmer distinction between these possibilities is sought on the basis of the dependence of the weight-average partition coefficient, determined by frontal gel chromatography, on total solute concentration (up to 56.6 g/liter). This analysis accounts for the composition-dependence of the ratio of the activity coefficients of partitioning monomer in mobile and stationary phases. It is concluded that all results are consistent with an indefinite self-association of lysozyme governed by a single association constant of 4.61 × 10² liter/mole.     

A thermodynamic model for the self-association of human spectrin

The self-association of human spectrin at 28.8 degrees C in 0.11 M salt (pH 7.5) has been studied by means of sedimentation equilibrium. Coincidence of omega function plots as a function of total spectrin concentration (0-2 g/L) indicated that equilibrium was achieved and that no significant concentration of solute was incapable of participating in the self-association reaction. On the basis of the root-mean-square deviation of the fits and the randomness of the residuals, the behavior can be described equally well, either by a cooperative isodesmic model, in which K12 approximately 2 x 10(6) M-1 and all other K approximately 10(6) M-1, or by an attenuated scheme in which K(i-1)i approximately (3.5 x 10(6)/i M-1. The returned values of the second virial coefficient, B, for both these models fall within the range calculated from the charge and Stokes radius of spectrin. A mechanism for spectrin self-association consistent with both schemes is proposed in which spectrin heterodimers undergo a reversible opening at the self-association interface. These open heterodimers then undergo indefinite self-association to form a series of open-chain oligomers in dynamic equilibrium with closed-loop oligomers.     

A Hilly path through the thermodynamics and statistical mechanics of protein solutions

The opus of Don Winzor in the fields of physical and analytical biochemistry is a major component of that certain antipodean approach to this broad area of research that blossomed in the second half of the twentieth century. The need to formulate problems in terms of thermodynamic nonideality posed the challenge of describing a clear route from molecular interactions to the parameters that biochemists routinely measure. Mapping out this route required delving into the statistical mechanics of solutions of macromolecules, and at every turn mathematically complex, rigorous, general results that had previously been derived previously, often by Terrell Hill, came to the fore. Central to this work were the definition of the “thermodynamic activity”, the pivotal position of the polynomial expansion of the osmotic pressure in terms of molar concentration and the relationship of virial coefficients to details of the forces between limited-size groups of interacting molecules. All of this was richly exploited in the task of taking account of excluded volume and electrostatic interactions, especially in the use of sedimentation equilibrium to determine values of constants for molecular association reactions. Such an approach has proved relevant to the study of molecular interactions generally, even those between the main macromolecular solute and components of the solvent, by using techniques such as exclusion and affinity chromatography as well as light scattering.     

Average quantities accessible to the analysis of sedimentation equilibrium data for self-association systems

This paper shows that analysis of sedimentation equilibrium data, searching for average molecular weights, gives quantities which are dependent on the total protein concentration of the samples, while knowledge of the molecular weight of the monomeric subunits allows a more meaningful search for the concentrations of the individual polymeric components of the system. From these the various average molecular weights can be construed, and the various dissociation equilibrium constants evaluated. Also, in this paper considerations are proposed on the meaning of nonideality terms in associating systems and possible ways for estimating them. As an example the proposed procedures have been applied to measurements of sedimentation equilibrium in carbonmonoxyhemoglobin.     

Sedimentation equilibrium in a solution containing an arbitrary number of solute species at arbitrary concentrations: theory and application to concentrated solutions of ribonuclease

Simple expressions are derived describing the equilibrium concentration gradient of each species in a solution containing an arbitrary number of solute species at arbitrary concentration, as a function of the concentration of all species. Quantitative relationships between the species gradients and experimentally observable signal gradients are presented. The expressions are model-free and take into account both attractive and repulsive interactions between all species. In order to analyze data obtained from strongly nonideal solutions, a statistical thermodynamic model for repulsive solute-solute interactions is required. The relations obtained are utilized to analyze the dependence of the equilibrium gradient of ribonuclease A in phosphate-buffered saline, pH 7.4, upon total protein concentration. Experimental results are interpreted in the context of a model for weak self-association leading to the formation of significant amounts of oligomers at total protein concentrations exceeding 25 g/l.