Determination of the sedimentation coefficient distribution by least-squares boundary modeling

A new method is presented for the calculation of apparent sedimentation coefficient distributions g*(s) for the size-distribution analysis of polymers in sedimentation velocity experiments. Direct linear least-squares boundary modeling by a superposition of sedimentation profiles of ideal nondiffusing particles is employed. It can be combined with algebraic noise decomposition techniques for the application to interference optical ultracentrifuge data at low loading concentrations with significant systematic noise components. Because of the use of direct boundary modeling, residuals are available for assessment of the quality of the fits and the consistency of the g*(s) distribution with the experimental data. The method can be combined with regularization techniques based on F statistics, such as used in the program CONTIN, or alternatively, the increment of s values can be adjusted empirically. The method is simple, has advantageous statistical properties, and reveals precise sedimentation coefficients. The new least-squares ls-g*(s) exhibits a very high robustness and resolution if data acquired over a large time interval are analyzed. This can result in a high resolution for large particles, and for samples with a high degree of heterogeneity. Because the method does not require a high frequency of scans, it can also be easily used in experiments with the absorbance optical scanning system. Published 2000 John Wiley & Sons, Inc.     

Tools for the Quantitative Analysis of Sedimentation Boundaries Detected by Fluorescence Optical Analytical Ultracentrifugation

Fluorescence optical detection in sedimentation velocity analytical ultracentrifugation allows the study of macromolecules at nanomolar concentrations and below. This has significant promise, for example, for the study of systems of high-affinity protein interactions. Here we describe adaptations of the direct boundary modeling analysis approach implemented in the software SEDFIT that were developed to accommodate unique characteristics of the confocal fluorescence detection system. These include spatial gradients of signal intensity due to scanner movements out of the plane of rotation, temporal intensity drifts due to instability of the laser and fluorophores, and masking of the finite excitation and detection cone by the sample holder. In an extensive series of experiments with enhanced green fluorescent protein ranging from low nanomolar to low micromolar concentrations, we show that the experimental data provide sufficient information to determine the parameters required for first-order approximation of the impact of these effects on the recorded data. Systematic deviations of fluorescence optical sedimentation velocity data analyzed using conventional sedimentation models developed for absorbance and interference optics are largely removed after these adaptations, resulting in excellent fits that highlight the high precision of fluorescence sedimentation velocity data, thus allowing a more detailed quantitative interpretation of the signal boundaries that is otherwise not possible for this system.     

Sedimentation equilibrium analysis of interference optical data by systematic noise decomposition.

An analytical method is described for removal of systematic signal offsets from interference optical data of sedimentation equilibrium gradients. It is demonstrated that the time-invariant signal contributions can be extracted from hydrodynamic modeling of interference profiles acquired during the approach to sedimentation equilibrium. This method is based on a technique for the explicit algebraic calculation of time-invariant noise components from sedimentation data, recently described for the direct modeling of sedimentation velocity experiments (P. Schuck and B. Demeler, Biophys. J. 76, 2288-2296, 1999). The calculated systematic signal offset is very well defined by the experimental data, stable over time, and its calculation is robust and to a large extent independent of the hydrodynamic model. The calculated time-invariant signal can be used to reduce the systematic errors in the measured sedimentation equilibrium profiles by more than an order of magnitude. It is shown that the resulting net equilibrium fringe profiles after subtraction of the time-invariant noise component allow equilibrium analyses consistent with those obtained from absorbance profiles. However, due to a higher dynamic range and the higher number of data points, the parameters derived from the net interference analysis can exhibit significantly improved precision. The presented study demonstrates the feasibility and potential of this analytical method for full exploitation of the remarkable precision of the interference optical data acquisition system, allowing sedimentation equilibrium experiments at loading concentrations below 0.05 mg/ml.     

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.     

The effective time of centrifugation for the analysis of boundary spreading in sedimentation velocity experiments

This investigation establishes a likely order of magnitude for the zero-time correction factor governing the effective time of centrifugation that is pertinent in the analysis of boundary spreading in sedimentation velocity experiments. This correction is shown to be too small to unduly affect the magnitudes of sedimentation and diffusion coefficients deduced from the application of computer software incorporating the printout value of ω2t and an effective position of the air-solution meniscus that is obtained as an additional parameter in the analysis involving nonlinear least-squares curve-fitting of sedimentation velocity distributions to the Lamm equation. Although this procedure slightly underestimates the actual meniscus position (r(a)), uncertainty about its exact location precludes adoption of the alternative approach with r(a) fixed and the correction to ω2t regarded as the additional curve-fitting parameter.     

Size-distribution analysis of proteins by analytical ultracentrifugation: strategies and application to model systems.

Strategies for the deconvolution of diffusion in the determination of size-distributions from sedimentation velocity experiments were examined and developed. On the basis of four different model systems, we studied the differential apparent sedimentation coefficient distributions by the time-derivative method, g(s*), and by least-squares direct boundary modeling, ls-g*(s), the integral sedimentation coefficient distribution by the van Holde-Weischet method, G(s), and the previously introduced differential distribution of Lamm equation solutions, c(s). It is shown that the least-squares approach ls-g*(s) can be extrapolated to infinite time by considering area divisions analogous to boundary divisions in the van Holde-Weischet method, thus allowing the transformation of interference optical data into an integral sedimentation coefficient distribution G(s). However, despite the model-free approach of G(s), for the systems considered, the direct boundary modeling with a distribution of Lamm equation solutions c(s) exhibited the highest resolution and sensitivity. The c(s) approach requires an estimate for the size-dependent diffusion coefficients D(s), which is usually incorporated in the form of a weight-average frictional ratio of all species, or in the form of prior knowledge of the molar mass of the main species. We studied the influence of the weight-average frictional ratio on the quality of the fit, and found that it is well-determined by the data. As a direct boundary model, the calculated c(s) distribution can be combined with a nonlinear regression to optimize distribution parameters, such as the exact meniscus position, and the weight-average frictional ratio. Although c(s) is computationally the most complex, it has the potential for the highest resolution and sensitivity of the methods described.     

Direct sedimentation analysis of interference optical data in analytical ultracentrifugation

Sedimentation data acquired with the interference optical scanning system of the Optima XL-I analytical ultracentrifuge can exhibit time-invariant noise components, as well as small radial-invariant baseline offsets, both superimposed onto the radial fringe shift data resulting from the macromolecular solute distribution. A well-established method for the interpretation of such ultracentrifugation data is based on the analysis of time-differences of the measured fringe profiles, such as employed in the g(s*) method. We demonstrate how the technique of separation of linear and nonlinear parameters can be used in the modeling of interference data by unraveling the time-invariant and radial-invariant noise components. This allows the direct application of the recently developed approximate analytical and numerical solutions of the Lamm equation to the analysis of interference optical fringe profiles. The presented method is statistically advantageous since it does not require the differentiation of the data and the model functions. The method is demonstrated on experimental data and compared with the results of a g(s*) analysis. It is also demonstrated that the calculation of time-invariant noise components can be useful in the analysis of absorbance optical data. They can be extracted from data acquired during the approach to equilibrium, and can be used to increase the reliability of the results obtained from a sedimentation equilibrium analysis.     

A method for measuring sedimentation and diffusion of macromolecules in capillary tubes by total intensity and quasi-elastic light-scattering techniques.

Light scattered from a macromolecular solution in a capillary tube is used to determine both the sedimentation and translational diffusion coefficients. The capillary tube is spun in a preparative centrifuge, removed, and placed in a light-scattering photometer equipped with a scanning mechanism. The intensity distribution of scattered light along the tube represents the concentration profile in the tube and provides the measure of boundary migration. The sedimentation coefficient is determined from this measure and the applied centrifugal field. The diffusion coefficient is obtained from a time-autocorrelation analysis of fluctuations in intensity of light scattered from any fixed point of the profile. These coefficients were obtained for two monodisperse systems, R17 bacteriophage and 28s ribosomal rat liver RNA. The molecular weights obtained from ratios of these coefficients are in good agreement with literature values. In the sedimentation analysis, deviations from linearity between boundary displacement and applied field were found to be less than 1%. This precision confirms that the boundary is stable for the capillary geometry even in the absence of a preformed density gradient. The sedimentation coefficients of identical samples were also measured with the Spinco Model E analytical ultracentrifuge; results of the two methods agree to within 4%. As a consequence of the capillary tube geometry and light-scattering detection, sedimentation coefficients can be obtained from sample volumes of less than 100 μl. This detection techniques is thus far demonstrated to be at least an order of magnitude more sensitive than Schlieren optics, thereby useful when uv absorption is not applicable. For diffusion measurements there are also several inherent advantages. The diffusion coefficient is obtained from the identical sample, and scanning provides the capability to measure D from various parts of the sedimentation profiles and thereby directly explore concentration dependence, homogeneity, and integrity of the sample. The capillary tube with a layer of silicone oil over the sample and centrifugation provides an effective method to cleanse the solution and trap all dust.     

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.     

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.     

Fluorescence detection for the XLI analytical ultracentrifuge

Analytical ultracentrifugation (AUC) provides first-principle hydrodynamic and thermodynamic information concerning the size, shape and interactions of macromolecules. The fundamental measurement needed in AUC is the macromolecular concentration as a function of radial position and time. Currently, the Beckman Coulter XLI analytical ultracentrifuge may be equipped with absorbance and refractive detectors, which provide complementary concentration determinations. For detecting trace quantities of materials, fluorescence detection offers unique advantages over either absorbance or interference detection. A prototype fluorescence detector for the XLI analytical ultracentrifuge has been developed and its characteristics determined. An Ar(+) laser provides a continuous 488-nm excitation beam. Radial resolution is achieved by scanning the focused beam along a radial axis. Detection of the fluorescence signal uses a co-axial, front-face optical configuration to reduce inaccuracies in the concentration caused by inner filter effects. A high-speed A/D data acquisition system allows the fluorescence intensity to be monitored continuously and at a sufficiently high angular resolution so that at any radial position the intensities from all of the samples may be acquired at each revolution. The fluorescence detector is capable of detecting concentrations as low as 300 pM for fluorescein-like labels. The radial resolution of the fluorescence detector is comparable to that of the absorbance system. Both sedimentation velocity and sedimentation equilibrium measurements may be made with the fluorescence detector. Results are presented comparing data acquired using the fluorescence with those acquired using the absorbance detector.     

Molecular mass determination by sedimentation velocity experiments and direct fitting of the concentration profiles.

A new method for the direct molecular mass determination from sedimentation velocity experiments is presented. It is based on a nonlinear least squares fitting procedure of the concentration profiles and simultaneous estimation of the sedimentation and diffusion coefficients using approximate solutions of the Lamm equation. A computer program, LAMM, was written by using five different model functions derived by Fujita (1962, 1975) to describe the sedimentation of macromolecules during centrifugation. To compare the usefulness of these equations for the analysis of hydrodynamic results, the approach was tested on data sets of Claverie simulations as well as experimental curves of some proteins. A modification for one of the model functions is suggested, leading to more reliable sedimentation and diffusion coefficients estimated by the fitting procedure. The method seems useful for the rapid molecular mass determination of proteins larger than 10 kDa. One of the equations of the Archibald type is also suitable for compounds of low molecular mass, probably less than 10 kDa, because this model function requires neither the plateau region nor a meniscus free of solute.     

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.     

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.     

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.     

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.     

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.     

Sedimentation equilibrium analysis of protein interactions with global implicit mass conservation constraints and systematic noise decomposition

Sedimentation equilibrium is a powerful tool for the characterization of protein self-association and heterogeneous protein interactions. Frequently, it is applied in a configuration with relatively long solution columns and with equilibrium profiles being acquired sequentially at several rotor speeds. The present study proposes computational tools, implemented in the software SEDPHAT, for the global analysis of equilibrium data at multiple rotor speeds with multiple concentrations and multiple optical detection methods. The detailed global modeling of such equilibrium data can be a nontrivial computational problem. It was shown previously that mass conservation constraints can significantly improve and extend the analysis of heterogeneous protein interactions. Here, a method for using conservation of mass constraints for the macromolecular redistribution is proposed in which the effective loading concentrations are calculated from the sedimentation equilibrium profiles. The approach is similar to that described by Roark (Biophys. Chem. 5 (1976) 185-196), but its utility is extended by determining the bottom position of the solution columns from the macromolecular redistribution. For analyzing heterogeneous associations at multiple protein concentrations, additional constraints that relate the effective loading concentrations of the different components or their molar ratio in the global analysis are introduced. Equilibrium profiles at multiple rotor speeds also permit the algebraic determination of radial-dependent baseline profiles, which can govern interference optical ultracentrifugation data, but usually also occur, to a smaller extent, in absorbance optical data. Finally, the global analysis of equilibrium profiles at multiple rotor speeds with implicit mass conservation and computation of the bottom of the solution column provides an unbiased scale for determining molar mass distributions of noninteracting species. The properties of these tools are studied with theoretical and experimental data sets.     

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.