The Thermal Structural Transition of α-Crystallin Inhibits the Heat Induced Self-Aggregation

-crystallin, the major constituent of human lens, is a member of the heat-shock proteins family and it is known to have a quaternary structural transition at . The presence of calcium ions and/or temperature changes induce supramolecular self-aggregation, a process of relevance in the cataractogenesis. Here we investigate the potential effect of the bovine -crystallin's structural transition on the self-aggregation process. Along all the temperatures investigated, aggregation proceeds by forming intermediate molecular assemblies that successively aggregate in clusters. The final morphology of the aggregates, above and below , is similar, but the aggregation kinetics are completely different. The size of the intermediate molecular assemblies, and their repulsive energy barrier show a marked increase while crossing . Our results highlight the key role of heat modified form of -crystallin in protecting from aggregation and preserving the transparency of the lens under hyperthermic conditions.


Introduction
Cataract, the opacity of the eye lens, is an age-onset pathology that affects nearly 50% of the world's population over the age of 65, and is the leading cause of blindness worldwide [1]. Pathological studies of cataractous lenses have revealed that cataracts are composed of protein aggregates that precipitate in eye lens cells. The prevalent proteins within the eye lens are the crystallins. Lens transparency is thought to be maintained by a liquid-like, short range order present in highly concentrated solutions of these proteins [2,3]. In mammals, there are three classes of crystallins denoted a, b, and c of which a-crystallin is the most abundant. acrystallin is an oligomeric protein, having a molecular mass of about 80071200kDa, composed of two types of subunits, aA and aB crystallins, each having a molecular mass of about 20 kDa and aA crystallins in a ratio of 3 to 1 with respect to aB [4]. The increase in light scattering in old and cataractous lenses can be ascribed to alterations in lens crystallins interactions due to age related posttranslational modification of a-crystallin [5][6][7][8]. The alterations are triggered by lens cells exposition to elevated temperatures or other stress factors like Ca 2z ions, that disrupt the liquid-like molecular order and promote the formation of large scattering particles [9,10], following pathways that include both changes in the secondary structure and in the state of assembly [11,12].
To monitor the heat and Ca 2z induced changes that occur in the structural domain of lens-crystallin different techniques like circular dichroism, fluorescence, Small Angle X-Ray scattering (SAXS) were used [13,14]. Differential scanning calorimetric studies on a-crystallin [14] show two endothermic transitions, a first ranging from 35 0 C to 51 0 C, peaked at T c~4 5 0 C and a second major transition peaked at T II C~6 0 0 C. The transition at T c~4 5 0 C has been found to be biologically relevant [15]. At this temperature, a-crystallin undergoes a minor change in its tertiary structure accompanying the exposure of its hydrophobic surfaces [16,17], whereas its secondary structure is relatively unchanged.
Here we focused on the effects of this structural transition on the a-crystallin self-aggregation. Several aggregations have been induced by changing temperature, and therefore, by generating different heat-modified a-crystallin forms [18]. At temperature larger than T c the kinetic pattern of the a-crystallin aggregation and the structural features of the clusters can be described according to the reaction limited cluster-cluster aggregation theory (RLCA) [19]. Aggregation occurs by initially forming the basic aggregation units, the high molecular weight forms of a-crystallin (HMW) [20], that successively continue to diffuse, collide and form rather compact fractal aggregates (d f~2 :15). Although the final morphology of the aggregates is similar, the aggregation kinetics are completely different below and above T c , together with the size of the HMW, and their repulsive energy barrier (DE b ). An abrupt increase in (DE b ) reveals a mechanism that markedly protects from aggregation preserving the transparency of the lens.

Materials and Methods
Preparation and aggregation of a-crystallin a-crystallin from bovine eye lens was prepared according to Andreasi et al. [19].
The a-crystallin fractions suspended in 10 mM Tris-HCl buffer, pH 7.4, were thoroughly mixed and pooled together. The purified protein was divided into aliquots and kept in the same buffer at {20 0 C until used. Just before the experiment, the samples were thawed and centrifuged at 5000g(Eppendorf 5418) for 30min at 4 0 C, and the supramolecular aggregates already formed were discarded. The supernatant was filtered through a 0:22mm Millipore low-retention filter directly into the measuring cuvette. Protein concentration was determined by using an absorption coefficient of A 0:1% 1cm~0 :81 at 280 nm [21]. Aggregation of acrystallin (1.0 mg/ml) was induced by quenching samples at the desiderated temperature and by the addition of 16mM CaCl 2 . Indeed heating provokes the generation of particularly reactive isoforms of a-crystallin [14], and calcium ions stabilize the aggregates while they are forming and allow their continuous growth [9]. The whole set of measurements have been performed on different aliquots of the same sample. Five aggregations process for each temperature have been followed.

Dynamic light scattering
Dynamic light scattering (DLS) (24) provides information on the aggregation kinetics and on the clusters dimension and evolution as the aggregation proceeds. DLS measurements were performed during aggregation by using a commercial computer-interfaced scattering instruments ALV/SLS-5000 system from ALV, Langen, Germany, equipped with a 50mW HeNe laser operating at 632.8 nm. The beam was focused to a spot (at 1=e 2 ) of *50mm and the detection of the scattered light was carried out with a monomode fiber coupled to a photomultiplier, both mounted on a stepping motor-controlled rotating arm. The sample was contained in a cylindrical cell (8mm inner diameter) immersed in a toluenefilled index-matching vat whose temperature was controlled with a resolution of 0:01 0 C. In this instrument, the usable wavevector range varies from *4 : 10 4 to *3 : 10 5 cm {1 . The autocorrelation function of the photopulses were performed by a 256-channel digital correlator (ALV-5000). Counts per second were used to measure the scattered intensity during the aggregation. DLS technique measures the intensity autocorrelation function g 2 (t)~SI(t)I(tzt)T=SIT 2 where t is the lag time and brackets represent the ensemble average [22]. The g 2 (t) can be related to the field autocorrelation function g 1 (t) through the Siegert relation g 2 (t)~1zbg 2 1 (t) where b is an instrumental constant (in our set-up b~1). The mathematical form of g 1 (t) depends on the physical properties of the investigated system. For monodisperse particles, the electric field autocorrelation function decays exponentially following g 1 (t)~e {Ct ,where the decay rate C depends on the particle translational diffusion coefficient according to C~Dq 2 . For a polydisperse sample, g 1 (t) is more complex than a single exponential. In this condition, the derivative of g 1 measures the intensity weighted average decay rate of the clusters: To determine C I experimentally, we fitted the logarithm of the measured autocorrelation function g 1 , to a second-order polynomial, according to the cumulant expansion [23]: Where we assumed C I~C1 . In aggregating systems, because of cluster-mass polydispersity, what we actually measure is an average effective diffusion coefficient that can be expressed as: where p I (r) is the normalized intensity-weighted radius distribution function describing the distribution of the fraction of the intensity scattered by a particle of hydrodynamic radius r and decay rate C(r), given by: where g is the water viscosity and k the Boltzmann constant. The intensity weighted average effective hydrodynamic radius r I can be obtained using Stokes-Einstein Relation The complete distribution of decay rates can also be recovered by introducing from the relation [22,24]: The recovery of the p I (r) distribution, a classical ill posed problem, can be obtained by performing a regularized Laplace inversion of the intensity autocorrelation function using the software CONTIN [25]. In the presence of highly polydisperse fractal clusters, we also have to account for the volume and for the inner structure of clusters, therefore we need a complete recovery of the normalized number-weighted radius distribution function p N (r) [24]: where M(r) is the mass of a cluster of hydrodynamic radius r and P(qr) is the scattering form factor of the particle. Finally once p N (r) is known, the mean hydrodynamic radius can be easily determined by:

Reaction limited cluster aggregation
A key to the understanding of proteins aggregation is the behavior of the energy of interaction between two approaching particles. It has been demonstrated that for a wide variety of proteins, this can be understood within the Derjaguin-Landau-Verwey-Overbeek (DLVO) model [26]. The key parameter is the repulsive energy barrier between two approaching particles. If the height of this energy barrier E b is sufficiently large compared to k B T, the particles will be unable to stick to one another when their diffusive motion causes them to collide, and the particle will be stable against aggregation. If E b is reduced to much less than k B T, every collision will result in the particles sticking together, leading to very rapid aggregation, limited only by the rate of diffusion induced collisions between the clusters. This regime is therefore called diffusion-limited colloid aggregation (DLCA) [27][28][29]. For DLCA, computer simulation and several different experimental techniques show that clusters, characterized by a typical fractal dimension d f~1 :8, are essentially monodisperse in that their mass distribution is bell shaped and peaked around an average mass value, which grows linearly with time. By contrast, if E b remains comparable to, or larger than, k B T, many collisions must occur before two particles can stick to one another. In this case, the aggregation rate is limited by the probability of overcoming the repulsive barrier P*exp({E b =k B T) (sticking probability), leading to much slower aggregation. In this regime, called reaction limited aggregation (RLCA) [30][31][32], clusters have a structure more dense than in the DLCA, and with the typical d f~2 :1. In each case, however, as particles stick together to become clusters, the clusters themselves continue to diffuse, collide, and aggregate. Each of these two regimes is characterized by a different time evolution of the average cluster mass, M, of the shape of clustermass distribution function, and of the fractal dimension of the resulting clusters. In RLCA regime the average cluster radius is an exponential function of time [19,30,33]: where r 0 is the basic aggregation units and the aggregation rate k agg is a constant that depends on the sticking probability and therefore on the repulsive energy barrier [30,34]: Clusters formed in the RLCA regime show an extremely high mass polydispersity, described by a power law, up to a cutoff mass M c , after which it decreases exponentially following: where t is a characteristic exponent (i.e t~1:5 in RLCA) [30][31][32].
The DLCA and RLCA regimes must be considered universal in that their features do not depend on the nature of interacting forces between particles. Spectroscopic and microscopy techniques usually detect the normalized number-weighted radius distribution function, p N (r), instead of the cluster-mass distribution function, N(M). In the case of RLCA aggregation we obtain: where r c is the cut-off radius of clusters with a mass M c .

Results
To characterize the extent of the aggregation process, we performed dynamic light scattering experiments by measuring the time evolution of the intensity weighted average hydrodynamic radius of the clusters, r I , determined according to Eq.2, Eq.3, and Eq.5. The results for samples at different temperatures above T c~4 5 0 C are reported in Fig. 1.
After an initial, fast, increase of r I a second, slower, exponential growth, is observed.
The first increase of r I is ascribed to the initial conversion of the protein from the native to the heat-and calcium-induced conformers that rapidly bind to form high molecular weight species (HMW) [19]. The second exponential growth, already investigated in detail [19], is consistent with an RLCA process where HMWs, of radius r 0 , after a large number of collisions can stick together. The time at which appears the crossing between these two steps it is called t lag .
Fits of Eq.9 to experimental data allow to recover k agg , r 0 for each aggregation process ( Table 1). The value of t lag is determined as the t value where r(t lag )~r 0 .
All the aggregations, carried out at different temperatures, show the same behavior. The initially formed basic aggregation units aggregate forming fractal clusters, accordingly to an RLCA process characterized by a temperature dependent rate constant (i.e. higher is the temperature faster is the aggregation rate). The size of basic aggregation units,instead, is independent on temperature with an average value of r 0(TwTc)~( 25:2+1:1)nm.
By decreasing temperature below T c the time evolution of the aggregation process undergoes to a dramatic modification (Fig. 2) and the fitting equation used is slightly modified according to the expression r~r 0 e kagg(t{t lag ) . Basic aggregation units are formed over a longer time and their average size is smaller (r 0(TvTc)~1 7:5+0:9) (Tab.1). At times larger than t lag an exponential increase of the hydrodynamic radius, r I , is still observed. In order to verify if aggregations below T c , although different with respect to those observed at TwT c , are still consistent with an RLCA process, we investigate the shape of the cluster-number distribution (p N (r)) that can be recovered directly from the intensity autocorrelation function [24,33]. In Fig. 3 we report representative distributions of p N (r) taken at different times during the aggregation occurred at T~39:61 0 C. Number weighted radius distributions shift toward higher values with time, and, more relevant, all distributions appear broad and highly asymmetric. Fits of Eq.12 to experimental data evidence that, in agreement with the RLCA aggregation theory, the p N (r) is well described by a power law up to the cut-off radius, r c , after which it becomes an exponential distribution. Contextually, the fractal dimension does not vary as aggregation proceeds, keeping an average value of =(2:10+0:11nm). This peculiar behavior imply an universal scaling of the number distribution, i.e. p N (r) depends only on the cut-off radius r c that increases exponentially with time as r c *e t=t0 (see inset of Fig. 3) [30], while is independent on the detailed nature of the aggregating particle. In Fig. 4 we report the whole aggregate's size distribution, recovered from aggregations occurred at two different temperatures, below (full dots) and above T c (open diamonds), when the average radius reach a specific value (r N~4 0nm). The overlap of the aggregates' size distribution curves before and after the structural transition is a further evidence that the shape of aggregates' size distribution is independent on the size of the aggregating particle (which has different dimension above and below T c ) as expected for RLCA aggregations [30].
Therefore above and below T c the aggregations are well characterized in the framework of RLCA theory: the final morphology of the aggregates is similar, but the aggregation kinetics are completely different.
A closer look of Tab.1, indeed, evidences that the aggregation rate k agg depends on T both above and below T c , with an abrupt change across T c . This jump of the k agg value is a consequence of the phase transition [35] and it can be easily visualized in Fig. 5 where the logarithm of k agg against the inverse temperature is reported. In this plot two different behaviors of k agg can be identified above and below T c (the 1/T c value is indicated by a dashed line). Fit of eq.10 to experimental data allows to recover two distinct energy barriers for the aggregation process: E b1~( 33,79+1,29) Kcal/mol and E b2~( 137,6+9,3) Kcal/mol (below and above T c respectively). The rate of formation of HMWs (k HMW ), estimated as t lag %k {1 HMW [36], reveals, instead, a single Arrhenius behavior along all the temperature range investigated (inset Fig. 5) showing that the rate of formation of intermediate aggregation units is independent on the alphacrystallin phase transition.
Discussion a-crystallin is the most abundant lens protein of the mammalian eye, and its aggregates are the main scattering elements strongly involved in the process of cataractogenesis. The increase in light scattering in old and cataractous lenses can be ascribed to Table 1. a-crystallin aggregation constants.   alterations in lens crystallins interactions due to age related posttranslational modification of a-crystallin [5][6][7][8]. The alterations are triggered by lens cells exposition to elevated temperatures or other stress factors like Ca 2z ions, that disrupts the liquid-like molecular order and promote the formation of large scattering particles [9,10]. Supramolecular structure of crystallins substantially varies both in lenses of different vertebrate species and in various parts of the same lens [37], therefore the understanding of the aggregation mechanism and the fractal clustering could be an important tool to characterize lens ageing and crystallin function.
Here we monitor changes in the a-crystallin's aggregation process induced by the thermal structural transition (T c~4 5 0 C).
At all the temperatures investigated supramolecular aggregation of a-crystallin could be described according to the reaction limited cluster-cluster aggregation theory. Aggregation of the Ca 2z and heat-modified proteins occurs initially by rapidly forming the first clusters or basic aggregating units, corresponding to HMW [20,38]. After that clusters themselves continue to diffuse, collide, and aggregate. As the aggregation proceeds, clusters with different masses are formed and stick each other. The average cluster size of aggregating HMWs increases exponentially in time and their fractal dimension *2:15 indicates that aggregates sample all the possible mutual configurations before they stick together.
The radius of the HMW is r 0 *18nm, and r 0 *25nm, below and above T c , respectively. This difference highlights a substantial alteration in the packing of a-crystallin subunits. Accordingly, HMW molecular weight increases 2:7 times and HMW concentration decreases by the same factor. A lower HMW concentration reduces the probability of collision decreasing the rate of formation of aggregates. Nevertheless, modifications in subunit's packaging do not affects the kinetic of formation of the HMWs (see inset of Fig. 5).
The aggregation rate, instead, undergoes to an overall abrupt change when crossing T c . Above T c , the energy barrier that must be crossed to create larger particles is *100Kcal=mol higher than below T c . Accordingly, the probability that an activated state occurs along the aggregation process results k agg (T~47,99)=k agg (T~41:29)^15 times lower than above T c . Therefore, the formation of large scattering particles in old and cataractous lenses is inhibited at high temperature by the structural transition that occurs at T c . As already pointed out, there is an effective inhibition only of the aggregation step: the rate of aggregation of HMW changes steeply at T c , whereas the rate of formation of HMW is unmodified at T c . Therefore, the change in tertiary structure occurring at the endothermic phase transition at T c~4 5 0 C [16,17], triggers a major reorganization of a-crystallin subunits in the HMWs. A smaller number of more stable HMW particlesis formed, and accordingly the aggregation is inhibited.
Lens crystallin is particularly recessive to deleterious effects from elecromagnetic radiations that are known to be a potential risk factor for cataract and other eyes diseases. Indeed, its aqueous content favors radiation absorption and the very weak vascularization makes difficult to stand fast temperature increases [39].
In this context, the natural self-protective mechanism that we report preserves the lens from premature opacification throughout the lifespan of the organism [40], abruptly reducing the formation of aggregates in the lens fiber cells under hyperthermic conditions, such as those determined by extended exposure to microwaves or other electromagnetic radiation emitted by cell and cordless phones, wireless communications, monitors and even high voltage lines [41][42][43][44][45]. (open diamonds) below and above T c , when the average radius reach a specific value (r N~4 0nm). The overlap of the aggregates size distribution curves before and after the structural transition is an evidence that the shape of size distribution of aggregates depends only on the cluster-cluster interaction potential (i.e DLVO in our case) while temperature simply modulate the extent of aggregation. doi:10.1371/journal.pone.0018906.g004 Figure 5. Logarithm of k agg against inverse temperature. Two different behaviors can be identified above and below T c (1/T c value is indicated by a dashed line). Fit of eq.10 to experimental data allows to determine two different values of the barrier energy: E b1~( 33,79+1,29)Kcal/mol and E b2~( 137,6+9,3) Kcal/mol for the two kinds of aggregation process (below and above T c respectively). The rate of formation k HMW of HMWs, estimated measuring t lag %k {1 HMW , reveal a normal Arrhenius behavior with temperature (inset), revealing that the same abrupt change doesn't involve the first aggregation step. doi:10.1371/journal.pone.0018906.g005