Novel size-independent modeling of the dilute solution conformation of the immunoglobulin IgG Fab' domain using SOLPRO and ELLIPS

Terpene polyacrylate TPA5 shows favorable molecular hydrodynamic properties as a potential bioinspired archaeological wood consolidant

There is currently a pressing need for the development of novel bioinspired consolidants for waterlogged, archaeological wood. Bioinspired materials possess many advantages, such as biocompatibility and sustainability, which makes them ideal to use in this capacity. Based on this, a polyhydroxylated monomer was synthesised from α-pinene, a sustainable terpene feedstock derived from pine trees, and used to prepare a low molar mass polymer TPA5 through free radical polymerisation. This polymer was extensively characterised by NMR spectroscopy (chemical composition) and molecular hydrodynamics, primarily using analytical ultracentrifugation reinforced by gel filtration chromatography and viscometry, in order to investigate whether it would be suitable for wood consolidation purposes. Sedimentation equilibrium indicated a weight average molar mass M_w of (4.3 ± 0.2) kDa, with minimal concentration dependence. Further analysis with MULTISIG revealed a broad distribution of molar masses and this heterogeneity was further confirmed by sedimentation velocity. Conformation analyses with the Perrin P and viscosity increment ν universal hydrodynamic parameters indicated that the polymer had an elongated shape, with both factors giving consistent results and a consensus axial ratio of ~ 4.5. These collective properties-hydrogen bonding potential enhanced by an elongated shape, together with a small injectable molar mass-suggest this polymer is worthy of further consideration as a potential consolidant.

Simulation and Theory of Antibody Binding to Crowded Antigen-Covered Surfaces

In this paper we introduce a fully flexible coarse-grained model of immunoglobulin G (IgG) antibodies parametrized directly on cryo-EM data and simulate the binding dynamics of many IgGs to antigens adsorbed on a surface at increasing densities. Moreover, we work out a theoretical model that allows to explain all the features observed in the simulations. Our combined computational and theoretical framework is in excellent agreement with surface-plasmon resonance data and allows us to establish a number of important results. (i) Internal flexibility is key to maximize bivalent binding, flexible IgGs being able to explore the surface with their second arm in search for an available hapten. This is made clear by the strongly reduced ability to bind with both arms displayed by artificial IgGs designed to rigidly keep a prescribed shape. (ii) The large size of IgGs is instrumental to keep neighboring molecules at a certain distance (surface repulsion), which essentially makes antigens within reach of the second Fab always unoccupied on average. (iii) One needs to account independently for the thermodynamic and geometric factors that regulate the binding equilibrium. The key geometrical parameters, besides excluded-volume repulsion, describe the screening of free haptens by neighboring bound antibodies. We prove that the thermodynamic parameters govern the low-antigen-concentration regime, while the surface screening and repulsion only affect the binding at high hapten densities. Importantly, we prove that screening effects are concealed in relative measures, such as the fraction of bivalently bound antibodies. Overall, our model provides a valuable, accurate theoretical paradigm beyond existing frameworks to interpret experimental profiles of antibodies binding to multi-valent surfaces of different sorts in many contexts.

Antibodies are the main working horses of the human immune system. Remarkably, no matter the size or the shape of the pathological intruders, these extremely flexible three-lobe molecules are able to form a complex, thus eliciting an immune response. What makes antibodies so effective? To answer this and other questions, we have developed a simplified computational scheme to simulate the dynamics of many antibodies interacting with each other and with antigens. Coarse-grained models are a great opportunity, as they give access to a true multi-scale approach to biologically relevant problems. In this work, our innovative method allowed us to simulate the binding process of many antibodies to surface-adsorbed antigens. This led us to elucidate and quantify many important physical aspects of their biological function in agreement with experiments, such as the role of their flexibility and crowding effects at the hapten-covered surface, which were shown to finely regulate the avidity.

Single-laser polarization FRET (polFRET) on the cell surface

A new method for the simultaneous detection of rotational mobility and proximity of cell surface receptors is presented based on cell-by-cell basis measurement of polarized fluorescence intensity components of the donor and acceptor of a FRET system. In addition to the FRET efficiency and the donor and acceptor concentrations, the method makes also possible the determination of the rotational characteristics and the associated fraction of the donors (FRET-fraction). The method is illustrated with flow cytometric and rFLIM measurements on donor-acceptor systems comprising fluorescently labeled whole antibodies and their Fab fragments against epitopes of the MHCI and MHCII cell surface receptors on human lymphoblast cells. Fluorescence anisotropy of donor and acceptor and FRET efficiency were measured for samples of different acceptor-to-donor concentration ratios. Acceptor anisotropy proved to be more sensitive than the donor anisotropy for sensing FRET. After determining the rotational constants of the donor-conjugated antibodies by measurements of FRET in the steady state, and by rFLIM as a reference, the associated fractions of the MHCI and MHCII molecules in their clusters were determined. Besides the flow cytometer and the wide-field rFLIM used in this study, the method can be applied also in other devices capable of dual-anisotropy detection.

Hydrodynamic modelling of protein conformation in solution: ELLIPS and HYDRO

The last three decades has seen some important advances in our ability to represent the conformation of proteins in solution on the basis of hydrodynamic measurements. Advances in theoretical modeling capabilities have been matched by commensurate advances in the precision of hydrodynamic measurements. We consider the advances in whole-body (simple ellipsoid-based) modeling—still useful for providing an overall idea of molecular shape, particularly for those systems where only a limited amount of data is available—and outline the ELLIPS suite of algorithms which facilitates the use of this approach. We then focus on bead modeling strategies, particularly the surface or shell–bead approaches and the HYDRO suite of algorithms. We demonstrate how these are providing great insights into complex issues such as the conformation of immunoglobulins and other multi-domain complexes.

Global fit and structure optimization of flexible and rigid macromolecules and nanoparticles from analytical ultracentrifugation and other dilute solution properties

The calculation of hydrodynamic and other solution properties from structural information (size and shape or flexibility) of macromolecules and nanoparticles is feasible thanks to existing theories and computational tools. Here we review our recent advances in the inverse problem of extracting structural information from those properties. The concepts of equivalent radii and ratios of radii are particularly useful in global-fitting structural analysis, when one has to treat simultaneously with various properties, eventually for a series of samples. Based on the equivalent radii or their ratios, we define target functions that measure the adequacy of a given structure to fit a set of experimental properties. Structural determination is carried out by minimization of those target functions. We review a variety of examples. Some of them refer to the simple, yet important models like ellipsoids, cylinders and wormlike chains, whose structure is determined by optimization of the model parameters. In other, more complex cases, properties are calculated with computational tools like programs in the HYDRO suite. We have devised other tools to make the structure optimization from the results of those calculations in a quite direct, simple and systematic manner.

Mechanism of GTPase-activity-induced self-assembly of human guanylate binding protein 1

Human guanylate binding protein 1 (hGBP1) belongs to the dynamin superfamily of large GTPases (LGs). In the course of GTP hydrolysis, the protein undergoes structural changes leading to self-assembly of the protein, which is a characteristic property of all family members. For self-assembly, the protein employs two distinct interaction sites, one of which is located within the LG domain of the protein located at the N-terminus, and the second is located in the C-terminal alpha-helical domain. Here, we identify intramolecular contacts between the LG domain and the helical part of hGBP1, which relay nucleotide-dependent structural changes from the N-terminus to the C-terminus and thereby mediate tetramer formation of the protein through a second contact site at the C-terminus. Furthermore, we demonstrate the impact of this intramolecular communication on the enzymatic activity of hGBP1 and on its cellular localization.

Effect of PEGylation on the Solution Conformation of Antibody Fragments

Covalent attachment of poly(ethylene glycol) (PEG) to therapeutic antibody fragments has been found effective in prolonging the half-life of the protein molecule in vivo. In this study analytical ultracentrifugation (AUC) in combination with small angle X-ray scattering (SAXS) has been applied to a number of antibody fragments and to their respective PEGylated conjugates. Despite the large increase in molecular weight due to the attachment of a 20-40 kDa PEG moiety, the PEGylated conjugates have smaller sedimentation coefficients, s, than their parent antibody fragments, due to a significant increase in frictional ratio f/f(o) (from approximately 1.3 to 2.3-2.8): the solution hydrodynamic properties of the conjugates are clearly dominated by the PEG moiety (f/f(o) approximately 3.0). This observation is reinforced by SAXS data at high values of r (separation of scattering centres within a particle) that appear dominated by the PEG part of the complex. By contrast, SAXS data at low values of r suggest that there are no significant conformational changes of the protein moiety itself after PEGylation The location of the PEGylation site within the conjugate was identified, and found to be consistent with expectation from the conjugation chemistry.

Solution conformation of wild-type and mutant IgG3 and IgG4 immunoglobulins using crystallohydrodynamics: possible implications for complement activation

We have employed the recently described crystallohydrodynamic approach to compare the time-averaged domain orientation of human chimeric IgG3wt (wild-type) and IgG4wt as well as two hinge mutants of IgG3 and an IgG4S331P (mutation from serine to proline at position 331, EU numbering) mutant of IgG4. The approach involves combination of the known shape of the Fab and Fc regions from crystallography with hydrodynamic data for the Fab and Fc fragments and hydrodynamic and small angle x-ray scattering data for the intact IgG structures. In this way, ad hoc assumptions over hydration can be avoided and model degeneracy (uniqueness problems) can be minimized. The best fit model for the solution structure of IgG3wt demonstrated that the Fab regions are directed away from the plane of the Fc region and with a long extended hinge region in between. The best fit model of the IgG3m15 mutant with a short hinge (and enhanced complement activation activity) showed a more open, but asymmetric structure. The IgG3HM5 mutant devoid of a hinge region (and also devoid of complement-activation activity) could not be distinguished at the low-resolution level from the structure of the enhanced complement-activating mutant IgG3m15. The lack of inter-heavy-chain disulphide bond rather than a significantly different domain orientation may be the reason for the lack of complement-activating activity of the IgG3HM5 mutant. With IgG4, there are significant and interesting conformational differences between the wild-type IgG4, which shows a symmetric structure, and the IgG4S331P mutant, which shows a highly asymmetric structure. This structural difference may explain the ability of the IgG4S331P mutant to activate complement in stark contrast to the wild-type IgG4 molecule which is devoid of this activity.

A dynamical study of antibody-antigen encounter reactions

The effects of internal dynamics in diffusion-driven encounters between macro-molecules represent a problem of broad relevance in molecular biology. In this view, we investigate a typical antigen-antibody reaction chain, based on a coarse-grained mechanical model parameterized directly upon results from single-molecule experiments. We demonstrate that the internal dynamics is a crucial factor in the encounter process. To describe our numerical results, we formulate a simple, intuitive theoretical framework, and we develop it analytically. This enables us to show that the inner dynamics of antibody molecules results in a cooperative behavior of their individual sub-units. Along the same lines, we also investigate the case of double binding to multi-valent antigens. Our results quantify the enhancement of avidity afforded by the double binding in excellent agreement with the available experimental data.

Oligomeric state of lipocalin-1 (LCN1) by multiangle laser light scattering and fluorescence anisotropy decay

Multiangle laser light scattering and fluorescence anisotropy decay measurements clarified the oligomeric states of native and recombinant tear lipocalin (lipocalin-1, TL). Native TL is monomeric. Recombinant TL (5-68 microM) with or without the histidine tag shows less than 7% dimer formation that is not in equilibrium with the monomeric form. Fluorescence anisotropy decay showed a correlation time of 9-10 ns for TL (10 microM-1 mM). Hydrodynamic calculations based on the crystallographic structure of a monomeric TL mutant closely concur with the observed correlation time. The solution properties calculated with HYDROPRO and SOLPRO programs from the available crystallographic structure of a monomeric TL mutant concur closely with the observed fluorescence anisotropy decay. The resulting model shows that protein topology is the major determinant of rotational correlation time and accounts for deviation from the Stokes-Einstein relation. The data challenge previous gel filtration studies to show that native TL exists predominantly as a monomer in solution rather than as a dimer. Delipidation of TL results in a formation of a complex oligomeric state (up to 25%). These findings are important as the dynamic processes in the tear film are limited by diffusional, translational as well as rotational, properties of the protein.

Equivalent radii and ratios of radii from solution properties as indicators of macromolecular conformation, shape, and flexibility

The equivalent radius for any solution property is the radius of a spherical particle having the same value of solution property as that of the macromolecule under consideration. Equivalent radii for different properties present a dependence on size and shape that are more similar than the values of the properties themselves. Furthermore, the ratios of equivalent radii of two properties depend on the conformation (shape or flexibility), but not on the absolute sizes. We define equivalent radii and their ratios, and describe their evaluation for some common models of rigid and flexible macromolecules. Using radii and ratios, we have devised procedures to fit macromolecular models to experimental properties, allowing the determination of the model parameters. Using these quantities, we can construct target functions for an equilibrated, unbiased optimization. The procedures, which have been implemented in public-domain computer programs, are illustrated for rigid, globular proteins, and the rodlike tobacco mosaic virus, and for semiflexible, wormlike heparin molecules.

Off-rate and concentration diversity in multidonor-derived dimers of immunoglobulin G

IgG-dimers in multidonor-derived preparations of IgG antibodies represent not only agents of therapeutic potential, but also molecules of basic immunological interest since their composition mirrors the currently unknown range of clonal concentrations and affinities. To analyze this fundamental type of diversity, a computational model is developed in agreement with a density functional theory and used to simulate the dissociation kinetics of dimers separated from a 5000 donor-derived IgG preparation (protein concentration: 0.74 mg/mL) via superimposition of 8100 arbitrary combinations of off-rates and initial concentrations. The Greedy algorithm-like procedure described requires iterative and consecutive changes of 8 from a total of 11 fitting parameters and allows to approximate the probability density distributions of either quantities within defined limits (apparent off-rates: approximately 4 x 10(-4) to 9 x 10(-17)s(-1); concentrations: approximately 3 x 10(-20) to 1 x 10(-11)M) by lognormal distributions of log-log(10)-type, each of them adapted with four particular parameters, as well as the number of different dimer populations ( approximately 2 x 10(13)). Moreover, reasonably dimensioned equilibrium constants involved in monovalent and bivalent random IgG dimerization are estimated by using a mean on-rate of 2.5 x 10(5)M(-1)s(-1) and interrelationships of molecular parameters derived from known models for antibody-antigen interaction.

Cooperative binding of the yeast Spt10p activator to the histone upstream activating sequences is mediated through an N-terminal dimerization domain

The yeast Spt10p activator is a putative histone acetyltransferase (HAT) possessing a sequence-specific DNA-binding domain (DBD) which binds to the upstream activation sequences (UAS elements) in the histone gene promoters. Spt10p binds to a pair of histone UAS elements with extreme positive cooperativity. The molecular basis of this cooperativity was addressed. Spt10p (640 residues) is an elongated dimer, but the isolated DBD (residues 283–396) is a monomer and binds non-cooperatively to DNA. A Spt10p fragment comprising the N-terminal domain (NTD), HAT domain and DBD (residues 1–396) binds cooperatively and is a dimer, whereas an overlapping Spt10p fragment comprising the DBD and C-terminal domains (residues 283–640) binds non-cooperatively and is a monomer. These observations imply that cooperative binding requires dimerization. The isolated NTD (residues 1–98) is a dimer and is responsible for dimerization. We propose that cooperativity involves a conformational change in the Spt10p dimer which facilitates the simultaneous recognition of two UAS elements. In vivo, deletion of the NTD results in poor growth, but does not prevent the binding at the HTA1 promoter, suggesting that dimerization is biologically important. Residues 1–396 are sufficient for normal growth, indicating that the critical functions of Spt10p reside in the N-terminal domains.

Crystallohydrodynamics of protein assemblies: Combining sedimentation, viscometry, and x-ray scattering

Crystallohydrodynamics describes the domain orientation in solution of antibodies and other multidomain protein assemblies where the crystal structures may be known for the domains but not the intact structure. The approach removes the necessity for an ad hoc assumed value for protein hydration. Previous studies have involved only the sedimentation coefficient leading to considerable degeneracy or multiplicity of possible models for the conformation of a given protein assembly, all agreeing with the experimental data. This degeneracy can be considerably reduced by using additional solution parameters. Conformation charts are generated for the three universal (i.e., size-independent) shape parameters P (obtained from the sedimentation coefficient or translational diffusion coefficient), nu (from the intrinsic viscosity), and G (from the radius of gyration), and calculated for a wide range of plausible orientations of the domains (represented as bead-shell ellipsoidal models derived from their crystal structures) and after allowance for any linker or hinge regions. Matches are then sought with the set of functions P, nu, and G calculated from experimental data (allowing for experimental error). The number of solutions can be further reduced by the employment of the D max parameter (maximum particle dimension) from x-ray scattering data. Using this approach we are able to reduce the degeneracy of possible solution models for IgG3 to a possible representative structure in which the Fab domains are directed away from the plane of the Fc domain, a structure in accord with the recognition that IgG3 is the most efficient complement activator among human IgG subclasses.

Numerical path integration technique for the calculation of transport properties of proteins

We present a new technique for the computation of both the translational diffusivity and the intrinsic viscosity of macromolecules, and apply it here to proteins. Traditional techniques employ finite element representations of the surface of the macromolecule, taking the surface to be a union of spheres or of polygons, and have computation times that are O(m(3)) where m is the number of finite elements. The new technique, a numerical path integration method, has computation times that are only O(m). We have applied the technique to approximately 1000 different protein structures. The computed translational diffusivities and intrinsic viscosities are, to lowest order, proportional respectively to N(-1/3)(R) and N(0)(R), where N(R) is the number of amino acid residues in the protein. Our calculations also show some correlation with the shape of the molecule, as represented by the ratio m(2)/m(3), where m(2) and m(3) are, respectively, the middle and the smallest of the three principal moments of inertia. Comparisons with a number of experimental results are also performed, with results generally consistent to within experimental error.

Estimating domain orientation of two human antibody IgG4 chimeras by crystallohydrodynamics

A modified crystallohydrodynamic approach introduced in 2001 is applied to two human IgG4 constructs from mouse IgG1. The constructs were point mutants of the chimeric antibody molecule cB72.3(gamma4): cB72.3(gamma4A), devoid of inter-chain disulfide bridging, and cB72.3(gamma4P), which has full inter-chain bridging. As before, the known crystallographic structures for the Fab and Fc domains were combined with the measured translational frictional ratios to obtain an estimate for the apparent time-averaged hydration of the domains and hence for that of the intact molecule. The original approach was modified with the hydrated dimensions of the domains being applied, rather than the anhydrous crystallographic dimensions, for assessing the inter-domain orientations using the algorithms HYDROSUB and SOLPRO. Both chimeric IgG4 molecules were found to have open, rather than compact, structures, in agreement with the previous study on wild-type human IgG4. The insertion of a frictionless connector between the domains was necessary, however, for representing the cB72.3(gamma4A) chimera. It therefore appears that the inter-chain disulfide bonds act as physical constraints in the cB72.3(gamma4P) chimera, forcing the antibody domains together and producing a less elongated structure than that of cB72.3(gamma4A). The open structures produced for the two IgG4 chimeras showed similarity to those structures identified for murine IgG1 and IgG2a molecules through X-ray crystallography.

O-glycan sialylation and the structure of the stalk-like region of the T cell co-receptor CD8

Studies of mucins suggest that the structural effects of O-glycans are restricted to steric interactions between peptide-linked GalNAc residues and adjacent polypeptide residues. It has been proposed, however, that differential O-glycan sialylation alters the structure of the stalk-like region of the T cell co-receptor, CD8, and that this, in turn, modulates ligand binding (Daniels, M. A., Devine, L., Miller, J. D., Moser, J. M., Lukacher, A. E., Altman, J. D., Kavathas, P., Hogquist, K. A., and Jameson, S. C. (2001) Immunity 15, 1051-1061; Moody, A. M., Chui, D., Reche, P. A., Priatel, J. J., Marth, J. D., and Reinherz, E. L. (2001) Cell 107, 501-512). We characterize the glycosylation of soluble, chimeric forms of the alphaalpha- and alphabeta-isoforms of murine CD8 containing the O-glycosylated stalk of rat CD8alphaalpha, and we show that the stalk O-glycans are differentially sialylated in CHO K1 versus Lec3.2.8.1 cells (82 versus approximately 6%, respectively). Sedimentation analysis indicates that the Perrin functions, Pexp, which reflect overall molecular shape, are very similar (1.61 versus 1.54), whereas the sedimentation coefficients (s) of the CHO K1- and Lec3.2.8.1-derived proteins differ considerably (3.73 versus 3.13 S). The hydrodynamic properties of molecular models also strongly imply that the sialylated and non-sialylated forms of the chimera have parallel, equally highly extended stalks ( approximately 2.6 A/residue). Our analysis indicates that, as in the case of mucins, the overall structure of O-glycosylated stalk-like peptides is sialylation-independent and that the functional effects of differential CD8 O-glycan sialylation need careful interpretation.

Orientation of a Y-shaped biomolecule adsorbed on a charged surface

The adsorption and orientation properties of a Y-shaped biomolecule, which models an immunoglobulin (Ig), on a charged surface are analyzed mesoscopically by Monte Carlo simulations. The orientation is a consequence of the interplay between van der Waals interactions and electrostatic interactions. For adsorption dominated by van der Waals attraction, the molecule prefers lying flat on the surface. For weak attraction, we observe a depletion zone in the concentration profile, which can result in a negative surface excess. A secondary peak is found for strong adsorption. For electrostatically dominated adsorption, the orientation is mainly determined by electric dipole and a vertically adsorbed molecule can be attained as it possesses strong electric dipole. Our study provides an explanation for experimental observations of preferential orientation.

Hydrodynamic properties of rigid macromolecules composed of ellipsoidal and cylindrical subunits

A procedure is devised for the calculation of hydrodynamic properties of rigid macromolecules composed subunits that are modeled as ellipsoids of revolution and cylinders. Owing to the axial symmetry of these shapes, smooth shell models can be constructured for the subunit structure. The bead shell model so constructed is employed for the calculation of the properties. A computer program, HYDROSUB, has been written implementing both the model building and the hydrodynamic calculation. A detailed example of the use of this methodology is presented for the case of the solution properties of the human antibody molecule immunoglobulin G3 (IgG3). Finally, hints are given on other uses and applications of the procedure.

Comparative investigations of biopolymer hydration by physicochemical and modeling techniques

The comparative investigation of biopolymer hydration by physicochemical techniques, particularly by small-angle X-ray scattering, has shown that the values obtained differ over a wide range, depending on the nature of the polymer and the environmental conditions. In the case of simple proteins, a large number of available data allow the derivation of a realistic average value for the hydration (0.35 g of water per gram of protein). As long as the average properties of proteins are considered, the use of such a default value is sufficient. Modeling approaches may be used advantageously, in order to differentiate between different assumptions and hydration contributions, and to correctly predict hydrodynamic properties of biopolymers on the basis of their three-dimensional structure. Problems of major concern are the positioning and the properties of the water molecules on the biopolymer surface. In this context, different approaches for calculating the molecular volume and surface of biopolymers have been applied, in addition to the development of appropriate hydration algorithms.

Crystallohydrodynamics for solving the hydration problem for multi-domain proteins: open physiological conformations for human IgG

Hydrodynamic methods provide a route for studying the low-resolution conformation--in terms of time-averaged spatial orientation of the Fab' and Fc domains relative to each other--of the human IgG subclasses, IgG1, IgG2, IgG3 and IgG4 in the environment in which many exist naturally---a solution. Representative modelling strategies are now available using 'shell-bead' or 'shell' modelling of the surface of the molecules with the size-independent programme SOLPRO [J. Garcia de la Torre, S.E. Harding, B. Carrasco, Eur. Biophys. J. 28 (1999) 119-132]. The shell model fits to the equivalent inertial surface ellipsoids of the published crystal structures for the Fab' and Fc domains of IgG are made and an apparent hydration delta(app) of 0.51g/g for Fab' and 0.70 g/g for the glycoprotein Fc are obtained, which yield an average value of (0.59+/-0.07) g/g for the intact antibody (2 Fab'+1 Fc). The relative orientations of these domains for each of the IgG subclasses is then found (using where appropriate a cylindrical hinge) from SOLPRO by modelling the Perrin function, P (i.e. 'frictional ratio due to shape') using this delta(app) and experimentally measured sedimentation coefficients. All the IgG subclasses appear as open, rather than compact structures with the degree of openness IgG3>IgG1>(IgG2, IgG4), with IgG3 and IgG1 non-coplanar. The hingeless mutant IgGMcg, with s degrees (20,w) approximately 6.8 S yields a coplanar structure rather similar to IgG2 and IgG4 and consistent with its crystallographic structure. The extension of this procedure for representing solution conformations of other antibody classes and other multi-domain proteins is indicated.

Signaling lymphocytic activation molecule (CDw150) is homophilic but self-associates with very low affinity

Signaling lymphocytic activating molecule ((SLAM) CDw150) is a glycoprotein that belongs to the CD2 subset of the immunoglobulin superfamily and is expressed on the surface of activated T- and B-cells. It has been proposed that SLAM is homophilic and required for bidirectional signaling during T- and B-cell activation. Previous work has suggested that the affinity of SLAM self-association might be unusually high, undermining the concept that protein interactions mediating transient cell-cell contacts, such as those involving leukocytes, have to be weak in order that such contacts are readily reversible. Using surface plasmon resonance-based methods and analytical ultracentrifugation (AUC), we confirm that SLAM is homophilic. However, we also establish a new theoretical treatment of surface plasmon resonance-derived homophilic binding data, which indicates that SLAM-SLAM interactions (solution K(d) approximately 200 micrometer) are in fact considerably weaker than most other well characterized protein-protein interactions at the cell surface (solution K(d) approximately 0.4-20 micrometer), a conclusion that is supported by the AUC analysis. Whereas further analysis of the AUC data imply that SLAM could form "head to head" dimers spanning adjacent cells, the very low affinity raises important questions regarding the physiological role and/or properties of such interactions.