Observation of liquid–liquid phase separation for eye lens γS-crystallin (original) (raw)

Abstract

γS-crystallin (γS) is an important human and bovine eye lens protein involved in maintaining the transparency of the eye. By adding small amounts of polyethylene glycol (PEG) to the binary aqueous bovine γS solutions, we have observed liquid–liquid phase separation (LLPS) at −8°C and revealed that, in the binary γS–water system, this phase transition would occur at −28°C. We have measured both the effect of PEG concentration on the LLPS temperature and protein/PEG partitioning between the two liquid coexisting phases. We use our measurements of protein/PEG partitioning to determine the nature and the magnitude of the γS-PEG interactions and to quantitatively assess the effectiveness of PEG as a crystallizing agent for γS. We use our measurements of LLPS temperature as a function of protein and PEG concentration to successfully determine the location of the critical point for the binary γS-water system. This phase transition cannot be observed in the absence of PEG because it is inaccessible due to the freezing of the system. Our findings indicate that the effective interactions between γS molecules in the binary γS–water system are attractive. We compare the magnitude of the attraction found for γS with the results obtained for the other γ-crystallins for which the critical temperature is located above the freezing point of the system. This work suggests that PEG can be used to reveal the existence of LLPS for a much wider range of binary protein–water systems than known previously.

Keywords: protein‖cataract‖polyethylene glycol‖phase transition‖partitioning


In protein condensation diseases the primary initial step in pathogenesis is the loss of solubility of the proteins, resulting in the formation of a condensed phase. Examples of such phases include a dense liquid phase, amorphous aggregates, or crystals that form in cataract (1, 2) or the formation of fibers that are responsible for sickle cell disease (3) and Alzheimer's disease (4). In all cases the attraction between proteins is the driving force for the onset of these pathologies. Thus, the characterization of the interactions between proteins is the first step for building strategies aimed at inhibiting such diseases (5).

Liquid–liquid phase separation (LLPS), which is induced by protein–protein attractive interactions, has been directly implicated in cataract. It has been shown that this is the mechanism for lens opacification in cold cataract and that it is also operative in galactosaemic and x-ray cataracts and in hereditary cataracts in mice (1, 6). The lens proteins primarily involved in phase separation are the γ-crystallins, a family of lens-specific monomeric proteins (1). Previous LLPS studies of binary aqueous solutions of four pure individual members of the bovine γ-crystallins have shown that these proteins fall into two groups: those that show an upper consolute point at ≈38°C (γE and γF) and those that show an upper consolute point at ≈5°C (γB and γD) (7, 8).

Another important member of the γ-crystallin family is γS-crystallin (γS). Although of approximately the same size (hydrodynamic radius ≈2.5 nm) and molecular mass (≈21 kDa), the expression of γS is different from that of the other proteins. γS is expressed after birth and its synthesis increases with age, whereas in bovine and human lenses the other γ-crystallins are expressed mainly prenatally and their synthesis sharply decreases after birth. It has also been found that the phase behavior of binary aqueous solutions of pure γS differs markedly from the other γ-crystallins. Specifically, γS solutions do not exhibit LLPS, even at temperatures as low as −10°C. Furthermore, the presence of γS in γB and γF aqueous solutions significantly lowers the phase separation temperatures of these proteins in the high concentration region. These facts suggest that the presence of γS in the lens helps to maintain transparency (9, 10).

We have experimentally and theoretically shown that polyethylene glycol (PEG) can be used to reveal LLPS in aqueous protein solutions (11). LLPS of binary protein–water systems may not be directly observable because it is inaccessible because of either the freezing of the system or the denaturation of the protein at high temperature (11, 12). The experiments with PEG can be used to deduce the location of the phase boundaries of the binary system. Here we show that LLPS of bovine γS aqueous solutions can be induced on adding a small amount of PEG. We have measured the phase boundaries for LLPS in the ternary γS–PEG–water system. These results are used to successfully locate the liquid–liquid phase boundary for the binary γS–water below the freezing point of the system. Therefore, we can determine that the interactions between γS molecules are attractive and compare the magnitude of the attraction with the other γ-crystallins.

The LLPS properties of the ternary protein–PEG–water system are described by a coexistence surface in the phase diagram. This coexistence surface represents the LLPS temperature, _T_ph, as a function of protein concentration, _c_1, and PEG concentration, _c_2. At a fixed temperature, the LLPS properties of the system are described by an isothermal coexistence curve, which gives the concentrations (cInline graphic, cInline graphic) and (cInline graphic, cInline graphic) of the coexisting phases, I and II. The partitioning of the components in the two coexisting phases is described by tie lines connecting the points (cInline graphic, cInline graphic) and (cInline graphic, cInline graphic) of the coexistence curve. The critical point, (cInline graphic, cInline graphic), is defined as the point on the coexistence curve where the condition (cInline graphic, cInline graphic) = (cInline graphic, cInline graphic) occurs. The location of this point as a function of the temperature is described by a critical line on the coexistence surface (11).

Materials and Methods

Bovine γS, which has a molecular weight of 20,800 g/mol, was isolated from 1- to 6-week-old calf lenses, obtained by overnight express from Antech (Tyler, TX). Pure γS was isolated and purified from these lenses by standard procedures as described (710). The purity of the native sample was at least 98%, based on both ion-exchange and size-exclusion HPLC. To suppress aggregation of γS, the purified protein was stored at 4°C in acetate buffer (0.2 M, pH 5) that contained DTT (0.02 M) and sodium azide (0.02%) (10, 13). PEG with average molecular weight of 1,450 g/mol (PEG1450) was purchased from Sigma-Aldrich and used without further purification.

Protein–PEG aqueous solutions were prepared as follows. The purified γS was dialyzed exhaustively into sodium phosphate buffer (0.1 M, pH 7.1) that contained sodium azide (0.02%). Solutions containing dilute γS in phosphate buffer were concentrated by ultrafiltration. When the desired protein concentration was reached, a known weight of PEG was added to the protein solution. The concentration of γS in the samples was determined by UV absorption at 280 nm, using the extinction coefficient value of 1.85 mg−1⋅ml⋅cm−1 calculated from the sequence (www.expasy.ch). The concentration of PEG in the samples was calculated by using the mass of PEG and the total volume of the solution. The volumetric contribution of each component was obtained by multiplying the mass of the component by the corresponding specific volume, i.e., 0.71 ml/g (14) for the protein, 0.84 ml/g for PEG (15), and 0.992 ml/g for the buffer (11).

The LLPS temperature for a given protein–PEG aqueous solution was determined by gradually lowering the temperature of the sample until clouding was observed. The cloud point was determined by the examination of the transmitted intensity-temperature profile as described by Liu et al. (9).

The coexisting phases were obtained by quenching the sample at fixed temperature below the coexistence surface. If after about 24 h the two coexisting phases had not separated by gravity, centrifugation was used for the separation. The protein concentration in each phase was determined by UV absorption. To determine the PEG concentration in each of the coexisting phases, an aliquot of known weight was taken from each phase, and the PEG was separated from the protein by ultrafiltration (Amicon Microcon YM-10). The concentration of PEG in the filtrate of known volume was determined by using a standardized refractive index detector (Perkin–Elmer LC-30 RI). Thus, the concentration of PEG in each aliquot could be calculated from its measured value in the filtrate. To account for the small effects of the two buffer salts (dibasic and monobasic sodium phosphate), the refractive index measurements were performed after isocratic elution of the filtrate solution on a size exclusion HPLC column (10 × 300 mm, Superdex 75 HR 10/30 from Amersham Pharmacia). The procedure was verified with protein–PEG aqueous solutions of known composition. In all cases, the measured protein and PEG concentrations in the two coexisting phases were consistent with the protein and PEG concentrations in the original, homogeneous samples.

Results

Our protein–PEG–buffer solutions can be described essentially as ternary protein–PEG–water systems. This approximation is justified by the insignificant difference in buffer concentration measured between the two coexisting phases.

We have previously shown that it is convenient to present our results in terms of the protein (component 1) concentration _c_1, and the PEG (component 2) concentration _c_2s inside the protein-free volume. If v is the protein specific volume and φ = c_1_v is the protein volume fraction, then _c_2s is equal to _c_2/(1 − φ). The choice of _c_2s (instead of _c_2) is convenient because of the following argument. We consider as a reference system one in which the particles of component 2 have a negligible size and do not introduce specific chemical interactions in the system. In this reference case, the different values of _c_2 between the two coexisting phases are entirely determined by the exclusion of component 2 from the volume occupied by component 1. Here, the condition of chemical equilibrium is expressed by the relation cInline graphic = cInline graphic. Thus, in our systems, any observed difference in _c_2s between two coexisting phases is related to the finite size and specific chemical properties of the PEG molecules (11).

In Fig. 1, we present our experimental data (solid circles) of the coexisting values (cInline graphic, cInline graphic) and (cInline graphic, cInline graphic) for the γS–PEG1450–water solutions at 265.1 K. The pairs of points representing the coexisting phases are connected by straight lines (tie lines). The dashed curve is an eye guide representing the coexistence curve. The coexisting phases were obtained by quenching samples with the same initial protein concentration of ≈300 mg/ml but with different initial PEG concentrations. Because the protein critical concentration, cInline graphic, is equal to the average concentration (cInline graphic + cInline graphic)/2 in the limit of cInline graphiccInline graphic = 0, we determined cInline graphic by plotting (cInline graphic + cInline graphic)/2 as a function of cInline graphiccInline graphic. The resulting value of the critical concentration (vertical bar in Fig. 1) is 260 ± 10 mg/ml, which corresponds to a critical protein volume fraction: φc = 0.18 ± 0.01. To quantify PEG partitioning, we consider the normalized slope of the tie lines at the critical point (∂_c_2s/∂_c_1)Inline graphic/cInline graphic. This quantity was obtained by plotting the incremental ratio (ln cInline graphic − ln cInline graphic)/(cInline graphiccInline graphic) as a function of cInline graphiccInline graphic and considering the limit of cInline graphiccInline graphic = 0. We found (∂_c_2s/∂_c_1)Inline graphic/cInline graphic = (−2.1 ± 0.1) 10−3⋅mg−1⋅ml.

Figure 1.

Figure 1

Coexisting surfaces at constant temperature for the ternary γS–PEG1450–water system. The pairs of points representing the coexisting phases (circles) are connected by the tie lines (solid lines). The dashed curve is a guide for the eye and the vertical bar locates the critical point at cInline graphic = 260 mg/ml. The normalized slope of the tie lines at the critical point is (∂_c_2s/∂_c_1)Inline graphic/cInline graphic = (−2.1 ± 0.1)10−3 mg−1⋅ml.

In Fig. 2, we present our measurements of the LLPS temperature, _T_ph, for the ternary γS–PEG1450–water solutions at constant PEG concentration of _c_2s = 74 mg/ml (circles), _c_2s = 66 mg/ml (squares), _c_2s = 60 mg/ml (diamonds), and _c_2s = 53 mg/ml (triangles). In all cases, the PEG concentration is smaller than the semidilute crossover concentration _c_2s = 80 mg/ml (16). In Fig. 2, we can see that _T_ph linearly increases as _c_1 increases. We can fit each experimental data set with the linear function _T_ph = _T_c(_c_2s) + (∂_T_ph/∂_c_1)Inline graphic(_c_1 − cInline graphic) (solid lines in Fig. 2). Using cInline graphic = 260 mg/ml we obtained the values for both the critical temperatures, _T_c(_c_2s), and the slopes (∂_T_ph/∂_c_1)Inline graphic. In Fig. 3, we plot _T_c(_c_2s) as a function of _c_2s. We observe that _T_c is approximately a linear function of _c_2s with (∂_T_ph/∂_c_2s)_c_1 = (0.27 ± 0.03) K⋅mg−1⋅ml. In Fig. 4, we plot (∂_T_ph/∂_c_1)Inline graphic as a function of _c_2s. We see that (∂_T_ph/∂_c_1)Inline graphic increases linearly with _c_2s and starts at zero when _c_2s = 0 as it is required for a binary protein–water system (9). We found (∂_T_ph/∂_c_1)Inline graphic/cInline graphic = (0.674 ± 0.004)10−3 K⋅mg−2⋅ml2.

Figure 2.

Figure 2

LLPS temperature, _T_ph, at constant PEG1450 concentration of _c_2s = 74 mg/ml (circles), _c_2 = 66 mg/ml (squares), _c_2s = 60 mg/ml (diamonds), and _c_2s = 53 mg/ml (triangles). The solid lines are linear fits to the experimental data. The protein critical concentration is cInline graphic = 260 mg/ml.

Figure 3.

Figure 3

Critical temperature, _T_c, for the ternary γS–PEG1450–water system as a function of _c_2s. The dashed line is a linear fit to the experimental data with slope (∂_T_ph/∂_c_2s)_c_1 = (0.27 ± 0.03) K⋅mg−1⋅ml.

Figure 4.

Figure 4

Slope of the critical temperature, (∂_T_ph/∂_c_1)Inline graphic/_c_2s, for the ternary γS–PEG1450–water system as a function of _c_2s. The dashed line is a linear fit to the experimental data with slope (∂_T_ph/∂_c_1)Inline graphic/_c_2s = (0.674 ± 0.004) 10−3 K⋅mg−2⋅ml2.

The behavior of the coexistence surface is described by the three slopes (∂_c_2s/∂_c_1)Tph, (∂_T_ph/∂_c_2s)_c_1, and (∂_T_ph/∂_c_1)_c_2s, which are related to each other by the mathematical relationship (9, 11):

graphic file with name M51.gif 1

We apply Eq. 1 on the critical line and use the measured values of (∂_c_2s/∂_c_1)Inline graphic/cInline graphic and (∂_T_ph/∂_c_1)Inline graphic/cInline graphic to determine that (∂_T_ph/∂_c_2s)Inline graphic = 0.320 ± 0.013 K⋅mg−1⋅ml. This value is in satisfactory agreement with the value obtained by the direct analysis of the dependence of _T_c on _c_2s shown in Fig. 3 [i.e., (∂_T_ph/∂_c_2s)c1 = (0.27 ± 0.03) K⋅mg−1⋅ml].

It is also convenient to describe the essential features of the coexistence surface within our experimental domain in terms of the variables _c_1 and _c_2 as shown in Fig. 5. All our data are taken in the vicinity of the point on the coexistence surface specified by the coordinates _c_1 = 260 mg/ml, _c_2 = 53 mg/ml, and _T_ph = 265 K. At this point, the values of the two principal slopes are (∂_T_ph/∂_c_2)Inline graphic = (0.38 ± 0.04) K⋅mg−1⋅ml and (∂_c_2/∂_c_1)Inline graphic/cInline graphic = (−3.0 ± 0.1)10−3 mg−1⋅ml.

Figure 5.

Figure 5

Coexistence surface, _T_ph(_c_1,_c_2), for the ternary γS–PEG1450–water system. Our measurements of _T_ph(_c_1,_c_2) reported in Fig. 2 cover the range 220 < _c_1 (mg/ml) < 360, 40 < _c_2 (mg/ml) < 60. A representative solid curve at constant PEG concentration is generated by the intersection of the coexistence surface and the plane at _c_2 = 53 mg/ml. Our concentration measurements on the coexisting phases reported in Fig. 1 are represented by the isothermal coexistence curve (solid curve) generated by the intersection of the coexistence surface with the plane at _T_ph = 265 K. The three tie lines at _T_ph = 265 K are calculated from Eq. 4 by using q = 0.36. The circle at _c_2 = 0 locates the critical point for the binary γS–water system. The location of the critical point as a function of temperature is described by the critical line (solid curve at _c_1 = 260 mg/ml).

We can estimate the critical temperature, TInline graphic, for the binary γS–water system by using the relationship TInline graphic = Tc − (∂_T_ph/∂_c_2)Inline graphic_c_2 to the critical line. This expression assumes that both the protein critical concentration, cInline graphic, and the slope, (∂_T_ph/∂_c_2)Inline graphic, are not functions of PEG concentration. This assumption was demonstrated to be correct for the ternary γD–PEG–water systems when PEG concentration is ≈50 mg/ml (11). We obtain TInline graphic ≈ 245 K for the binary γS–water system. In Fig. 5, the critical point for the binary γS–water system is located by the circle at _c_2 = 0.

In Table 1, we report the value of the critical temperature, TInline graphic, and the critical volume fraction, φ_c_, for the binary γS–water system. We also report the corresponding values for the high TInline graphic γ-crystallins (γE and γF) and the low TInline graphic γ-crystallins (γB and γD) (8). All of the proteins in Table 1 have essentially the same critical volume fraction (φ_c_ ≈ 0.2), indicating that the range of the inter-protein interactions is essentially the same for all five γ-crystallins (17). However, the difference in critical temperature of ≈66 K between the high TInline graphic proteins and γS, and the difference of ≈33 K between the low TInline graphic proteins and γS is the result of a corresponding difference between the magnitudes of the attraction energy. We therefore conclude that the magnitude of the attraction between γS molecules is, respectively, ≈20% and ≈10% smaller than the molecular attraction found for the high TInline graphic and the low TInline graphic proteins.

Table 1.

Critical points for the bovine γ-crystallins

γ-Crystallins φc TInline graphic, K
γE 0.22 ± 0.01 313
γF 0.22 ± 0.01 309
γB 0.19 ± 0.01 278
γD 0.20 ± 0.01 278
γS 0.18 ± 0.01 245

Discussion

The experimentally observable feature of the coexistence surface, (∂_c_2/∂_c_1)Inline graphic/cInline graphic, is fully determined by the effective protein and PEG chemical potentials μ1 and μ2 (11). If PEG concentration is relatively small, these two chemical potentials can be expressed in the following way (11):

graphic file with name M75.gif 2a
graphic file with name M76.gif 2b

where μ̂i ≡ (μ_i_ − μInline graphic)/RT (with i = 1,2), μInline graphic is the standard part of the chemical potential, T is the temperature, and R is the ideal gas constant. The quantity μ̂′1(_c_1, T) ≡ μ̂1(_c_1, 0, T) corresponds to the binary protein–water system and ξ(_c_1, T) is the first term in a series expansion describing the contribution to μ̂2 over and above the ideal mixing entropy. The two cross-derivatives of the chemical potentials, (∂μ̂1/∂_c_2)_c_1,T and (∂μ̂2/∂_c_1)_c_2,T, are both equal to (∂ξ/∂_c_1)T. This quantity characterizes the protein–PEG interactions (11).

The Slope of the Tie Lines.

If we differentiate Eq. 2a with respect to _c_2 at constant _c_1 and T, and Eq. 2b with respect to _c_1 at constant μ̂2 and T, we obtain

graphic file with name M79.gif 3

The normalized derivative, (∂_c_2/∂_c_1)μ̂2,T coincides with the limiting value of the slope of the tie lines at the critical point and has been experimentally measured. This quantity determines the effectiveness of PEG as a crystallizing agent for γS. Indeed, if protein crystals (the solid phase) are in thermodynamic equilibrium with the liquid phase, μ̂1 must be equal to its value in the solid phase, μ̂Inline graphic. For relatively low protein concentrations, μ̂′1(_c_1, T) ≈ ln _c_1. Thus, using the equality of chemical potentials and Eq. 2a, we find ln _c_1 ≈ μ̂Inline graphic − (∂ξ/∂_c_1)T c2. This relationship gives the protein solubility, _c_1, as a function of PEG concentration, _c_2. The dependence of μ̂Inline graphic on _c_1 and _c_2 is expected to be very weak. Thus the fractional decrease of protein solubility with PEG concentration is fully determined by the value of (∂ξ/∂_c_1)T, which is provided by the slope of the tie line (see Eq. 3) (11). Thus, PEG1450 reduces the solubility of γS. We therefore believe that PEG is a useful tool not only for inducing LLPS of γS solutions but also for favoring the formation of γS crystals by solubility reduction.

The Excluded Volume Model of Protein–PEG Interactions.

We now analyze the experimental coexisting values (φI, cInline graphic) and (φII, cInline graphic) to gain insight into the physical factors that are responsible for the γS–PEG1450 interactions. We apply a simple excluded volume model that directly relates PEG partitioning to the difference in free volume fractions between the two coexisting phases. If we assume that PEG molecules can be described as ideal polymer coils and that protein molecules have a spherical shape, then due to steric hindrance, each protein will be surrounded by an adjacent region where the centers of mass of the coils are excluded, and the width, δ, of the resulting depletion layer will be proportional to the gyration radius, _R_s, of the polymer coil (1823). If α is the volume fraction available to the centers of mass of the coils, the condition of chemical equilibrium becomes cInline graphic/αI = cInline graphic/αII, where _c_2/α is the polymer concentration in the free volume (21). An approximate expression for the free volume fraction, α, as a function of the hard sphere volume fraction, φ, is given by the well-established scaled particle theory (2426):

graphic file with name M87.gif 4

where η = φ/(1 − φ), A = 3_q_ + 3_q_2 + _q_3, B = 9_q_2/2 + 3_q_3, C = 3_q_3, and q = δ/R is the depletion layer thickness normalized with respect to the hard sphere radius, R. The prefactor, 1 − φ, in Eq. 4, is the volume fraction not occupied by the spheres, whereas the exponential factor describes the effect of the depletion layers.

Because cInline graphic/cInline graphic = α(q, φII)/α(q, φI), we substitute the expression for α provided by Eq. 4 into the ratio of the free volume fractions and thereby can determine an apparent q value for each experimental tie line. We find that, within the experimental error, this quantity does not depend on the tie-line position in the coexistence curve and report the average apparent value of = 0.36 ± 0.01. On the other hand, theoretically we expect that q = k(_R_g/R) where k, in the case of hard sphere-ideal chain interactions, is a known function of _R_g/R (23). For PEG1450, _R_g = 1.56 nm (14), whereas R = 1.80 nm is estimated from the molecular weight and specific volume of γS. Because _R_g/R = 0.87, the theoretical value of k is 1.03 (23) and the theoretical value of q is 0.89. Therefore, the apparent value of q is numerically ≈2.5 times smaller than the corresponding theoretical values. This significant discrepancy between the apparent and the theoretical q was also observed in the case of ternary γD–PEG water systems. It appears that weak attractive interactions between PEG and γS induce a flattening of the polymer coils near the protein surface and a corresponding reduction of the depletion layer thickness (11).

Summary and Conclusions

We have experimentally observed the liquid–liquid phase transition in ternary γS–PEG–water ternary systems. We have used the experimental coexisting values (φI, cInline graphic) and (φII, cInline graphic) together with the slope, (∂c_2/∂_c_1)Inline graphic/cInline graphic, to determine the γS–PEG interactions. We have used the experimental observable feature of the coexistence surface, (∂_T_ph/∂_c_2_s)_c_1, to successfully determine the location of the upper consolute point for the binary γS–water system. This phase transition could not previously be observed in absence of PEG because it is inaccessible due to the freezing of the system.

It is possible to argue that liquid–liquid transitions have not been observed for many binary protein–water systems because they fall outside the accessible experimental temperature range (from ≈−10°C to ≈40°C). If this is the case, PEG can be used to reveal the existence of LLPS for a much wider range of proteins than known previously.

Acknowledgments

We acknowledge Neer Asherie, Ajay Pande, Jayanti Pande, Alexander Chernov, and Seth Fraden for useful comments. This work was supported by National Aeronautics and Space Administration Grant NAG8-1659 and National Institutes of Health Grant EY05217.

Abbreviations

PEG

polyethylene glycol

LLPS

liquid–liquid phase separation

γS

γS-crystallin

References