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.     

Determination of small differences in buoyant density of macromolecules by analytical density gradient equilibrium centrifugation.

A sensitive method is proposed for the determination of small differences between the buoyant densities of different species of monodisperse macromolecules by analytical density gradient equilibrium centrifugation. The procedure involves the measurement at sedimentation equilibrium of the bandwidths of the concentration distribution of the separate macromolecules and of a mixture of the different species. The difference in buoyant densities can then be estimated from the difference between the bandwidths.     

The functional binding epitope of a high affinity variant of human growth hormone mapped by shotgun alanine-scanning mutagenesis: insights into the mechanisms responsible for improved affinity

A high-affinity variant of human growth hormone (hGH(v)) contains 15 mutations within site 1 and binds to the hGH receptor (hGHR) approximately 400-fold tighter than does wild-type (wt) hGH (hGH(wt)). We used shotgun scanning combinatorial mutagenesis to dissect the energetic contributions of individual residues within the hGH(v) binding epitope and placed them in context with previously determined structural information. In all, the effects of alanine substitutions were determined for 35 hGH(v) residues that are directly contained in or closely border the binding interface. We found that the distribution of binding energy in the functional epitope of hGH(v) differs significantly from that of hGH(wt). The residues that contributed the majority of the binding energy in the wt interaction (the so-called binding "hot spot") remain important, but their contributions are attenuated in the hGH(v) interaction, and additional binding energy is acquired from residues on the periphery of the original hotspot. Many interactions that inhibited the binding of hGH(wt) are replaced by interactions that make positive contributions to the binding of hGH(v). These changes produce an expanded and diffused hot spot in which improved affinity results from numerous small contributions distributed broadly over the interface. The mutagenesis results are consistent with previous structural studies, which revealed widespread structural differences between the wt and variant hormone-receptor interfaces. Thus, it appears that the improved binding affinity of hGH(v) site 1 was not achieved through minor adjustments to the wt interface, but rather, results from a wholesale reconfiguration of many of the original binding elements.     

Isolation and characterization of a DnaJ-like protein in rats: the C-terminal 10-kDa domain of hsc70 is not essential for stimulating the ATP-hydrolytic activity of hsc70 by a DnaJ-like protein.

A DnaJ-like protein, RDJ1, was isolated from a rat brain cDNA library. The protein is predicted to have 397 amino acid residues and shares 99% identity to that of HDJ2, a human DnaJ-like protein. RDJ1 was also shown to rescue the temperature-sensitive lethality of a strain containing a mutated cytosolic DnaJ in yeast, ydj1-151. Fragments containing the J-domain of RDJ1 either with or without the G/F motif were expressed in Escherichia coli. The purified proteins stimulated the ATPase activity of hsc70 and of the 60-kDa N-terminal fragment of hsc70. These results imply that RDJ1 can interact with the N-terminal 60-kDa fragment of hsc70 to activate ATP hydrolysis by hsc70.     

Dissecting the energetics of an antibody-antigen interface by alanine shaving and molecular grafting.

Alanine-scanning mutagenesis on human growth hormone (hGH) identified 5 primary determinants (Arg 8, Asn 12, Arg 16, Asp 112, and Asp 116) for binding to a monoclonal antibody (MAb 3) (Jin L, Fendly BM, Wells JA, 1992, J Mol Biol 226:851-865). To further analyze the energetic importance of residues surrounding these five, we mutated all neighboring residues to alanine in groups of 7-16 (a procedure we call alanine shaving). Even the most extremely mutated variant, with 16 alanine substitutions, caused less than a 10-fold reduction in binding affinity to MAb3. By comparison, mutating any 1 of the 5 primary determinants to alanine caused a 6- to > 500-fold reduction in affinity. Replacing any of the 4 charged residues (Arg 8, Arg 16, Asp 112, and Asp 116) with a homologous residue (i.e., Arg to Lys or Asp to Glu) caused nearly as large a reduction in affinity as the corresponding alanine replacement. It was possible to graft the 5 primary binding determinants onto a nonbinding homologue of hGH, human placental lactogen (hPL), which has 86% sequence identity to hGH. The grafted hPL mutant bound 10-fold less tightly than hGH to MAb3 but bound as well as hGH when 2 additional framework mutations were introduced. Attempts to recover binding affinity by grafting the MAb3 epitope onto more distantly related scaffolds having a similar 4-helix bundle motif, such as human prolactin (23% sequence identity) or granulocyte colony-stimulating factor, were unsuccessful.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hemoglobin equilibrium analysis by the multiangle laser light-scattering method

Dimer-tetramer and monomer-dimer-tetramer equilibria of tetrameric hemoglobins and their single chains in the CO form, respectively, were evaluated using the microbatch multiangle light-scattering (MALS) analysis system. The molecular weights of human Hb A and Hb F in the CO form were dependent on concentration. The dissociation constants to dimers of Hb A and Hb F were 2.58 x 10(-6) and 0.66 x 10(-6), respectively. Equilibration of single globin chains, including alpha, beta, and gamma chains, was also evaluated by the same method. The dissociation constants of alpha-chain dimers to monomers, of beta-chain tetramers to monomers, and of gamma-chain tetramers to dimers were 14 x 10(-6), 25 x 10(-17), and 6.86 x 10(-6) M, respectively. These results indicate that the MALS analysis system can not only determine molecular weight but also characterize protein-protein interactions of multi-subunit proteins.     

The extended multidomain solution structures of the complement protein Crry and its chimeric conjugate Crry-Ig by scattering,analytical ultracentrifugation and constrained modelling: implications for function and therapy

Complement receptor-related gene/protein y (Crry) is a cell membrane-bound regulator of complement activation found in mouse and rat. Crry contains only short complement/consensus repeat (SCR) domains. X-ray and neutron scattering was performed on recombinant rat Crry containing the first five SCR domains (rCrry) and mouse Crry with five SCR domains conjugated to the Fc fragment of mouse IgG1 (mCrry-Ig) in order to determine their solution structures at medium resolution. The radius of gyration R(G) of rCrry was determined to be 4.9-5.0 nm, and the R(G) of the cross-section was 1.2-1.5 nm as determined by X-ray and neutron scattering. The R(G) of mCrry-Ig was 6.6-6.7 nm, and the R(G) of the cross-section were 2.3-2.4 nm and 1.3 nm. The maximum dimension of rCrry was 18 nm and that for mCrry-Ig was 26 nm. The neutron data indicated that rCrry and mCrry-Ig have molecular mass values of 45,000 Da and 140,000 Da, respectively, in agreement with their sequences, and sedimentation equilibrium data supported these determinations. Time-derivative velocity experiments gave sedimentation coefficients of 2.4S for rCrry and 5.4S for mCrry-Ig. A medium-resolution model of rCrry was determined using homology models that were constructed for the first five SCR domains of Crry from known crystal and NMR structures, and linked by randomly generated linker peptide conformations. These trial-and-error calculations revealed a small family of extended rCrry structures that best accounted for the scattering and ultracentrifugation data. These were shorter than the most extended rCrry models as the result of minor bends in the inter-SCR orientations. The mCrry-Ig solution data were modelled starting from a fixed structure for rCrry and the crystal structure of mouse IgG1, and was based on conformational searches of the hinge peptide joining the mCrry and Fc fragments. The best-fit models showed that the two mCrry antennae in mCrry-Ig were extended from the Fc fragment. No preferred orientation of the antennae was identified, and this indicated that the accessibility of the antennae for the molecular targets C4b and C3b was not affected by the covalent link to Fc. A structural comparison between Crry and complement receptor type 1 indicated that the domain arrangement of Crry SCR 1-3 is as extended as that of the CR1 SCR 15-17 NMR structure.     

Mutational analysis of human DNase I at the DNA binding interface: implications for DNA recognition,catalysis, and metal ion dependence.

Human deoxyribonuclease I (DNase I), an enzyme used to treat cystic fibrosis patients, has been systematically analyzed by site-directed mutagenesis of residues at the DNA binding interface. Crystal structures of bovine DNase I complexed with two different oligonucleotides have implicated the participation of over 20 amino acids in catalysis or DNA recognition. These residues have been classified into four groups based on the characterization of over 80 human DNase I variants. Mutations at any of the four catalytic amino acids His 134, His 252, Glu 78, and Asp 212 drastically reduced the hydrolytic activity of DNase I. Replacing the three putative divalent metal ion-coordinating residues Glu 39, Asp 168, or Asp 251 led to inactive variants. Amino acids Gln 9, Arg 41, Tyr 76, Arg 111, Asn 170, Tyr 175, and Tyr 211 were also critical for activity, presumably because of their close proximity to the active site, while more peripheral DNA interactions stemming from 13 other positions were of minimal significance. The relative importance of these 27 positions is consistent with evolutionary relationships among DNase I across different species, DNase I-like proteins, and bacterial sphingomyelinases, suggesting a fingerprint for a family of DNase I-like proteins. Furthermore, we found no evidence for a second active site that had been previously implicated in Mn2+-dependent DNA degradation. Finally, we correlated our mutational analysis of human DNase I to that of bovine DNase I with respect to their specific activity and dependence on divalent metal ions.     

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).     

Quantitative evaluation of the chicken lysozyme epitope in the HyHEL-10 Fab complex: free energies and kinetics.

The hen (chicken) egg-white lysozyme (HEWL) epitope for the monoclonal antibody HyHEL-10 Fab (Fab-10) was investigated by alanine scan mutagenesis. The association rate constants (k(on)) for the HEWL Fab-10 complexes were obtained from the homogenous solution method described in the preceding paper (Taylor et al., 1998). A new method for determining the dissociation rate constant (k(off)) for the complex, by trapping nascent free antibody with an inactive HEWL mutant is described. The values of k(on) fall within a factor of 2 of the wild-type (WT) HEWL value (1.43+/-0.13 X 10(6)M(-1)s(-1)), while the increases in k(off)more nearly reflect the total change in free energies of the complex (deltadeltaG(D)). The dissociation constants (K(D)) were measured directly in those cases where satisfactory kinetic data could not be obtained. The Y20A, K96A, and K97A HEWL.Fab-10 complexes are destabilized by more than 4 kcal/mol compared to the WT complex. The R21A, L75A, and D101A antibody complexes are moderately destabilized (0.7 < deltadeltaG(D)< or = 1.0 kcal/mol). Additional mutations of the "hotspot" residues (Tyr20, Lys96, Lys97) were constructed to probe, more precisely, the nature of their contributions to complex formation. The results show that the entire hydrocarbon side chains of Tyr20 and Lys97, and only the epsilon-amino group of Lys96, contribute to the stability of the complex. The value of deltadeltaG(D) for the R21A mutant complex is a distinct outlier in the Arg21 replacement series demonstrating the importance of supplementing alanine scan mutagenesis with additional mutations.     

Energetic analysis of an antigen/antibody interface: alanine scanning mutagenesis and double mutant cycles on the HyHEL-10/lysozyme interaction.

Alanine scanning mutagenesis of the HyHEL-10 paratope of the HyHEL-10/HEWL complex demonstrates that the energetically important side chains (hot spots) of both partners are in contact. A plot of deltadeltaG(HyHEL-10_mutant) vs. deltadeltaG(HEWL_mutant) for the five of six interacting side-chain hydrogen bonds is linear (Slope = 1). Only 3 of the 13 residues in the HEWL epitope contribute >4 kcal/mol to the free energy of formation of the complex when replaced by alanine, but 6 of the 12 HyHEL-10 paratope amino acids do. Double mutant cycle analysis of the single crystallographically identified salt bridge, D32H/K97, shows that there is a significant energetic penalty when either partner is replaced with a neutral side-chain amino acid, but the D32(H)N/K97M complex is as stable as the WT. The role of the disproportionately high number of Tyr residues in the CDR was evaluated by comparing the deltadeltaG values of the Tyr --> Phe vs. the corresponding Tyr --> Ala mutations. The nonpolar contacts in the light chain contribute only about one-half of the total deltadeltaG observed for the Tyr --> Ala mutation, while they are significantly more important in the heavy chain. Replacement of the N31L/K96 hydrogen bond with a salt bridge, N31D(L)/K96, destabilizes the complex by 1.4 kcal/mol. The free energy of interaction, deltadeltaG(int), obtained from double mutant cycle analysis showed that deltadeltaG(int) for any complex for which the HEWL residue probed is a major immunodeterminant is very close to the loss of free energy observed for the HyHEL-10 single mutant. Error propagation analysis of double mutant cycles shows that data of atypically high precision are required to use this method meaningfully, except where large deltadeltaG values are analyzed.     

Purification, stabilization, and concentration of very weak protein-protein complexes: Shifting the association equilibrium via complex selective adsorption on lowly activated supports

Very weak protein-protein interactions are very difficult to detect because these complexes could be under the detection limit or they tend to dissociate. Here, using as a model the antibody-antigen interaction weaken by the presence of dioxane, we have shown a strategy for the protein complexes purification by selective adsorption of the associated proteins. This strategy is based on the use of poorly activated anionic exchanger supports to selectively adsorb large complexes. This selective adsorption of the associated proteins shifted the association equilibrium of the soluble proteins toward the associated form. Thus, in the presence of 15% v/v dioxane, a concentration that is able to almost fully break the immunocomplex (less that 3% of the immunocomplex appeared associated when soluble antigen-antibody mixture was cross-linked with aldehyde-dextran), we can obtain more than 90% of the fully pure immunocomplex from the non-associated protein, adsorbed on anionic exchanger supports having a very low activation. This simple strategy may be a very useful tool to solve one of the most relevant challenges in the modern proteomics, the detection of very weak protein-protein interactions.     

The structure analysis of Humanin analog, AGA-(C8R)HNG17, by circular dichroism and sedimentation equilibrium: comparison with the parent molecule

A 24-amino acid peptide, Humanin (HN), is a novel peptide that protects neuronal cells in vitro and in vivo from Alzheimer's disease-related toxicities. We have shown before that the structures of HN and a 1000-fold more active analog, HNG, with a Ser14Gly mutation are largely disordered. During additional mutational analysis, a shorter 17-amino acid form, AGA-(C8R)HNG17, was accidentally discovered to have a 100-fold higher activity than HNG. Here we have characterized the structural properties of the AGA-(C8R)HNG17 analog by circular dichroism (CD) and sedimentation equilibrium analysis. First, the structure in water was characterized, since these peptides have been dissolved in water prior to biological analysis. The AGA-(C8R)HNG17 peptide exhibited extensive beta-sheet structure in water, completely different from the aqueous HN and HNG structures. The beta-sheet structure was converted to a disordered structure upon dilution into phosphate-buffered saline (PBS) at low peptide concentration (e.g., below 0.2mg/ml), which was similar to the structure of HN and HNG, observed under similar conditions. Sedimentation equilibrium analysis showed that the AGA-(C8R)HNG17 analog was essentially monomeric in PBS, while HNG showed extensive aggregation. Such aggregation of HNG was observed when the peptide was added to the serum-containing cell culture media. Thus, the mutations introduced into the AGA-(C8R)HNG17 analog generated a peptide different from the parent HNG and HN peptides in the self-association properties and hence the solubility, which most likely contributed to the increased biological activity of the AGA-(C8R)HNG17 analog.     

The Switch that Does Not Flip: The Blue-Light Receptor YtvA from Bacillus subtilis Adopts an Elongated Dimer Conformation Independent of the Activation State as Revealed by a Combined AUC and SAXS Study

Photoreceptors play an important role in plants and bacteria by converting extracellular stimuli into intracellular signals. One distinct class are the blue-light-sensitive phototropins harboring a light-oxygen-voltage (LOV) domain coupled to various effector domains. Photon absorption by the chromophore within the LOV domain results in an activation of the output domain via mechanisms that are hitherto not well understood. The photoreceptor YtvA from Bacillus subtilis is a bacterial analog of phototropins, consists of an LOV and a sulfate transporter/anti-sigma factor antagonist domain, and is involved in the response of the bacterium to environmental stress. We present here analytical ultracentrifugation studies and small-angle X-ray scattering experiments, showing that YtvA is a dimer. On the basis of these results, we present a low-resolution model of the dimer in the dark and the lit state of the protein. In addition, we show that YtvA does not change its oligomerization state or its overall shape upon light activation.