Short-time transport properties in dense suspensions: From neutral to charge-stabilized colloidal spheres

Testing mixing rules for structural and dynamical quantities in multi-component crowded protein solutions

Crowding effects significantly influence the phase behavior and the structural and dynamic properties of the concentrated protein mixtures present in the cytoplasm of cells or in the blood serum. This poses enormous difficulties for our theoretical understanding and our ability to predict the behavior of these systems. While the use of course grained colloid-inspired models allows us to reproduce the key physical solution properties of concentrated monodisperse solutions of individual proteins, we lack corresponding theories for complex polydisperse mixtures. Here, we test the applicability of simple mixing rules in order to predict solution properties of protein mixtures. We use binary mixtures of the well-characterized bovine eye lens proteins α and γB crystallin as model systems. Combining microrheology with static and dynamic scattering techniques and observations of the phase diagram for liquid-liquid phase separation, we show that reasonably accurate descriptions are possible for macroscopic and mesoscopic signatures, while information on the length scale of the individual protein size requires more information on cross-component interaction.

Wavenumber-dependent dynamic light scattering optical coherence tomography measurements of collective and self-diffusion

We demonstrate wavenumber-dependent DLS-OCT measurements of collective and self-diffusion coefficients in concentrated silica suspensions across a broad q-range, utilizing a custom home-built OCT system. Depending on the sample polydispersity, either the collective or self-diffusion is measured. The measured collective-diffusion coefficient shows excellent agreement with hard-sphere theory and serves as an effective tool for accurately determining particle sizes. We employ the decoupling approximation for simultaneously measuring collective and self-diffusion coefficients, even in sufficiently monodisperse suspensions, using a high-speed Thorlabs OCT system. This enables particle size and volume fraction determination without the necessity of wavenumber-dependent measurements. We derive a relationship between the particle number-based polydispersity index and the ratio of self and collective mode amplitudes in the autocorrelation function and utilize it to measure the particle number-based polydispersity index. Notably, the polydispersity determined in this manner demonstrates improved sensitivity to smaller particle sizes compared to the standard intensity-based DLS cumulant analysis performed on dilute samples.

Combining Scattering Experiments and Colloid Theory to Characterize Charge Effects in Concentrated Antibody Solutions

Charges and their contribution to protein-protein interactions are essential for the key structural and dynamic properties of monoclonal antibody (mAb) solutions. In fact, they influence the apparent molecular weight, the static structure factor, the collective diffusion coefficient, or the relative viscosity, and their concentration dependence. Further, charges play an important role in the colloidal stability of mAbs. There exist standard experimental tools to characterize mAb net charges, such as the measurement of the electrophoretic mobility, the second virial coefficient, or the diffusion interaction parameter. However, the resulting values are difficult to directly relate to the actual overall net charge of the antibody and to theoretical predictions based on its known molecular structure. Here, we report the results of a systematic investigation of the solution properties of a charged IgG1 mAb as a function of concentration and ionic strength using a combination of electrophoretic measurements, static and dynamic light scattering, small-angle X-ray scattering, and tracer particle-based microrheology. We analyze and interpret the experimental results using established colloid theory and coarse-grained computer simulations. We discuss the potential and limits of colloidal models for the description of the interaction effects of charged mAbs, in particular pointing out the importance of incorporating shape and charge anisotropy when attempting to predict structural and dynamic solution properties at high concentrations.

Diffusing Wave Spectroscopy Measurements of Colloidal Suspension Dynamics

Diffusing wave spectroscopy (DWS) is used to measure the dynamics of charged silica particles between the volume fractions 0.065 ≤ ϕ ≤ 0.352 (weight percentages from 12.7 to 55.8 wt %). The short-time diffusivity averaged over the scattering vectors sampled by DWS D ¯ (ϕ) decreases with an increasing concentration. An effective hard-sphere model that accounts for hydrodynamic interactions and a double-layer repulsion fits the values up to an effective volume fraction ϕ e f f = ϕ b ^ 3 ≈ 0.6 , where b ^ is the excluded shell radius normalized by the particle radius b ^ = b/a = 1.3. While DWS measurements of diffusivity are sensitive to repulsive interactions, we show that they are relatively insensitive to attraction, such as those due to secondary minima in the interaction potential or weak depletion interaction.

Structure and short-time diffusion of concentrated suspensions consisting of silicone-stabilised PMMA particles: a quantitative analysis taking polydispersity effects into account

We characterise structure and dynamics of concentrated suspensions of silicone-stabilised PMMA particles immersed in index-matching decalin-tetralin mixtures by means of static and quasielastic light scattering experiments. These particles can reproducibly be prepared via a comparatively easy route and are thus promising model systems with hard-sphere interaction. We demonstrate the hard-sphere behaviour of dense suspensions of these systems rigorously taking polydispersity effects into account. Structure factors S(Q) can in the entire range of volume fractions with liquid-like structure quantitatively be modelled using a multi-component Percus-Yevick ansatz regarding the particle size distribution and the form factor assuming a core-shell model with a scattering length density gradient in the PMMA core. Herewith, hydrodynamic functions H(Q) are in the whole accessible Q-range beyond the second maximum of H(Q) quantitatively modelled using a rescaled δγ-approach for all investigated volume fractions. With these data, previously provided characterisation of dilute systems is extended: the excellent agreement of structural and dynamic properties with theoretical predictions for hard spheres demonstrates the suitability of these particles as a model system for hard spheres.

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.

Approaching the hard sphere limit in colloids suitable for confocal microscopy - the end of a decade lasting quest

PMMA-PHSA particles serve as the hard sphere model system since the 1980s. We investigate the fluid structure of fluorescent ones in three different solvents by laser scanning confocal microscopy: a decalin-tetrachloroethylene (TCE)-mixture and a decalin-cyclohexylbromide (CHB)-mixture with and without tetrabutylammoniumbromide (TBAB). The experimental 3D radial distribution functions are modeled by analytical theory and computer simulations taking polydispersity and the experimental position uncertainty into account. The quantitative comparison between experiment and simulation/theory establishes hard sphere like behavior for particles in decalin-TCE for a wide range of particle packing fractions. To the best of our knowledge, we present the first experimental dataset of a fluid structure that agrees convincingly with Percus-Yevick over a wide concentration range. Furthermore, charged sphere behavior is confirmed both for the decalin-CHB and the decalin-CHB-TBAB solvents, and it is demonstrated that a finite particle concentration reduces screening in the decalin-CHB-TBAB system compared to the bulk solvent.

Nanoscale Bending Dynamics in Mixed-Chain Lipid Membranes

Lipids that have two tails of different lengths are found throughout biomembranes in nature, yet the effects of this asymmetry on the membrane properties are not well understood, especially when it comes to the membrane dynamics. Here we study the nanoscale bending fluctuations in model mixed-chain 14:0–18:0 PC (MSPC) and 18:0–14:0 PC (SMPC) lipid bilayers using neutron spin echo (NSE) spectroscopy. We find that despite the partial interdigitation that is known to persist in the fluid phase of these membranes, the collective fluctuations are enhanced on timescales of tens of nanoseconds, and the chain-asymmetric lipid bilayers are softer than an analogous chain-symmetric lipid bilayer with the same average number of carbons in the acyl tails, di-16:0 PC (DPPC). Quantitative comparison of the NSE results suggests that the enhanced bending fluctuations at the nanosecond timescales are consistent with experimental and computational studies that showed the compressibility moduli of chain-asymmetric lipid membranes are 20% to 40% lower than chain-symmetric lipid membranes. These studies add to growing evidence that the partial interdigitation in mixed-chain lipid membranes is highly dynamic in the fluid phase and impacts membrane dynamic processes from the molecular to mesoscopic length scales without significantly changing the bilayer thickness or area per lipid.

Light scattering from mixtures of interacting, nonionic micelles with hydrophobic solutes

Model equations for the Rayleigh ratio and the electric field autocorrelation function are derived using thermodynamic fluctuation theory applied to crowded solute-containing micellar solutions and microemulsions with negligible molecular species and polydispersity. This theory invokes non-equilibrium thermodynamics and enforces local equilibrium between molecular solute, surfactant, and the various micellar species, in order to elucidate the influence of self-assembly on light scattering correlation functions. We find that self-assembly driven variations in the average micelle radius and aggregation number along gradients in concentration, which were previously shown to drive strong multicomponent diffusion effects expressed via the ternary diffusivity matrix [D], do not affect the scattering functions in the limit of zero local polydispersity. Hence, theoretical predictions for the Rayleigh ratio and the field autocorrelation function for ternary mixtures of solute-containing, locally monodisperse micellar solutions are identical to those developed for binary mixtures of monodisperse, colloidal hard spheres. However, self-assembly driven multicomponent diffusion phenomena are predicted to influence the thermodynamic driving forces for diffusion in these mixtures. In support of our theoretical results, measurements for the Rayleigh ratio and the field autocorrelation function for ternary aqueous solutions of decaethylene glycol monododecyl ether (C₁₂E₁₀) with either decane or limonene solute were performed for several molar ratios and volume fractions up to ϕ ≈ 0.25, and for binary mixtures of C₁₂E₁₀/water up to ϕ ≈ 0.5. Excellent agreement between our light scattering theory and experimental data is achieved for low to moderate volume fractions (ϕ < 0.3), and at higher concentrations when our theoretical results are corrected to account for micelle dehydration.

Short-Time Transport Properties of Bidisperse Suspensions of Immunoglobulins and Serum Albumins Consistent with a Colloid Physics Picture

The crowded environment of biological systems such as the interior of living cells is occupied by macromolecules with a broad size distribution. This situation of polydispersity might influence the dependence of the diffusive dynamics of a given tracer macromolecule in a monodisperse solution on its hydrodynamic size and on the volume fraction. The resulting size dependence of diffusive transport crucially influences the function of a living cell. Here, we investigate a simplified model system consisting of two constituents in aqueous solution, namely, of the proteins bovine serum albumin (BSA) and bovine polyclonal gamma-globulin (Ig), systematically depending on the total volume fraction and ratio of these constituents. From high-resolution quasi-elastic neutron spectroscopy, the separate apparent short-time diffusion coefficients for BSA and Ig in the mixture are extracted, which show substantial deviations from the diffusion coefficients measured in monodisperse solutions at the same total volume fraction. These deviations can be modeled quantitatively using results from the short-time rotational and translational diffusion in a two-component hard sphere system with two distinct, effective hydrodynamic radii. Thus, we find that a simple colloid picture well describes short-time diffusion in binary mixtures as a function of the mixing ratio and the total volume fraction. Notably, the self-diffusion of the smaller protein BSA in the mixture is faster than the diffusion in a pure BSA solution, whereas the self-diffusion of Ig in the mixture is slower than in the pure Ig solution.

Modeling photoacoustic pressure generation in colloidal suspensions at different volume fractions based on a multi-scale approach

Further development of quantitative photoacoustic tomography requires understanding the photoacoustic pressure generation by modeling the generation process. This study modeled the initial photoacoustic pressure in colloidal suspensions, used as tissue phantoms, at different volume fractions on a multi-scale approach. We modeled the thermodynamic and light scattering properties on a microscopic scale with/without treating the hard-sphere interaction between colloidal particles. Meanwhile, we did the light energy density on a macroscopic scale. We showed that the hard-sphere interaction significantly influences the initial pressure and related quantities at a high volume fraction except for the thermodynamic properties. We also showed the initial pressure at the absorber inside the medium logarithmically decreases with increasing the volume fractions. This result is mainly due to the decay of the light energy density with light scattering. Our numerical results suggest that modeling light scattering and propagation is crucial over modeling thermal expansion.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the d-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that d-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns–ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Neutron spectroscopy, dynamic light scattering, X-ray diffraction, and MD-simulations were used to investigate the effect of ligand binding on the structure and diffusive dynamics of Escherichia coli aspartate alpha-decarboxylase.

Dynamical Behavior of Disordered Regions in Disease-Related Proteins Revealed by Quasielastic Neutron Scattering

Background and Objectives: Intrinsically disordered proteins (IDPs) and proteins containing intrinsically disordered regions (IDRs) are known to be involved in various human diseases. Since the IDPs/IDRs are fluctuating between many structural substrates, the dynamical behavior of the disease-related IDPs/IDRs needs to be characterized to elucidate the mechanisms of the pathogenesis of the diseases. As protein motions have a hierarchy ranging from local side-chain motions, through segmental motions of loops or disordered regions, to diffusive motions of entire molecules, segmental motions, as well as local motions, need to be characterized. Materials and Methods: Combined analysis of quasielastic neutron scattering (QENS) spectra with the structural data provides information on both the segmental motions and the local motions of the IDPs/IDRs. Here, this method is applied to re-analyze the QENS spectra of the troponin core domain (Tn-CD), various mutants of which cause the pathogenesis of familial cardiomyopathy (FCM), and α-synuclein (αSyn), amyloid fibril formation of which is closely related to the pathogenesis of Parkinson's disease, collected in the previous studies. The dynamical behavior of wild-type Tn-CD, FCM-related mutant Tn-CD, and αSyn in the different propensity states for fibril formation is characterized. Results: In the Tn-CD, the behavior of the segmental motions is shown to be different between the wild type and the mutant. This difference is likely to arise from changes in the intramolecular interactions, which are suggested to be related to the functional aberration of the mutant Tn-CD. In αSyn, concerted enhancement of the segmental motions and the local motions is observed with an increased propensity for fibril formation, suggesting the importance of these motions in fibril formation. Conclusions: Characterization of the segmental motions as well as the local motions is thus useful for discussing how the changes in dynamical behavior caused by the disease-related mutations and/or environmental changes could be related to the functional and/or behavioral aberrations of these proteins.

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions

We present guidelines to estimate the effect of electrostatic repulsion in sedimenting dilute particle suspensions. Our results are based on combined Langevin dynamics and lattice Boltzmann simulations for a range of particle radii, Debye lengths and particle concentrations. They show a simple relationship between the slope K of the concentration-dependent sedimentation velocity and the range χ of the electrostatic repulsion normalized by the average particle-particle distance. When χ → 0, the particles are too far away from each other to interact electrostatically and K = 6.55 as predicted by the theory of Batchelor. As χ increases, K likewise increases as if the particle radius increased in proportion to χ up to a maximum around χ = 0.4. Over the range χ = 0.4-1, K relaxes exponentially to a concentration-dependent constant consistent with known results for ordered particle distributions. Meanwhile the radial distribution function transitions from a disordered gas-like to a liquid-like form. Power law fits to the concentration-dependent sedimentation velocity similarly yield a simple master curve for the exponent as a function of χ, with a step-like transition from 1 to 1/3 centered around χ = 0.6.

Self-Diffusive Properties of the Intrinsically Disordered Protein Histatin 5 and the Impact of Crowding Thereon: A Combined Neutron Spectroscopy and Molecular Dynamics Simulation Study

Intrinsically disordered proteins (IDPs) are proteins that, in comparison with globular/structured proteins, lack a distinct tertiary structure. Here, we use the model IDP, Histatin 5, for studying its dynamical properties under self-crowding conditions with quasi-elastic neutron scattering in combination with full atomistic molecular dynamics (MD) simulations. The aim is to determine the effects of crowding on the center-of-mass diffusion as well as the internal diffusive behavior. The diffusion was found to decrease significantly, which we hypothesize can be attributed to some degree of aggregation at higher protein concentrations, (≥100 mg/mL), as indicated by recent small-angle X-ray scattering studies. Temperature effects are also considered and found to, largely, follow Stokes–Einstein behavior. Simple geometric considerations fail to accurately predict the rates of diffusion, while simulations show semiquantitative agreement with experiments, dependent on assumptions of the ratio between translational and rotational diffusion. A scaling law that previously was found to successfully describe the behavior of globular proteins was found to be inadequate for the IDP, Histatin 5. Analysis of the MD simulations show that the width of the distribution with respect to diffusion is not a simplistic mirroring of the distribution of radius of gyration, hence, displaying the particular features of IDPs that need to be accounted for.

Non-Gaussian, transiently anomalous, and ergodic self-diffusion of flexible dumbbells in crowded two-dimensional environments: Coupled translational and rotational motions

We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction ϕ and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on ϕ, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing ϕ. We also analyze the relative motion in the dimers, finding that larger ϕ suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(ϕ) of translational and rotational motion of the dumbbells an exponential decay with ϕ for weak and a power-law variation D(ϕ)∝(ϕ-ϕ^{★})^{2.4} for strong crowding is found. A comparison of simulation results with theoretical predictions for D(ϕ) is discussed and some relevant experimental systems are overviewed.

Diffusion and sedimentation in colloidal suspensions using multiparticle collision dynamics with a discrete particle model

We study self-diffusion and sedimentation in colloidal suspensions of nearly hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to experimental data and Brownian dynamics simulations with free-draining hydrodynamics (BD) as well as pairwise far-field hydrodynamics described using the Rotne-Prager-Yamakawa mobility tensor (BD+RPY). The dynamics in MD+MPCD suggest that the colloidal particles are only partially coupled to the solvent at short times. However, the long-time self-diffusion coefficient in MD+MPCD is comparable to that in experiments, and the sedimentation coefficient in MD+MPCD is in good agreement with that in experiments and BD+RPY, suggesting that MD+MPCD gives a reasonable description of hydrodynamic interactions in colloidal suspensions. The discrete-particle MD+MPCD approach is convenient and readily extended to more complex shapes, and we determine the long-time self-diffusion coefficient in suspensions of nearly hard cubes to demonstrate its generality.

Geometrical Influence on Particle Transport in Cross-Flow Ultrafiltration: Cylindrical and Flat Sheet Membranes

Cross-flow membrane ultrafiltration (UF) is used for the enrichment and purification of small colloidal particles and proteins. We explore the influence of different membrane geometries on the particle transport in, and the efficiency of, inside-out cross-flow UF. For this purpose, we generalize the accurate and numerically efficient modified boundary layer approximation (mBLA) method, developed in recent work by us for a hollow cylindrical membrane, to parallel flat sheet geometries with one or two solvent-permeable membrane sheets. Considering a reference dispersion of Brownian hard spheres where accurate expressions for its transport properties are available, the generalized mBLA method is used to analyze how particle transport and global UF process indicators are affected by varying operating parameters and the membrane geometry. We show that global process indicators including the mean permeate flux, the solvent recovery indicator, and the concentration factor are strongly dependent on the membrane geometry. A key finding is that irrespective of the many input parameters characterizing an UF experiment and its membrane geometry, the process indicators are determined by three independent dimensionless variables only. This finding can be very useful in the design, optimization, and scale-up of UF processes.

Microfluidic osmotic compression of a charge-stabilized colloidal dispersion: Equation of state and collective diffusion coefficient

We show, using a model coupling mass transport and liquid theory calculations for a charge-stabilized colloidal dispersion, that diffusion significantly limits measurement times of its equation of state (EOS), osmotic pressure vs composition, using the osmotic compression technique. Following this result, we present a microfluidic chip allowing one to measure the entire EOS of a charged dispersion at the nanoliter scale in a few hours. We also show that time-resolved analyses of relaxation to equilibrium in this microfluidic experiment lead to direct estimates of the collective diffusion coefficient of the dispersion in Donnan equilibrium with a salt reservoir.

Roughness induced rotational slowdown near the colloidal glass transition

HYPOTHESIS

In concentrated suspensions, the dynamics of colloids are strongly influenced by the shape and topographical surface characteristics of the particles. As the particles get into close proximity, surface roughness alters the translational and rotational Brownian motions in different ways. Eventually, the rotations will get frustrated due to geometric hindrance from interacting asperities.

EXPERIMENTS

We use model raspberry-like colloids to study the effect of roughness on the translational and rotational dynamics. Using Confocal Scanning Laser Microscopy and particle tracking, we simultaneously resolve the two types of Brownian motion and obtain the corresponding Mean Squared Displacements for varying concentrations up to the maximum packing fraction.

FINDINGS

Roughness not only lowers the concentration of the translational colloidal glass transition, but also generates a broad concentration range in which the rotational Brownian motion changes signature from high-amplitude diffusive to low-amplitude rattling. This hitherto not reported second glass transition for rough spherical colloids emerges when the particle intersurface distance becomes comparable to the roughness length scale. Our work provides a unifying understanding of the surface characteristics' effect on the rotational dynamics during glass formation and provides a microscopic foundation for many roughness-related macroscale phenomena in nature and technology.

Temperature and salt controlled tuning of protein clusters

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (c_p) and the concentrations of the trivalent salts YCl₃ and LaCl₃ (c_s). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on c_p, c_s and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(c_p, c_s, T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

Microsecond hydrodynamic interactions in dense colloidal dispersions probed at the European XFEL

Many soft-matter systems are composed of macromolecules or nanoparticles suspended in water. The characteristic times at intrinsic length scales of a few nanometres fall therefore in the microsecond and sub-microsecond time regimes. With the development of free-electron lasers (FELs) and fourth-generation synchrotron light-sources, time-resolved experiments in such time and length ranges will become routinely accessible in the near future. In the present work we report our findings on prototypical soft-matter systems, composed of charge-stabilized silica nanoparticles dispersed in water, with radii between 12 and 15 nm and volume fractions between 0.005 and 0.2. The sample dynamics were probed by means of X-ray photon correlation spectroscopy, employing the megahertz pulse repetition rate of the European XFEL and the Adaptive Gain Integrating Pixel Detector. We show that it is possible to correctly identify the dynamical properties that determine the diffusion constant, both for stationary samples and for systems driven by XFEL pulses. Remarkably, despite the high photon density the only observable induced effect is the heating of the scattering volume, meaning that all other X-ray induced effects do not influence the structure and the dynamics on the probed timescales. This work also illustrates the potential to control such induced heating and it can be predicted with thermodynamic models.

Machine learning of lubrication correction based on GPR for the coupled DPD-DEM simulation of colloidal suspensions

Hydrodynamic interactions have a major impact on the suspension properties, but they are absent in atomic and molecular fluids due to a lack of intervening medium at close range. To reproduce the correct hydrodynamic interactions, lubrication correction is essential to compensate the missing short-range hydrodynamics from the fluids. However, lubrication correction requires many simulations in particle-based simulations of colloidal suspensions. To address the problem, we employ an active learning strategy based on Gaussian process regression (GPR) for normal and tangential lubrication corrections to significantly reduce the number of necessary simulations and apply the correction to the coupled multiscale simulation of monodisperse hard-sphere colloidal suspensions. In particular, a single-particle dissipative particle dynamics (DPD) model with parameter correction is used to describe the solvent-solvent and colloid-solvent interactions, and a discrete element method (DEM) model to depict the colloid-colloid frictional contacts. The lubrication correction results demonstrate that only six and four independent simulations (observation points for GPR training) are required to achieve accurate normal and tangential lubrication corrections, respectively. To validate the machine learning of lubrication correction based on GPR, we investigate the self-diffusion coefficients of colloids, suspension rheology and microstructure using the coupled DPD-DEM model with GPR lubrication correction. Our simulation results show that the machine learning of lubrication correction based on GPR is effective and the lubrication corrected DPD-DEM model is indeed capable of accurately capturing hydrodynamic interactions and correctly reproducing dynamical and rheological properties of colloidal suspensions. Moreover, the machine learning of lubrication correction based on GPR is not limited to the coupled DPD-DEM simulation of colloidal suspensions presented here, but can be easily applied to other particle-based simulations of particulate suspensions.

Probing the balance of Telechelic Polymers Bridging and Screened-Coulomb Interactions in Microemulsion System

The structural properties of neutral oil-in-water (O/W) microemulsions nanodroplets in the presence of telechelic polymers PEO modified with two hydrophobic end functions (PEO- 2m) were studied by Small-Angle Neutron Scattering (SANS). We show that the unknown attractive interaction introduced by the addition of PEO-2m, can could be balanced by a known repulsion interaction due to the addition of a cationic surfactant Cp +, which causes a progressive transformation of the bridges connecting the nanodroplets into loops. For this, we compute the structure factor S(q) with an effective pair potential, using the Ornstein Zernicke (OZ) integral equation approach with the Hypernetted Chain (HNC) closure relation. Generally, the agreement between the neutron scattering spectra and the numerical study is reasonable and allows a detailed description for each sample.

Self-diffusion of nonspherical particles fundamentally conflicts with effective sphere models

Modeling diffusion of nonspherical particles presents an unsolved and considerable challenge, despite its importance for the understanding of crowding effects in biology, food technology and formulation science. A common approach in experiment and simulation is to map nonspherical objects on effective spheres to subsequently use the established predictions for spheres to approximate phenomena for nonspherical particles. Using numerical evaluation of the hydrodynamic mobility tensor, we show that this so-called effective sphere model fundamentally fails to represent the self-diffusion in solutions of ellipsoids as well as rod-like assemblies of spherical beads. The effective sphere model drastically overestimates the slowing down of self-diffusion down to volume fractions below 0.01. Furthermore, even the linear term relevant at lower volume fraction is inaccurate, linked to a fundamental misconception of effective sphere models. To overcome the severe problems related with the use of effective sphere models, we suggest a protocol to predict the short-time self-diffusion of rod-like systems, based on simulations with hydrodynamic interactions that become feasible even for more complex molecules as the essential observable shows a negligible system-size effect.

Multicomponent diffusion of interacting, nonionic micelles with hydrophobic solutes

Ternary diffusion coefficient matrices [D] were measured using the Taylor dispersion method, for crowded aqueous solutions of decaethylene glycol monododecyl ether (C₁₂E₁₀) with either decane or limonene solute. The matrix [D], for both systems, was found to be highly non-diagonal, and concentration dependent, over a broad domain of solute to surfactant molar ratios and micelle volume fractions. A recently developed theoretical model, based on Batchelor's theory for gradient diffusion in dilute, polydisperse mixtures of interacting spheres, was simplified by neglecting local polydispersity, and effectively used to predict [D] with no adjustable parameters. Even though the model originates from dilute theory, the theoretical results were in surprisingly good agreement with experimental data for concentrated mixtures, with volume fractions up to φ≈ 0.47. In addition, the theory predicts eigenvalues D_- and D₊ that correspond to long-time self and gradient diffusion coefficients, respectively, for monodisperse spheres, in reasonable agreement with experimental data.

Collective diffusion coefficient of a charged colloidal dispersion: interferometric measurements in a drying drop

In the present work, we use Mach-Zehnder interferometry to thoroughly investigate the drying dynamics of a 2D confined drop of a charged colloidal dispersion. This technique makes it possible to measure the colloid concentration field during the drying of the drop at a high accuracy (about 0.5%) and with a high temporal and spatial resolution (about 1 frame per s and 5 μm per pixel). These features allow us to probe mass transport of the charged dispersion in this out-of-equilibrium situation. In particular, our experiments provide the evidence that mass transport within the drop can be described by a purely diffusive process for some range of parameters for which the buoyancy-driven convection is negligible. We are then able to extract from these experiments the collective diffusion coefficient of the dispersion D(φ) over a wide concentration range φ = 0.24-0.5, i.e. from the liquid dispersed state to the solid glass regime, with a high accuracy. The measured values of D(φ) ≃ 5-12D0 are significantly larger than the simple estimate D0 given by the Stokes-Einstein relation, thus highlighting the important role played by the colloidal interactions in such dispersions.

A neutron scattering perspective on the structure, softness and dynamics of the ligand shell of PbS nanocrystals in solution

Small-angle neutron and X-ray scattering, neutron backscattering and neutron time-of-flight spectroscopy are applied to reveal the structure of the ligand shell, the temperature-dependent diffusion properties and the phonon spectrum of PbS nanocrystals functionalized with oleic acid in deuterated hexane. The nanocrystals decorated with oleic acid as well as the desorbed ligand molecules exhibit simple Brownian diffusion with a Stokes-Einstein temperature-dependence and inhibited freezing. Ligand molecules desorbed from the surface show strong spatial confinement. The phonon spectrum of oleic acid adsorbed to the nanocrystal surface exhibits hybrid modes with a predominant Pb-character. Low-energy surface modes of the NCs are prominent and indicate a large mechanical softness in solution. This work provides comprehensive insights into the ligand-particle interaction of colloidal nanocrystals in solution and highlights its effect on the diffusion and vibrational properties as well as their mechanical softness.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.

Anisotropic dynamics and kinetic arrest of dense colloidal ellipsoids in the presence of an external field studied by differential dynamic microscopy

Anisotropic dynamics on the colloidal length scale is ubiquitous in nature. Of particular interest is the dynamics of systems approaching a kinetically arrested state. The failure of classical techniques for investigating the dynamics of highly turbid suspensions has contributed toward the limited experimental information available up until now. Exploiting the recent developments in the technique of differential dynamic microscopy (DDM), we report the first experimental study of the anisotropic collective dynamics of colloidal ellipsoids with a magnetic hematite core over a wide concentration range approaching kinetic arrest. In addition, we have investigated the effect of an external magnetic field on the resulting anisotropic collective diffusion. We combine DDM with small-angle x-ray scattering and rheological measurements to locate the glass transition and to relate the collective short- and long-time diffusion coefficients to the structural correlations and the evolution of the zero shear viscosity as the system approaches an arrested state.

Self-diffusion of a highly concentrated monoclonal antibody by fluorescence correlation spectroscopy: insight into protein-protein interactions and self-association

The dynamic behavior of monoclonal antibodies (mAbs) at high concentration provides insight into protein microstructure and protein-protein interactions (PPI) that influence solution viscosity and protein stability. At high concentration, interpretation of the collective-diffusion coefficient Dc, as determined by dynamic light scattering (DLS), is highly challenging given the complex hydrodynamics and PPI at close spacings. In contrast, self-diffusion of a tracer particle by Brownian motion is simpler to understand. Herein, we develop fluorescence correlation spectroscopy (FCS) for the measurement of the long-time self-diffusion of mAb2 over a wide range of concentrations and viscosities in multiple co-solute formulations with varying PPI. The normalized self-diffusion coefficient D0/Ds (equal to the microscopic relative viscosity ηeff/η0) was found to be smaller than η/η0. Smaller ratios of the microscopic to macroscopic viscosity (ηeff/η) are attributed to a combination of weaker PPI and less self-association. The interaction parameters extracted from fits of D0/Ds with a length scale dependent viscosity model agree with previous measurements of PPI by SLS and SAXS. Trends in the degree of self-association, estimated from ηeff/η with a microviscosity model, are consistent with oligomer sizes measured by SLS. Finally, measurements of collective diffusion and osmotic compressibility were combined with FCS data to demonstrate that the changes in self-diffusion between formulations are due primarily to changes in the protein-protein friction in these systems, and not to protein-solvent friction. Thus, FCS is a robust and accessible technique for measuring mAb self-diffusion, and, by extension, microviscosity, PPI and self-association that govern mAb solution dynamics.

Structure and dynamics in suspensions of soft core-shell colloids in the fluid regime

We report on a detailed experimental study of the structure and short-time dynamics in fluid-regime suspensions of soft core-shell spherical particles with different molecular weights of the chains forming the soft outer shell, and therefore different degrees of particle softness, using 3D dynamic light scattering (3D-DLS). Owing to the particle softness, the liquid-crystal coexistence regime is found to be broader than that of hard-sphere (HS) suspensions. Static light scattering in the dilute regime yields form factors that can be described using a spherical core-shell model and second virial coefficients A₂ > 0 indicative of purely repulsive interactions. The particle-particle interactions are longer ranged for all considered systems except those of the smaller molecular weight chain grafted particles which show a HS-like behavior. 3D-DLS experiments in the concentrated regime up to the liquid-crystal transition provide the short-time diffusion function, D(q), in a broad range of scattering wavenumbers, q, from which the structural (cage) and short-time self-diffusion coefficients D(q_m) and D_S = D(q ≫ q_m), respectively, are deduced as functions of the effective particle volume fraction, ϕ = c/c^*, where c^* is the overlap concentration, calculated using the hydrodynamic particle radius, R_H. The size of the nearest-neighbor cage of particles is characterized by 2π/q_m, with D(q) and the static structure factor S(q) attaining at q_m the smallest and largest values, respectively. Experimental data of D(q_m) and D_S are contrasted with analytic theoretical predictions based on a simplifying hydrodynamic radius model where the internal hydrodynamic structure of the core-shell particles is mapped on a single hydrodynamic radius parameter γ = R_H/R_eff, for constant direct interactions characterized by an (effective) hard-core radius R_eff. The particle softness is reflected, in particular, in the corresponding shape of the static structure factor, while the mean solvent (Darcy) permeability of the particles related to γ is reflected in the dynamic properties only. For grafted particles with longer polymer chains, D(q_m) and D_S are indicative of larger permeability values while particles with shorter chains are practically nonpermeable. The particle softness is also evident in the effective random close packing fraction estimated from the extrapolated zero-value limit of the cage diffusion coefficient D(q_m).

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

A Colloid Approach to Self-Assembling Antibodies

Concentrated solutions of monoclonal antibodies have attracted considerable attention due to their importance in pharmaceutical formulations; yet, their tendency to aggregate and the resulting high viscosity pose considerable problems. Here we tackle this problem by a soft condensed matter physics approach, which combines a variety of experimental measurements with a patchy colloid model, amenable of analytical solution. We thus report results of structural antibodies and dynamic properties obtained through scattering methods and microrheological experiments. We model the data using a colloid-inspired approach, explicitly taking into account both the anisotropic shape of the molecule and its charge distribution. Our simple patchy model is able to disentangle self-assembly and intermolecular interactions and to quantitatively describe the concentration-dependence of the osmotic compressibility, collective diffusion coefficient, and zero shear viscosity. Our results offer new insights on the key problem of antibody formulations, providing a theoretical and experimental framework for a quantitative assessment of the effects of additional excipients or chemical modifications and a prediction of the resulting viscosity.

Study of network composition in interpenetrating polymer networks of poly(N isopropylacrylamide) microgels: The role of poly(acrylic acid)

HYPOTHESIS

The peculiar swelling behaviour of poly(N-isopropylacrylamide) (PNIPAM)⁴-based responsive microgels provides the possibility to tune both softness and volume fraction with temperature, making these systems of great interest for technological applications and theoretical implications. Their intriguing phase diagram can be even more complex if poly(acrylic acid) (PAAc)⁵ is interpenetrated within PNIPAM network to form Interpenetrating Polymer Network (IPN)⁶ microgels that exhibit an additional pH-sensitivity. The effect of the PAAc/PNIPAM polymeric ratio on both swelling capability and dynamics is still matter of investigation.

EXPERIMENTS

Here we investigate the role of PAAc in the behaviour of IPN microgels across the volume phase transition through dynamic light scattering (DLS),⁷ transmission electron microscopy (TEM)⁸ and electrophoretic measurements as a function of microgel concentration and pH.

FINDINGS

Our results highlight that aggregation is favored at increasing weight concentration, PAAc content and pH and that a crossover PAAc content C_PAAc^∗9 exists above which the ionic charges on the microgel become relevant. Moreover we show that the softness of IPN microgels can be tuned ad hoc by changing the PAAc/PNIPAM ratio. These findings provide new insights into the possibility to control experimentally aggregation properties, charge and softness of IPN microgels by varying PAAc content.

Protein-Protein Interactions of Highly Concentrated Monoclonal Antibody Solutions via Static Light Scattering and Influence on the Viscosity

The ability to design and formulate mAbs to minimize attractive interactions at high concentrations is important for protein processing, stability, and administration, particularly in subcutaneous delivery, where high viscosities are often challenging. The strength of protein-protein interactions (PPIs) of an IgG1 and IgG4 monoclonal antibody (mAb) from low to high concentration was determined by static light scattering (SLS) and used to understand viscosity data. The PPI were tuned using NaCl and five organic ionic co-solutes. The PPI strength was quantified by the normalized structure factor S(0)/ S(0)_HS and Kirkwood-Buff integral G₂₂/ G_22,HS (HS = hard sphere) determined from the SLS data and also by fits with (1) a spherical Yukawa potential and (2) an interacting hard sphere (IHS) model, which describes attraction in terms of hypothetical oligomers. The IHS model was better able to capture the scattering behavior of the more strongly interacting systems (mAb and/or co-solute) than the spherical Yukawa potential. For each descriptor of PPI, linear correlations were obtained between the viscosity at high concentration (200 mg/mL) and the interaction strengths evaluated both at low (20 mg/mL) and high concentrations (200 mg/mL) for a given mAb. However, the only parameter that provided a correlation across both mAbs was the oligomer mass ratio ( m_oligomer/ m_{monomer+dimer}) from the IHS model, indicating the importance of self-association (in addition to the direct influence of the attractive PPI) on the viscosity.

Weak Shape Anisotropy Leads to a Nonmonotonic Contribution to Crowding, Impacting Protein Dynamics under Physiologically Relevant Conditions

The effect of a nonspherical particle shape on the dynamics in crowded solutions presents a significant challenge for a comprehensive understanding of interaction and structural relaxation in biological and soft matter. We report that small deviations from a spherical shape induce a nonmonotonic contribution to the crowding effect on the short-time cage diffusion compared with spherical systems, using molecular dynamics simulations with mesoscale hydrodynamics of a multiparticle collision dynamics fluid in semidilute systems with volume fractions smaller than 0.35. We show that the nonmonotonic effect due to anisotropy is caused by the combination of a reduced relative mobility over the entire concentration range and a looser and less homogeneous cage packing of nonspherical particles. Our finding stresses that nonsphericity induces new complexity, which cannot be accounted for in effective sphere models, and is of great interest in applications such as formulations as well as for the fundamental understanding of soft matter in general and crowding effects in living cells in particular.

Phase Behavior of Bowl-Shaped Colloids: Order and Dynamics in Plastic Crystals and Glasses

Charged fluorescent bowl-shaped colloids consisting of a polystyrene core surrounded by a poly(N-isopropylmethacrylamide) shell are obtained by nanoengineering spherical composite microgels. The phase diagram of these soft bowl-shaped colloids interacting through long-range Yukawa-type interactions is investigated using confocal laser scanning microscopy. The bowl-shaped structure leads to marked differences in phase-behavior compared to their spherical counterpart. With increasing number density, a transition from a fluid to a plastic crystal phase, with freely rotating particles, followed by a glass-like state is observed. It is found that the anisotropic bowl shape frustrates crystallization and slows down crystallization kinetics and causes the glass-like transition to shift to a significantly lower volume fraction than for the spheres. Quantitative analysis of the positional and orientational order demonstrates that the plastic crystal phase exhibits quasi-long range translational order and orientational disorder, while in the disordered glass-like phase the long-range translational order vanishes and short-range rotational order appears, dictated by the specific bowl shape. It is further shown that the different structural transitions are characterized by decoupling of the translational and orientational dynamics.

Nanosecond Tracer Diffusion as a Probe of the Solution Structure and Molecular Mobility of Protein Assemblies: The Case of Ovalbumin

Protein diffusion is not only an important process ensuring biological function but can also be used as a probe to obtain information on structural properties of protein assemblies in liquid solutions. Here, we explore the oligomerization state of ovalbumin at high protein concentrations by means of its short-time self-diffusion. We employ high-resolution incoherent quasielastic neutron scattering to access the self-diffusion on nanosecond timescales, on which interparticle contacts are not altered. Our results indicate that ovalbumin in aqueous (D₂O) solutions occurs in increasingly large assemblies of its monomeric subunits with rising protein concentration. It changes from nearly monomeric toward dimeric and ultimately larger than tetrameric complexes. Simultaneously, we access information on the internal molecular mobility of ovalbumin on the nanometer length scale and compare it with results obtained for bovine serum albumin, immunoglobulin, and β-lactoglobulin.

Short- and long-time diffusion and dynamic scaling in suspensions of charged colloidal particles

We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in concentrated suspensions of charge-stabilized colloidal spheres. In theory and simulation, the spheres are assumed to interact directly by a hard-core plus screened Coulomb effective pair potential. The intermediate scattering function, f_c(q, t), is calculated by elaborate accelerated Stokesian dynamics (ASD) simulations for Brownian systems where many-particle hydrodynamic interactions (HIs) are fully accounted for, using a novel extrapolation scheme to a macroscopically large system size valid for all correlation times. The study spans the correlation time range from the colloidal short-time to the long-time regime. Additionally, Brownian Dynamics (BD) simulation and mode-coupling theory (MCT) results of f_c(q, t) are generated where HIs are neglected. Using these results, the influence of HIs on collective and self-diffusion and the accuracy of the MCT method are quantified. It is shown that HIs enhance collective and self-diffusion at intermediate and long times. At short times self-diffusion, and for wavenumbers outside the structure factor peak region also collective diffusion, are slowed down by HIs. MCT significantly overestimates the slowing influence of dynamic particle caging. The dynamic scattering functions obtained in the ASD simulations are in overall good agreement with our dynamic light scattering (DLS) results for a concentration series of charged silica spheres in an organic solvent mixture, in the experimental time window and wavenumber range. From the simulation data for the time derivative of the width function associated with f_c(q, t), there is indication of long-time exponential decay of f_c(q, t), for wavenumbers around the location of the static structure factor principal peak. The experimental scattering functions in the probed time range are consistent with a time-wavenumber factorization scaling behavior of f_c(q, t) that was first reported by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)] for suspensions of hard spheres. Our BD simulation and MCT results predict a significant violation of exact factorization scaling which, however, is approximately restored according to the ASD results when HIs are accounted for, consistent with the experimental findings for f_c(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, f_s(q, t), and its non-Gaussian parameter α₂(t) and for the particle mean squared displacement W(t) and its time derivative. Since self-diffusion properties are not assessed in standard DLS measurements, a method to deduce W(t) approximately from f_c(q, t) is theoretically validated.

Hard-sphere-like dynamics in highly concentrated alpha-crystallin suspensions

The dynamics of concentrated suspensions of the eye-lens protein alpha crystallin have been measured using x-ray photon correlation spectroscopy. Measurements were made at wave vectors corresponding to the first peak in the hard-sphere structure factor and volume fractions close to the critical volume fraction for the glass transition. Langevin dynamics simulations were also performed in parallel to the experiments. The intermediate scattering function f(q,τ) could be fit using a stretched exponential decay for both experiments and numerical simulations. The measured relaxation times show good agreement with simulations for polydisperse hard-sphere colloids.