Self-association of human beta-galactocerebrosidase: Dependence on pH, salt, and surfactant.

Krabbe disease, also known as globoid cell leukodystrophy, is a rare genetic neurodegenerative disease caused by a deficiency of the galactocerebrosidase enzyme. To understand the association status of human beta-galactocerebrosidase (hGALC) in solution, we employed analytical ultracentrifugation (AUC). Our AUC results show that hGALC has a tendency for reversible self-association. Self-association decreases as the concentration of sodium chloride increases from 50 to 500 mM. This indicates that ionic interactions are involved in the association. The association is also dependent on pH, and high order oligomerization decreases as the pH increases from 4.5 to 7.5. Taken together, our results indicate that hGALC has the highest tendency for oligomerization at physiological ionic strength and pH (lysosomal lumen). This is the first report describing the self-associating property of hGALC in solution.

A Comprehensive Brownian Dynamics Approach for the Determination of Non-ideality Parameters from Analytical Ultracentrifugation.

Brownian dynamics (BD) has been applied as a comprehensive tool to model sedimentation and diffusion of nanoparticles in analytical ultracentrifugation (AUC) experiments. In this article, we extend the BD algorithm by considering space-dependent diffusion and solvent compressibility. With this, the changes in the sedimentation and diffusion coefficient from altered solvent properties at increased pressures are accurately taken into account. Moreover, it is demonstrated how the concept of space-dependent diffusion is employed to describe concentration-dependent sedimentation and diffusion coefficients, in particular, through the Gralen coefficient and the second virial coefficient. The influence of thermodynamic nonideality on diffusional properties can be accurately simulated and agree with well-known evaluation tools. BD simulations for sedimentation equilibrium and sedimentation velocity (SV) AUC experiments including effects of hydrodynamic and thermodynamic nonideality are validated by global evaluation in SEDANAL. The interplay of solvent compressibility and retrieved nonideality parameters can be studied utilizing BD. Finally, the second virial coefficient is determined for lysozyme from SV AUC experiments and BD simulations and compared to membrane osmometry. These results are in line with DLVO theory. In summary, BD simulations are established for the validation of nonideal sedimentation in AUC providing a sound basis for the evaluation of complex interactions even in polydisperse systems.

AUC measurements of diffusion coefficients of monoclonal antibodies in the presence of human serum proteins.

The goal of this work is to develop a preclinical method for quantitative hydrodynamic and thermodynamic analysis of therapeutic proteins in crowded environments like human serum. The method utilizes tracer amounts of fluorescently labeled monoclonal antibodies and the Aviv AU-FDS optical system. We have performed sedimentation velocity experiments as a function of mAb, human serum albumin and human IgG concentration to extract self- and cross-term hydrodynamic nonideality effects. SV measurements are consistently complicated by weak mAb-mAb and mAb-IgG interactions (Wright et al. in Anal Biochem 550:72-83, 2018). In an attempt to explore different approaches we have investigated measurements of diffusion coefficients by traditional synthetic boundary experiments. Here we present a new technique incorporated into SEDANAL that can globally analyze the full time course of synthetic boundary experiments. This approach also utilizes F-mAb against a high concentration of unlabeled carrier protein (HSA or IgG). In principle both diffusion and sedimentation coefficient information can be extracted including hydrodynamic and thermodynamic nonideality. The method can be performed at a traditional low speed (5-7K rpm) or at high speeds. The high speed method can also be used to measure D and s for small molecules like fluorescein (often contaminants of F-HSA and F-mAb). The advantage of synthetic boundary over the standard sedimentation velocity method is that it allows for higher precision determination of diffusion coefficients. The concentration dependence of D can be corrected for hydrodynamic nonideality effects by plotting D * (1 + kijcj) vs total carrier concentration. The slope of the fitted data allows an alternate approach to determine self- and cross-term thermodynamic nonideality. This method can also explore cross-term diffusion coefficient effects. These results are compared to dynamic light scattering approaches which are limited to kD determinations for solutions of pure protein.

Characterization of therapeutic antibodies in the presence of human serum proteins by AU-FDS analytical ultracentrifugation.

The preclinical characterization of biopharmaceuticals seeks to determine the stability, state of aggregation, and interaction of the antibody/drug with other macromolecules in serum. Analytical ultracentrifugation is the best experimental method to understand these factors. Sedimentation velocity experiments using the AU-FDS system were performed in order to quantitatively characterize the nonideality of fluorescently labeled therapeutic antibodies in high concentrations of human serum proteins. The two most ubiquitous serum proteins are human serum albumin, HSA, and γ-globulins, predominantly IgG. Tracer experiments were done pairwise as a function of HSA, IgG, and therapeutic antibody concentration. The sedimentation coefficient for each fluorescently labeled component as a function of the concentration of the unlabeled component yields the hydrodynamic nonideality (ks). This generates a 3x3 matrix of ks values that describe the nonideality of each pairwise interaction. The ks matrix is validated by fitting both 2:1 mixtures of HSA (1-40 mg/ml) and IgG (0.5-20 mg/ml) as serum mimics, and human serum dilutions (10-100%). The data are well described by SEDANAL global fitting with the ks nonideality matrix. The ks values for antibodies are smaller than expected and appear to be masked by weak association. Global fitting to a ks and K2 model significantly improves the fits.

Studying polyglutamine aggregation in Caenorhabditis elegans using an analytical ultracentrifuge equipped with fluorescence detection.

This work explores the heterogeneity of aggregation of polyglutamine fusion constructs in crude extracts of transgenic Caenorhabditis elegans animals. The work takes advantage of the recent technical advances in fluorescence detection for the analytical ultracentrifuge. Further, new sedimentation velocity methods, such as the multi-speed method for data capture and wide distribution analysis for data analysis, are applied to improve the resolution of the measures of heterogeneity over a wide range of sizes. The focus here is to test the ability to measure sedimentation of polyglutamine aggregates in complex mixtures as a prelude to future studies that will explore the effects of genetic manipulation and environment on aggregation and toxicity. Using sedimentation velocity methods, we can detect a wide range of aggregates, ranging from robust analysis of the monomer species through an intermediate and quite heterogeneous population of oligomeric species, and all the way up to detecting species that likely represent intact inclusion bodies based on comparison to an analysis of fluorescent puncta in living worms by confocal microscopy. Our results support the hypothesis that misfolding of expanded polyglutamine tracts into insoluble aggregates involves transitions through a number of stable intermediate structures, a model that accounts for how an aggregation pathway can lead to intermediates that can have varying toxic or protective attributes. An understanding of the details of intermediate and large-scale aggregation for polyglutamine sequences, as found in neurodegenerative diseases such as Huntington's Disease, will help to more precisely identify which aggregated species may be involved in toxicity and disease.

Interaction of Myosin Phosphatase Target Subunit (MYPT1) with Myosin Phosphatase-RhoA Interacting Protein (MRIP): A Role of Glutamic Acids in the Interaction.

Scaffold proteins bind to and functionally link protein members of signaling pathways. Interaction of the scaffold proteins, myosin phosphatase target subunit (MYPT1) and myosin phosphatase-RhoA interacting protein (MRIP), causes co-localization of myosin phosphatase and RhoA to actomyosin. To examine biophysical properties of interaction of MYPT1 with MRIP, we employed analytical ultracentrifugation and surface plasmon resonance. In regard to MRIP, its residues 724-837 are sufficient for the MYPT1/MRIP interaction. Moreover, MRIP binds to MYPT1 as either a monomer or a dimer. With respect to MYPT1, its leucine repeat region, LR (residues 991-1030) is sufficient to account for the MYPT1/MRIP interaction. Furthermore, point mutations that replace glutamic acids 998-1000 within LR reduced the binding affinity toward MRIP. This suggests that the glutamic acids of MYPT1 play an important role in the interaction.

Sedimentation Velocity: A Classical Perspective.

Here we give an overview of the history of sedimentation velocity analysis focusing on seminal and fundamental contributions that derived from early ultracentrifugation studies. We introduce the concepts of nonequilibrium thermodynamics and outline the derivation of the Svedberg and the Lamm equations and the requirements for including both hydrodynamic and thermodynamic nonideality. We introduce the phenomenological equations for coupled flows as developed from the principles of nonequilibrium or irreversible thermodynamics and derive a form of the Lamm equation that incorporates cross-diffusion coefficients and coupled gradient terms. We give an historical overview of solutions to the Lamm equation including Fujita-MacCosham solutions and Claverie finite-element numerical solutions and discuss the software that have implemented these solutions. We discuss the three major optical systems (absorbance, interference, and fluorescence) and recently developed multiwavelength systems. We also suggest a number of experimental practices and guidelines for optimizing the determination of s and D and discuss the appropriate centerpiece components and their utility. This chapter complements other recent reviews submitted by the authors (Correia, Lyons, Sherwood, & Stafford, 2015; Stafford, 2015) and should be considered an effort to revive the importance of irreversible thermodynamics in the understanding and analysis of sedimentation velocity ultracentrifugation data.

Simultaneous analysis of hydrodynamic and optical properties using analytical ultracentrifugation equipped with multiwavelength detection.

Analytical ultracentrifugation (AUC) has proven to be a powerful tool for the study of particle size distributions, particle shapes, and interactions with high accuracy and unrevealed resolution. In this work we show how the analysis of sedimentation velocity data from the AUC equipped with a multiwavelength detector (MWL) can be used to gain an even deeper understanding of colloidal and macromolecular mixtures. New data evaluation routines have been integrated in the software SEDANAL to allow for the handling of MWL data. This opens up a variety of new possibilities because spectroscopic information becomes available for individual components in mixtures at the same time using MWL-AUC. For systems of known optical properties information on the hydrodynamic properties of the individual components in a mixture becomes accessible. For the first time, the determination of individual extinction spectra of components in mixtures is demonstrated via MWL evaluation of sedimentation velocity data. In our paper we first provide the informational background for the data analysis and expose the accessible parameters of our methodology. We further demonstrate the data evaluation by means of simulated data. Finally, we give two examples which are highly relevant in the field of nanotechnology using colored silica and gold nanoparticles of different size and extinction properties.

Assembly of the SLIP1-SLBP complex on histone mRNA requires heterodimerization and sequential binding of SLBP followed by SLIP1.

The SLIP1-SLBP complex activates translation of replication-dependent histone mRNAs. In this report, we describe how the activity of the SLIP1-SLBP complex is modulated by phosphorylation and oligomerization. Biophysical characterization of the free proteins shows that whereas SLIP1 is a homodimer that does not bind RNA, human SLBP is an intrinsically disordered protein that is phosphorylated at 23 Ser/Thr sites when expressed in a eukaryotic expression system such as baculovirus. The bacterially expressed unphosphorylated SLIP1-SLBP complex forms a 2:2 high-affinity (K(D) < 0.9 nM) heterotetramer that is also incapable of binding histone mRNA. In contrast, phosphorylated SLBP from baculovirus has a weak affinity (K(D) ~3 µM) for SLIP1. Sequential binding of phosphorylated SLBP to the histone mRNA stem-loop motif followed by association with SLIP1 is required to form an "active" ternary complex. Phosphorylation of SLBP at Thr171 promotes dissociation of the heterotetramer to the SLIP1-SLBP heterodimer. Using alanine scanning mutagenesis, we demonstrate that the binding site on SLIP1 for SLBP lies close to the dimer interface. A single-point mutant near the SLIP1 homodimer interface abolished interaction with SLBP in vitro and reduced the abundance of histone mRNA in vivo. On the basis of these biophysical studies, we propose that oligomerization and SLBP phosphorylation may regulate the SLBP-SLIP1 complex in vivo. SLIP1 may act to sequester SLBP in vivo, protecting it from proteolytic degradation as an inactive heterotetramer, or alternatively, formation of the SLIP1-SLBP heterotetramer may facilitate removal of SLBP from the histone mRNA prior to histone mRNA degradation.

Disulfide cross-linked antiparallel actin dimer.

Oxidation of actin monomer (G-actin) with copper o-phenanthroline resulted in a rapid, high yield of disulfide cross-linked dimer. The cross-link is due to an intermolecular disulfide bond between actin Cys374 of each molecule, resulting in a tail-to-tail, i.e., antiparallel, actin dimer. Analytical ultracentrifugation profiles of G-actin can be ascribed to the existence of actin monomers with very little, if any, dimer. Thus, actin dimers are not energetically favorable, indicating that cross-linked dimers are formed during random diffusional collisions. On the other hand, a similar oxidation of actin polymer (F-actin) resulted in a much lower yield of the cross-linked actin dimer that showed no sign of leveling off. Therefore, it is proposed that the cross-linked dimer from actin polymer is due to collisional complexes of actin monomers that are in equilibrium with the polymer during actin treadmilling. These results account for the reported observation that during the early stages of actin polymerization (where the actin monomer concentration is high) cross-linked antiparallel actin dimers are formed in relatively high yield whereas none are formed at later stages of polymerization. These findings raise questions concerning the validity of the antiparallel actin dimer model of in vitro actin polymerization that is based on the assumption that the ability to form cross-linked actin dimers implies the existence of stable dimers.

Equilibrium self-association of tropomyosin.

It has recently been reported that tropomyosin exists exclusively as a dimer in physiological salt conditions. It is shown in the present work using analytical ultracentrifugation that, on the contrary, tropomyosin is in equilibrium between monomer, dimer and tetramer with a weak tendency to dimerize and tetramerize. Such a finding has consequences for the assembly of the tropomyosin-actin complex.

Characterization of the self-association of human interferon-α2b, albinterferon-α2b, and pegasys.

The self-association of human interferon-α2b (hIFN-α2b), albinterferon-α2b (a recombinant protein with human serum albumin and hIFN-α2b peptides fused together in a single polypeptide chain), and Pegasys (PEGylated hIFN-α2a) was characterized by analytical ultracentrifugation analyses. By examining the apparent sedimentation coefficient distribution profiles of each protein at different concentrations, it was concluded that the above three proteins are self-associating in albinterferon-α2b formulation buffer. By model fitting of sedimentation data using SEDANAL software, the stoichiometry and equilibrium constants of the self-association of these proteins were characterized. The self-association of hIFN-α2b results in the formation of stable dimers, fast-reversible tetramers, octamers, and hexadecamers. In contrast, although both albinterferon-α2b and Pegasys are self-associated, their self-association stoichiometries are significantly different from that of hIFN-α2b. The self-association of albinterferon-α2b results in the formation of reversible dimers and trimers, whereas the self-association of Pegasys gives only reversible dimers. The self-association behaviors of hIFN-α2b and albinterferon-α2b involves attractive electrostatic forces, which can be suppressed to a negligible level in low pH (pH 4.0-4.5) and high salt concentration (400 mM NaCl) buffer, allowing quantification of their size variant contents by sedimentation velocity analysis.

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.

The kinetic mechanism of mouse myosin VIIA.

Myosin VIIa is crucial in hearing and visual processes. We examined the kinetic and association properties of the baculovirus expressed, truncated mouse myosin VIIa construct containing the head, all 5IQ motifs and the putative coiled coil domain (myosin VIIa-5IQ). The construct appears to be monomeric as determined by analytical ultracentrifugation experiments, and only single headed molecules were detected by negative stain electron microscopy. The relatively high basal steady-state rate of 0.18 s(-1) is activated by actin only by ∼3.5-fold resulting in a V(max) of 0.7 s(-1) and a K(ATPase) of 11.5 µM. There is no single rate-limiting step of the ATP hydrolysis cycle. The ATP hydrolysis step (M·T M·D·P) is slow (12 s(-1)) and the equilibrium constant (K(H)) of 1 suggests significant reversal of hydrolysis. In the presence of actin ADP dissociates with a rate constant of 1.2 s(-1). Phosphate dissociation is relatively fast (>12 s(-1)), but the maximal rate could not be experimentally obtained at actin concentrations ≤ 50 µM because of the weak binding of the myosin VIIa-ADP-P(i) complex to actin. At higher actin concentrations the rate of attached hydrolysis (0.4 s(-1)) becomes significant and partially rate-limiting. Our findings suggest that the myosin VIIa is a "slow", monomeric molecular motor with a duty ratio of 0.6.

SEDVIEW, real-time sedimentation analysis.

The ability to obtain a sedimentation coefficient distribution as the run proceeds, and to get an early idea of the quality of a particular sample, has not been made available in real-time during the run in any of the existing software packages. It is desirable on many occasions to be able to see the number of components present in a sample at an early stage of the run. The ability to ascertain the extent of heterogeneity of sample would help enormously to reduce the amount of time that is necessary to obtain that information. Most software packages currently available require that the run be completed before analysis is carried out or at least some of the early scans analyzed off-line to determine if the run should continue. A software package called SEDVIEW has been developed by us to allow early analysis in real-time.

Modular structure of smooth muscle Myosin light chain kinase: hydrodynamic modeling and functional implications.

Smooth muscle myosin light chain kinase (smMLCK) is a calcium-calmodulin complex-dependent enzyme that activates contraction of smooth muscle. The polypeptide chain of rabbit uterine smMLCK (Swiss-Prot entry P29294) contains the catalytic/regulatory domain, three immunoglobulin-related motifs (Ig), one fibronectin-related motif (Fn3), a repetitive, proline-rich segment (PEVK), and, at the N-terminus, a unique F-actin-binding domain. We have evaluated the spatial arrangement of these domains in a recombinant 125 kDa full-length smMLCK and its two catalytically active C-terminal fragments (77 kDa, residues 461-1147, and 61 kDa, residues 461-1002). Electron microscopic images of smMLCK cross-linked to F-actin show particles at variable distances (11-55 nm) from the filament, suggesting that a well-structured C-terminal segment of smMLCK is connected to the actin-binding domain by a long, flexible tether. We have used structural homology and molecular dynamics methods to construct various all-atom representation models of smMLCK and its two fragments. The theoretical sedimentation coefficients computed with HYDROPRO were compared with those determined by sedimentation velocity. We found agreement between the predicted and observed sedimentation coefficients for models in which the independently folded catalytic domain, Fn3, and Ig domains are aligned consecutively on the long axis of the molecule. The PEVK segment is modeled as an extensible linker that enables smMLCK to remain bound to F-actin and simultaneously activate the myosin heads of adjacent myosin filaments at a distance of >or=40 nm. The structural properties of smMLCK may contribute to the elasticity of smooth muscle cells.

Hydrodynamic and mass spectrometry analysis of nearly-intact human fibrinogen, chicken fibrinogen, and of a substantially monodisperse human fibrinogen fragment X.

The shape and solution properties of fibrinogen are affected by the location of the C-terminal portion of the Aalpha chains, which is presently still controversial. We have measured the hydrodynamic properties of a human fibrinogen fraction with these appendages mostly intact, of chicken fibrinogen, where they lack 11 characteristic 13-amino acids repeats, and of human fragment X, a plasmin early degradation product in which they have been removed. The human fibrinogen/fragment X samples were extensively characterized by SDS-PAGE/Western blotting and mass spectrometry, allowing their composition to be precisely determined. The solution properties of all samples were then investigated by analytical ultracentrifugation and size-exclusion HPLC coupled with multi-angle light scattering and differential pressure viscometry detectors. The measured parameters suggest that the extra repeats have little influence on the overall fibrinogen conformation, while a significant change is brought about by the removal of the C-terminal portion of the Aalpha chains beyond residue Aalpha200.

Effect of kinetics on sedimentation velocity profiles and the role of intermediates.

We have previously presented a tutorial on direct boundary fitting of sedimentation velocity data for kinetically mediated monomer-dimer systems [Correia and Stafford, 2009]. We emphasized the ability of Sedanal to fit for the k(off) values and measure their uncertainty at the 95% confidence interval. We concluded for a monomer-dimer system the range of well-determined k(off) values is limited to 0.005-10(-5) s(-1) corresponding to relaxation times of approximately 70 to approximately 33,000 s. More complicated reaction schemes introduce the potential complexity of low concentrations of an intermediate that may also influence the kinetic behavior during sedimentation. This can be seen in a cooperative ABCD system (A+B --> C; B+C --> D) where C, the 1:1 complex, is sparsely populated (K(1)=10(4) M(-1), K(2)=10(8) M(-1)). Under these conditions a k(1,off)<0.01 s(-1) produces slow kinetic features. The low concentration of species C contributes to this effect while still allowing the accurate estimation of k(1,off) (although k(2,off) can readily compensate and contribute to the kinetics). More complex reactions involving concerted assembly or cooperative ring formation with low concentrations of intermediate species also display kinetic effects due to a slow flux of material through the sparsely populated intermediate states. This produces a kinetically limited reaction boundary that produces partial resolution of individual species during sedimentation. Cooperativity of ring formation drives the reaction and thus separation of these two effects, kinetics and energetics, can be challenging. This situation is experimentally exhibited by systems that form large oligomers or rings and may especially contribute to formation of micelles and various protein aggregation diseases including formation of beta-amyloid and tau aggregates. Simulations, quantitative parameter estimation by direct boundary fitting and diagnostic features for these systems are presented with an emphasis on the features available in Sedanal to simulate and analyze kinetically mediated systems.

Phosphorylation dependence of hsp27 multimeric size and molecular chaperone function.

The molecular chaperone Hsp27 exists as a distribution of large oligomers that are disassembled by phosphorylation at Ser-15, -78, and -82. It is controversial whether the unphosphorylated Hsp27 or the widely used triple Ser-to-Asp phospho-mimic mutant is the more active molecular chaperone in vitro. This question was investigated here by correlating chaperone activity, as measured by the aggregation of reduced insulin or alpha-lactalbumin, with Hsp27 self-association as monitored by analytical ultracentrifugation. Furthermore, because the phospho-mimic is generally assumed to reproduce the phosphorylated molecule, the size and chaperone activity of phosphorylated Hsp27 were compared with that of the phospho-mimic. Hsp27 was triply phosphorylated by MAPKAP-2 kinase, and phosphorylation was tracked by urea-PAGE. An increasing degree of suppression of insulin or alpha-lactalbumin aggregation correlated with a decreasing Hsp27 self-association, which was the least for phosphorylated Hsp27 followed by the mimic followed by the unphosphorylated protein. It was also found that Hsp27 added to pre-aggregated insulin did not reverse aggregation but did inhibit these aggregates from assembling into even larger aggregates. This chaperone activity appears to be independent of Hsp27 phosphorylation. In conclusion, the most active chaperone of insulin and alpha-lactalbumin was the Hsp27 (elongated) dimer, the smallest Hsp27 subunit observed under physiological conditions. Next, the Hsp27 phospho-mimic is only a partial mimic of phosphorylated Hsp27, both in self-association and in chaperone function. Finally, the efficient inhibition of insulin aggregation by Hsp27 dimer led to the proposal of two models for this chaperone activity.

Extracting equilibrium constants from kinetically limited reacting systems.

It has been known for some time that slow kinetics will distort the shape of a reversible reaction boundary. Here we present a tutorial on direct boundary fitting of sedimentation velocity data for a monomer-dimer system that exhibits kinetic effects. Previous analysis of a monomer-dimer system suggested that rapid reaction behavior will persist until the relaxation time of the system exceeds 100 s (reviewed in Kegeles and Cann, 1978). Utilizing a kinetic integrator feature in Sedanal (Stafford and Sherwood, 2004), we can now fit for the k(off) values and measure the uncertainty at the 95% confidence interval. For the monomer-dimer system the range of well determined k(off) values is limited to 0.005 to 10(-5) s(-1) corresponding to relaxation times (at a loading concentration of the Kd) of approximately 70 to approximately 33,000 s. For shorter relaxation times the system is fast and only the equilibrium constant K but not k(off) can be uniquely determined. For longer relaxation times the system is irreversibly slow, and assuming the system was at initial equilibrium before the start of the run, only the equilibrium constant K but not k(off) can be uniquely determined.