Increased surface charge in the protein chaperone Spy enhances its anti-aggregation activity

Chaperones are essential components of the protein homeo-stasis network. There is a growing interest in optimizing chaperone function, but exactly how to achieve this aim is unclear. Here, using a model chaperone, the bacterial protein Spy, we dem-onstrate that substitutions that alter the electrostatic potential of Spy ’ s concave, client-binding surface enhance Spy ’ s anti-aggrega-tion activity. We show that this strategy is more efficient than one that enhances the hydrophobicity of Spy ’ s surface. Our findings thus challenge the traditional notion that hydrophobic interactions are the major driving forces that guide chaperone – substrate binding. Kinetic data revealed that both charge- and hydrophobicity-enhanced Spy variants release clients more slowly, resulting in a greater “ holdase ” activity. However, increasing short-range hydrophobic interactions deleteriously affected Spy ’ s ability to capture substrates, thus reducing its in vitro chaperone activity toward fast-aggregating substrates. Our strategy in chaperone surface engineering therefore sought to fine-tune the different molecular forces involved in chaperone – substrate interactions rather than focusing on enhancing hydrophobic interactions. These results improve our understanding of the mechanistic basis of chaperone – client interactions and illustrate how protein surface – based mutational strategies can facilitate the rational improvement of molecular chaperones.

Chaperones are essential components of the protein homeostasis network. There is a growing interest in optimizing chaperone function, but exactly how to achieve this aim is unclear.
Here, using a model chaperone, the bacterial protein Spy, we demonstrate that substitutions that alter the electrostatic potential of Spy's concave, client-binding surface enhance Spy's anti-aggregation activity. We show that this strategy is more efficient than one that enhances the hydrophobicity of Spy's surface. Our findings thus challenge the traditional notion that hydrophobic interactions are the major driving forces that guide chaperonesubstrate binding. Kinetic data revealed that both charge-and hydrophobicity-enhanced Spy variants release clients more slowly, resulting in a greater "holdase" activity. However, increasing short-range hydrophobic interactions deleteriously affected Spy's ability to capture substrates, thus reducing its in vitro chaperone activity toward fast-aggregating substrates. Our strategy in chaperone surface engineering therefore sought to fine-tune the different molecular forces involved in chaperonesubstrate interactions rather than focusing on enhancing hydrophobic interactions. These results improve our understanding of the mechanistic basis of chaperone-client interactions and illustrate how protein surface-based mutational strategies can facilitate the rational improvement of molecular chaperones.
Molecular chaperones are evolutionarily conserved, intricate molecular machines designed to combat protein folding problems (1,2). Many are up-regulated under stress conditions and act not only to prevent protein aggregation, but also to assist in protein refolding. Some even facilitate the disaggregation of damaged proteins to allow them to refold (3). Chaperones also function during de novo protein folding, in part by decreasing nonnative interactions between folding intermediates and thus acting to minimize aggregation and promote productive folding (4,5). Although chaperone action is spatially and temporally complex, the fundamental mechanism by which chaperones recognize their substrates must ultimately be controlled by basic intermolecular forces. Because the exposure of hydrophobic segments on the surface of unfolded or partially folded proteins is critical for their aggregation (6)(7)(8), it has been assumed that molecular chaperones recognize and bind these hydrophobic regions on their substrates to protect them from aggregation (9)(10)(11). One very widely cited review for instance makes the statement that chaperones "typically recognize hydrophobic amino acid side chains exposed by nonnative proteins" (9). Very similar statements are present in many other reviews; sheer repetition has made the vital nature of hydrophobicity in chaperone action almost axiomatic.
However, in addition to hydrophobic interactions, it is important to consider that electrostatic interactions can also play an important role in the initial binding of chaperones to unfolded substrates (12)(13)(14). Electrostatic forces can speed association reactions and have been recognized as a dominant long-range force in mediating protein-protein associations (15). High-speed binding events are often correlated with favorable electrostatic interactions (16,17). Kinetic studies of the GroEL-barnase interaction, for instance, revealed that substrate association is strongly dependent on ionic strength (18). This dependence is evident only in the early stages of complex formation. Thus, for this chaperone-substrate interaction, an initial, transient ionic interaction exists before any hydrophobic binding occurs (18). Similarly, client binding to the ATP-independent chaperone Spy is also extremely sensitive to ionic strength; as the ionic strength increases, the rate of association decreases (19). When assayed in vivo, the activity of Spy also shows a strong dependence on salt concentration, further demonstrating the important role that electrostatic interactions play in Spy chaperone activity (20).
Detailed knowledge of how chaperones interact with substrates can guide efforts to optimize the performance of molecular chaperones. However, unlike with enzymes, where substrate specificity has frequently been altered by directed evolution and protein design, efforts to alter the specificity or efficiency of molecular chaperones have rarely been successfully undertaken (21). This is in part due to the lack of detailed information about substrate-chaperone interactions, because chaperones are often promiscuous and because precise substrate-chaperone binding sites are often poorly defined (12,22). Fortunately, directed evolution techniques offer an opportunity for altering a chaperone's binding properties even in the absence of detailed knowledge of the chaperone's binding sites and mechanism. These types of experiments can provide a unique perspective in efforts to understand the working mechanisms of chaperones, for example, through the analysis of how beneficial mutations contribute to function (23,24). Previously, using genetic selections to stabilize unstable variants of immunity protein 7 (Im7, a 10-kDa Escherichia coli protein with a well-characterized folding pathway (25)), we succeeded in isolating Spy variants that were able to increase the stability of these destabilized Im7 mutants. Most of these appear to work by expanding hydrophobic regions on the concave surfaces of Spy, where the client protein is believed to bind (26). We wondered, however, if increasing surface hydrophobicity is the only avenue for improving a chaperone's function or if nature can provide other solutions.
Here, using an improved in vivo selection, we discovered new point mutations that enhance Spy's activity. Some of these new variants act not by increasing hydrophobicity but instead by increasing the positive charges on Spy's concave surface, allowing Spy to capture its substrates more quickly. Combining these positively charged point mutations further enhanced Spy's anti-aggregation activity; however, combining the previously discovered hydrophobicity-enhancing mutations did not. By performing kinetic measurements on these hydrophobically enhanced mutants, we were able to attribute this inability to further increase Spy's anti-aggregation activity to a considerable reduction in the rates at which they bind to the substrate. Additionally, we demonstrated that rapid capture of the substrate is more important than longer substrate binding in enhancing Spy's ability to suppress rapid aggregation. Despite these mechanistic differences, we found that the performance of Spy variants in vivo is largely determined by their expression levels, suggesting that instead of relying on a single trait, nature appears to balance multiple activity determinants to optimize Spy performance.

Isolation and initial characterization of new Spy variants with enhanced chaperone activities
To identify the determinants of Spy chaperone activity, we performed directed evolution to enhance Spy activity using our previously established approach (26) but with an important modification. Instead of subjecting the spy gene to random point mutagenesis, the Spy-encoding library was prepared with a DNA shuffling-based method called "random insertionaldeletional strand exchange mutagenesis" (RAISE) (Fig. S1A) (27). This was done, in part, to overcome the codon bias of the error-prone PCR approach that was used to generate the previous spy mutant library (28,29). The selection relies on Spy's ability to stabilize the L53A I54A mutant of Im7 inserted into the b-lactamase stability biosensor (30). Improved folding of Im7 L53A I54A allows the two split fragments of b-lactamase to re-associate and restore its enzymatic activity. Selection for penicillin resistance of the cells co-expressing the library and the b-lactamase (bla)-Im7 L53A I54A stability biosensor allowed us to isolate 42 Spy variants (Fig. S1B). Although 27 (64%) of these Spy variants contained known activity-enhancing mutations (26) : (Q100L, H96L, Q49L, L32P, and Q25R), we obtained six new gain-of-function mutations, including four isolated as single substitutions (F29V, L32Q, G48A, and Q52R) and two as high-frequency mutations (H65Y and H96R) that were identified in combination with other amino acid substitutions ( Fig. 1A and Fig. S2).
To verify that the high-frequency mutations (H65Y and H96R) are able to exert a positive effect on their own and that the six new mutations (F29V, L32Q, G48A, Q52R, H65Y, and H96R) act to enhance the activity of Spy as a chaperone, we reintroduced the individual mutations into the spy gene on pCDFTrc-Spy and assayed for penicillin resistance of cells coexpressing the bla-Im7 L53A I54A biosensor (26). As expected, we found that the relative minimal inhibitory concentrations (MICs) of cells expressing the Spy variants were higher than the MIC of cells expressing WT Spy (Fig. 1B), verifying that all these mutations are on their own sufficient to enhance the antibiotic resistance of the biosensor. Next, to test whether the newly selected Spy variants are generally more active toward other substrates, we purified these Spy variant proteins and assayed for their abilities to prevent the aggregation of DTTreduced a-lactalbumin (a-LA) and heat-denatured malate dehydrogenase (MDH). All six of these variants increased Spy activity on a-LA and four showed increased activity toward MDH ( Fig. 1C and Fig. S3A), indicating that beneficial mutations selected in vivo on one substrate (Im7 L53A I54A) also function better in vitro toward at least one other substrate. Overall, activity-enhancing variants obtained in the current and previous selection were localized to three regions of Spy: the disordered N terminus, the flexible linker between helix 1 and 2, and the regions near the junction between helix 2 and 3, and helix 3 and 4 ( Fig. 1A). Most of these mutants were found in regions that, when changed, presumably minimally perturb Spy's overall structure. These mutations may affect local flexibility, which if true would be consistent with our previous observations that flexibility is important for Spy's action (26).
Our previous result revealed that tighter binding to substrates, slower substrate release rates, and increased structural flexibility all contributed to enhancing Spy's activity (26). To explore whether the new Spy variants used similar strategies to improve their function, we conducted thermodynamic and kinetic characterization of these new variants. Consistent with our previous findings, we observed a correlation between the dissociation constants (K d s) of the six new Spy variants toward the substrate Im7 H40W L53A I54A and their abilities to stabilize Im7 L53A I54A in vivo (Fig. 1D). However, unlike with our previous observations in which 70% of the mutants decreased the thermodynamic stability of Spy, only half of our new Spy variants showed dramatically decreased thermodynamic stabilities (F29V, L32Q, and H96R); the stabilities of the other three variants essentially remained unaltered (Fig. 1E). In addition, all the new Spy variants had similar secondary structures to the WT protein as measured by circular dichroism (CD) at 25°C (Fig. S3B). Therefore, the activity-increasing mechanisms for half of the new variants (G48A, Q52R, and H65Y) cannot be easily explained by a decrease in stability or a change in overall secondary structure. These results motivated us to conduct a more detailed investigation in the hope of uncovering additional factor(s) that help determine Spy's anti-aggregation activity.
Increased hydrophobicity at the concave surface contributes to the increased chaperone activities of some Spy variants Previously, we proposed that increasing surface hydrophobicity is beneficial because three of the activity-enhancing Charges regulate chaperone activity variants that we had previously isolated (Q49L, H96L, and Q100L) are located immediately adjacent to the two predominantly hydrophobic patches (P1 and P2) on the concave surface of Spy (26). The His-65 and Gly-48 residues in the newly isolated mutants are also located near these patches ( Fig. 2A), suggesting that the activity of the H65Y and G48A variants might be explained by a similar mechanism. To evaluate the extent to which these mutations change Spy's hydrophobicity, we performed 4,49-bis-1-anilinonaphthalene-8-sulfonate (bis-ANS) titration on these two variants and WT Spy. Bis-ANS is an environmentally sensitive fluorescent probe, the quantum yield of which increases upon binding to hydrophobic regions of proteins (31). In addition to the Q49L, H96L, and Q100L "super Spy" variants, our assay also intended to include a "nonsuper Spy" control with one of its hydrophobic residues located near Spy's hydrophobic patches replaced by a hydrophilic serine residue. To identify the hydrophobic residue that tolerates the serine substitution, we surveyed all the hydrophobic residues located in the two main hydrophobic patches (P1 and P2) and one hydrophobic residue near P1, which are L34S, I42S, M46S, M53S, M64S, I68S, M85S, M93S, M97S, and I103S (Fig. S4A). We found that the M53S mutation had less impact on the overall protein stability than all the other substitutions so we used it as the nonsuper Spy control (Fig. S4B). When co-expressed with the bla-Im7 L53A I54A biosensor, the M53S mutant decreased the MIC to ;70% of the antibiotic resistance of the strain expressing WT Spy (Fig. S5A). At each bis-ANS concentration tested, all variants except M53S produced fluorescent signals higher than that of WT, indicating that the surface hydrophobicity of all variants except the control M53S was increased (Fig. 2B). Using reverse titration, we determined the stoichiometry of bis-ANS binding and found that all variants except M53S bound more bis-ANS molecules than the WT ( Fig. 2B and Fig. S6), a further indication that these super Spy mutants had enhanced surface hydrophobicity. In addition, we observed a good correlation between the bis-ANS binding capacities and the anti-aggregation activities of the variants (Fig. 2C). toward penicillin in strains expressing the Spy variants and the bla-Im7 L53A I54A folding biosensor containing a destabilized Im7 variant. C, all newly obtained Spy variants showed higher activity in preventing the aggregation of DTT-reduced a-LA (a-LA agg. Prev) compared with WT (n = 3, Student's t test). **p , 0.01; ***p , 0.001. Assay buffer contained 100 mM NaCl. D, variants with higher affinity (lower K d ) toward Im7 H40W L53A I54A as measured by bio-layer interferometry exhibit higher in vivo chaperone activity toward Im7 L53A I54A. E, thermodynamic stabilities of the Spy variants determined by urea-induced denaturation at 25°C indicate that only three new variants (F29V, L32Q, and H96R) and the Q100L control are destabilized compared with WT. DG NU is the free energy difference between the unfolded and folded states. Error bars in (B), (C), and (D) (for MIC normal ), S.D. of three independent measurements whose data points are shown; error bars in (D) (for K d) and (E) are fitting errors.
For the M53S control mutant, we found that its chaperone activity also decreased in vitro, although its stability and secondary structure were not fundamentally altered compared with the WT protein ( Fig. S4B and Fig. S5, B and C). Thus, the decrease in M53S's activity might be attributed to the reduction in surface hydrophobicity. Together, these results suggest that Spy's activity can be correlated to its accessible hydrophobic surface area.
Optimizing surface charge confers Spy with higher activity than maximizing surface hydrophobicity It is unlikely that the enhanced activities of the other two variants, Q52R and H96R, are attributable to increased surface hydrophobicity because the mutant amino acid is charged. Spy has an overall positively charged concave surface (Fig. 3A). This feature allows Spy to use long-range electrostatic interactions to form the initial complex with its substrates (19). Because Q52R and H96R both increase the positive charge on the mutated sites, we suspected that they might function by enhancing these initial electrostatic interactions between Spy and its substrates. To test this hypothesis, we increased the concentration of salt in the a-LA aggregation assay. We found that Q52R and H96R remained substantially more active than WT Spy in 100 mM salt but no longer showed substantial enhancement in 150 mM salt ( Fig. 1C and Fig. 3, B and C). This result suggests that the activity enhancement of these two variants emanates from enhanced electrostatic interactions with substrates. We wondered if two other positions, which are adjacent to Spy's main hydrophobic areas and also to sites of super Spy substitutions, namely, Gln-49 and Gln-100, could also tolerate positively charged mutations. To address this question, we constructed the Q49R and Q100R variants and assayed their antiaggregation activity toward a-LA, which is negatively charged overall ( Fig. 3A) but presumably presents both hydrophobic and charged surfaces upon unfolding. Consistent with the idea that electrostatic interactions are important for Spy-substrate interaction, we found that the chaperone activities of these variants toward a-LA were also improved relative to WT Spy ( Fig. 3D). Additionally, whereas Q49R was almost as active as the previously isolated super Spy variant Q49L, the activities of the Q100R and H96R variants were even higher than their corresponding previously isolated leucine-substituted super Spy variants (Fig. 3D).
We then wondered if combinations of activity-enhancing mutations, particularly those with the same type of substitution, could result in a chaperone with further enhanced activity. Indeed, the combination of H96R and Q100R resulted in a roughly 2-fold increase in Spy activity compared with the single variants, indicating that the effects of these two arginine substitutions are additive (Fig. 3D). However, the incorporation of a third arginine substitution (Q49R) to form the triple variant Spy Q49R H96R Q100R had negligible effects. We next tested whether activity enhancement was specific to the positive charge on the arginine residues by mutating Gln-49, His-96, and Gln-100 to glutamic acid, whose side chain is similar to that of glutamine except that it bears a negative charge. We found that the chaperone activities of these more negatively charged variants were substantially decreased compared with those of WT (Fig. 3E), even though the in vitro stabilities of Q49E and H96E were almost unchanged from that of the WT (Fig. S7, A and B). Q100E, despite having a 9°C decrease in the temperature at the midpoint of the unfolding curve, is 94% folded at 25°C, the temperature at which the activity assay was conducted (Fig. S7C). This result indicates that the observed activity enhancement for the arginine-substituted variants depends on the positive nature of the charge.
In contrast, combining leucine substitutions produced very different results: the combination of H96L and Q100L, with or without the additional Q49L, decreased activities to WT Spy levels (Fig. 3D). The effect of these combinations is somewhat surprising given the common assumption that chaperones are primarily attracted to their substrates via hydrophobic interactions, an assumption that might lead one to naively think that increasing surface hydrophobicity would always be advantageous for a chaperone. The nonadditive effect of the leucine substitutions instead suggests that surface hydrophobicity has Charges regulate chaperone activity been optimized rather than maximized by evolution, as is common for many traits that are subject to evolutionary tradeoffs.

Decreased substrate association rates hinder chaperone activities for hydrophobically enhanced variants
Previous work revealed that electrostatic interactions dominate the early stage of Spy-substrate complex formation, whereas hydrophobic interactions dominate the later stage to further stabilize the complex (19). Our observations about the effects of various combinations of mutations on activity could be interpreted to indicate that optimizing chaperone-client interactions at early stages of the interaction is more important for aggregation prevention than optimizing interactions that occur at later stages. To test whether the behavior of our newly isolated Spy variants is consistent with these findings, we determined the thermodynamic and kinetic parameters of their interaction with reduced, carboxymethylated a-LA using a biolayer interferometry-based association/dissociation assay at the physiological salt concentration (150 mM NaCl). We found that the arginine-substituted and leucine-substituted Spy variants all bound to a-LA more tightly than did WT Spy, and the combination of either type of mutation resulted in even higher affinities. The leucine-substituted variants had ;2-fold higher affinities (K d s) than the corresponding arginine-substituted variants, with H96L being the only exception in that it bound a-LA as tightly as did H96R (Fig. 4A). . The concave surface of Spy is dominated by positive charges (blue), whereas the overall surface of native a-LA is negatively charged (red). The electrostatic surface potential was calculated in PyMol with the Adaptive Poisson-Boltzmann Solver tools 2.1 plugin. A color scale for the charge distribution from -8 to 8 was used. B and C, anti-aggregation activities of Spy variants Q52R and H96R in the presence of 100 mM (B) or 150 mM NaCl (C). Representative curves from duplicate measurements are shown (B-C). D, both leucine-substituted (yellow) and arginine-substituted (blue) Spy variants increase their abilities to prevent a-LA aggregation (a-LA agg. Prev). However, the charge-enhanced variants are in general more active than their corresponding hydrophobically enhanced counterparts (n = 3, Student's t test). ***p , 0.001; ns, not significant. E, introducing negative charges on the concave surface of Spy results in decreased chaperone activity. Chaperone activities of glutamate-substituted Spy variants at residues Q49, H96, and Q100 were assayed for their abilities to prevent aggregation of DTT-reduced a-LA (a-LA agg. Prev). Glutamate-substituted Spy variants showed decreased chaperone activities compared with Spy WT (n = 3, Student's t test). ***p , 0.001. Error bars, S.D. of three independent measurements whose data points are shown.
We observed a similar trend with the substrate release rate (i.e. the dissociation rate, k off ); all Spy variants had decreased k off values relative to WT, with the leucine-substituted variants exhibiting lower values than the corresponding arginine-substituted ones (Fig. 4B). This trend was also observed when assaying two other substrates, Im7 H40W L53A I54A and a-casein (Fig. S8). Previously, based on the observation that the substrate release rate was less affected by ionic strength than the binding process, it was concluded that substrate release is controlled mainly by hydrophobic forces, whereas electrostatic interactions contribute more to the initial stages of complex formation than to substrate release (19). Our new observation that the leucine-substituted variants hold the substrate longer than their argininesubstituted counterparts do is consistent with this conclusion. Furthermore, when increasing the temperature of kinetic assays from 25°C to 32°C, we found an increase in the k off rates for all the variants tested (Fig. S9). As temperature increases, electrostatic interactions become strengthened due to reduced solvent screening (32,33), which then will tend to decrease the k off values if the substrate release process is also primarily governed by electrostatic interactions. Our observed k off increase instead argues the involvement of hydrophobic forces, the strength of which peaks at around 16°C and diminishes as the temperature is raised or lowered relative to this temperature (34). Thus, the increased k off rates from 25°C to 32°C is consistent with the conclusion that the client release process is dominated by hydrophobic interactions.
To effectively prevent aggregation, Spy must recognize the aggregation-prone species quickly before they can form aggregation nuclei. We found that the three arginine-substituted variants (Q49R, H96R, and Q100R) had slightly increased association rates (k on ), but when these mutations were combined, the k on rates decreased to that of the WT protein (for H96R Q100R (2R)) or even slightly below WT (for Q49R H96R Q100R (3R)) (Fig. 4C). In contrast, only one leucine-substituted variant (Q49L) acted to increase k on , whereas the H96L Q100L (2L) or Q49L H96L Q100L (3L) substitutions decreased k on by over 2fold, which might account for the reduced anti-aggregation activities of these variants. These results also suggest that the k on value must be kept above a specific threshold to allow for effective substrate capture. The 2L and 3L variants with rate constants of ;10 5 M 21 s 21 are at the lower end of the diffusion-controlled regime (the "basal" value for diffusion-controlled association rate constants ranges from 10 5 to 10 6 M 21 s 21 ) (15). To go beyond the basal diffusion-controlled rate constant, intermolecular forces must be active during the diffusion process leading to the formation of a transient complex (15). Long-range electrostatic interactions are thought to be the main type of interaction capable of speeding up the formation of the initial complex and may therefore help explain why the charge-enhanced single arginine-substituted variants have accelerated substrate association rates and enhanced activities. In contrast, short-range forces, including hydrophobic and van der Waals interactions, are relatively ineffective in affecting diffusion-limited interactions and contribute much less to rate enhancements compared with long-range electrostatic interactions (15). Our kinetic results thus provide a plausible explanation as to why the more hydrophobic 2L and 3L variants showed no enhanced anti-aggregation activity toward a-LA.
To further examine the involvement of electrostatic interactions in the binding reaction of these Spy variants, we repeated the kinetic measurements at different ionic strengths. We saw a decrease in the k on rates of Spy WT and the Q49R, H96R, and Q100R mutants as the salt concentration increased from 150 mM to 250 mM (Fig. S10A). The Q49L, H96L, 2R, and 3R mutants were less sensitive to changes in salt concentration from 150 mM to 200 mM, but all showed a decrease in k on rates at 250 mM salt. Interestingly, the Q100L, 2L, and 3L mutants showed a different trend in k on rates: their k on rates increased with increasing ionic strength (Fig. S10A). This trend is not surprising given that higher ionic strength is known to contribute to stronger hydrophobic interactions (35,36). Thus, our results suggest that higher ionic strength increases the k on rates of the leucine-substituted Spy variants, possibly by enhancing their hydrophobic interactions with a-LA. On the other hand, consistent with previous reports that electrostatic interactions also contribute to the maintenance of Spy-substrate complexes (19,20), we found that the k off rates of all Spy variants increased with increasing salt concentrations, which resulted in their overall lower affinities at higher ionic strengths (Fig. S10, B and C). These kinetic measurements confirmed the involvement of both electrostatic and hydrophobic interactions in Spy's substrate binding and release and also revealed the interplay of different molecular forces as ionic strength changes. , and association rate constants k on (C) of leucine-substituted (yellow) and arginine-substituted (blue) Spy variants toward reduced, carboxymethylated a-LA (unfolded) were determined by bio-layer interferometry. Generally, leucine-substituted variants bind substrate tighter (A), recognize substrate slower (B), and release substrate slower (C) than their corresponding arginine-substituted counterparts (n = 2, Student's t test). *p , 0.05; **p , 0.01; ns, not significant. Error bars, S.D. of two independent measurements whose data points are shown.

Charges regulate chaperone activity
Enabling rapid substrate association is likely to be more effective than attempting to decrease substrate dissociation rates in preventing aggregation To determine whether the apparently important role of fast association in preventing aggregation can be generalized to other substrate proteins, we measured the anti-aggregation activity of Spy toward another substrate, MDH.
We used two methods to denature MDH: 8 M urea and heat. We did this to test different types of unfolding intermediates that aggregate at different rates. We reasoned that a rapidly aggregating MDH unfolding intermediate would be more dependent on chaperone capture to delay the formation of the aggregation nuclei than a slowly aggregating intermediate. If fast association is indeed a key factor in preventing aggregation, given the marked difference in the k on rates, we would expect to observe a more pronounced difference in the activities of arginine-or leucine-substituted Spy variants in preventing aggregation of fast-aggregating MDH unfolding intermediates compared with their activities toward slow-aggregating MDH intermediates. We found that MDH that had been chemically denatured using 8 M urea aggregated very rapidly at both 32°C and 25°C because the T 50 (times to reach half of the maximum aggregation signal) of MDH alone at these temperatures were short: 8 s at 32°C and 15 s at 25°C (Fig. 5, A and B). However, when the incubation temperature was further decreased to 22°C, the aggregation of urea-denatured MDH was significantly slowed; the T 50 was extended to 56 s. This indicated to us that reduction in temperature could change the aggregating kinetics of urea-denatured MDH and could change the rapidaggregating nature of urea-denatured MDH at high temperatures to a more slowly aggregating behavior (Fig. 5C). Heat-denatured MDH in the absence of Spy aggregated at 43°C (T 50 = 201 s) even slower than did urea-denatured MDH at 22°C (T 50 = 56 s) (Fig. 5D). Using these methods, we could observe two rapidly aggregating species (urea-denatured MDH at 32°C and 25°C) and two slowly aggregating models (urea-denatured MDH at 22°C and heat-denatured MDH at 43°C).
We found that for rapidly aggregating MDH, all leucine-substituted Spy variants, including the 2L and 3L variants, displayed equal or slightly worse activities than did WT Spy, whereas all arginine-substituted Spy variants, including the 2R and 3R variants, performed better than did WT Spy (Fig. 5, A and B and Fig. S11, A and B). This result is consistent with the observations made with a-LA aggregation, which is also rapid (Fig. 3D). For slowly aggregating MDH, almost all leucine-substituted variants were more active than WT, and the Q49L, Q100L, and 3L variants were nearly as active as their corresponding arginine-substituted variants (Fig. 5, C and D and Fig.  S11, C and D). Taking all these results into account, we conclude that for relatively slow-aggregating MDH, the formation of a more stable Spy-MDH complex still allows for the productive inhibition of aggregation. However, if one is seeking to prevent problematically fast aggregation, then promoting rapid substrate association rates may be more impactful than attempting to stabilize the complex.

Activities of the charge-enhanced Spy variants are limited by decreased protein levels in vivo
Our in vitro data showed that the charge-enhanced variants at Gln-49, His-96, and Gln-100 were more effective than the WT or even the hydrophobicity-enhanced variants at inhibiting the aggregation of fast-aggregating substrates. This result led us to wonder why the Q49R and Q100R variants or their combinations did not emerge during our selections. To answer this question, to systematically analyze the effects of chaperone action of these two residues that are clearly vital for chaperone action, and to better tie together our in vivo and in vitro results, we mutated Gln-100 to each of the other 19 amino acids and assayed the chaperone activities of these variants in vivo using the bla-Im7 L53A I54A stability biosensor. Unexpectedly given our current interpretation based on in vitro assays, but consistent with our previous interpretations, we found that the antibiotic resistance of the strains expressing these Spy variants was reasonably well-correlated with the hydropathicity of the substituting residue (Fig. 6A). There was also a correlation, although weaker, between the hydropathicity of the substituting residue and the amount of bla-Im7 L53A I54A (Fig. 6B). These results indicate that there is a relationship between the hydropathicity of the substituting residue and the stabilization it confers. We then examined the steady-state levels of each Spy variant and found that many of the variants with hydrophilic substitutions had overall lower levels compared with the variants substituted with hydrophobic amino acids (Fig. 6, C  and D). In particular, the Q100R mutant, which had the lowest MIC value, was also expressed in very low amounts. Therefore, the in vivo performance of charge-enhanced Spy variants was probably impaired by their steady-state levels. This indicates that these variants may have in vivo folding problems and thus be more prone to proteolysis in vivo.

Discussion
A single protein can undergo an almost bewildering variety of conformational changes during the folding process. This, combined with the fact that most chaperones assist in the folding of a broad variety of substrates, makes understanding the nature of chaperone-substrate interaction difficult from both a theoretical and practical standpoint, but understanding this interaction is vital because it represents one of the core problems in chaperone biology. It is understandable that scientists have attempted to simplify this problem by using the following, quite attractive, logic: folding and unfolding proteins display hydrophobic regions, the hydrophobic effect drives such regions together, which should cause proteins to aggregate, and chaperones prevent aggregation; ergo, chaperones recognize hydrophobic regions. Early attempts at solving this problem by mapping chaperone specificity using peptide display approaches provided some evidence for this model, but close examination of even this early data revealed a more complex picture: the recognition site for Hsp70, though it does contain hydrophobic residues, also contains positive charge (37). The commonly held view that molecular chaperones bind their substrates primarily via hydrophobic interactions was thus too simple a picture even in the early days of chaperone biology (5,38). Another somewhat neglected point concerning protein aggregation is that in addition to being triggered by exposure of sticky hydrophobic surface patches, self-association of native-like structures can also result in aggregation (39). Both theoretical and experimental evidences Figure 6. Lower steady-state levels of Spy variants harboring hydrophilic mutations at Gln-100 may explain their decreased in vivo activities. A, relative MICs of Spy variants compared with WT Spy were plotted against the amino acids' hydropathy index. B, amounts of bla-Im7 L53A I54A in strains expressing different Spy variants (normalized to the bla-Im7 L53A I54A amount in the strain expressing WT Spy) were plotted against the amino acids' hydropathy index. C, steady-state levels of each Spy variant and the bla-Im7 L53A I54A fusion protein co-expressed in the same strain were analyzed by Western blotting with anti-Spy and anti-TEM-1 b-lactamase antibodies, respectively. Endogenous DnaK was used as a loading control. D, amounts of different Spy variants (normalized to the amount of WT Spy) were plotted against the amino acids' hydropathy index. A, B, and D, hydropathy index is according to Zviling (57); yellow dots, amino acids with a hydropathy index above 0; gray dots, amino acids with a hydropathy index below 0. Error bars in (A), S.D. of three independent measurements; error bars in (D), S.D. of two technical repeats.
Charges regulate chaperone activity also suggest that electrostatic interactions can drive aggregation (40,41). Though electrostatic repulsion between protein molecules disfavors self-association, asymmetric distribution of surface charge can allow for alignment of patches with complementary charges, which increases a protein's aggregation propensity (39). Recapitulating these overlooked facts will help us to achieve our long-term goal, i.e. to understand protein folding in sufficient detail so that it can be manipulated in vivo.
Our mechanistic characterization of activity-enhanced Spy variants enabled us to delineate the particular contributions that various intermolecular forces play in contributing to Spy's anti-aggregation activity (Fig. 7). This has allowed us to gain a much more nuanced view of chaperone-substrate interaction than the current view, where hydrophobic interactions are thought to dominate. We found that Spy variants possessing discrete surface properties have very different effects on their ability to prevent aggregation; specifically, Spy variants with more positively charged concave surfaces are able to interact with aggregation-prone species quickly and retain these substrates on their surfaces for a relatively longer period of time Figure 7. Mechanism of the charge-enhanced or hydrophobicity-enhanced Spy variants in preventing client aggregation. For fast-aggregating clients (upper panel), both WT Spy and the charge-enhanced variants are capable of rapidly binding unfolded clients, but the latter also releases the substrate more slowly, resulting in very few aggregates. In contrast, although the hydrophobicity-enhanced variants are capable of forming tighter complexes with clients, their slow association rates limit their ability to prevent rapid aggregation. For slow-aggregating clients (lower panel), the very slow dissociation rates of the hydrophobicity-enhanced variants make up for their shortage in recognizing clients, resulting in almost equally holding activities as the charge-enhanced variants. Boxed regions indicate the processes more impactful in determining the aggregation of different types of clients.
because of slower dissociation rates. This results in less aggregation-dependent precipitation. These results contradict the widely held assumption that chaperones and their clients generally interact through hydrophobic interactions. Hydrophobically enhanced variants were also able to reduce the formation of aggregates because of their property of only allowing for the very slow release of their substrates. Hydrophobic substitutions may be effective as long as their rate of binding is above some threshold value that allows for effective substrate capture and these substitutions do not substantially change the distribution of the original electrostatic potential on the concave surface of Spy.
There are additional indications that chaperone function is not solely dependent on hydrophobic interactions. For instance, it is known that GroEL preferentially recognizes amphipathic regions within its substrate proteins (42). The related eukaryotic chaperonin TRiC was also found to rely on a combinatorial substrate-recognition mechanism that involves a distinct, evolutionarily conserved pattern of polar and hydrophobic residues. Similar to Spy, charged and polar residues on TRiC contribute to enhancing the association rates, and hydrophobic residues are responsible for decreasing the dissociation rates (43). Likewise, the binding sites of the intrinsically disordered protein Tau on the chaperone Hsp90 are composed of scattered patches that have both hydrophobic and negative charge characteristics, allowing the chaperone to make a large number of low-affinity contacts with various substrates (44). The electrostatic interactions between Skp and its outer membrane protein substrates were found to be important for substrate delivery and folding into the liposomes in vitro (45). It has been shown that trigger factor recognizes the peptide segments enriched in aromatic and basic residues within its substrate proteins (46). In addition, SecB also appears to recognize both hydrophobic and positively charged regions, allowing it to deliver substrates with different hydrophobicity and electrostatic potential to SecA for translocation (47,48). Thus, many molecular chaperones appear to have a bipartite nature associated with their substrate binding sites. This may enable their adaptation to the different structural properties present in their diverse substrates and thus facilitate promiscuous substrate recognition. Our previous and current work suggests that Spy also possesses this bipartite nature and optimizing either hydrophobic or electrostatic properties can convert it into a more efficient chaperone toward different in vitro substrates.
Engineering molecular chaperones, especially improving their anti-aggregation or disaggregation activities, is important because it may provide crucial insights that could be helpful in treating several age-related neurodegenerative diseases whose hallmarks include toxic protein misfolding and aggregation (21,49). One successful example of engineering to modify aggregation activity is the optimization of the Hsp104 chaperone, which disaggregates and eradicates preformed fibrils of many amyloid proteins (24,(50)(51)(52). Improved variants identified in an in vivo screen against a-synuclein and other amyloid proteins in yeast revealed a hot spot: Ala-503; substitution with any residue except proline at this site made Hsp104 more active (24). However, directed evolution of chaperones to improve their general activity is not entirely straightforward. For exam-ple, the GroEL chaperonin appears to possess a narrow mutational window; mutations that made it better at folding GFP decreased its ability to fold its native substrates (23). Our work with Spy also revealed that mutants obtained in an in vivo screen against one type of substrate might not work efficiently toward another type of substrate. We saw this with a-LA and chemically denatured MDH whose in vivo aggregation is difficult to model because these substrates aggregate very rapidly in vitro, so rapidly that chaperones are unable to act on them. Multiple factors can contribute to chaperone activity. These factors may already have been optimized for each chaperone by evolution, and improving on this evolutionary balance may be difficult. However, knowledge is power, and more detailed knowledge about chaperone-substrate interactions will undoubtedly be helpful in designing chaperones.
In this study, we have presented a genetic means of modulating chaperone activity. The observation that chaperones can act by modifying the surface electrostatic potential and hydrophobic properties suggests that even evolutionarily optimized chaperones can be further improved to more effectively coordinate substrate binding and release. Although our findings here are based on the optimization of Spy, the binary nature of many chaperone surfaces supports similar optimization strategies for these chaperones, including ATP-dependent ones such as Hsp90. To this end, we believe that our approach, along with insights gained from other successful engineering endeavors, such as improving overall structural flexibility (26), enhancing ATPase activity (23), and combining functional domains derived from different molecular chaperones (53), will provide a powerful means of enhancing chaperone activity. We expect that as more key mechanistic and functional insights into the details of chaperone-substrate interaction emerge, the field of rational chaperone design will advance and will eventually be able to support the development of therapeutic chaperones targeting neurodegenerative diseases.

RAISE library construction and selection
Two independent RAISE libraries, each with ;5 3 10 4 clones, were constructed as previously described with the following modifications (27). 10 mg of the PCR product of the spy gene was digested by 0.1 units of DNaseI in a 20-ml reaction at 37°C for either 25 s or 18 s to generate a random fragment mixture, with a range of 75 to 150 bp or 100 to 300 bp, respectively. Subsequently, 2 pmol of the digested fragments were 39-tailed with random nucleotides by 2.5 units of terminal deoxynucleotidyl transferase enzyme in a 20-ml reaction at 37°C for 1 h. Then, 1 pmol of 39-tailed fragments were assembled into fulllength spy via self-priming PCR mediated by 2 units of Deep VentR (exo-) polymerase in a 20-ml reaction. Finally, the mutated spy libraries were amplified with primers complementary to the end of the spy gene, ligated in the pCDF vector, and transformed into E. coli strain Trans1T1.
Library selection was conducted by transforming the pCDFTrc-Spy library into a spy and ampC null strain HW41, which also contained the pBR322 bla::GSlinker Im7 L53A I54A vector. Subsequently, the transformants were plated on agar plates supplemented with 3 mg/ml penicillin and 0.1 mM isopropyl b-D-thiogalactopyranoside, and resistant colonies were picked for sequencing after overnight incubation at 37°C.

Mutagenesis and protein purification
Mutations to spy were generated via the QuikChange sitedirected mutagenesis method (Agilent, catalog no. 200518). All Spy variants were purified according to Quan et al. (26) with the following modifications. After initial purification with immobilized nickel-nitrilotriacetic acid affinity chromatography, His 6 -SUMO-Spy was eluted with 300 mM imidazole and digested by ubiquitin-like-specific protease 1 in solution while dialyzing against buffer A (20 mM HEPES and 0.5 mM EDTA, pH 7.5) before undergoing cation exchange chromatography. Protein concentrations were determined by absorbance at 280 nm, and all proteins were aliquoted into 1.5-ml Eppendorf tubes and stored at -80°C prior to use.

Chaperone activity assay in vivo
Spot titer experiments were performed to quantify the relative in vivo chaperone activities of Spy variants. Generally, cells expressing each Spy variant and the bla-Im7 L53A I54A sensor were grown to mid-log phase, serially diluted, and plated onto agar plates containing increasing concentrations of penicillin V (0-3500 mg/ml) and 0.1 mM isopropyl b-D-thiogalactopyranoside. Then, the relative MIC (MIC mutant /MIC WT ) of Spy variants toward penicillin V was calculated as described previously (26,30).

a-LA aggregation assay
Bovine a-LA, type III (Sigma-Aldrich) was used as a model substrate to quantify the specific in vitro holding activities of Spy variants as described previously (26). Aggregation of a-LA (50 mM) was initiated by the addition of 20 mM DTT in a buffer containing 50 mM sodium phosphate, 100 mM sodium chloride, and 5 mM EDTA at pH 6.9. The light scattering of a-LA in the absence or presence of Spy was monitored at 360 nm using a Synergy HTX Multi-Mode Microplate Reader with a 5-min detection period for 10 h at room temperature. For relative activity quantification, a standard curve was produced by plotting the slopes calculated from light scattering data between 120 and 180 min against 22 different Spy WT:a-LA ratios. Then, slopes of light scattering determined at the same time range for different Spy variants:a-LA ratios (0.075, 0.1, 0.15, and 0.2) were calculated to quantify the relative activity according to the standard curve. The final average relative activities of Spy variants were generated from eight relative activity values (four concentrations in duplicate).

MDH aggregation assay
For thermal-denatured MDH aggregation, 0.5 mM MDH was added into 43°C pre-warmed 40 mM HEPES-KOH, pH 7.5 buffer in the presence or absence of 0.25 mM Spy or Spy variants under stirring. MDH aggregation was monitored by light scattering at 360 nm using a Lumina fluorescence spectrometer supplied with a Peltier temperature controller (Thermo Fisher Scientific).
For chemically denatured MDH aggregation, MDH (denatured in 8 M urea for 1 h at 25°C) was diluted into 40 mM HEPES-KOH, pH 7.5 buffer with a final concentration of 0.5 mM in the presence or absence of 0.05 mM Spy or Spy variants under stirring. Light scattering at 360 nm was recorded with a Lumina fluorescence spectrometer at 22, 25, and 32°C.

Bio-layer interferometry
Kinetics of Spy variants toward biotinylated and carboxymethylated a-LA and biotinylated Im7 H40W L53A I54A were measured with an Octet RED96 system (ForteBio, Menlo Park, CA). Carboxymethylated a-LA was prepared to obtain a stable, totally unfolded substrate as previously described (54). The carboxymethylated a-LA and Im7 H40W L53A I54A were biotinylated with EZ-link NHS-biotin (Thermo Fisher Scientific) with a 1:1 molar ratio of biotin to protein. 20 mg/ml biotinylated and carboxymethylated a-LA or biotinylated Im7 H40W L53A I54A was immobilized to the streptavidin sensor, whereas different concentrations (0.25, 0.5, 1, 1.5, and 2 mM) of Spy WT and variants in 40 mM HEPES and 150 mM NaCl, pH 7.5 were added to a 96-well plate (for carboxymethylated a-LA, the assay buffer contains 150-250 mM NaCl). One kinetic cycle at a defined Spy concentration included the following steps: ligand binding for 25 min, wash for 5 min, equilibrium for 1 min, association for 10 min, and dissociation for 10 min. Curve fitting and calculation of association rate constants (k on ), dissociation rate constants (k off ), and dissociation constants (K d ) were performed as described previously (26).

Thermostability measurement
For urea-induced unfolding of WT Spy and variants, ellipticity at 226 nm was monitored for each protein over increasing concentrations of urea (0-9 M) with a Chirascan circular dichroism (CD) spectrometer (Applied Photophysics Corp., Surrey, UK). All proteins were brought to a concentration of 0.1 mg/ml and incubated in denaturation buffer (30 mM diglycine, pH 8.0, and 0.1 M NaCl) supplemented with different concentrations of urea for 30 min before measurement. The Gibbs free energy of unfolding of each protein was determined via the linear extrapolation method; entire unfolding curves were fitted by the method of Santoro and Bolen using Origin 8.0 software (OriginLab) as described previously (55).
The thermal unfolding of WT Spy and variants was followed by monitoring the ellipticity at 222 nm over the temperature range of 20-75°C, with an interval of 1°C and a heating rate of 1.5°C/min. The thermal unfolding curves were fitted with the Gibbs-Helmholtz equation using Origin 8.0 software (Origin-Lab) and the transition temperature at the midpoint of the unfolding curve for each protein was calculated as previously described (56).

Bis-ANS fluorescence
Bis-ANS was dissolved in HN buffer (40 mM HEPES-NaOH and 150 mM NaCl, pH 7.5) to obtain a 7 mM stock. To estimate the number of bis-ANS bound to each Spy molecule, 0.5 mM Spy variants in HN buffer were titrated with bis-ANS over a range of concentrations (0-6 mM). The emission of bis-ANS at 493.4 nm was recorded with an excitation wavelength of 395 nm. The observed fluorescence intensity was plotted as a function of bis-ANS concentration, and the curves were fitted with the hyperbolic function shown in Equation 1, where F obs is the recorded bis-ANS fluorescence signal, F max is the maximal fluorescence reached upon saturation of the complex, [A] is the concentration of bis-ANS, and K is the apparent binding constant of bis-ANS. To obtain the stoichiometry of ANS binding, a reverse titration was performed using a fixed bis-ANS concentration of 1 mM and increasing concentrations of Spy or Spy variants (0-3 mM) (Fig. S6). The reverse titration curves were fitted with Equation 2, where f obs is the recorded bis-ANS fluorescence signal in the reverse titration, f max is the maximal ANS fluorescence when saturated with Spy, [S] is the concentration of Spy, and K9 is the apparent binding constant of Spy. Then, the stoichiometry of bis-ANS was calculated according to Equation 3, where n is the stoichiometry of bis-ANS bound to one Spy and [S 0 ] is the Spy concentration used in the first titration experiment (0.5 mM).

Statistical analysis
Appropriate two-tailed unpaired student's t test was employed to compare two independent groups. p values less than 0.05 were considered significant. Statistical analysis was conducted with GraphPad Prism software (version 5.0).

Data availability
All data are contained with this article and in the supporting information.