A model for sedimentation in inhomogeneous media. II. Compressibility of aqueous and organic solvents

The effects of solvent compressibility on the sedimentation behavior of macromolecules as observed in analytical ultracentrifugation are examined. Expressions for the density and pressure distributions in the solution column are derived and combined with the finite element solution of the Lamm equation in inhomogeneous media to predict the macromolecular concentration distributions under different conditions. Independently, analytical expressions are derived for the sedimentation of non-diffusing particles in the limit of low compressibility. Both models are quantitatively consistent and predict solvent compressibility to result in a reduction of the sedimentation rate along the solution column and a continuous accumulation of solutes in the plateau region. For both organic and aqueous solvents, the calculated deviations from the sedimentation in incompressible media can be very large and substantially above the measurement error. Assuming conventional configurations used for sedimentation velocity experiments in analytical ultracentrifugation, neglect of the compressibility of water leads to systematic errors underestimating sedimentation coefficients by approximately 1% at a rotor speeds of 45000 rpm, but increasing to 2-5% with increasing rotor speeds and decreasing macromolecular size. The proposed finite element solution of the Lamm equation can be used to take solvent compressibility quantitatively into account in direct boundary models for discrete species, sedimentation coefficient distributions or molar mass distributions. Using the analytical expressions for the sedimentation of non-diffusing particles, the ls-g*(s) distribution of apparent sedimentation coefficients is extended to the analysis of sedimentation in compressible solvents. The consideration of solvent compressibility is highly relevant not only when using organic solvents, but also in aqueous solvents when precise sedimentation coefficients are needed, for example, for hydrodynamic modeling.     

Size-distribution analysis of macromolecules by sedimentation velocity ultracentrifugation and lamm equation modeling

A new method for the size-distribution analysis of polymers by sedimentation velocity analytical ultracentrifugation is described. It exploits the ability of Lamm equation modeling to discriminate between the spreading of the sedimentation boundary arising from sample heterogeneity and from diffusion. Finite element solutions of the Lamm equation for a large number of discrete noninteracting species are combined with maximum entropy regularization to represent a continuous size-distribution. As in the program CONTIN, the parameter governing the regularization constraint is adjusted by variance analysis to a predefined confidence level. Estimates of the partial specific volume and the frictional ratio of the macromolecules are used to calculate the diffusion coefficients, resulting in relatively high-resolution sedimentation coefficient distributions c(s) or molar mass distributions c(M). It can be applied to interference optical data that exhibit systematic noise components, and it does not require solution or solvent plateaus to be established. More details on the size-distribution can be obtained than from van Holde-Weischet analysis. The sensitivity to the values of the regularization parameter and to the shape parameters is explored with the help of simulated sedimentation data of discrete and continuous model size distributions, and by applications to experimental data of continuous and discrete protein mixtures.     

Determination of molecular parameters by fitting sedimentation data to finite-element solutions of the Lamm equation.

A method for fitting experimental sedimentation velocity data to finite-element solutions of various models based on the Lamm equation is presented. The method provides initial parameter estimates and guides the user in choosing an appropriate model for the analysis by preprocessing the data with the G(s) method by van Holde and Weischet. For a mixture of multiple solutes in a sample, the method returns the concentrations, the sedimentation (s) and diffusion coefficients (D), and thus the molecular weights (MW) for all solutes, provided the partial specific volumes (v) are known. For nonideal samples displaying concentration-dependent solution behavior, concentration dependency parameters for s(sigma) and D(delta) can be determined. The finite-element solution of the Lamm equation used for this study provides a numerical solution to the differential equation, and does not require empirically adjusted correction terms or any assumptions such as infinitely long cells. Consequently, experimental data from samples that neither clear the meniscus nor exhibit clearly defined plateau absorbances, as well as data from approach-to-equilibrium experiments, can be analyzed with this method with enhanced accuracy when compared to other available methods. The nonlinear least-squares fitting process was accomplished by the use of an adapted version of the "Doesn't Use Derivatives" nonlinear least-squares fitting routine. The effectiveness of the approach is illustrated with experimental data obtained from protein and DNA samples. Where applicable, results are compared to methods utilizing analytical solutions of approximated Lamm equations.     

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.     

A new approximate whole boundary solution of the Lamm differential equation for the analysis of sedimentation velocity experiments

Sedimentation velocity is one of the best-suited physical methods for determining the size and shape of macromolecular substances or their complexes in the range from 1 to several thousand kDa. The moving boundary in sedimentation velocity runs can be described by the Lamm differential equation. Fitting of suitable model functions or solutions of the Lamm equation to the moving boundary is used to obtain directly sedimentation and diffusion coefficients, thus allowing quick determination of size, shape and other parameters of macromolecules. Here we present a new approximate whole boundary solution of the Lamm equation that simultaneously allows the specification of sedimentation and diffusion coefficients with deviations smaller than 1% from the expected values.     

An improved function for fitting sedimentation velocity data for low-molecular-weight solutes.

Many traditional approaches to the analysis of sedimentation velocity data work poorly with data for low-molecular-weight solutes, which have sedimentation boundaries that are severely broadened by diffusion. An approach that has previously had some success is to directly fit these broad boundaries to approximate solutions of the Lamm equation that directly account for the high diffusion. However, none of the available approximate solutions work well at times both early and late in the run, or give boundary shapes that are highly accurate, especially for species of molecular weight < 10,000. An improved fitting function has been developed to overcome some of these limitations. The new function adds two correction terms to the Fujita-MacCosham solution. The optimum coefficients for these new correction terms were determined by a least-squares approach. The accuracy and limitations of fitting with this new function were tested against synthetic data sets obtained by finite-element methods, for analysis of samples containing either single species or several noninteracting species. We also compare the strengths and weaknesses of this method of analysis, and its ability to work with noisy data, relative to recently developed time-derivative methodologies.     

Determination of sedimentation coefficients for small peptides.

Direct fitting of sedimentation velocity data with numerical solutions of the Lamm equations has been exploited to obtain sedimentation coefficients for single solutes under conditions where solvent and solution plateaus are either not available or are transient. The calculated evolution was initialized with the first experimental scan and nonlinear regression was employed to obtain best-fit values for the sedimentation and diffusion coefficients. General properties of the Lamm equations as data analysis tools were examined. This method was applied to study a set of small peptides containing amphipathic heptad repeats with the general structure Ac-YS-(AKEAAKE)nGAR-NH2, n = 2, 3, or 4. Sedimentation velocity analysis indicated single sedimenting species with sedimentation coefficients (s(20,w) values) of 0.37, 0.45, and 0.52 S, respectively, in good agreement with sedimentation coefficients predicted by hydrodynamic theory. The described approach can be applied to synthetic boundary and conventional loading experiments, and can be extended to analyze sedimentation data for both large and small macromolecules in order to define shape, heterogeneity, and state of association.     

Identification and interpretation of complexity in sedimentation velocity boundaries.

Synthetic sedimentation velocity boundaries were generated using finite-element solutions to the original and modified forms of the Lamm equation. Situations modeled included ideal single- and multicomponent samples, concentration-dependent samples, noninteracting multicomponent samples, and reversibly self-associating samples. Synthetic boundaries subsequently were analyzed using the method of van Holde and Weischet, and results were compared against known input parameters. Results indicate that this analytical method provides rigorous diagnostics for virtually every type of sample complexity encountered experimentally. Accordingly, both the power and utility of sedimentation velocity experiments have been significantly expanded.     

Analyses of sedimentation equilibrium data

A numerical procedure is presented which can quite adequately compute the molecular weight averages as a function of solute concentration from sedimentation equilibrium data for homogeneous systems and for monomer-dimer associating systems with a possible extension to heterogeneous systems where monotonic variation in the weight average molecular weight is observed such as in weakly associating or dissociating systems. The procedure utilizes the method of orthogonal polynomials for curve fitting which allows for a rapid determination of best fit with minimal round off error. The procedure is particularly applicable in cases where the concentration of solute at the meniscus can be considered to be neither appreciable and reasonably well determined as in low speed sedimentation equilibrium experiments, nor essentially zero as in high speed sedimentation equilibrium experiments where the calculations become somewhat more simplified. The use of moderate speed sedimentation equilibrium has the advantage of providing a more broad concentration distribution in the centrifuge cell which yields more extensive information concerning dissociating systems yet still provides results at low solute concentrations where most solutes can be considered to be behaving ideally.     

Self-association studies on the EphB2 receptor SAM domain using analytical ultracentrifugation

The self-association behavior of the Eph-kinases SAM domain has been studied in phosphate buffer, pH 7.4, containing 0.14 M NaCl using concentration-dependent sedimentation equilibrium experiments. Only weak interactions typical for a monomer-dimer equilibrium up to at least 12 mg/mL were observed. Such concentrated solutions require a consideration of the non-ideality expressed by virial coefficients. A special centrifuge equation was used for the global analysis to estimate equilibrium constants based on the thermodynamic activities of the reactants. When neglecting this, the parameters deviate by about 20%. Association constants for dimerization of the EphB2-SAM domain vary between 163 M(-1) at 10 degrees C and 395 M(-1) at 32 degrees C, indicating hydrophobic forces are involved in the dimerization process. In solutions of about 12 mg/mL, less than 50% dimers are in solution and higher oligomers can be excluded.     

Molecular weights from approach-to-sedimentation equilibrium data using nonlinear regression analysis

Data obtained from early times during the transient period of sedimentation equilibrium experiments are analyzed using an approximate solution to the Lamm equation to estimate s/D. The Cr versus r data obtained at several times during approach-to-equilibrium are analyzed using a nonlinear least squares algorithm and Fujita's approximate solution. This procedure was tested using D-Ser13-somatostatin, ribonuclease, and ovalbumin. The results obtained demonstrate that for monodisperse samples s/D may be rapidly and reliably estimated using this method.     

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.     

Band centrifugation of macromolecules in self-generating density gradients. II. Sedimentation and diffusion of macromolecules in bands

Three approaches to the simultaneous sedimentation and diffusion of hands or zones of noninteracting homogeneous macromolecules are examined: (1) The authors' method of moments: (2) the transport me of Sehumaker and Rosenbloom; and (3) the stochastic solution of the Lamm equation due to Gehatia and Katehalski. All three methods indicate that the motion of the maximum of the hand may be used to evaluate the sedimentation coefficient. The moment, method provides relations which appear to be useful for measuring diffusion coefficients. Relations are given for the analysis of resolved components. The problem of measuring sedimentation coefficients of macromolecules with concentration-dependent sedimentation coefficients is examined. Methods are described for evaluating the sedimentation coefficient in these systems and for obtaining the sedimentation coefficient at infinite dilution. Methods are described for determining the weight-average sedimentation coefficient in Multi-component systems, and the differential and integral distribution of sedimentation coefficients of macromolecules with low-diffusion coefficients.     

Simultaneous determination of the individual concentration gradients of two solute species in a centrifuged mixture: application to analytical ultracentrifugation

A procedure for determining the absolute activity of 14C-labeled and 3H-labeled solutes in a mixture from the measured counts per minute in two scintillation energy windows is described. It is shown that the method described here provides a substantially more accurate determination of 3H activity in the presence of a larger 14C activity, and a more accurate determination of 14C activity in the presence of a larger 3H activity, than does the standard dual label analysis implemented in a Beckman LS 3801 scintillation counter. The new dual label procedure is combined with the automated fractionation procedure of Attri and Minton [(1986) Anal. Biochem. 152, 319-328] to permit the gradients of each of two differently radiolabeled solute species in a mixture to be individually determined following centrifugation. It is shown that the sedimentation coefficients of each of two differently labeled noninteracting proteins in a mixture may be readily determined in a sedimentation velocity experiment, and that the molecular weights of each of two such proteins in a mixture may be readily determined in a sedimentation equilibrium experiment.     

An approximate solution to the Lamm equation

An approximate solution to the Lamm equation subject to the initial and boundary conditions for conventional sedimentation velocity experiments is derived and compared with the approximate solution of Fujita and MacCosham. Calculations with this solution demonstrate that the half-height method of estimating sedimentation coefficients yields correct values for epsilon < 0.02.     

Non-ideality by sedimentation velocity of halophilic malate dehydrogenase in complex solvents

We have investigated the potential of sedimentation velocity analytical ultracentrifugation for the measurement of the second virial coefficients of proteins, with the goal of developing a method that allows efficient screening of different solvent conditions. This may be useful for the study of protein crystallization. Macromolecular concentration distributions were modeled using the Lamm equation with the approximation of linear concentration dependencies of the diffusion constant, D = D(o) (1 + k(D)c), and the reciprocal sedimentation coefficient s = s(o)/(1 + k(s)c). We have studied model distributions for their information content with respect to the particle and its non-ideal behavior, developed a strategy for their analysis by direct boundary modeling, and applied it to data from sedimentation velocity experiments on halophilic malate dehydrogenase in complex aqueous solvents containing sodium chloride and 2-methyl-2,4-pentanediol, including conditions near phase separation. Using global modeling for three sets of data obtained at three different protein concentrations, very good estimates for k(s) and s degrees and also for D degrees and the buoyant molar mass were obtained. It was also possible to obtain good estimates for k(D) and the second virial coefficients. Modeling of sedimentation velocity profiles with the non-ideal Lamm equation appears as a good technique to investigate weak inter-particle interactions in complex solvents and also to extrapolate the ideal behavior of the particle.     

Analytical ultracentrifugation as a tool for studying protein interactions

The last two decades have led to significant progress in the field of analytical ultracentrifugation driven by instrumental, theoretical, and computational methods. This review will highlight key developments in sedimentation equilibrium (SE) and sedimentation velocity (SV) analysis. For SE, this includes the analysis of tracer sedimentation equilibrium at high concentrations with strong thermodynamic non-ideality, and for ideally interacting systems, the development of strategies for the analysis of heterogeneous interactions towards global multi-signal and multi-speed SE analysis with implicit mass conservation. For SV, this includes the development and applications of numerical solutions of the Lamm equation, noise decomposition techniques enabling direct boundary fitting, diffusion deconvoluted sedimentation coefficient distributions, and multi-signal sedimentation coefficient distributions. Recently, effective particle theory has uncovered simple physical rules for the co-migration of rapidly exchanging systems of interacting components in SV. This has opened new possibilities for the robust interpretation of the boundary patterns of heterogeneous interacting systems. Together, these SE and SV techniques have led to new approaches to study macromolecular interactions across the entire spectrum of affinities, including both attractive and repulsive interactions, in both dilute and highly concentrated solutions, which can be applied to single-component solutions of self-associating proteins as well as the study of multi-protein complex formation in multi-component solutions.     

The determination of density and molecular weight distributions of lipoproteins by sedimentation equilibrium.

A method is presented by which an experimental record of total concentration as a function of radial distance, obtained in a sedimentation equilibrium experiment conducted with a noninteracting mixture in the absence of a density gradient, may be analyzed to obtain the unimodal distributions of molecular weight and of partial molar volume when these vary concomitantly and continuously. Particular attention is given to the caracterization of classes of lipoproteins exhibiting Gaussian distributions of these quantities, although the analysis is applicable to other types of unimodal distribution. Equations are also formulated permitting the definition of the corresponding distributions of partial specific volume and of density. The analysis procedure is based on a method (employing Laplace transforms) developed previously, but differs from it in that it avoids the necessity of differentiating experimental results, which introduces error. The method offers certain advantages over other procedures used to characterize and compare lipoprotein samples (exhibiting unimodal distributions) with regard to the duration of the experiment, economy of the sample, and, particularly, the ability to define in principle all of the relevant distributions from one sedimentation equilibrium experiment and an external measurement of the weight average partial specific volume. These points and the steps in the analysis procedure are illustrated with experimental results obtained in the sedimentation equilibrium of a sample of human serum low density lipoprotein. The experimental parameters (such as solution density, column height, and angular velocity) used in the conduction of these experiments were selected on the basis of computer-simulated examples, which are also presented. These provide a guide for other workers interested in characterizing lipoproteins of this class.     

Sedimentation velocity analysis of heterogeneous protein-protein interactions: Lamm equation modeling and sedimentation coefficient distributions c(s)

We describe algorithms for solving the Lamm equations for the reaction-diffusion-sedimentation process in analytical ultracentrifugation, and examine the potential and limitations for fitting experimental data. The theoretical limiting case of a small, uniformly distributed ligand rapidly reacting with a larger protein in a "constant bath" of the ligand is recapitulated, which predicts the reaction boundary to sediment with a single sedimentation and diffusion coefficient. As a consequence, it is possible to express the sedimentation profiles of reacting systems as c(s) distribution of noninteracting Lamm equation solutions, deconvoluting the effects of diffusion. For rapid reactions, the results are quantitatively consistent with the "constant bath" approximation, showing c(s) peaks at concentration-dependent positions. For slower reactions, the deconvolution of diffusion is still partially successful, with c(s) resolving peaks that reflect the populations of sedimenting species. The transition between c(s) peaks describing reaction boundaries of moderately strong interactions (K(D) approximately 10(-6) M) or resolving sedimenting species was found to occur in a narrow range of dissociation rate constant between 10(-3) and 10(-4) s(-1). The integration of the c(s) peaks can lead to isotherms of species populations or s-value of the reaction boundary, respectively, which can be used for the determination of the equilibrium binding constant.     

Plant annexins form calcium-independent oligomers in solution

The oligomeric state in solution of four plant annexins, namely Anx23(Ca38), Anx24(Ca32), Anx(Gh1), and Anx(Gh2), was characterized by sedimentation equilibrium analysis and gel filtration. All proteins were expressed and purified as amino-terminal His(n) fusions. Sequencing of the Anx(Gh1) construct revealed distinct differences with the published sequence. Sedimentation equilibrium analysis of Anx23(Ca38), Anx24(Ca32), and Anx(Gh1) suggests monomer-trimer equilibria for each protein with association constants in the range of 0.9 x 10(10)-1.7 x 10(11) M(-2). All four proteins were subjected to analytical gel filtration under different buffer conditions. Observations from this experiment series agree quantitatively with the ultracentrifugation results, and strongly suggest calcium independence of the annexin oligomerization behavior; moreover, binding of calcium ions to the proteins seems to require disassembly of the oligomers. Anx(Gh2) showed a different elution profile than the other plant annexins; while having only a very small trimer content, this annexin seems to exist in a monomer-dimer equilibrium in solution.