Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Exploring the Evolutionary History of Kinetic Stability in the α-Lytic Protease Family

In addition to encoding the tertiary fold and stability, the primary sequence of a protein encodes the folding trajectory and kinetic barriers that determine the speed of folding. How these kinetic barriers are encoded is not well understood. Here, we use evolutionary sequence variation in the α-lytic protease (αLP) protein family to probe the relationship between sequence and energy landscape. αLP has an unusual energy landscape: the native state of αLP is not the most thermodynamically favored conformation and, instead, remains folded due to a large kinetic barrier preventing unfolding. To fold, αLP utilizes an N-terminal pro region similar in size to the protease itself that functions as a folding catalyst. Once folded, the pro region is removed, and the native state does not unfold on a biologically relevant time scale. Without the pro region, αLP folds on the order of millennia. A phylogenetic search uncovers αLP homologs with a wide range of pro region sizes, including some with no pro region at all. In the resulting phylogenetic tree, these homologs cluster by pro region size. By studying homologs naturally lacking a pro region, we demonstrate they can be thermodynamically stable, fold much faster than αLP, yet retain the same fold as αLP. Key amino acids thought to contribute to αLP's extreme kinetic stability are lost in these homologs, supporting their role in kinetic stability. This study highlights how the entire energy landscape plays an important role in determining the evolutionary pressures on the protein sequence.

Ssy5 is a signaling serine protease that exhibits atypical biogenesis and marked S1 specificity

Ssy5 is a signaling endoprotease that plays a key role in regulating central metabolism, cellular aging, and morphological transitions important for growth and survival of yeast (Saccharomyces cerevisiae) cells. In response to extracellular amino acids, Ssy5 proteolytically activates the transcription factors Stp1 and Stp2, leading to enhanced Ssy1-Ptr3-Ssy5 (SPS) sensor-regulated gene expression. Ssy5 comprises a catalytic (Cat) domain and an extensive regulatory prodomain. Ssy5 is refractory to both broad-spectrum and serine protease-specific inhibitors, confounding its classification as a protease, and no information about Ssy5's cleavage-site preferences and its mechanism of substrate selection is available. Here, using mutational and inhibition experiments, we investigated the biogenesis and catalytic properties of Ssy5 and conclusively show that it is a serine protease. Atypical for the majority of serine proteases, Ssy5's prodomain was obligatorily required in cis during biogenesis for the maturation of the proteolytic activity of the Cat domain. Autolysis and Stp1 and Stp2 cleavage occurred between a cysteine (at the P1 site) and a serine or alanine (at the P'1 site) and required residues with short side chains at the P1 site. Substitutions in the Cat domain affecting substrate specificity revealed that residues Phe-634, His-661, and Gly-671 in the S1-binding pocket of this domain are important for Ssy5 catalytic function. This study confirms that the signaling protease Ssy5 is a serine protease and provides a detailed understanding of the biogenesis and intrinsic properties of this key enzyme in yeast.

Protein folding at extreme temperatures: Current issues

The range of temperatures compatible with life is currently estimated from -25°C, as exemplified by metabolically active bacteria between sea ice crystals, and up to 122°C in hydrothermal vents as exemplified by the archaeon Methanopyrus kandleri. In the context of protein folding, as soon as a polypeptide emerges from the ribosome, it is exposed to the effects of environmental temperatures. Recent investigations have shown that the rate of protein folding is not adapted to extreme temperatures and should be very fast at high temperature and low in cold environments. This lack of adaptation is driven by kinetic constraints on protein stability. To counteract the deleterious effects of fast protein folding in hyperthermophiles, chaperones such as the Trigger Factor hold and slow down the rate of folding intermediates. Prolyl isomerization, a rate-limiting step in the folding of many proteins, is strongly temperature-dependent and impairs folding of psychrophilic proteins in the cold. This is compensated by reduction of the proline content in cold-adapted proteins, by an increased number of prolyl isomerases encoded in the genome of psychrophilic microorganisms and by overexpression of prolyl isomerases under low temperature cultivation. After folding, the native state is reached and although extremophilic proteins share the same fold, dramatic differences in stability have been recorded by differential scanning calorimetry.

Evolutionary trend toward kinetic stability in the folding trajectory of RNases H

Proper folding of proteins is critical to producing the biological machinery essential for cellular function. The rates and energetics of a protein's folding process, which is described by its energy landscape, are encoded in the amino acid sequence. Over the course of evolution, this landscape must be maintained such that the protein folds and remains folded over a biologically relevant time scale. How exactly a protein's energy landscape is maintained or altered throughout evolution is unclear. To study how a protein's energy landscape changed over time, we characterized the folding trajectories of ancestral proteins of the ribonuclease H (RNase H) family using ancestral sequence reconstruction to access the evolutionary history between RNases H from mesophilic and thermophilic bacteria. We found that despite large sequence divergence, the overall folding pathway is conserved over billions of years of evolution. There are robust trends in the rates of protein folding and unfolding; both modern RNases H evolved to be more kinetically stable than their most recent common ancestor. Finally, our study demonstrates how a partially folded intermediate provides a readily adaptable folding landscape by allowing the independent tuning of kinetics and thermodynamics.

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

Propeptides as modulators of functional activity of proteases

Most proteases are synthesized in the cell as precursor-containing propeptides. These structural elements can determine the folding of the cognate protein, function as an inhibitor/activator peptide, mediate enzyme sorting, and mediate the protease interaction with other molecules and supramolecular structures. The data presented in this review demonstrate modulatory activity of propeptides irrespective of the specific mechanism of action. Changes in propeptide structure, sometimes minor, can crucially alter protein function in the living organism. Modulatory activity coupled with high variation allows us to consider propeptides as specific evolutionary modules that can transform biological properties of proteases without significant changes in the highly conserved catalytic domains. As the considered properties of propeptides are not unique to proteases, propeptide-mediated evolution seems to be a universal biological mechanism.

Interactions with the bifunctional interface of the transcriptional coactivator DCoH1 are kinetically regulated

Pterin-4a-carbinolamine dehydratase (PCD) is a highly conserved enzyme that evolved a second, unrelated function in mammals, as a transcriptional coactivator. As a coactivator, PCD is known as DCoH or dimerization cofactor of the transcription factor HNF-1. These two activities are associated with a change in oligomeric state: from two dimers interacting as an enzyme in the cytoplasm to a dimer interacting with a dimer of HNF-1 in the nucleus. The same interface of DCoH forms both complexes. To determine how DCoH partitions between its two functions, we studied the folding and stability of the DCoH homotetramer. We show that the DCoH1 homotetramer is kinetically trapped, meaning once it forms it will not dissociate to interact with HNF-1. In contrast, DCoH2, a paralog of DCoH1, unfolds within hours. A simple mutation in the interface of DCoH2 from Ser-51 to Thr, as found in DCoH1, increases the kinetic stability by 9 orders of magnitude, to τ(½) ∼ 2 million years. This suggests that the DCoH1·HNF-1 complex must co-fold to interact. We conclude that simple mutations can dramatically affect the dissociation kinetics of a complex. Residue 51 represents a "kinetic hot spot" instead of a "thermodynamic hot spot." Kinetic regulation allows PCD to adopt two distinct functions. Mutations in DCoH1 associated with diabetes affect both functions of DCoH1, perhaps by disrupting the balance between the two DCoH complexes.

Increasing the reaction rate of hydroxynitrile lyase from Hevea brasiliensis toward mandelonitrile by copying active site residues from an esterase that accepts aromatic esters

The natural substrate of hydroxynitrile lyase from rubber tree (HbHNL, Hevea brasiliensis) is acetone cyanohydrin, but synthetic applications usually involve aromatic cyanohydrins such as mandelonitrile. To increase the activity of HbHNL toward this unnatural substrate, we replaced active site residues in HbHNL with the corresponding ones from esterase SABP2 (salicylic acid binding protein 2). Although this enzyme catalyzes a different reaction (hydrolysis of esters), its natural substrate (methyl salicylate) contains an aromatic ring. Three of the eleven single-amino-acid-substitution variants of HbHNL reacted more rapidly with mandelonitrile. The best was HbHNL-L121Y, with a kcat 4.2 times higher and high enantioselectivity. Site-saturation mutagenesis at position 121 identified three other improved variants. We hypothesize that the smaller active site orients the aromatic substrate more productively.

Functional modulation of a protein folding landscape via side-chain distortion

Ultrahigh-resolution (< 1.0 Å) structures have revealed unprecedented and unexpected details of molecular geometry, such as the deformation of aromatic rings from planarity. However, the functional utility of such energetically costly strain is unknown. The 0.83 Å structure of α-lytic protease (αLP) indicated that residues surrounding a conserved Phe side-chain dictate a rotamer which results in a ~6° distortion along the side-chain, estimated to cost 4 kcal/mol. By contrast, in the closely related protease Streptomyces griseus Protease B (SGPB), the equivalent Phe adopts a different rotamer and is undistorted. Here, we report that the αLP Phe side-chain distortion is both functional and conserved in proteases with large pro regions. Sequence analysis of the αLP serine protease family reveals a bifurcation separating those sequences expected to induce distortion and those that would not, which correlates with the extent of kinetic stability. Structural and folding kinetics analyses of family members suggest that distortion of this side-chain plays a role in increasing kinetic stability within the αLP family members that use a large Pro region. Additionally, structural and kinetic folding studies of mutants demonstrate that strain alters the folding free energy landscape by destabilizing the transition state (TS) relative to the native state (N). Although side-chain distortion comes at a cost of foldability, it suppresses the rate of unfolding, thereby enhancing kinetic stability and increasing protein longevity under harsh extracellular conditions. This ability of a structural distortion to enhance function is unlikely to be unique to αLP family members and may be relevant in other proteins exhibiting side-chain distortions.

The influence of active site loop mutations on the thermal stability of azurin from Pseudomonas aeruginosa

The copper site and overall structures of azurin (AZ) variants in which the amicyanin (AMI) and plastocyanin (PC) metal binding loops have been introduced, AZAMI and AZPC, respectively, are similar to that of AZ, whereas the loop conformations resemble those in the native proteins. To assess the influence of these loop mutations on stability, the thermal unfolding of AZAMI and AZPC has been investigated by differential scanning calorimetry, absorption and fluorescence spectroscopy. The calorimetric profiles of both variants exhibit a complex shape consisting of two endothermic peaks and an exothermic peak. The temperature of the maximum heat of absorption for the single endothermic peak is 82.7°C for AZ, whereas for AZAMI and AZPC the most intense endothermic peaks are at 74.9 and 68.1°C comparable to values for AMI and PC, respectively. Denaturation investigated using the temperature dependence of the absorbance at ∼600nm and Trp emission, also demonstrates decreased stability for both loop mutants. The thermal transition between the native and the denaturated states is irreversible, scan rate dependent and consistent with the two-state irreversible model. The structure of the active-site loop has a dramatic effect on the kinetic stability and the unfolding pathway of cupredoxins.

Thermodynamic and kinetic stability of a large multi-domain enzyme from the hyperthermophile Aeropyrum pernix

The multi-domain enzyme isocitrate dehydrogenase from the hyperthermophile Aeropyrum pernix was studied by denaturant-induced unfolding. At pH 7.5, changes in circular dichroism ellipticity and intrinsic fluorescence showed a complex unfolding transition, whereas at pH 3.0, an apparently two-state and highly reversible unfolding occurred. Analytical ultracentrifugation revealed the dissociation from dimer to monomer at pH 3.0. The thermodynamic and kinetic stability were studied at pH 3.0 to explore the role of inter-domain interactions independently of inter-subunit interplay on the wild type and R211M, a mutant where a seven-membered inter-domain ionic network has been disrupted. The unfolding and folding transitions occurred at slightly different denaturant concentrations even after prolonged equilibration time. The difference between the folding and the unfolding profiles was decreased in the mutant R211M. The apparent Gibbs free energy decreased approximately 2 kcal/mol and the unfolding rate increased 4.3-fold in the mutant protein, corresponding to a decrease in activation free energy of unfolding of 0.86 kcal/mol. These results suggest that the inter-domain ionic network might be responsible for additional stabilization through a significant kinetic barrier in the unfolding pathway that could also explain the larger difference observed between the folding and unfolding transitions of the wild type.

Thermo- and mesostabilizing protein interactions identified by temperature-dependent statistical potentials

The goal of controlling protein thermostability is tackled here through establishing, by in silico analyses, the relative weight of residue-residue interactions in proteins as a function of temperature. We have designed for that purpose a (melting-) temperature-dependent, statistical distance potential, where the interresidue distances are computed between the side-chain geometric centers or their functional centers. Their separate derivation from proteins of either high or low thermal resistance reveals the interactions that contribute most to stability in different temperature ranges. Thermostabilizing interactions include salt bridges and cation-pi interactions (especially those involving arginine), aromatic interactions, and H-bonds between negatively charged and some aromatic residues. In contrast, H-bonds between two polar noncharged residues or between a polar noncharged residue and a negatively charged residue are relatively less stabilizing at high temperatures. An important observation is that it is necessary to consider both repulsive and attractive interactions in overall thermostabilization, as the degree of repulsion may also vary with temperature. These temperature-dependent potentials are not only useful for the identification of meso- and thermostabilizing pair interactions, but also exhibit predictive power, as illustrated by their ability to predict the melting temperature of a protein based on the melting temperature of homologous proteins.

Unfolding simulations reveal the mechanism of extreme unfolding cooperativity in the kinetically stable alpha-lytic protease

Kinetically stable proteins, those whose stability is derived from their slow unfolding kinetics and not thermodynamics, are examples of evolution's best attempts at suppressing unfolding. Especially in highly proteolytic environments, both partially and fully unfolded proteins face potential inactivation through degradation and/or aggregation, hence, slowing unfolding can greatly extend a protein's functional lifetime. The prokaryotic serine protease α-lytic protease (αLP) has done just that, as its unfolding is both very slow (t_1/2 ∼1 year) and so cooperative that partial unfolding is negligible, providing a functional advantage over its thermodynamically stable homologs, such as trypsin. Previous studies have identified regions of the domain interface as critical to αLP unfolding, though a complete description of the unfolding pathway is missing. In order to identify the αLP unfolding pathway and the mechanism for its extreme cooperativity, we performed high temperature molecular dynamics unfolding simulations of both αLP and trypsin. The simulated αLP unfolding pathway produces a robust transition state ensemble consistent with prior biochemical experiments and clearly shows that unfolding proceeds through a preferential disruption of the domain interface. Through a novel method of calculating unfolding cooperativity, we show that αLP unfolds extremely cooperatively while trypsin unfolds gradually. Finally, by examining the behavior of both domain interfaces, we propose a model for the differential unfolding cooperativity of αLP and trypsin involving three key regions that differ between the kinetically stable and thermodynamically stable classes of serine proteases.

Proteins, synthesized as linear polymers of amino acids, fold up into compact native states, burying their hydrophobic amino acids into their interiors. Protein folding minimizes the non-specific interactions that unfolded protein chains can make, which include aggregation with other proteins and degradation by proteases. Unfortunately, even in the native state, proteins can partially unfold, opening up regions of their structure and making these adverse events possible. Some proteins, particularly those in harsh environments full of proteases, have evolved to virtually eliminate partial unfolding, significantly reducing their rate of degradation. This elimination of partial unfolding is termed “cooperative,” because unfolding is an all-or-none process. One class of proteins has diverged into two families, one bacterial and highly cooperative and the other animal and non-cooperative. We have used detailed simulations of unfolding for members of each family, α-lytic protease (bacterial) and trypsin (animal) to understand the unfolding pathways of each and the mechanism for the differential unfolding cooperativity. Our results explain prior biochemical experiments, reproduce the large difference in unfolding cooperativity between the families, and point to the interface between α-lytic protease's two domains as essential to establishing unfolding cooperativity. As seen in an unrelated protein family, generation of a cooperative domain interface may be a common evolutionary response for ensuring the highest protein stability.

Intra-molecular chaperone: the role of the N-terminal in conformational selection and kinetic control

The vast majority of the proteins in nature are under thermodynamic control, consistent with the universally accepted notion that proteins exist in their thermodynamically most stable state. Yet, recently a number of examples of proteins whose fold is under kinetic control have come to light. Their functions and environments vary. The first among these are some proteases, discovered in the early 1990s. There, an N-terminal proregion is self-cleaved after the protein folded, leaving the remainder of the chain in a kinetically trapped state. A related scenario was observed for microcin J25, an antibacterial peptide. This peptide presents a trapped covalently knotted conformation. The third and the most recently discovered case is the multidrug-resistant transporter protein, P-glycoprotein. There, a synonymous 'silent' mutation leads to ribosome stalling with a consequent altered kinetically trapped state. Here we argue that in all three examples, the N-terminal plays the role of an intra-molecular chaperone, that is, the N-terminal conformation selects among all competing local conformations of a downstream segment. By providing a pattern, the N-terminal chaperone segment assists the protein folding process. If the N-terminal is subsequently cleaved, the protein can be under kinetic control, since it is trapped in a thermodynamically less-stable state.

Effect of proline substitutions on stability and kinetic properties of a cold adapted subtilase

A cold adapted subtilisin-like serine proteinase from a Vibrio species is two amino acids shorter at the N-terminus than related enzymes adapted to higher temperatures and has a 15 residues' C-terminal extension relative to the highly homologous thermophilic enzyme aqualysin I from Thermus aquaticus. These enzymes are produced as pro-enzymes with an N-terminal chaperone sequence for correct folding and a C-terminal signal peptide for secretion, which are subsequently cleaved off by autocatalysis to give the mature enzyme. A truncated form of the Vibrio proteinase where the C-terminal extension was removed and two residues near the N-terminus were substituted with proline, to resemble the N- and C-terminal regions in aqualysin I, resulted in increased thermostability and diminished catalytic efficiency. The proline substitutions shift the site of autocatalytic cleavage at the N-terminus by two amino acids, apparently by rigidifying the terminal residues and support the formation of a beta-sheet that fixes the N-terminus to the main body of the protein.

Between-species variation in the kinetic stability of TIM proteins linked to solvation-barrier free energies

Theoretical, computational, and experimental studies have suggested the existence of solvation barriers in protein unfolding and denaturation processes. These barriers are related to the finite size of water molecules and can be envisioned as arising from the asynchrony between water penetration and breakup of internal interactions. Solvation barriers have been proposed to play roles in protein cooperativity and kinetic stability; therefore, they may be expected to be subject to natural selection. We study the thermal denaturation, in the presence and in the absence of chemical denaturants, of triosephosphate isomerases (TIMs) from three different species: Trypanosoma cruzi, Trypanosoma brucei, and Leishmania mexicana. In all cases, denaturation was irreversible and kinetically controlled. Surprisingly, however, we found large differences between the kinetic denaturation parameters, with T. cruzi TIM showing a much larger activation energy value (and, consequently, much lower room-temperature, extrapolated denaturation rates). This disparity cannot be accounted for by variations in the degree of exposure to solvent in transition states (as measured by kinetic urea m values) and is, therefore, to be attributed mainly to differences in solvation-barrier contributions. This was supported by structure-energetics analyses of the transition states and by application of a novel procedure to estimate from experimental data the solvation-barrier impact at the entropy and free-energy levels. These analyses were actually performed with an extended protein set (including six small proteins plus seven variants of lipase from Thermomyces lanuginosus and spanning a wide range of activation parameters), allowing us to delineate the general trends of the solvation-barrier contributions. Overall, this work supports that proteins sharing the same structure and function but belonging to different organisms may show widely different solvation barriers, possibly as a result of different levels of the selection pressure associated with cooperativity, kinetic stability, and related factors.

Thermal stabilization of the protozoan Entamoeba histolytica alcohol dehydrogenase by a single proline substitution

Analysis of the three-dimensional structures of two closely related thermophilic and hyperthermophilic alcohol dehydrogenases (ADHs) from the respective microorganisms Entamoeba histolytica (EhADH1) and Thermoanaerobacter brockii (TbADH) suggested that a unique, strategically located proline residue (Pro275) at the center of the dimerization interface might be crucial for maintaining the thermal stability of TbADH. To assess the contribution of Pro275 to the thermal stability of the ADHs, we applied site-directed mutagenesis to replace Asp275 of EhADH1 with Pro (D275P-EhADH1) and conversely Pro275 of TbADH with Asp (P275D-TbADH). The results indicate that replacing Asp275 with Pro significantly enhances the thermal stability of EhADH1 (DeltaT(1/2) <or= +10 degrees C), whereas the reverse mutation in the thermophilic TbADH (P275D-TbADH) reduces the thermostability of the enzyme (DeltaT(1/2) <or= -18.8 degrees C). Analysis of the crystal structures of the thermostabilized mutant D275P-EhADH1 and the thermocompromised mutant P275D-TbADH suggest that a proline residue at position 275 thermostabilized the enzymes by reducing flexibility and by reinforcing hydrophobic interactions at the dimer-dimer interface of the tetrameric ADHs.

Yeast cytosine deaminase mutants with increased thermostability impart sensitivity to 5-fluorocytosine

Prodrug gene therapy (PGT) is a treatment strategy in which tumor cells are transfected with a 'suicide' gene that encodes a metabolic enzyme capable of converting a nontoxic prodrug into a potent cytotoxin. One of the most promising PGT enzymes is cytosine deaminase (CD), a microbial salvage enzyme that converts cytosine to uracil. CD also converts 5-fluorocytosine (5FC) to 5-fluorouracil, an inhibitor of DNA synthesis and RNA function. Over 150 studies of CD-mediated PGT applications have been reported since 2000, all using wild-type enzymes. However, various forms of CD are limited by inefficient turnover of 5FC and/or limited thermostability. In a previous study, we stabilized and extended the half-life of yeast CD (yCD) by repacking of its hydrophobic core at several positions distant from the active site. Here we report that random mutagenesis of residues selected based on alignment with similar enzymes, followed by selection for enhanced sensitization to 5FC, also produces an enzyme variant (yCD-D92E) with elevated T(m) values and increased activity half-life. The new mutation is located at the enzyme's dimer interface, indicating that independent mutational pathways can lead to an increase in stability, as well as a more subtle effect on enzyme kinetics. Each independently derived set of mutations significantly improves the enzyme's performance in PGT assays both in cell culture and in animal models.

Kinetic folding of Haloferax volcanii and Escherichia coli dihydrofolate reductases: haloadaptation by unfolded state destabilization at high ionic strength

Salts affect protein stability by multiple mechanisms (e.g., the Hofmeister effect, preferential hydration, electrostatic effects and weak ion binding). These mechanisms can affect the stability of both the native state and the unfolded state. Previous equilibrium stability studies demonstrated that KCl stabilizes dihydrofolate reductases (DHFRs) from Escherichia coli (ecDHFR, E. coli DHFR) and Haloferax volcanii (hvDHFR1, H. volcanii DHFR encoded by the hdrA gene) with similar efficacies, despite adaptation to disparate physiological ionic strengths (0.2 M versus 2 M). Kinetic studies can provide insights on whether equilibrium effects reflect native state stabilization or unfolded state destabilization. Similar kinetic mechanisms describe the folding of urea-denatured ecDHFR and hvDHFR1: a 5-ms stopped-flow burst-phase species that folds to the native state through two sequential intermediates with relaxation times of 0.1-3 s and 25-100 s. The latter kinetic step is very similar to that observed for the refolding of hvDHFR1 from low ionic strength. The unfolding of hvDHFR1 at low ionic strength is relatively slow, suggesting kinetic stabilization as observed for some thermophilic enzymes. Increased KCl concentrations slow the urea-induced unfolding of ecDHFR and hvDHFR1, but much less than expected from equilibrium studies. Unfolding rates extrapolated to 0 M denaturant, k(unf)(H(2)O), are relatively independent of ionic strength, demonstrating that the KCl-induced stabilization of ecDHFR and hvDHFR1 results predominantly from destabilization of the unfolded state. This supports the hypothesis from previous equilibrium studies that haloadaptation harnesses the effects of elevated salt concentrations on the properties of the aqueous solvent to enhance protein stability.