Tyrosine hydrogen bonds make a large contribution to protein stability

The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.     

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

Increasing protein stability: Importance of ΔCp and the denatured state

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site‐directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the β‐turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28°C higher than the wild‐type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, ΔC_p, by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.     

Contribution of hydrogen bonding to the conformational stability of ribonuclease T1.

For 30 years, the prevailing view has been that the hydrophobic effect contributes considerably more than hydrogen bonding to the conformational stability of globular proteins. The results and reasoning presented here suggest that hydrogen bonding and the hydrophobic effect make comparable contributions to the conformational stability of ribonuclease T1 (RNase T1). When RNase T1 folds, 86 intramolecular hydrogen bonds with an average length of 2.95 A are formed. Twelve mutants of RNase T1 [Tyr----Phe (5), Ser----Ala (3), and Asn----Ala (4)] have been prepared that remove 17 of the hydrogen bonds with an average length of 2.93 A. On the basis of urea and thermal unfolding studies of these mutants, the average decrease in conformational stability due to hydrogen bonding is 1.3 kcal/mol per hydrogen bond. This estimate is in good agreement with results from several related systems. Thus, we estimate that hydrogen bonding contributes about 110 kcal/mol to the conformational stability of RNase T1 and that this is comparable to the contribution of the hydrophobic effect. Accepting the idea that intramolecular hydrogen bonds contribute 1.3 +/- 0.6 kcal/mol to the stability of systems in an aqueous environment makes it easier to understand the stability of the "molten globule" states of proteins, and the alpha-helical conformations of small peptides.     

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.     

Contribution of single tryptophan residues to the fluorescence and stability of ribonuclease Sa

Ribonuclease Sa (RNase Sa) contains no tryptophan (Trp) residues. We have added single Trp residues to RNase Sa at sites where Trp is found in four other microbial ribonucleases, yielding the following variants of RNase Sa: Y52W, Y55W, T76W, and Y81W. We have determined crystal structures of T76W and Y81W at 1.1 and 1.0 A resolution, respectively. We have studied the fluorescence properties and stabilities of the four variants and compared them to wild-type RNase Sa and the other ribonucleases on which they were based. Our results should help others in selecting sites for adding Trp residues to proteins. The most interesting findings are: 1), Y52W is 2.9 kcal/mol less stable than RNase Sa and the fluorescence intensity emission maximum is blue-shifted to 309 nm. Only a Trp in azurin is blue-shifted to a greater extent (308 nm). This blue shift is considerably greater than observed for Trp71 in barnase, the Trp on which Y52W is based. 2), Y55W is 2.1 kcal/mol less stable than RNase Sa and the tryptophan fluorescence is almost completely quenched. In contrast, Trp59 in RNase T1, on which Y55W is based, has a 10-fold greater fluorescence emission intensity. 3), T76W is 0.7 kcal/mol more stable than RNase Sa, indicating that the Trp side chain has more favorable interactions with the protein than the threonine side chain. The fluorescence properties of folded Y76W are similar to those of the unfolded protein, showing that the tryptophan side chain in the folded protein is largely exposed to solvent. This is confirmed by the crystal structure of the T76W which shows that the side chain of the Trp is only approximately 7% buried. 4), Y81W is 0.4 kcal/mol less stable than RNase Sa. Based on the crystal structure of Y81W, the side chain of the Trp is 87% buried. Although all of the Trp side chains in the variants contribute to the unusual positive circular dichroism band observed near 235 nm for RNase Sa, the contribution is greatest for Y81W.     

Mutational and structural-based analyses of the osmolyte effect on protein stability

It is known that several naturally occurring substances known as osmolytes increase the conformational stability of proteins. Bolen and co-worker proposed the osmophobic theory, which asserts the osmolyte effect occurs because of an unfavorable interaction of osmolytes mainly with the protein backbone, based on the results on the transfer Gibbs energy of amino acids (Deltag) [Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963]. In this paper, we report the effect of sarcosine on the conformational stability (DeltaG) of RNase Sa (96 residues and one disulfide bond) and four mutant proteins. The thermal denaturation curves for RNase Sa in sarcosine fitted a two-state model on nonlinear least-squares analysis. All the RNase Sa proteins were stabilized by sarcosine. For example, the increase in stability of the wild-type protein in 4 M sarcosine due to the osmolyte effect (Delta(o)DeltaG) is 3.2 kcal/mol. Mutational analysis of the osmolyte effect indicated that the changed Delta(o)DeltaG values upon mutation (Delta(m)Delta(o)DeltaG), as estimated from the Deltag values, are similar to the experimental values. Structural-based analysis of the osmolyte effect was also performed using model denatured structures: (a) a fully extended model (single chain) with no disulfide bond, (b) two-part, unfolded models (two chains) with a disulfide bond constructed through molecular dynamic (MD) simulation, and (c) a two-part, folded model (two chains). The two-part, unfolded models were expected to be more suitable as denatured structures. The Delta(o)DeltaG values calculated using the two-part, unfolded models were more consistent with experimental values than those calculated using the fully extended and two-part, folded models. This suggests that MD simulation is useful for testing denatured structures. These results indicate that the osmophobic theory can explain the osmolyte effect on protein stability.     

Comparison of local and global stability of an analogue of a disulfide-folding intermediate with those of the wild-type protein in bovine pancreatic ribonuclease A: identification of specific regions of stable structure along the oxidative folding pathway

We have identified specific regions of the polypeptide chain of bovine pancreatic ribonuclease A (RNase A) that are critical for stabilizing the oxidative folding intermediate des-[40-95] (with three native disulfide bonds but lacking the fourth native Cys40-Cys95 disulfide bond) in an ensemble of largely disordered three-disulfide precursors (3S if des-[40-95]). A stable analogue of des-[40-95], viz., [C40A, C95A] RNase A, which contains three out of four native disulfide pairings, was previously found to have a three-dimensional structure very similar to that of the wild-type protein. However, it is determined here from GdnHCl denaturation experiments to have significantly reduced global stability, i.e., = 4.5 kcal /mol at 20 degrees C and pH 4.6. The local stability of [C40A, C95A] RNase A was also examined using site-specific amide (2)H/(1)H exchange measurements at pD 5.0 to determine the individual unfolding free energy of specific residues under both strongly native (12 degrees C) and more destabilizing (20 degrees C) conditions. Comparison of the relative stabilities at specific amide sites of [C40A, C95A] RNase A at both temperatures with the corresponding values for the wild-type protein at 35 degrees C corroborates previous experimental evidence that unidentified intramolecular contacts in the vicinity of the preferentially formed native one-disulfide (C65-C72) loop are crucial for stabilizing early folding intermediates, leading to des-[40-95]. Moreover, values of for residues at or near the third alpha-helix, and in part of the second beta-sheet of [C40A, C95A] RNase A, indicate that these two regions of regular backbone structure contribute to stabilizing the global chain fold of the des-[40-95] disulfide-folding intermediate in the wild-type protein. More significantly, we have identified numerous specific residues in the first alpha-helix and the first beta-sheet of the protein that are stabilized in the final step of the major oxidative regeneration pathway of RNase A (des-[40-95] --> N).     

Ribonuclease Sa conformational stability studied by NMR-monitored hydrogen exchange

The conformational stability of ribonuclease Sa (RNase Sa) has been measured at the per-residue level by NMR-monitored hydrogen exchange at pH* 5.5 and 30 degrees C. In these conditions, the exchange mechanism was found to be EXII. The conformational stability calculated from the slowest exchanging amide groups was found to be 8.8 kcal/mol, in close agreement with values determined by spectroscopic methods. RNase Sa is curiously rich in acidic residues (pI = 3.5) with most basic residues being concentrated in the active-site cleft. The effects of dissolved salts on the stability of RNase Sa was studied by thermal denaturation experiments in NaCl and GdmCl and by comparing hydrogen exchange rates in 0.25 M NaCl to water. The protein was found to be stabilized by salt, with the magnitude of the stabilization being influenced by the solvent exposure and local charge environment at individual amide groups. Amide hydrogen exchange was also measured in 0.25, 0.50, 0.75, and 1.00 M GdmCl to characterize the unfolding events that permit exchange. In contrast to other microbial ribonucleases studied to date, the most protected, globally exchanging amides in RNase Sa lie not chiefly in the central beta strands but in the 3/10 helix and an exterior beta strand. These structural elements are near the Cys7-Cys96 disulfide bond.     

The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins

We recently found that Asp-His-Ser/Thr-Trp hydrogen-bonded tetrads are widely and uniquely present in the WD40-repeat proteins. WDR5 protein is a seven WD40-repeat propeller with five such tetrads. To explore the effect of the tetrad on the structure and stability of WD40-repeat proteins, the wild-type WDR5 and its seven mutants involving the substitutions of tetrad residues have been isolated. The crystal structures of the wild-type WDR5 and its three WDR5 mutants have been determined by X-ray diffraction method. The mutations of the tetrad residues are found not to change the basic structural features. The denaturing profiles of the wild type and the seven mutants with the use of denaturant guanidine hydrochloride have been studied by circular dichroism spectroscopy to determine the folding free energies of these proteins. The folding free energies of the wild type and the S62A, S146A, S188A, D192E, W330F, W330Y, and D324E mutants are measured to be about -11.6, -2.7, -3.1, -2.9, -3.6, -7.1, -7.0, and -7.5 kcal/mol, respectively. These suggest that (1) the hydrogen bonds in these hydrogen bond networks are unusually strong; (2) each hydrogen-bonded tetrad provides over 12 kcal/mol stability to the protein; thus, the removal of any single tetrad would cause unfolding of the protein; (3) since there are five tetrads, the protein must be in a highly unstable state without the tetrads, which might be related to its biological functions.     

Increasing protein stability by altering long-range coulombic interactions.

It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.     

Entropy drives integrin alphaIIbbeta3:echistatin binding--evidence from surface plasmon resonance spectroscopy.

This investigation examined the molecular mechanisms that enable the alphaIIbbeta3 integrin to bind efficiently, tightly, and selectively to echistatin, an RGD disintegrin. We used surface plasmon resonance spectroscopy to measure the rate, extent, and stability of complexes formed between micellar alphaIIbbeta3 and recombinant echistatin (rEch) mutants, immobilized on the surface of a biosensor chip. alphaIIbbeta3 bound readily and tightly to wild-type RGD-rEch and RGDF-rEch but not to RGA-rEch or AGD-rEch, demonstrating that both of those charged moieties contribute to integrin recognition. van't Hoff analysis of the temperature dependence of the RGD-rEch K d data yielded an unfavorable enthalpy change, Delta H degrees = 14 +/- 3 kcal/mol, offset by a favorable entropy term, TDelta S degrees = 23 +/- 3 kcal/mol. Eyring analysis of the temperature dependence of the kinetic parameters yielded Delta H a degrees (++) = 9 +/- 2 kcal/mol and TDelta S a degrees (++) = -4 +/- 2 kcal/mol, indicating that a substantial activation enthalpy barrier and a smaller activation entropy hinder assembly of the encounter complex. Thus, equilibrium thermodynamic data demonstrate that entropy is the dominant factor stabilizing integrin:echistatin binding, while transition-state thermodynamic parameters indicate that the rate of complex formation is enthalpy-limited. When electrostatic contacts are the major source of receptor:ligand stability, theory and experiment indicate that entropy-favorable ion-pair desolvation often provides the driving force for molecular recognition.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

Stabilization of human RNase 1 by introduction of a disulfide bond between residues 4 and 118

In order to stabilize human RNase 1 by introduction of an intramolecular cross-link, a mutant protein (4-118CL RNase 1), in which Arg4 and Val118 are replaced with cysteine residues and linked by a disulfide bond, was designed and expressed in Escherichia coli as inclusion bodies. The 4-118CL RNase 1 that refolded under redox conditions was a monomer without free SH groups and retained 11% of the activity of the wild-type recombinant RNase 1, indicating that the mutant enzyme was correctly folded with the formation of an additional disulfide bond between Cys4 and Cys118. From guanidium chloride denaturation experiments based on the assumption of a two-state transition for unfolding, it was demonstrated that the introduction of the present cross-link increased the thermodynamic stability of RNase 1 by 2.0 kcal/mol. This value was lower than that, 5.4 kcal/mol, theoretically calculated from the reduction of chain entropy of the unfolded state due to the introduction of the cross-link. These results suggest that the present cross-link also destabilized the folded state of RNase 1 by 3.4 kcal/mol. Along with the increase in the thermodynamic stability, the stability of RNase 1 against trypsin digestion was also significantly increased by the introduction of this cross-link. It is likely, although not proven, that stabilized human RNases are favorable for clinical use, because human RNase-based immunotoxins should have long half-lives as to proteolytic degradation after endocytosis.     

Protein structural plasticity exemplified by insertion and deletion mutants in T4 lysozyme.

To further investigate the ways in which proteins respond to changes in the length of the polypeptide chain, a series of 32 insertions and five deletions were made within nine different alpha-helices of T4 lysozyme. In most cases, the inserted amino acid was a single alanine, although in some instances up to four residues, not necessarily alanine, were used. Different insertions destabilized the protein by different amounts, ranging from approximately 1 to 6 kcal/mol. In one case, no protein could be obtained. An "extension" mutant in which the carboxy terminus of the molecule was extended by four alanines increased stability by 0.3 kcal/mol. For the deletions, the loss in stability ranged from approximately 3 to 5 kcal/mol. The structures of six insertion mutants, as well as one deletion mutant and the extension mutant, were determined, three in crystal forms nonisomorphous with wild type. In all cases, including previously described insertion mutants within a single alpha-helix, there appears to be a strong tendency to preserve the helix by translocating residues so that the effects of the insertion are propagated into a bend or loop at one end or the other of the helix. In three mutants, even the hydrophobic core was disrupted so as to permit the preservation of the alpha-helix containing the insertion. Translocation (or "register shift") was also observed for the deletion mutant, in this case a loop at the end of the helix being shortened. In general, when translocation occurs, the reduction in stability is only moderate, averaging 2.5 kcal/mol. Only in the most extreme cases does "bulging" or "looping-out" occur within the body of an alpha-helix, in which case the destabilization is substantial, averaging 4.9 kcal/mol. Looping-out can occur for insertions close to the end of a helix, in which case the destabilization is less severe, averaging 2.6 kcal/mol. Mutant A73-[AAA] as well as mutants R119-[A] and V131-[A], include shifts in the backbone of 3-6 A, extending over 20 residues or more. As a result, residues 114-142, which form a "cap" on the carboxy-terminal domain, undergo substantial reorganizations such that the interface between this "cap" and the rest of the protein is altered substantially. In the case of mutant A73-[AAA], two nearby alpha-helices, which form a bend of approximately 105 degrees in the wild-type structure, reorganize in the mutant structure to form a single, essentially straight helix. These structural responses to mutation demonstrate the plasticity of protein structures and illustrate ways in which their three-dimensional structures might changes during evolution.     

Are turns required for the folding of ribonuclease T1?

Ribonuclease T1 (RNase T1) is a small, globular protein of 104 amino acids for which extensive thermodynamic and structural information is known. To assess the specific influence of variations in amino acid sequence on the mechanism for protein folding, circularly permuted variants of RNase T1 were constructed and characterized in terms of catalytic activity and thermodynamic stability. The disulfide bond connecting Cys-2 and Cys-10 was removed by mutation of these residues to alanine (C2, 10A) to avoid potential steric problems imposed by the circular permutations. The original amino-terminus and carboxyl-terminus of the mutant (C2, 10A) were subsequently joined with a tripeptide linker to accommodate a reverse turn and new termini were introduced throughout the primary sequence in regions of solvent-exposed loops at Ser-35 (cp35S1), Asp-49 (cp49D1), Gly-70 (cp70G1), and Ser-96 (cp96S1). These circularly permuted RNase T1 mutants retained 35-100% of the original catalytic activity for the hydrolysis of guanylyl(3'-->5')cytidine, suggesting that the overall tertiary fold of these mutants is very similar to that of wild-type protein. Chemical denaturation curves indicated thermodynamic stabilities at pH 5.0 of 5.7, 2.9, 2.6, and 4.6 kcal/mol for cp35S1, cp49D1, cp70G1, and cp96S1, respectively, compared to a value of 10.1 kcal/mol for wild-type RNase T1 and 6.4 kcal/mol for (C2, 10A) T1. A fifth set of circularly permuted variants was attempted with new termini positioned in a tight beta-turn between Glu-82 and Gln-85. New termini were inserted at Asn-83 (cp83N1), Asn-84 (cp84N1), and Gln-85 (cp85Q1). No detectable amount of protein was ever produced for any of the mutations in this region, suggesting that this turn may be critical for the proper folding and/or thermodynamic stability of RNase T1.     

Conformational stability and activity of ribonuclease T1 with zero, one, and two intact disulfide bonds

Ribonuclease T1 has two disulfide bonds linking cysteine residues 2-10 and 6-103. We have prepared a derivative of ribonuclease T1 in which one disulfide bond is broken and the cysteine residues carboxymethylated, (2-10)-RCM-T1, and three derivatives in which both disulfides are broken and the cysteine residues reduced, R-T1, carboxamidomethylated, RCAM-T1, or carboxymethylated, RCM-T1. The RNA hydrolyzing activity of these proteins has been measured, and urea and thermal denaturation studies have been used to determine conformational stability. The activity, melting temperature, and conformational stability of the proteins are: ribonuclease T1 (100%, 59.3 degrees C, 10.2 kcal/mol), (2-10)-RCM-T1 (86%, 53.3 degrees C, 6.8 kcal/mol), R-T1 (53%, 27.2 degrees C, 3.0 kcal/mol), RCAM-T1 (43%, 21.2 degrees C, 1.5 kcal/mol), and RCM-T1 (35%, 16.6 degrees C, 0.9 kcal/mol). Thus, the conformational stability is decreased by 3.4 kcal/mol when one disulfide bond is broken and by 7.2-9.3 kcal/mol when both disulfide bonds are broken. It is quite remarkable that RNase T1 can fold and function with both disulfide bonds broken and the cysteine residues carboxymethylated. The large decrease in the stability is due mainly to an increase in the conformational entropy of the unfolded protein which results when the constraints of the disulfide bonds on the flexibility are removed. We propose a new equation for predicting the effect of a cross-link on the conformational entropy of a protein: delta Sconf = -2.1 - (3/2)R 1n n, where n is the number of residues between the side chains which are cross-linked. This equation gives much better agreement with experimental results than other forms of this equation which have been used previously.     

Flipping the switch from monomeric to dimeric CV-N has little effect on antiviral activity

Cyanovirin-N can exist in solution in monomeric and domain-swapped dimeric forms, with HIV-antiviral activity being reported for both. Here we present results for CV-N variants that form stable solution dimers: the obligate dimer [DeltaQ50]CV-N and the preferential dimer [S52P]CV-N. These variants exhibit comparable DeltaG values (10.6 +/- 0.5 and 9.4 +/- 0.5 kcal.mol(-1), respectively), similar to that of stabilized, monomeric [P51G]CV-N (9.8 +/- 0.5 kcal.mol(-1)), but significantly higher than wild-type CV-N (4.1 +/- 0.2 kcal.mol(-1)). During folding/unfolding, no stably folded monomer was observed under any condition for the obligate dimer [DeltaQ50]CV-N, whereas two monomeric, metastable species were detected for [S52P]CV-N at low concentrations. This is in contrast to our previous results for [P51G]CV-N and wild-type CV-N, for which the dimeric forms were found to be the metastable species. The dimeric mutants exhibit comparable antiviral activity against HIV and Ebola, similar to that of wild-type CV-N and the stabilized [P51G]CV-N variant.     

Increasing protein conformational stability by optimizing beta-turn sequence

Protein conformational stability is an important concern in many fields. Here, we describe a strategy for significantly increasing conformational stability by optimizing beta-turn sequence. Proline and glycine residues are statistically preferred at several beta-turn positions, presumably because their unique side-chains contribute favorably to conformational stability in certain beta-turn positions. However, beta-turn sequences often deviate from preferred proline or preferred glycine. Therefore, our strategy involves replacing non-proline and non-glycine beta-turn residues with preferred proline or preferred glycine residues. Here, we develop guidelines for selecting appropriate mutations, and present results for five mutations (S31P, S42G, S48P, T76P, and Q77G) that significantly increase the conformational stability of RNase Sa. The increases in stability ranged from 0.7 kcal/mol to 1.3 kcal/mol. The strategy was successful in overlapping or isolated beta-turns, at buried (up to 50%) or completely exposed sites, and at relatively flexible or inflexible sites. Considering the significant number of beta-turn residues in every globular protein and the frequent deviation of beta-turn sequences from preferred proline and preferred glycine residues, this simple, efficient strategy will be useful for increasing the conformational stability of proteins.