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.


I. INTRODUCTION
Macromolecular crowding is a ubiquitous and relevant factor in biological fluids such as the cytoplasm and blood serum, expressing the seemingly trivial finding that macromolecules in natural and biotechnological solutions are affected by the presence of other molecules. 1,2The kinetic and thermodynamic consequences such as varied reaction rates and equilibria 3 and lowered solution stability 4,5 govern largely protein diffusion and interaction, as well as the phase behavior and viscosity of protein solutions.For globular proteins, a colloidal picture has proven highly useful for an in-depth and quantitative understanding of concentrated solutions of proteins, which is of high relevance for understanding of function and misfunction in protein solutions, related in particular, to diseases caused by protein condensation. 6,7Particular examples include the liquid-liquid phase behavior, [8][9][10][11][12] cluster formation [13][14][15][16][17] as well as the slowing down of selfand cage diffusion, [18][19][20][21][22][23] eventually leading to dynamical arrest. 16,24,25or pure solutions of one protein, the physical picture and modeling approaches up to rather concentrated solutions are well developed and can be used also to predict properties with qualitative and even reasonably quantitative accuracy.Here, pure solutions also include systems where a protein is dissolved in a complex mixture of different additives such as salts and co-solvents that is then treated as an effective solvent. 12,268][29][30] In addition, computational approaches are more involved and costly, as interactions between larger particles in a bath of smaller ones require the simulation of a large number of particles. 24][35] In this paper, we systematically and comprehensively explore a well-established binary protein system of a and c B crystallin protein from the bovine eye lens fluid as a first step toward an approximate understanding of multi-component systems.
The class of crystallin proteins present in the eye lens has been intensively investigated in the past, due to their link to the eye lens function and conditions such as cataract.7][38][39] Briefly, three distinct protein families are found in the vertebrate eye lens, called a, b, and c crystallin. 40In this work, we focus on two crystallin proteins extracted from the bovine lens: a crystallin and c B crystallin.
Bovine a crystallin, a multi-subunit protein with a molecular weight Mw % 800 kDa and an average diameter d % 15 nm, has been successfully described with a polydisperse hard sphere model. 18,24,41,42ombining x-ray photon correlation spectroscopy (XPCS), neutron spin echo spectroscopy (NSE), dynamic (DLS) and static (SLS) light scattering, zero-shear viscosity measurements, and molecular dynamics simulations of highly concentrated a crystallin solutions, the dynamical arrest transition (i.e., glass transition) has been shown to occur at volume fractions around / u 0:58, 18,19,24 driven by the socalled caging effect for hard sphere systems. 43,44 fundamentally different scenario is observed for the monomeric c B crystallin fraction (Mw u 21 kDa, d u 4 nm), a subclass of the bovine c crystallins.It has been shown to interact via weak, shortrange attractions that lead to a liquid-liquid phase separation of c B crystallin solutions upon lowering the temperature. 25,45,46Dynamic properties, such as collective diffusion coefficients measured via DLS, thus show a strong dependence on the temperature close to the liquidliquid coexistence curve (or binodal) due to critical phenomena. 25The arrest transition for c B crystallin occurs at volume fractions significantly lower than hard spheres 43,47 due to short-range attractions and patchy interactions. 19In fact, interprotein interactions induce the formation of transient networks and/or clusters that slow down the dynamics at relatively low volume fractions, which is further enhanced by the nonspherical monomeric shape. 48According to a recent study performed via microrheology measurements, the glass transition of c B crystallin occurs at a volume fraction between / % 0.27 and 0.31. 49nterestingly, the glass transition for the c B crystallin shows no temperature dependence in the investigated temperature-range.
A comprehensive characterization of mixtures of crystallin proteins is so far missing.While phase behavior and the underlying protein interactions have been addressed, 29,31,50,51 dynamical and rheological properties remained uncharacterized, despite their high relevance for a complete understanding of eye function and eye conditions such as presbyopia.
This study attempts to explore simple mixing rules based on the detailed knowledge of the individual crystallins and test whether such models are capable of reproducing the concentration and temperature dependence of key dynamic parameters such as the collective diffusion coefficient and the relative viscosity of crystallin mixtures.For binary mixtures of a and c B crystallin, we provide a characterization of their (i) phase behavior, (ii) the osmotic compressibility measured by SLS, (iii) structural properties via small angle x-ray scattering (SAXS) (iv) zero-shear viscosities utilizing both DLS-and confocal laser scanning microscopy (CLSM)-based microrheology measurements, and (v) collective diffusion coefficients from DLS and NSE experiments.We then compare the results of the ac B binary mixtures to the theoretical predictions, calculated by appropriately weighted combinations of the two properties from the individual or pure crystallin solutions.Given the large number of individual parameters discussed below, we provide a glossary of the used variables for pure and binary lens crystallin solutions in Table I.

II. MIXING RULES
Before discussing specific results, we present the basic idea of mixing rules for the prediction of relevant structural and dynamic properties of multi-component protein solutions by calculating a weighted average of the corresponding quantities of the different single-component systems.While the general idea of mixing rules is simple, the choice of the specific form invokes an underlying physical picture.Previous theoretical and experimental studies have provided evidence that for complex protein mixtures as for example present in cells with their thousands of different proteins, important features such as the location of phase boundaries or the diffusion of individual proteins are primarily determined by the effective volume fraction of the particular protein species. 52,53It has been pointed out that overall, specific strong attractions between different protein species in the cell cytosol are largely absent, and they thus interact not unlike an effective hard sphere liquid.As a result, self-diffusion of individual proteins is then for example primarily dependent on the overall protein volume fraction.Phase boundaries for liquid-liquid phase separation (LLPS) or crystallization, on the other hand, depend on the interactions between a particular protein type and thus largely on the effective volume fraction of this protein species.
For our study, we build on these findings and on the known properties of the two lens crystallin proteins investigated.We assume that a crystallin interacts with itself and with c B crystallin through excluded volume interactions only, whereas c B crystallin also experiences a short-range attractive interactions with itself, resulting in LLPS at lower temperatures.As a result, we then construct a series of mixing rules based on the assumption that the behavior of a and c B crystallin in ac B mixtures can be described as effective single-component solutions with an effective volume fraction that also takes into account the reduced accessible volume due to the presence of another hard spherelike component.These effective volume fractions / a;eff and / c B ;eff can then be expressed as and where / a and / c B are the actual volume fractions of the two proteins in solution.

A. Static scattering intensity
In order to illustrate how this picture translates into a mixing rule for static quantities, we start with the description of the static scattering intensity I(q) as a function of the magnitude of the scattering vector q as measured for example in a SAXS experiment.For a binary mixture of hard spheres with in total N particles, I(q) is given by where we have omitted some prefactors relating specific molecular properties to the absolute scattering intensity, and where F i ðqÞ and F j ðqÞ are the scattering amplitudes of an individual particle for the two particle species, and S ij ðqÞ the so-called partial structure factors, respectively. 54ithin the underlying assumption of describing the contributions of each species to a solution property as those from an effective one-component system defined by an effective volume fraction, we thus rewrite the full scattered intensity I ac B ðqÞ of the binary system as a sum of two effective self-terms as follows: The scattering of the mixed solution is thus described only by an intensity-weighted average of the single-component solutions.Here, c a and c c B are the concentration of the pure crystallins in mg/mL and the total protein concentration c ac B ¼ c a þ c c B .I FF a ðqÞ and I FF c B ðqÞ are the concentration-normalized scattered intensities of a low-concentration, non-interacting solution (i.e., form factor).
S / a;eff ðqÞ and S / c B ;eff ðqÞ are the structure factors for the pure components taken at the effective volume fraction, as calculated using Eqs.( 1) and (2).With this approach, we thus assume that the effect of mutual interactions between a and c B crystallins can be represented by an increased excluded volume effect in the interactions between the individual crystallins (i.e., aa and c B -c B ).

B. Diffusion coefficients
While mixing rules for static scattering can thus be assumed as simple additive quantities of individual components, dynamic quantities show an additional complexity: in addition to the direct excluded volume interaction, hydrodynamic interactions that depend mainly on the total volume fraction also affect the dynamics.To take the hydrodynamic slowing down into account, we base the mixing rule on the well-known relation for collective colloidal diffusion, 55 DðqÞ ¼ D 0 HðqÞ SðqÞ ; where the free diffusion coefficient D 0 is modulated by both the static structure factor S(q) and the hydrodynamic function H(q) to obtain the lengthscale-dependent diffusion function D(q).

Gradient diffusion from dynamic light scattering
In DLS experiments, the characteristic length scale 2p=q probed is much larger than the protein diameter d.DLS therefore probes the so-called gradient diffusion coefficient given by Dðqd !0Þ.Thus, we use the established relation for H(0) in the limit of qd !0 to correct for hydrodynamic effects where U s =U 0 represents the concentration-dependent short-time sedimentation velocity. 56For D 0 , we calculate the intensity average from the individual diffusion coefficients D 0;a and D 0;c B at their effective volume fractions.We thus obtain the effective average free diffusion coefficient D 0;ac B from The corresponding effective average S eff ðqÞ in the limit of qd !0 is calculated using

Local collective diffusion from neutron spin echo spectroscopy
For cases where q % 2p=d such as for neutron spin echo measurements, the above approach is not possible anymore, as H(q) is dependent on the structure and at the same time on the total volume fraction.We thus use the weighted average of the measured q-dependent collective diffusion coefficients D a ðqÞ and D b ðqÞ, where I a ðqÞ ¼ c a I FF a ðqÞS / a;eff ðqÞ and . D a ðqÞ and D c B ðqÞ denote the measured diffusion coefficient at a volume fraction corresponding to the effective volume fraction / a;eff and / c B ;eff , respectively.We remark that the measured diffusion coefficient thus do not contain the full hydrodynamic slowing down, so that a formal analogy to Eq. ( 5) is only preserved to a certain degree.While some shortcomings due to the neglected hydrodynamic slowing down at elevated volume fractions are to be expected, a thorough test how far these mixing rules can work is nevertheless of general interest.

C. Viscosity
For viscosity measurements, we apply a simple Arrhenius relation for the relative viscosity g r of mixtures 57,58 to theoretically predict g r ð/ ac B Þ of the binary mixture from the relative viscosities of the individual protein solutions at their effective volume fraction g r ð/ a;eff Þ and g r ð/ c b ;eff Þ, Here, we assume that a particular protein species i diffuses in a viscous medium formed by the other protein species j and the solvent.The viscosity of this effective solvent is given by that of the solution of the protein j at its effective volume fraction / j;eff , i.e., by gð/ j;eff Þ.The contribution of protein i to the total viscosity is then described by the corresponding relative viscosity of a solution of protein i at a volume fraction / i;eff , i.e., by gð/ i;j Þ=gð/ j;eff Þ ¼ g r ð/ i;eff Þ.The individual values of g r ð/ a;eff Þ and g r ð/ c B ;eff Þ were obtained from the measured concentration dependence of g r ð/ a Þ and g r ð/ c B Þ using a parameterization as described in supplementary material.We stress that the simplicity of the mixing rules presented should not be seen as a serious drawback, but rather as a first test of how far simple concepts could be used to describe and predict the behavior of complex mixtures of different proteins under highly crowded conditions as often found in the cytosol of living cells or in blood serum.Section III presents different experimental quantities over a large range of concentrations, which are then compared with the predictions from these mixing rules.Importantly, these experimental quantities also encompass macro-, meso-as well as microscopic properties, and we will discuss the important consequences of these different length scales in the applicability of the mixing rules.

III. RESULTS AND DISCUSSION A. Phase diagram
We first test our approach with a macroscopic property and investigate the location of the coexistence curve for liquid-liquid phase separation in a À c B protein mixtures.LLPS in the eye lens has been investigated intensively for its link to cataract formation and found to be directly related to c B À c B protein-protein interactions. 6,8,46The results from our characterization of the location of the binodal for LLPS in these mixtures as determined from cloud point measurements for three different a : c B volume ratios (3:1, 1:1, and 1:3) are shown in Fig. 1(a).
When compared to the binodal in pure solutions of c B crystallin from an earlier study 25 (gray triangles), we see that an increasing amount of a crystallin has a significant effect, shifting the binodal to higher overall protein concentrations.It is interesting to look at this in view of the evolutionary pressure on the eye lens to maintain a cytosol composition that allows it to achieve high transparency at the very high protein concentrations required to generate the necessary index of refraction profile across the lens. 6,59As we will discuss below, the different crystallins allow indeed for a densely packed cell interior that optimizes index of refraction, transparency as well as viscosity.
If our reasoning in Sec.II is correct, we would expect to find the location of the binodal primarily determined by the relative volume fraction of c B corrected by the reduced available volume caused by the a crystallins present, i.e., by / c B ;eff calculated using Eq. ( 2).This is demonstrated in Fig. 1(b), and we see that indeed all data points collapse around the binodal for the pure c B system when using / c B ;eff instead of / acB .
The use of such a simple mixing rule based on the picture of a protein embedded in an effective solvent, where the additional other protein species present only result in a reduced available sample volume and thus an effectively increased volume fraction, is of course thermodynamically incorrect.We remark that a complete modeling of the phase diagram of ac B binary mixtures has in fact been performed previously. 29,51This work was also extended to a complete description of the resulting light scattering intensity based on a model of sticky spheres, allowing for a detailed investigation of the effects of proteinprotein interactions (a À a; a À c B and c B À c B ) on the intensity I(q) at q ¼ 0. 60 However, the resulting multidimensional phase diagrams become far too complex to be useful in a simple framework aiming at a qualitative understanding of the stability, the effect of temperature and composition, and the existence and location of an arrest line in crystallin mixtures with compositions similar to those in eye lens cells.Moreover, these calculations are difficult to extend to dynamic quantities such as the relative viscosity or the collective diffusion coefficient.We therefore aim to develop and test simple models to describe phenomenological observations in a semi-quantitative way, amenable for numerical calculations of thermodynamic, structural and dynamic properties that could not be achieved otherwise and would require extensive computer simulations instead.

B. Viscosity
Next, we extend our approach to a macroscopic dynamic quantity and investigate the relative zero-shear viscosity g r of ac B binary mixtures.We employ measurements using both DLS-based and MPT-based microrheology, as described in the Materials and Methods section.To ensure that experiments are made on mixtures in the onephase region of the phase diagram, all measurements are performed at temperatures well above T cloud .
Figures 2(a The volume fraction dependence of g r of the mixture is well reproduced by the mixing rule in Eq. (10) [Figs.2(a) and 2(b)].(see supplementary material for a detailed description).The combination of the g r of a and c B crystallins taken at their effective volume fractions / c B ;eff and / a;eff indeed appears to provide a good prediction of the binary mixture.This is quite remarkable given the very different mechanisms behind dynamical arrest for the two different lens proteins.While a crystallin exhibits a hard sphere-like behavior with an arrest transition at around / g;a % 0:58, 24 c B crystallin has been found to show a behavior described by a model of particles with a weak short range attraction combined with additional attractive patches.This results in a dynamical arrest transition at significantly lower volume fractions of / g;c B % 0:3 À 0:35, presumably caused by the presence of transient clusters. 19,25This good agreement suggests that on macroscopic length scales probed by the zero shear viscosity, and given strongly differing viscosities, viscoelastic relaxations arising from protein-protein interactions can be decoupled for the two proteins.At the effective volume fractions / a;eff ( / g;a measured for our binary mixtures, a crystallin solutions have rather low viscosities and relax quickly, allowing for a treatment of the a crystallin as an effective solvent for c B crystallin, which shows significantly higher viscosity and slower relaxation times at these effective volume fractions due to the formation of transient clusters. Due to the attractive interactions between c B crystallins, their solutions undergo LLPS at lower temperatures as demonstrated in Fig. 1.This also has consequences for scattering experiments, as there will be a significant contribution from critical scattering, resulting in a strongly enhanced intensity in SLS and a slowing down of the measured correlation function obtained in DLS as the temperature approaches T cloud . 25While these additional contributions from critical opalescence and critical slowing down make microrheology measurements using DLS-based tracer microrheology much more difficult and induce additional errors (see supplementary material), this can be overcome by using MPT-based microrheology instead.However, while temperature has a strong effect on the phase behavior and the results from scattering experiments, the relative viscosity shows no measurable temperature dependence.This indicates that there are no or only weak critical contributions to the g r of binary ac B mixtures.
An interesting feature is observed when plotting and comparing relative viscosity data for pure crystallins and the binary mixture as a function of mass concentration c in mg/mL.As shown in Figs.2(c) and 2(d), g r from the different systems collapse onto a master curve that extends to rather high concentrations.This coincidence suggests that for these specific proteins, g r is determined by the actual weight concentration for low to intermediate protein concentrations, independent of the ratio between the pure components.It is important to point out that this g r master curve holds up to concentrations c of around 300 mg/mL.For higher c, g r diverges at the arrest value c g of the particular protein species or mixture, which is different for the two proteins.In particular, for a crystallin the arrest transition occurs at % 341 mg /mL [red dashed lines in Figs.2(c) and 2(d)], while dynamical arrest for c B crystallin occurs at % 395 mg/mL [green dashed lines in Figs.2(c) and 2(d)].However, it is important to point out that we do not expect to find such a behavior for arbitrary proteins and consider this an effect of the very different voluminosities of a and c B crystallin that compensate for the different underlying arrest mechanisms.Instead, we expect that a mixing rule such as given in Eq. (10), where effective volume fractions are considered, should have a much wider applicability.

C. Collective gradient diffusion
Next, we move on to mesoscopic length scales of many protein diameters and investigate collective gradient diffusion of ac B binary mixtures using dynamic light scattering.The protein solution is no longer considered to behave as an effective homogeneous liquid, but as a concentration field, in which concentration fluctuations or inhomogeneous density distributions relax with a rate expressed through the gradient diffusion coefficient D ¼ Dð/; q !0Þ.
Figure 3 presents the measured diffusion coefficients D ac B ð/ ac B Þ for the 1:1 mixture at two temperatures T ¼ 25 and 35 C for a series of volume fractions / ac B .For both temperatures, we observe a decrease in Dð/ ac B Þ with increasing volume fraction.Also shown is the concentration dependence of the pure components D a ð/ a Þ and D c B ð/ c B Þ, respectively. 24,25hile a crystallin solutions exhibit a concentration dependence of Dð/Þ that is typical for hard spheres with an almost linear increase in Dð/Þ with increasing volume fraction, 24 the /-dependence of Dð/Þ is much more complex for c B crystallin.For this protein, we find a non-monotonic /-dependence due to the combination of contributions from a strongly T-dependent critical slowing down in the vicinity of the critical point and a regular non-critical contribution leading to dynamical arrest at / % 0:35. 25For the 1:1 ac B mixture, we observe a behavior that is roughly intermediate between the two individual protein solutions.
We now attempt to use Eq. ( 5) together with the simple mixing rules given in Eqs. ( 6)-( 8) in order to calculate the theoretical predictions for the mixed system.To determine D 0 for the simple mixing rule in Eq. ( 5), we use D 0;a ¼ 2.2 Â 10 À11 and D 0;c B ¼ 8.74 Â 10 À11 m 2 s À1 at T ¼ 20 C, which we then adapt for the correct solvent viscosities at the temperatures measured (see supplementary material for details). 24,25For S(0) in the case of a crystallin, we employ the Carnahan-Starling equation for monodisperse hard spheres S CS ð0Þ for each temperature are taken from Ref. 25 and parameterized with a simple polynomial function.
The agreement between the measured concentration dependence of the gradient diffusion coefficient Dð/ ac B Þ and the predictions from Eq. ( 5) together with the simple mixing rules given in Eqs. ( 6)-( 8) is quite remarkable, especially when considering that the individual proteins exhibit such a vastly different behavior (Fig. 3).We believe that the hybrid approach used here is essential for this, where we calculate the thermodynamic driving force for the decay of concentration fluctuations given by the osmotic compressibility S(0) based on the particular a À a and c B À c B interactions using the effective volume fractions of both protein components, while we estimate the contributions from hydrodynamic interactions from the overall volume fraction.This allows us to also incorporate contributions from critical fluctuations for proteins that undergo LLPS, while the frictional resistance that the different proteins experience should primarily depend on the total volume fraction of proteins in solution.

D. Solution structure
After having successfully reproduced macro-and mesoscopic thermodynamic and dynamic properties, we next investigate whether the proposed mixing rules can also describe structural and dynamic solution properties on more local or molecular length scales.We first investigate the solution structure on length scales from a fraction of the protein diameter up to a few protein molecules as obtained from SAXS measurements.
We first verify that Eq. ( 4) correctly reproduces the measured SAXS intensity I(q) for protein mixtures in the absence of contributions from protein-protein interactions, i.e., in the limit of infinite dilution where S / a;eff ðqÞ ¼ S / c B ;eff ðqÞ ¼ 1.Under these conditions, we measure an effective form factor I FF ac B ðqÞ that corresponds to an intensity-weighted average of the two individual form factors of a and c B crystallin given by Figure 4 shows the concentration-normalized form factor I FF ac B ðqÞ from SAXS for dilute solutions of ac B binary mixtures at different a:c B mixing ratios, i.e., 1:3 (red circles), 1:1 (green squares), and 3:1 (blue triangles).At this dilute volume fraction of / ac B ¼ 0:003, effects from protein-protein interactions can be neglected.This is indeed demonstrated in Fig. 4, where we see that the scattering from the mixed solutions is quantitatively reproduced by a linear combination (colormatched lines) of the individual form factors for the pure a (purple line) and c B (orange line) crystallins also shown in Fig. 4.This also clearly indicates the absence of any self-association under these conditions.
We next characterize the concentration and temperature dependence of the solution structure for 1:1 mixtures using SAXS.Figures 5(a The theoretical predictions are based on the scattering intensities measured for the pure a and c B crystallin solutions using the mixing rule in Eq. ( 4) and their effective volume fractions / a;eff and / c B ;eff [Eqs.( 1) and ( 2)], respectively.This requires theoretical expressions for the structure factors of the pure components S / a;eff ðqÞ and S / c B ;eff ðqÞ.
For the structure factor of a crystallin, we use the Percus-Yevick approximation for polydisperse hard spheres, where the particle size distribution is described with a Schulz distribution with an average particle diameter of 15 nm and a polydispersity of 20%, which has been shown to describe S a ðqÞ with very good accuracy on the entire accessible volume fraction range up to the arrest transition at / g;a % 0:58. 24For pure a crystallin solutions, we moreover expect no measurable temperature dependence of S / a ðqÞ for the temperatures investigated.
For c B crystallin solutions, the situation is more complex, as the additional attraction between proteins results in a very strong and non-monotonic temperature and concentration dependence of S / c B ðqÞ FIG. 3. Gradient diffusion coefficients of 1:1 ac B mixtures [D ac B ð/; q !0Þ] and pure components [D a ð/; q !0Þ; D c B ð/; q !0Þ)] in the q !0 limit at two temperatures and as a function of volume fraction.Lines: theoretical diffusion coefficients for binary mixtures calculated as described in the text (solid lines), and interpolated values from DLS experiments with the pure crystallins (dashed lines).Solid symbols: collective diffusion coefficients obtained via DLS measurements at T ¼ 25 C (red triangles) and 35 C (blue circles).

FIG. 4.
Scattering intensity from SAXS for ac B binary mixtures at different a (or c B ) ratios and a total / ¼ 0.003.Red, green, and blue open symbols identify the 1:3, 1:1, and 3:1 ratios, respectively.Form factors for pure crystallins are indicated as purple lines for a crystallin and orange lines for c B crystallin.Colored solid lines are the resultant fits from Eq. (11).
in the vicinity of the critical point.As previously reported in the literature, 25 the structure factor of c B crystallin is well described by a combination of two contributions from critical and non-critical scattering given by at q Շ 2p d , with the protein diameter d c B % 4 nm.S crit describes the contributions from critical concentration fluctuations and thus strongly depends on temperature and concentration, primarily given by the distance of the particular state point from the critical point or spinodal line.The second term S non crit represents the qdependent background, related to the athermal excluded volume interactions.
We calculate S c B ðq; T; /Þ with the following expression: The first term is a simple Ornstein-Zernike relation to describe the contributions from critical scattering in Eq. ( 13), where S crit ð0Þ is given by the osmotic compressibility and n is the static correlation length. 25or the non-critical contribution, we consider exclude volume interactions between the ellipsoid-shaped c B crystallins only, which we describe using a so-called random-phase approximation, where P c B ðqÞ is the form factor of c B , normalized to Iðq; 0Þ ¼ 1.The actual values used are taken from a detailed investigation of the structural and dynamic properties of c B -crystallin described in Ref. 25.More details about the procedure and an example of the different contributions to the scattering intensity from a c B crystallin solution can be found in supplementary material.
The comparison between the experimental data and the calculated curves shown in Fig. 5 reveals good agreement for higher q-values at q տ 0:5 nm À1 , and also the low-q limit Ið0Þ=c ac B obtained from SLS is quite well reproduced at both temperatures.However, at intermediate q-values around 0.07 nm À1 Շ q Շ 0:5 nm À1 we observe systematic differences for both temperatures.There appears to be a more prominent local structure on length scales of the a-crystallin diameter than predicted by our mixing rule model.In order to better understand the origins of this discrepancy, it is quite instructive to look at the individual contributions to the scattering intensity arising from the two proteins also shown in Fig. 5.The high-q part of I ac B ðqÞ=c ac B obtained from SAXS is dominated by the scattering from c B -crystallin, while the scattering from a-crystallin dominates for 0.1 nm À1 Շ q Շ 0:4 nm À1 , i.e., close to the nearest neighbor peak of a-crystallin solutions at comparable protein volume fractions.Finally, at low-q values it depends on temperature whether the scattering for a-crystallin or c B -crystallin dominates due to the strong dependence of S crit ðq; T; /Þ in Eq. (12).
In order to better understand the limits of our mixing rule approach, we also calculate the measured effective structure factor using where I ac B ðqÞ corresponds to the measured scattering intensity of the mixture, c ac B is the total weight concentrations and I FF ac B ðqÞ is the measured concentration-normalized form factor of the mixture.The measured values of S eff ðqÞ are shown in Fig. 6, together with the calculated values based on the mixing rule approach.Here we use to calculate S eff ðqÞ, where I FF a ðqÞ and I FF c B ðqÞ are the measured form factors of the two proteins, while S / a;eff ðqÞ and S / c B ;eff ðqÞ are calculated as described above.
Figure 6 indeed shows that the local structure is well reproduced for q տ 0:5 nm À1 , where the scattered intensity from a crystallins has strongly decreased due to the form factor of the large protein, and where the data thus primarily describe the local structural correlations between c B crystallins.Here, S / c B ;eff ðqÞ is also almost independent of temperature, as these local correlations are dominated by excluded volume effects.At low q-values, S eff ðqÞ exhibits a clear temperature dependence, reflecting the enhanced contribution form critical scattering from c B with decreasing temperature.Moreover, on length scales comparable to the a crystallin diameter and larger, the structural correlations are much more pronounced than predicted.Furthermore, the additional low-q value from SLS indicates the existence of additional weak attractions that result in a significant upturn of S / c B ;eff ðqÞ at low q.
The c B contributions from critical fluctuations and the resulting static correlation length at these temperatures are too small to cause this considerable correlation hole and subsequent strong upturn at low q.However, previous results from SANS experiments and computer simulations of a À c B mixtures have provided compelling evidence for the presence of mutual weak attractions between these two proteins, resulting effectively in a higher stability against phase separation when compared to a situation where the two protein species interact via excluded volume interactions only. 28,50Computer simulations that allow for the calculation of the different partial structure factors [S aa ðqÞ; S ac B ðqÞ, and S c B c B ðqÞ in Eq. ( 3)] indeed show that S aa ðqÞ exhibits a behavior that is strongly reminiscent of weakly attractive hard spheres, where the c B proteins act as temporary bonds between two a crystallins. 28t is quite intriguing that the mixing rule approach reproduces the osmotic compressibility of the mixture quite well, despite the fact that local correlations are not well described and strongly underestimated on length scales of the a crystallin diameter.We have characterized S(0) for different mixing ratios, temperatures and total protein volume fraction using SLS, and the results are summarized in Fig. 7.
Here, we see that Eq. ( 8) indeed reproduces the measured compressibility quite well for the higher total volume fractions / ¼ 0:16 and / ¼ 0:22, but fails to do so at the lowest volume fraction / ¼ 0:11.Moreover, the agreement is also less good for the mixing ratio a:c B ¼ 1 : 3 and the lower temperature T ¼ 25 C, where the measured values are much higher than predicted by the mixing rule approach.At high concentrations and higher temperatures where critical fluctuations are not important, the compressibility is primarily determined by strong excluded volume effects between the different proteins, which is reasonably described by our approach.On the other hand, addition of a smaller amount of a crystallin to c B crystallins enhances the criticality of c B solutions in a non-monotonic way, with a maximum at a mixing ratio of 1 : 3.Under these conditions, the mixture phase separates at a higher temperature, and critical scattering is thus significantly enhanced already at T ¼ 25 C. 28 Moreover, the effects of an additional weak attraction on S(0) become much more visible at lower volume fractions, where S(0) shows also a nonmonotonic behavior for attractive hard sphere systems with a maximum that depends on the nature and strength of the attraction. 62ith c B crystallins acting as temporary bonds between a crystallins, this can lead to enhanced scattering even far away from phase separation, and thus qualitatively explain the failure of the mixing rule calculations at lower total protein concentrations.
The analysis of the SAXS data has clearly demonstrated the limits of the simple mixing rules in describing the properties of protein mixtures at high concentrations, i.e., under crowded conditions.While they were able to reproduce the experimental data for the macroscopic and mesoscopic thermodynamic and dynamic properties almost quantitatively, they fail to describe the microscopic structural features for higher concentrations in the 1:1 mixture.This clearly indicates that we would need to explicitly consider the different partial structure factors, which would require us to properly treat also mutual interactions between a and c B crystallins.This is also in line with earlier simulation work, which had pointed out that a weak mutual attraction between a and c B crystallins resulted in an enhanced stability of the mixture, suppressing LLPS and increasing transparency of the solution.However, it is also interesting to note that these mixing rules are able to reproduce the osmotic compressibility quite well at high concentrations, i.e., under conditions existing in living cells, as long as the solution conditions are not close to a boundary for LLPS for some of its proteins.8) where we use the full q-dependence instead of only the q ¼ 0 limit.Open (SLS) and solid (SAXS) symbols are for the experimental results, solid lines for the predictions given by Eq. ( 8).

E. Collective diffusion on the local scale
As a final test, we explored collective diffusion on the microscopic length scale of individual molecules using neutron spin echo spectroscopy, i.e., the relaxation of density correlations such as cage formation in concentrated protein solutions.
We collected intermediate scattering functions I(q, t) for mixed solutions of a and c B crystallin at volume fractions of 0.11 and 0.22.We used several q values covering both the correlation hole below 0.1 nm À1 and the nearest-neighbor scale of c B crystallin beyond 2 nm À1 .From the rate constant CðqÞ of the initial decay of the measured correlation functions I(q, t), we obtained q-dependent diffusion coefficients DðqÞ ¼ C=q 2 (cf.supplementary material for details).
Figure 8 displays the results as a function of the fraction of c B crystallin (symbols).We observe clear non-monotonous trends of the diffusion coefficients with increasing c fraction.In all cases, the 3:1 a : c B mixtures relax faster than pure a crystallin.From this maximum value, D(q) decays to the value for pure c B crystallin.While the signatures at different q values differ mainly by an offset, the observed trend is significantly varied between the two explored volume fractions.
It has previously been reported that for proteins such as c B crystallin experiencing patchy attractions, collective diffusion at length scales of the nearest neighbor distance can be dramatically altered when compared to the predictions for simple colloid models based on an isotropic potential due to the formation of transient clusters. 19hile collective diffusion on mesoscopic length scales measured by DLS could be well described by a simple mixing rule given by Eq. ( 5), we can no longer use such an approach.Instead we calculate a weighted average of the diffusion coefficients of the pure systems D a ðqÞ and D c B ðqÞ at their respective effective volume fractions / a;eff and / c B ;eff as given by Eq. ( 9).The thus predicted values of D(q) indeed show a qualitatively similar behavior as observed in the NSE experiments.
This particular dependence on composition is a result of the combination of the very different scattering intensities and diffusion coefficients for the two proteins.At both q-values, the scattering intensity from c B crystallin is much higher than that of a crystallin, resulting in ratios of I a ðqÞ=I c B ðqÞ as summarized in Table II.At the same time, the collective diffusion coefficient D c B ðqÞ exhibits a much stronger concentration dependence than D a ðqÞ (see supplementary material), resulting in a non-monotonic dependence of the intensity weighted D (q) as a function of the mixing ratio.
However, while the qualitative trend is correctly reproduced, the actual values calculated are systematically overestimating the measured D(q).This is particularly pronounced for the 3:1 mixture.Moreover, the measured q-dependence of D(q) is much more pronounced than predicted by Eq. ( 9), with significantly lower values of D(q) at q ¼ 2.2 nm À1 than at q ¼ 0.08 nm À1 .
The diffusion coefficients measured by NSE represent short-time collective relaxations on the nearest neighbor distance and are governed by the local neighbor structure as well as by hydrodynamic interactions.The effective static structure factor of the mixture S eff ðqÞ is quite well reproduced by our mixing rule approach at these q-values (see Fig. 6), indicating that the local structure on these rather short length scales of order the size of c B crystallin are not too different from those underlying the mixing rules used, i.e., assuming that they are determined primarily by the c B -c B and a-a interactions at their corresponding effective volume fractions.On the other hand, in our FIG.8. Diffusion coefficients D(q) from neutron spin echo spectroscopy in mixtures of a and c crystallin as a function of total volume fraction and composition for two qvalues (black symbols: q ¼ 0.8 nm À1 ; red symbols: q ¼ 2.    While this explains some of the overestimation of D(q) by our model, it can of course not explain the significant q-dependence in the measured data.
For a puristic theoretical understanding, one would need all pairwise partial structure factors, which are not available for the given experimental systems.While a complete understanding of these signatures would thus require comprehensive simulations or complicated parametrizations of complex colloid theory that also takes into account possible patchy interactions between different proteins, this deviation from the simplistic mixing rule implies that mutual interactions between the two proteins actually strongly affect the dynamics on local length scales, as intuitively expected.

F. Conclusions
Driven by the need to understand crowding effects in protein mixtures as present for example in the living cell, we have embarked on a systematic test of the application of simple mixing rules in order to reproduce a number of key properties.We focus, in particular, on macroscopic quantities such as the location of phase boundaries for liquid-liquid phase separation, the osmotic compressibility and the zero shear viscosity, properties on mesoscopic length scales such as the gradient diffusion coefficient, and finally structural and dynamic quantities on local length scales such as the static structure factor and the (local) short-time collective diffusion coefficient.There is no applicable theory capable of incorporating different sizes, interaction potentials, and anisotropy in shape and interactions common for crowded protein mixtures.While one could of course use computer simulations, the resulting computational costs still become unaffordable, in particular, when aiming to reproduce dynamic properties on large length and time scales.Therefore, we have defined different mixing rules for these properties that are all based on the same underlying principle: There is no specific interaction between different species other than excluded volume effects, and we estimate the measured quantities from a weighted sum of the properties of the individual proteins in crowded solutions at an effective volume fraction that takes into account the accessible volume due to the presence of other protein species.
Despite their simplicity, we obtain rather good agreement with experimental data on macroscopic and mesoscopic length scales, i.e., phase behavior, osmotic compressibility, zero shear viscosity, and gradient diffusion.It is only at low volume fractions and low temperatures close to the binodal for LLPS that deviations become quite significant.The agreement is less convincing at intermediate length scales of order the protein size of the larger protein a-crystallin, where the lack of a mutual attraction between c B and a-crystallins results in significant disagreement between calculated and measured values for both S(q) and D(q).At even shorter length scales of the size of the smaller protein c B crystallin, the finding is mixed, with very good agreement for S (q), but significant deviations for D(q).This clearly shows the sensitivity of the local short-time collective diffusion coefficient on patchy attractions that lead to the formation of transient clusters with life times comparable or larger to the decay time probed with NSE. 19This is quite in contrast to the mesoscopic gradient diffusion coefficient measured with DLS, where characteristic relaxation times are long compared to the cluster lifetime, and where thus the mixing rule relying on weighted individual properties become quite powerful.
While our approach will require further validation with more complex mixtures, the success of such simple rules for the present system of a binary protein solution, where the individual proteins exhibit vastly different solution properties, appears promising for future applications in real life systems such as cells or other highly crowded protein mixtures.

IV. METHODS A. Sample preparation
Eye lens proteins were extracted from calf eye lenses and purified according to previous work. 25,51In brief, a crystallin was extracted from the cortical extract solution via size-exclusion chromatography (SEC, Superdex 200 column), using a 52.4 mM phosphate buffer at pH 7.1 in H 2 O.The solvent of the collected fractions was exchanged with a D 2 O phosphate buffer with the addition of 20 mM of 1,4-dithiothreitol (DTT, Sigma-Aldrich, SE), 1 mM of ethylenediaminetetraacetic acid (EDTA, Sigma-Aldrich, SE) and 0.02wt.% sodium azide (Sigma-Aldrich, SE), pH 7.1.Finally, the new solution was concentrated via Amicon Ultra centrifugal filters of 10 kDa (Sigma-Aldrich, SE).The concentration determination was achieved via UV absorption spectroscopy measurements at k ¼ 280 nm and the specific absorption coefficient E 0:1%;280 nm a;1 cm ¼ 0.845 mL Á mg À1 Á cm. 51We converted the weight concentrations to volume fractions, according to the relation / ¼ c Á , where c is the concentration in mg/mL and is the voluminosity of the protein.For a crystallin, ¼ 1.7 mL/g. 24or c B crystallin, the filtered nuclear extract was purified through SEC (same as above), using a 275 mM sodium acetate buffer (Sigma-Aldrich, SE) in H 2 O, pH 4.5 as a mobile phase.The collected fractions were then further purified via ion-exchange chromatography (IEX, SP Sepharose Fast Flow, GE Healthcare, USA), using a sodium acetate buffer in H 2 O, pH 4.8 with a 0-325 mM sodium chloride (Sigma-Aldrich, SE) gradient.This allowed the collection of the pure c B protein according to its isoelectric point (pI). 63The following steps (solvent exchange, sample concentration, protein concentration determination using a specific absorption coefficient E 0:1%;280 nm c B ;1 cm ¼ 2.18 mL mg À1 cm, and conversion to volume fractions) were achieved as described above for a-crystallin.In this case, the voluminosity of c B corresponds to ¼ 0.71 mL/g. 46e investigated a : c B mixing ratios of 3:1, 1:1, and 1:3 per volume for the phase diagram, SLS, SAXS, and NSE, while we focused on the 1:1 mixing ratio for DLS and microrheology.

B. Static and dynamic light scattering
Static light scattering (SLS) experiments and transmission measurements for the determination of S(0) were performed on a homebuilt multi-angle light scattering (MA-DLS) instrument.This spectrometer not only allows for suppression of multiple scattering using the 3D cross correlation technology, but has the additional advantage of simultaneous measurements at 4 scattering angles h separated by 30 .The q-range covered by LS is 0.003-0.03nm À1 . 64dditional dynamic (DLS) and static (SLS) light scattering measurements were performed with a goniometer light scattering setup (3D LS Spectrometer, LS Instruments, AG), implemented with a 3D cross correlation unit to suppress contributions from multiple scattering. 65,66nless stated, all measurements were performed at a h range from 30 to 130 (steps of 5 ), corresponding to a scattering vector range q ¼ ð4pn=kÞ sinðh=2Þ % 0.006-0.023nm À1 .
For DLS, intensity auto-correlation functions g 2 ðq; tÞ À 1 were converted to the intermediate scattering functions , where t is the lag-time and b is the spatial coherence factor, determined by the intercept of g 2 ðq; t Þ À 1.The intermediate scattering functions from q-dependent DLS were then described by fitting the initial part of the correlation function g 1 ðq; tÞ until it has decayed to 80% of its initial amplitude, using a single exponential function g 1 ðq; tÞ ¼ expðÀC tÞ; where C is the relaxation rate and corresponds to the inverse of the relaxation time s.The collective diffusion coefficients D c were then determined according to the relation D c ¼ C=q 2 , which corresponds to the slope of the linear fit of C vs q 2 .

C. Microrheology
Microrheology was performed via DLS 61 and multiple particle tracking (MPT) using confocal laser scanning microscopy (CLSM). 67racer particles were prepared according to, 68 using polystyrene particles stabilized with covalently bound 20 kDa poly(ethylene) glycol chains.For DLS-based microrheology, we use a particle diameter d ¼ 300 nm.On the other hand, for CLMS-based MPT we rely on fluorescent-labeled polystyrene particles, d ¼ 1 lm.The colloidal stability of the tracer particles was tested a priori according to previous work. 61,68or DLS measurements, a volume of 1 lL of the tracer particle solution (concentration % 2 wt.%, number density % 10 12 particles Á cm À3 ) was added to 100 ll of protein solution.All DLS experiments were carried out at a single h ¼ 90 at T ¼ 25 and 35 C. The addition of tracer particles to the protein solutions results in a double-step relaxation process in the g 2 ðq; tÞ À 1 function, originating from the scattering contributions from the proteins and the tracer particles.
In our measurements, the slow relaxation mode resulting from the tracer particles (i.e., at large t) has an amplitude that is > 95% of the total g 2 ðq; tÞ À B function, making the additional fast relaxation contribution from the proteins negligible. 61he auto-correlation functions were described with a double exponential function, using an iterative non-linear fitting procedure based on Ref. 69: where B is the baseline, b is the spatial coherence factor, A is the amplitude of the first relaxation mode C 1 , i.e., the protein contribution, and C 2 is the decay from the tracer particle contribution.The diffusivity of the tracer particle in the protein sample was then calculated as D Sample ¼ Cq 2 .Using the Stokes-Einstein relation, we then calculated the relative viscosity (g r ) value, where g Sample corresponds to the zero shear viscosity of the sample and D Ref to the diffusion coefficient of the tracer particle dispersed in pure solvent with the solvent viscosity g Ref .
For MPT, a volume of 1 lL of the stock solution of the fluorescent particles (d ¼ 1 lm (ThermoFisher, F8820), concentration % 2 wt.%) was added to 100 ll of protein solution.5 ll protein solution with particles were transferred in a hermetically closed microscope glass slide.After their preparation, the samples were thermally equilibrated (T ¼ 25 and 35 C) in a Leica SP5 confocal laser scanning microscope equipped with a temperature chamber.Videos of the samples were performed by acquiring 6000 frames, with a rate of 0.13 frame s À1 in a sample depth of 10 lm (i.e., distance from the upper glass slide).Videos were then processed via IDL software (L3 Harris Geospatial) and dedicated scripts. 67,702][73] Finally, the viscosity was obtained via the Stokes-Einstein relation and normalized by the viscosity of D 2 O at the given temperature. 74In order to reach high / and determine the arrest transition of the sample, we rely on the procedure described in a previous work. 67Here, the sample is left to evaporate for a controlled amount of time (from % 1 to 15 min) before closing it hermetically.The difference between the weight of the sample at the moment of its transfer to the glass slide and the weight after the solvent evaporation, allow the determination of the final volume fraction.
The motivation for MPT-based microrheology measurements in addition to the DLS-based microrheology is twofold.First, with MPT microrheology, we bypass the contribution of critical fluctuations around the critical volume fraction that create interference with the signal of the tracer particles in DLS.Second, MPT allows to determine the point where the glass transition occurs. 67

D. Small angle x-ray scattering
Small angle x-ray scattering (SAXS) measurements were performed with a pinhole camera system (Ganesha 300 XL, SAXSLAB) equipped with a high brillance microfocus sealed tube and thermostated capillary stage.The accessible q-range for these measurements was from 0.06 nm À1 q 6 nm À1 .All measurements were corrected for background radiation, transmission, buffer, capillary, and renormalized for the protein concentration.

E. Neutron spin echo spectroscopy
Neutron spin echo (NSE) experiments were performed at the J-NSE instrument 75 at the Heinz Maier-Leibnitz Zentrum (Munich, DE).The probed q-range spanned from 0.1 to 2.3 nm À1 and a Fourier time ( ti ) between 0.03 and 140 ns.Intermediate scattering functions for NSE were fitted with a single exponential decay for the initial decay (t < 50 ns): Iðq; ti Þ ¼ expðÀC ti Þ and the diffusion coefficients were calculated via D NSE ¼ Cq 2 .For pure a and c B crystallin solutions, the diffusion coefficients were taken from Ref. 19 and interpolated to match the right volume fraction range (for details see supplementary material).

SUPPLEMENTARY MATERIAL
See the supplementary material for the following additional information: Examples for microrheology results from 3D-DLS and MPTbased measurements using confocal laser scanning microscopy (Fig. S1); a table with the values of the parameters used for the parameterization of the relative viscosities of the pure crystallin solutions (Table S1); concentration and temperature dependence of the gradient diffusion coefficient and the osmotic compressibility of c B crystallin determined by DLS and SLS (Fig. S2); visualization of compilation of c B crystallin structure factors as a function of volume fraction and temperature (Fig. S3); interpolation of diffusion coefficients from NSE as a function of volume fraction for the pure crystallin solutions (Fig. S4); q-dependent diffusion coefficients (Fig. S5), and examples of intermediate scattering functions (Fig. S6) from NSE for different binary crystallin mixtures.
) and 2(b) report g r for the pure a crystallin (red triangles) and c B crystallin (green squares) solutions as well as for the 1:1 binary mixture as a function of volume fraction / at T ¼ 25 and 35 C. The results for the ac B binary mixture are shown as orange (DLS-based microrheology) and yellow (MPT-based microrheology) circles.Representative g 2 ðq; t Þ À 1 functions from DLS-based microrheology, as well as mean-squared displacements and Van Hove functions from multiple particle tracking (MPT)-based microrheology are reported in Fig.S1in supplementary material.

FIG. 1 .
FIG. 1.(a) Phase diagram for three ac B binary mixtures with volume fraction ratios of / a : / c B ¼ 3 : 1 (blue diamonds), 1 : 1 (red squares), and 1 : 3 (black circles) in D 2 O phosphate buffer at pH ¼ 7.1.Shown are the cloud point temperatures T cloud as a function of the total volume fraction / acB .Also shown is the binodal for pure c B (gray triangles and gray solid line) (data for pure c B crystallin taken from Ref. 25).Data for / a : / c B ¼ 1 : 1 are from two different protein preparation batches (solid red squares batch 1, open red squares batch 2), and the data from batch 2 are shifted by 2.6 K in order to accommodate age and batch related shifts in the critical temperature of c B crystallin.(b) Same data shown as a function of / c B ;eff .

FIG. 2 .
FIG. 2. (a) and (b).g r vs volume fraction / for the pure crystallins (taken from Refs.49 and 61) and the ac B binary mixture in D 2 O, taken at T ¼ 25 and 35 C. The colored dashed lines represent the parameterized results for the pure crystallins (see supplementary material for details).The black dashed line shows the results of the theoretical prediction for the binary mixtures [Eq.(10)].(c) and (d).Relative viscosity g r vs protein total concentration for pure a and c B crystallins and ac B binary mixture at T ¼ 25 and 35 C, respectively.Also shown are the locations of the arrest lines for pure a and c B crystallin solutions.
Figure4shows the concentration-normalized form factor I FF ac B ðqÞ from SAXS for dilute solutions of ac B binary mixtures at different a:c B mixing ratios, i.e., 1:3 (red circles), 1:1 (green squares), and 3:1 (blue triangles).At this dilute volume fraction of / ac B ¼ 0:003, effects from protein-protein interactions can be neglected.This is indeed demonstrated in Fig.4, where we see that the scattering from the mixed solutions is quantitatively reproduced by a linear combination (colormatched lines) of the individual form factors for the pure a (purple line) and c B (orange line) crystallins also shown in Fig.4.This also clearly indicates the absence of any self-association under these conditions.We next characterize the concentration and temperature dependence of the solution structure for 1:1 mixtures using SAXS.Figures5(a) and 5(b) show the experimental SAXS curves and their theoretical predictions from Eq. (4) for a : c B ¼ 1 : 1 binary mixtures at / ac B ¼ 0:22 and two temperatures T ¼ 25 (a) and T ¼ 35 C (b), respectively.The theoretical predictions are based on the scattering intensities measured for the pure a and c B crystallin solutions using the mixing rule in Eq. (4) and their effective volume fractions / a;eff and / c B ;eff [Eqs.(1) and (2)], respectively.This requires theoretical expressions for the structure factors of the pure components S / a;eff ðqÞ and S / c B ;eff ðqÞ.For the structure factor of a crystallin, we use the Percus-Yevick approximation for polydisperse hard spheres, where the particle size distribution is described with a Schulz distribution with an average particle diameter of 15 nm and a polydispersity of 20%, which has been shown to describe S a ðqÞ with very good accuracy on the entire accessible volume fraction range up to the arrest transition at / g;a % 0:58.24For pure a crystallin solutions, we moreover expect no measurable temperature dependence of S / a ðqÞ for the temperatures investigated.For c B crystallin solutions, the situation is more complex, as the additional attraction between proteins results in a very strong and non-monotonic temperature and concentration dependence of S / c B ðqÞ

FIG. 5 .
FIG. 5. Concentration-normalized scattering intensity I ac B ðqÞ=c ac B from SAXS (solid black circles) and SLS (open black circle) for a : c B ¼ 1 : 1 binary mixtures at a total volume fraction / ¼ 0:22 for T ¼ 35 (a) and T ¼ 25 C (b).Also shown are the predicted qÀ dependence from Eq. (4) (solid black line), the case for no interactions given by Eq. (11) (dashed red line), and the contributions from a (dashed green line) and c B (dashed blue line) given by I a ðqÞ=c ac B ¼ c a Á I FF a ðqÞ Á S / a;eff ðqÞ= c ac B and I c B ðqÞ=c ac B ¼ c c B Á I FF c B ðqÞ Á S / c B ;eff ðqÞ=c ac B .

FIG. 6 .
FIG. 6.The measured effective structure factors S eff ðqÞ ¼ ½I ac B ðqÞ=c ac B =I FF ac B ðqÞ together with the predictions based on Eq. (8) where we use the full q-dependence instead of only the q ¼ 0 limit.Open (SLS) and solid (SAXS) symbols are for the experimental results, solid lines for the predictions given by Eq. (8).T ¼ 35 C: red color and T ¼ 25 C black color.
FIG. 6.The measured effective structure factors S eff ðqÞ ¼ ½I ac B ðqÞ=c ac B =I FF ac B ðqÞ together with the predictions based on Eq. (8) where we use the full q-dependence instead of only the q ¼ 0 limit.Open (SLS) and solid (SAXS) symbols are for the experimental results, solid lines for the predictions given by Eq. (8).T ¼ 35 C: red color and T ¼ 25 C black color.
FIG.8.Diffusion coefficients D(q) from neutron spin echo spectroscopy in mixtures of a and c crystallin as a function of total volume fraction and composition for two qvalues (black symbols: q ¼ 0.8 nm À1 ; red symbols: q ¼ 2.2 nm À1 .Solid symbols are measured and open symbols calculated values.The dashed lines connecting the calculated values are drawn as guides to the eye only.All data were measured at 35 C. (a) / ¼ 0:11 and (b): / ¼ 0:22.
/ a : / c B I a ðqÞ=I c B ðqÞ I a ðqÞ=I c B ðqÞ I a ðqÞ=I c B ðqÞ I a ðqÞ=I c B ðqÞ / ¼ 0:11 / ¼ 0:11

TABLE I .
Glossary of used variables for pure and mixed solutions.

TABLE II .
Ratio of the total scattering contribution from a crystallin to that of c B crystallin, I a ðqÞ=I c B ðqÞ, calculated using I a ðqÞ ¼ cðaÞI FF a ðqÞS / a;eff ðqÞ andI c B ðqÞ ¼ cðc B ÞI FF c B ðqÞS / c B ;effðqÞ, for different mixing ratios, at two total volume fractions and q-values.
2 nm À1 effective volume fractions only, i.e., we treat the other protein species as immobile obstacles that reduce the available free volume but do not contribute to hydrodynamic interactions.