The analysis of macromolecular interactions by sedimentation equilibrium

The study of macromolecular interactions by sedimentation equilibrium is a highly technical method that requires great care in both the experimental design and data analysis. The complexity of the interacting system that can be analyzed is only limited by the ability to deconvolute the exponential contributions of each of the species to the overall concentration gradient. This is achieved in part through the use of multi-signal data collection and the implementation of soft mass conservation. We illustrate the use of these constraints in SEDPHAT through the study of an A+B+B?AB+B?ABB system and highlight some of the technical challenges that arise. We show that both the multi-signal analysis and mass conservation result in a precise and robust data analysis and discuss improvements that can be obtained through the inclusion of data from other methods such as sedimentation velocity and isothermal titration calorimetry.     

Modern analytical ultracentrifugation in protein science: a tutorial review

Analytical ultracentrifugation (AU) is reemerging as a versatile tool for the study of proteins. Monitoring the sedimentation of macromolecules in the centrifugal field allows their hydrodynamic and thermodynamic characterization in solution, without interaction with any matrix or surface. The combination of new instrumentation and powerful computational software for data analysis has led to major advances in the characterization of proteins and protein complexes. The pace of new advancements makes it difficult for protein scientists to gain sufficient expertise to apply modern AU to their research problems. To address this problem, this review builds from the basic concepts to advanced approaches for the characterization of protein systems, and key computational and internet resources are provided. We will first explore the characterization of proteins by sedimentation velocity (SV). Determination of sedimentation coefficients allows for the modeling of the hydrodynamic shape of proteins and protein complexes. The computational treatment of SV data to resolve sedimenting components has been achieved. Hence, SV can be very useful in the identification of the oligomeric state and the stoichiometry of heterogeneous interactions. The second major part of the review covers sedimentation equilibrium (SE) of proteins, including membrane proteins and glycoproteins. This is the method of choice for molar mass determinations and the study of self-association and heterogeneous interactions, such as protein-protein, protein-nucleic acid, and protein-small molecule binding.     

Using Lamm-Equation modeling of sedimentation velocity data to determine the kinetic and thermodynamic properties of macromolecular interactions

The interaction of macromolecules with themselves and with other macromolecules is fundamental to the functioning of living systems. Recent advances in the analysis of sedimentation velocity (SV) data obtained by analytical ultracentrifugation allow the experimenter to determine important features of such interactions, including the equilibrium association constant and information about the kinetic off-rate of the interaction. The determination of these parameters is made possible by the ability of modern software to fit numerical solutions of the Lamm Equation with kinetic considerations directly to SV data. Herein, the SV analytical advances implemented in the software package SEDPHAT are summarized. Detailed analyses of SV data using these strategies are presented. Finally, a few highlights of recent literature reports that feature this type of SV data analysis are surveyed.     

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.     

The use of analytical sedimentation velocity to extract thermodynamic linkage

For 25 years, the Gibbs Conference on Biothermodynamics has focused on the use of thermodynamics to extract information about the mechanism and regulation of biological processes. This includes the determination of equilibrium constants for macromolecular interactions by high precision physical measurements. These approaches further reveal thermodynamic linkages to ligand binding events. Analytical ultracentrifugation has been a fundamental technique in the determination of macromolecular reaction stoichiometry and energetics for 85 years. This approach is highly amenable to the extraction of thermodynamic couplings to small molecule binding in the overall reaction pathway. In the 1980s this approach was extended to the use of sedimentation velocity techniques, primarily by the analysis of tubulin-drug interactions by Na and Timasheff. This transport method necessarily incorporates the complexity of both hydrodynamic and thermodynamic nonideality. The advent of modern computational methods in the last 20 years has subsequently made the analysis of sedimentation velocity data for interacting systems more robust and rigorous. Here we review three examples where sedimentation velocity has been useful at extracting thermodynamic information about reaction stoichiometry and energetics. Approaches to extract linkage to small molecule binding and the influence of hydrodynamic nonideality are emphasized. These methods are shown to also apply to the collection of fluorescence data with the new Aviv FDS.     

Beyond the second virial coefficient: Sedimentation equilibrium in highly non-ideal solutions

The general theory of sedimentation equilibrium (SE), applicable to mixtures of interacting sedimentable solutes at arbitrary concentration, is summarized. Practical techniques for the acquisition of SE data suitable for analysis are described. Experimental measurements and analyses of SE in concentrated protein solutions are reviewed. The method of non-ideal tracer sedimentation equilibrium (NITSE) is described. Experimental studies using NITSE to detect and quantitatively characterize intermolecular interactions in mixtures of dilute tracer species and concentrated proteins or polymers are reviewed.     

Diffusion of the Reaction Boundary of Rapidly Interacting Macromolecules in Sedimentation Velocity

Sedimentation velocity analytical ultracentrifugation combines relatively high hydrodynamic resolution of macromolecular species with the ability to study macromolecular interactions, which has great potential for studying dynamically assembled multiprotein complexes. Complicated sedimentation boundary shapes appear in multicomponent mixtures when the timescale of the chemical reaction is short relative to the timescale of sedimentation. Although the Lamm partial differential equation rigorously predicts the evolution of concentration profiles for given reaction schemes and parameter sets, this approach is often not directly applicable to data analysis due to experimental and sample imperfections, and/or due to unknown reaction pathways. Recently, we have introduced the effective particle theory, which explains quantitatively and in a simple physical picture the sedimentation boundary patterns arising in the sedimentation of rapidly interacting systems. However, it does not address the diffusional spread of the reaction boundary from the cosedimentation of interacting macromolecules, which also has been of long-standing interest in the theory of sedimentation velocity analytical ultracentrifugation. Here, effective particle theory is exploited to approximate the concentration gradients during the sedimentation process, and to predict the overall, gradient-average diffusion coefficient of the reaction boundary. The analysis of the heterogeneity of the sedimentation and diffusion coefficients across the reaction boundary shows that both are relatively uniform. These results support the application of diffusion-deconvoluting sedimentation coefficient distributions c(s) to the analysis of rapidly interacting systems, and provide a framework for the quantitative interpretation of the diffusional broadening and the apparent molar mass values of the effective sedimenting particle in dynamically associating systems.     

On the analysis of protein self-association by sedimentation velocity analytical ultracentrifugation

Analytical ultracentrifugation is one of the classical techniques for the study of protein interactions and protein self-association. Recent instrumental and computational developments have significantly enhanced this methodology. In this paper, new tools for the analysis of protein self-association by sedimentation velocity are developed, their statistical properties are examined, and considerations for optimal experimental design are discussed. A traditional strategy is the analysis of the isotherm of weight-average sedimentation coefficients s(w) as a function of protein concentration. From theoretical considerations, it is shown that integration of any differential sedimentation coefficient distribution c(s), ls-g(*)(s), or g(s(*)) can give a thermodynamically well-defined isotherm, as long as it provides a good model for the sedimentation profiles. To test this condition for the g(s(*)) distribution, a back-transform into the original data space is proposed. Deconvoluting diffusion in the sedimentation coefficient distribution c(s) can be advantageous to identify species that do not participate in the association. Because of the large number of scans that can be analyzed in the c(s) approach, its s(w) values are very precise and allow extension of the isotherm to very low concentrations. For all differential sedimentation coefficients, corrections are derived for the slowing of the sedimentation boundaries caused by radial dilution. As an alternative to the interpretation of the isotherm of the weight-average s value, direct global modeling of several sedimentation experiments with Lamm equation solutions was studied. For this purpose, a new software SEDPHAT is introduced, allowing the global analysis of several sedimentation velocity and equilibrium experiments. In this approach, information from the shape of the sedimentation profiles is exploited, which permits the identification of the association scheme and requires fewer experiments to precisely characterize the association. Further, under suitable conditions, fractions of incompetent material that are not part of the reversible equilibrium can be detected.     

Sedimentation velocity analysis of heterogeneous protein-protein interactions: sedimentation coefficient distributions c(s) and asymptotic boundary profiles from Gilbert-Jenkins theory

Interacting proteins in rapid association equilibrium exhibit coupled migration under the influence of an external force. In sedimentation, two-component systems can exhibit bimodal boundaries, consisting of the undisturbed sedimentation of a fraction of the population of one component, and the coupled sedimentation of a mixture of both free and complex species in the reaction boundary. For the theoretical limit of diffusion-free sedimentation after infinite time, the shapes of the reaction boundaries and the sedimentation velocity gradients have been predicted by Gilbert and Jenkins. We compare these asymptotic gradients with sedimentation coefficient distributions, c(s), extracted from experimental sedimentation profiles by direct modeling with superpositions of Lamm equation solutions. The overall shapes are qualitatively consistent and the amplitudes and weight-average s-values of the different boundary components are quantitatively in good agreement. We propose that the concentration dependence of the area and weight-average s-value of the c(s) peaks can be modeled by isotherms based on Gilbert-Jenkins theory, providing a robust approach to exploit the bimodal structure of the reaction boundary for the analysis of experimental data. This can significantly improve the estimates for the determination of binding constants and hydrodynamic parameters of the complexes.     

Characterization of weak protein dimerization by direct analysis of sedimentation equilibrium distributions: the INVEQ approach

Closer scrutiny has been accorded a recently reported procedure for characterizing weak protein dimerization by sedimentation equilibrium (INVEQ) in which the equilibrium distribution is analyzed as a dependence of radial distance on solute concentration rather than of solute concentration on radial distance. By demonstrating theoretically that the fundamental parameter derived from the analysis is simply the difference between the dimerization constant and the osmotic second virial coefficient for monomer-monomer interaction, this investigation refutes the original claim that independent estimates of these two parameters can be obtained by nonlinear curve fitting of the sedimentation equilibrium distribution. This criticism also applies to conventional analyses of sedimentation distributions by the commonly employed Beckman Origin and NONLIN software. Numerically simulated distributions are then analyzed to demonstrate limitations of the procedure and also to indicate a means of improving the reliability of the returned estimate of the dimerization constant. These features are illustrated by applying the original and revised analytical procedures to a sedimentation equilibrium distribution for alpha-chymotrypsin (pH 4.0, I 0.05 M).     

Analytical ultracentrifugation for the study of protein association and assembly

Analytical ultracentrifugation remains pre-eminent among the methods used to study the interactions of macromolecules under physiological conditions. Recent developments in analytical procedures allow the high resolving power of sedimentation velocity methods to be coupled to sedimentation equilibrium approaches and applied to both static and dynamic associations. Improvements in global modeling based on numerical solutions of the Lamm equation have generated new sedimentation velocity applications with an emphasis on data interpretation using sedimentation coefficient or molar mass distributions. Procedures based on the use of multiple optical signals from absorption and interference optics for the analysis of the sedimentation velocity and equilibrium behavior of more complex interactions have now been developed. New applications of tracer sedimentation equilibrium experiments and the development of a fluorescence optical system for the analytical ultracentrifuge extend the accessible concentration range over several orders of magnitude and, coupled with the new analytical procedures, provide powerful new tools for studies of both weak and strong macromolecular interactions in solution.     

Variable-Field Analytical Ultracentrifugation: I. Time-Optimized Sedimentation Equilibrium

Sedimentation equilibrium (SE) analytical ultracentrifugation (AUC) is a gold standard for the rigorous determination of macromolecular buoyant molar masses and the thermodynamic study of reversible interactions in solution. A significant experimental drawback is the long time required to attain SE, which is usually on the order of days. We have developed a method for time-optimized SE (toSE) with defined time-varying centrifugal fields that allow SE to be attained in a significantly (up to 10-fold) shorter time than is usually required. To achieve this, numerical Lamm equation solutions for sedimentation in time-varying fields are computed based on initial estimates of macromolecular transport properties. A parameterized rotor-speed schedule is optimized with the goal of achieving a minimal time to equilibrium while limiting transient sample preconcentration at the base of the solution column. The resulting rotor-speed schedule may include multiple over- and underspeeding phases, balancing the formation of gradients from strong sedimentation fluxes with periods of high diffusional transport. The computation is carried out in a new software program called TOSE, which also facilitates convenient experimental implementation. Further, we extend AUC data analysis to sedimentation processes in such time-varying centrifugal fields. Due to the initially high centrifugal fields in toSE and the resulting strong migration, it is possible to extract sedimentation coefficient distributions from the early data. This can provide better estimates of the size of macromolecular complexes and report on sample homogeneity early on, which may be used to further refine the prediction of the rotor-speed schedule. In this manner, the toSE experiment can be adapted in real time to the system under study, maximizing both the information content and the time efficiency of SE experiments.     

Analysis of high-affinity assembly for AMPA receptor amino-terminal domains

Analytical ultracentrifugation (AUC) and steady-state fluorescence anisotropy were used to measure the equilibrium dissociation constant (Kd) for formation of dimers by the amino-terminal domains (ATDs) of the GluA2 and GluA3 subtypes of AMPA receptor. Previous reports on GluA2 dimerization differed in their estimate of the monomer-dimer Kd by a 2,400-fold range, with no consensus on whether the ATD forms tetramers in solution. We find by sedimentation velocity (SV) analysis performed using absorbance detection a narrow range of monomer-dimer Kd values for GluA2, from 5 to 11 nM for six independent experiments, with no detectable formation of tetramers and no effect of glycosylation or the polypeptide linker connecting the ATD and ligand-binding domains; for GluA3, the monomer-dimer Kd was 5.6 μM, again with no detectable tetramer formation. For sedimentation equilibrium (SE) experiments, a wide range of Kd values was obtained for GluA2, from 13 to 284 nM, whereas for GluA3, the Kd of 3.1 μM was less than twofold different from the SV value. Analysis of cell contents after the ～1-week centrifuge run by silver-stained gels revealed low molecular weight GluA2 breakdown products. Simulated data for SE runs demonstrate that the apparent Kd for GluA2 varies with the extent of proteolysis, leading to artificially high Kd values. SV experiments with fluorescence detection for GluA2 labeled with 5,6-carboxyfluorescein, and fluorescence anisotropy measurements for GluA2 labeled with DyLight405, yielded Kd values of 5 and 11 nM, consistent with those from SV with absorbance detection. However, the sedimentation coefficients measured by AUC using absorbance and fluorescence systems were strikingly different, and for the latter are not consistent with hydrodynamic protein models. Thus, for unknown reasons, the concentration dependence of sedimentation coefficients obtained with fluorescence detection SV may be unreliable, limiting the usefulness of this technique for quantitative analysis.     

Evaluation of diffusion coefficients by means of an approximate steady-state condition in sedimentation velocity distributions

This investigation examined the feasibility of manipulating the rotor speed in sedimentation velocity experiments to spontaneously generate an approximate steady-state condition where the extent of diffusional spreading is matched exactly by the boundary sharpening arising from negative s–c dependence. Simulated sedimentation velocity distributions based on the sedimentation characteristics for a purified mucin preparation were used to illustrate a simple procedure for determining the diffusion coefficient from such steady-state distributions in situations where the concentration dependence of the sedimentation coefficient, s = s⁰/(1 + Kc), was quantified in terms of the limiting sedimentation coefficient as c → 0 (s⁰) and the concentration coefficient (K). Those simulations established that spontaneous generation of the approximate steady state could well be a feature of sedimentation velocity distributions for many unstructured polymer systems because the requirement that Kc_oω²s⁰/D be between 46 and 183 cm⁻² is not unduly restrictive. Although spontaneous generation of the approximate steady state is also a theoretical prediction for structured macromolecular solutes exhibiting linear concentration dependence of the sedimentation coefficient, s = s⁰(1 − kc), the required value of k is far too large for any practical advantage to be taken of this approach with globular proteins.     

MultiSig: a new high-precision approach to the analysis of complex biomolecular systems

MultiSig is a newly developed mode of analysis of sedimentation equilibrium (SE) experiments in the analytical ultracentrifuge, having the capability of taking advantage of the remarkable precision (~0.1 % of signal) of the principal optical (fringe) system employed, thus supplanting existing methods of analysis through reducing the ‘noise’ level of certain important parameter estimates by up to orders of magnitude. Long-known limitations of the SE method, arising from lack of knowledge of the true fringe number in fringe optics and from the use of unstable numerical algorithms such as numerical differentiation, have been transcended. An approach to data analysis, akin to ‘spatial filtering’, has been developed, and shown by both simulation and practical application to be a powerful aid to the precision with which near-monodisperse systems can be analysed, potentially yielding information on protein-solvent interaction. For oligo- and poly-disperse systems the information returned includes precise average mass distributions over both cell radial and concentration ranges and mass-frequency histograms at fixed radial positions. The application of MultiSig analysis to various complex heterogenous systems and potentially multiply-interacting carbohydrate oligomers is described.     

Extended Fujita approach to the molecular weight distribution of polysaccharides and other polymeric systems

In 1962 H. Fujita (H. Fujita, Mathematical Theory of Sedimentation Analysis, Academic Press, New York, 1962) examined the possibility of transforming a quasi-continuous distribution g(s) of sedimentation coefficient s into a distribution f(M) of molecular weight M for linear polymers using the relation f(M)=g(s)·(ds/dM) and showed that this could be done if information about the relation between s and M is available from other sources. Fujita provided the transformation based on the scaling relation s=κ(s)M(0.5), where κ(s) is taken as a constant for that particular polymer and the exponent 0.5 essentially corresponds to a randomly coiled polymer under ideal conditions. This method has been successfully applied to mucus glycoproteins (S.E. Harding, Adv. Carbohyd. Chem. Biochem. 47 (1989) 345-381). We now describe an extension of the method to general conformation types via the scaling relation s=κM(b), where b=0.4-0.5 for a coil, ～0.15-0.2 for a rod and ～0.67 for a sphere. We give examples of distributions f(M) versus M obtained for polysaccharides from SEDFIT derived least squares g(s) versus s profiles (P. Schuck, Biophys. J. 78 (2000) 1606-1619) and the analytical derivative for ds/dM performed with Microcal ORIGIN. We also describe a more direct route from a direct numerical solution of the integral equation describing the molecular weight distribution problem. Both routes give identical distributions although the latter offers the advantage of being incorporated completely within SEDFIT. The method currently assumes that solutions behave ideally: sedimentation velocity has the major advantage over sedimentation equilibrium in that concentrations less than 0.2mg/ml can be employed, and for many systems non-ideality effects can be reasonably ignored. For large, non-globular polymer systems, diffusive contributions are also likely to be small.     

Sedimentation velocity analysis of highly heterogeneous systems

This article discusses several improvements to the van Holde-Weischet (vHW) method [Biopolymers 17 (1978) 1387] that address its capability to deal with sedimentation coefficient distributions spanning a large range of s values. The method presented here allows the inclusion of scans early and late in the experiment that ordinarily would need to be excluded from the analysis due to ultracentrifuge cell end effects. Scans late in the experiment are compromised by the loss of a defined plateau region and by back-diffusion from the bottom of the cell. Early scans involve partial boundaries that have not fully cleared the meniscus. In addition, a major refinement of the algorithm for determining the boundary fractions is introduced, taking into account different degrees of radial dilution for different species in the system. The method retains its desirable model-independent properties (the analysis of sedimentation data does not require prior knowledge of a user-imposed model or range of sedimentation coefficients) and reports diffusion-corrected s value distributions, which can be presented either in a histogram format or the traditional integral distribution format. Data analyzed with the traditional vHW method are compared with those of the improved method to demonstrate the benefit from the added information in the analysis.     

Analysis of protein self-association under conditions of the thermodynamic non-ideality

Analysis of protein-protein interactions in highly concentrated solutions requires a consideration of the non-ideality in such solutions which is expressed by the virial coefficients. Different equations are presented to estimate effects of the thermodynamic non-ideality on the macromolecular interaction of self-associating proteins in sedimentation equilibrium experiments. Usually the influence of thermodynamic non-ideal behavior are described by concentration power series. The convergence of such power series is limited at high solute concentration. When expressing the thermodynamic non-ideality by an activity power series this disadvantage can be minimized. The developed centrifuge equations are the basis for a global analysis to estimate equilibrium constants and the corresponding thermodynamic activities of the reactants. Based on fit analysis of synthetic concentration profiles it was established that marked deviations from the expected association constants are observed for proteins with strong association forces between solute molecules. Considerable differences were also observed in weakly interacting systems. This was due to the excluded volume of the protein which is similar in magnitude to the binding constant. For interactions with moderate affinities values extremely close to the true binding values were obtained, as confirmed by experimental results with concanavalin A.