Energetics of side chain packing in staphylococcal nuclease assessed by exchange of valines, isoleucines, and leucines

To examine the importance of side chain packing to protein stability, each of the 11 leucines in staphylococcal nuclease was substituted with isoleucine and valine. The nine valines were substituted with leucine and isoleucine, while the five isoleucines, previously substituted with valine, were substituted with leucine and methionine. These substitutions conserve the hydrophobic character of these side chains but alter side chain geometry and, in some cases, size. In addition, eight threonine residues, previously substituted with valine, were substituted with isoleucine to test the importance of packing at sites normally not occupied by a hydrophobic residue. The stabilities of these 58 mutant proteins were measured by guanidine hydrochloride denaturation. To the best of our knowledge, this is the largest library of single packing mutants yet characterized. As expected, repacking stability effects are tied to the degree of side chain burial. The average energetic cost of moving a single buried methyl group was 0.9 kcal/mol, albeit with a standard deviation of 0.8 kcal/mol. This average is actually slightly greater than the value of 0.7-0.8 kcal/mol estimated for the hydrophobic transfer energy of a methylene from octanol to water. These results appear to indicate that van der Waals interactions gained from optimal packing are at least as important in stabilizing the native state of proteins as hydrophobic transfer effects.     

The effect of pressure and guanidine hydrochloride on azurins mutated in the hydrophobic core.

The unfolding of the blue-copper protein azurin from Pseudomonas aeruginosa by guanidine hydrochloride, under nonreducing conditions, has been studied by fluorescence techniques and circular dichroism. The denaturation transition may be fitted by a simple two-state model. The total free energy change from the native to the unfolded state was 9.4 +/- 0.4 kcal.mol-1, while a lower value (6.4 +/- 0.4 kcal.mol-1) was obtained for the metal depleted enzyme (apo-azurin) suggesting that the copper atom plays an important stabilization role. Azurin and apo-azurin were practically unaffected by hydrostatic pressure up to 3000 bar. Site-directed mutagenesis has been used to destabilize the hydrophobic core of azurin. In particular either hydrophobic residue Ile7 or Phe110 has been substituted with a serine. The free energy change of unfolding by guanidinium hydrochloride, resulted to be 5.8 +/- 0.3 kcal.mol-1 and 4.8 +/- 0.3 kcal.mol-1 for Ile7Ser and Phe110Ser, respectively, showing that both mutants are much less stable than the wild-type protein. The mutated apoproteins could be reversible denatured even by high pressure, as demonstrated by steady-state fluorescence measurements. The change in volume associated to the pressure-induced unfolding was estimated to be -24 mL.mol-1 for Ile7Ser and -55 mL.mol-1 for Phe110Ser. These results show that the tight packing of the hydrophobic residues that characterize the inner structure of azurin is fundamental for the protein stability. This suggests that the proper assembly of the hydrophobic core is one of the earliest and most crucial event in the folding process, bearing important implication for de novo design of proteins.     

Stabilization of GroEL minichaperones by core and surface mutations

We report the crystal structures of two hexa-substituted mutants of a GroEL minichaperone that are more stable than wild-type by 7.0 and 6.1 kcal mol(-1). Their structures imply that the increased stability results from multiple factors including improved hydrophobic packing, optimised hydrogen bonding and favourable structural rearrangements. It is commonly believed that protein core residues are immutable and generally optimized for energy, while on the contrary, surface residues are variable and hence unimportant for stability. But, it is now becoming clear that mutations of both core and surface residues can increase protein stability, and that protein cores are more flexible and thus more tolerant to mutation than expected. Sequence comparison of homologous proteins has provided a way to pinpoint the residues that contribute constructively to stability and to guide the engineering of protein stability. Stabilizing mutations identified by this approach are most frequently located at protein surfaces but with a few found in protein cores. In the latter case, local flexibility in the hydrophobic core is the key factor that allows the energetically favourable burial of larger hydrophobic side-chains without undue energetic penalties from steric clashes.     

Nine hydrophobic side chains are key determinants of the thermodynamic stability and oligomerization status of tumour suppressor p53 tetramerization domain.

The contribution of almost each amino acid side chain to the thermodynamic stability of the tetramerization domain (residues 326-353) of human p53 has been quantitated using 25 mutants with single-residue truncations to alanine (or glycine). Truncation of either Leu344 or Leu348 buried at the tetramer interface, but not of any other residue, led to the formation of dimers of moderate stability (8-9 kcal/mol of dimer) instead of tetramers. One-third of the substitutions were moderately destabilizing (<3.9 kcal/mol of tetramer). Truncations of Arg333, Asn345 or Glu349 involved in intermonomer hydrogen bonds, Ala347 at the tetramer interface or Thr329 were more destabilizing (4.1-5.7 kcal/mol). Strongly destabilizing (8.8- 11.7 kcal/mol) substitutions included those of Met340 at the tetramer interface and Phe328, Arg337 and Phe338 involved peripherally in the hydrophobic core. Truncation of any of the three residues involved centrally in the hydrophobic core of each primary dimer either prevented folding (Ile332) or allowed folding only at high protein concentration or low temperature (Leu330 and Phe341). Nine hydrophobic residues per monomer constitute critical determinants for the stability and oligomerization status of this p53 domain.     

Rotamer strain energy in protein helices - quantification of a major force opposing protein folding

It is widely believed that the dominant force opposing protein folding is the entropic cost of restricting internal rotations. The energetic changes from restricting side-chain torsional motion are more complex than simply a loss of conformational entropy, however. A second force opposing protein folding arises when a side-chain in the folded state is not in its lowest-energy rotamer, giving rotameric strain. chi strain energy results from a dihedral angle being shifted from the most stable conformation of a rotamer when a protein folds. We calculated the energy of a side-chain as a function of its dihedral angles in a poly(Ala) helix. Using these energy profiles, we quantify conformational entropy, rotameric strain energy and chi strain energy for all 17 amino acid residues with side-chains in alpha-helices. We can calculate these terms for any amino acid in a helix interior in a protein, as a function of its side-chain dihedral angles, and have implemented this algorithm on a web page. The mean change in rotameric strain energy on folding is 0.42 kcal mol-1 per residue and the mean chi strain energy is 0.64 kcal mol-1 per residue. Loss of conformational entropy opposes folding by a mean of 1.1 kcal mol-1 per residue, and the mean total force opposing restricting a side-chain into a helix is 2.2 kcal mol-1. Conformational entropy estimates alone therefore greatly underestimate the forces opposing protein folding. The introduction of strain when a protein folds should not be neglected when attempting to quantify the balance of forces affecting protein stability. Consideration of rotameric strain energy may help the use of rotamer libraries in protein design and rationalise the effects of mutations where side-chain conformations change.     

Acetyl coenzyme A binding by chloramphenicol acetyltransferase: long-range electrostatic determinants of coenzyme A recognition.

The possible involvement of arginyl and lysyl side chains of chloramphenicol acetyltransferase (CAT) in binding coenzyme A (CoA) was studied by means of chemical modification, site-directed mutagenesis, variation in ionic strength, use of competitive inhibitors or substrate analogues, and X-ray crystallography. Unlike a number of enzymes, including citrate synthase, CAT does not employ specific ion pairs with the phosphoanionic centers of CoA to bind the acetyl donor, and arginyl residues play no role in recognition of the coenzyme. Although phenylglyoxal inactivates CAT reversibly, it does so by the formation of an unstable adduct with a thiol group, that of Cys-31 in the chloramphenicol binding site. The inhibitory effect of increasing ionic strength on kcat/Km(acetyl-CoA) can be explained by long-range electrostatic interactions between CoA and the epsilon-amino groups of Lys-54 and Lys-177, both of which are solvent-accessible. The epsilon-amino group of Lys-54 contributes 1.3 kcal.mol-1 to the binding of acetyl-CoA via interactions with both the 3'- and 5'-phosphoanions of CoA. Lys-177 contributes only 0.4 kcal.mol-1 to the productive binding of acetyl-CoA, mediated by long-range (approximately 14 A) interactions with the 5'-alpha- and -beta-phosphoanions of CoA. The combined energetic contribution of Lys-54 and Lys-177 to acetyl-CoA binding (1.7 kcal.mol-1) is less than that previously demonstrated (2.4 kcal.mol-1) for a simple hydrophobic interaction between Tyr-178 and the adenine ring of CoA (Day & Shaw, 1992).(ABSTRACT TRUNCATED AT 250 WORDS)     

Core side-chain packing and backbone conformation in Lpp-56 coiled-coil mutants

Native proteins exhibit precise geometric packing of atoms in their hydrophobic interiors. Nonetheless, controversy remains about the role of core side-chain packing in specifying and stabilizing the folded structures of proteins. Here we investigate the role of core packing in determining the conformation and stability of the Lpp-56 trimerization domain. The X-ray crystal structures of Lpp-56 mutants with alanine substitutions at two and four interior core positions reveal trimeric coiled coils in which the twist of individual helices and the helix-helix spacing vary significantly to achieve the most favored superhelical packing arrangement. Introduction of each alanine "layer" into the hydrophobic core destabilizes the superhelix by 1.4 kcal mol(-1). Although the methyl groups of the alanine residues pack at their optimum van der Waals contacts in the coiled-coil trimer, they provide a smaller component of hydrophobic interactions than bulky hydrophobic side-chains to the thermodynamic stability. Thus, specific side-chain packing in the hydrophobic core of coiled coils are important determinants of protein main-chain conformation and stability.     

Sequence context modulates the stability of a GxxxG-mediated transmembrane helix-helix dimer

To quantify the relationship between sequence and transmembrane dimer stability, a systematic mutagenesis and thermodynamic study of the protein-protein interaction residues in the glycophorin A transmembrane helix-helix dimer was carried out. The results demonstrate that the glycophorin A transmembrane sequence dimerizes when its GxxxG motif is abolished by mutation to large aliphatic residues, suggesting that the sequence encodes an intrinsic propensity to self-associate independent of a GxxxG motif. In the presence of an intact GxxxG motif, the glycophorin A dimer stability can be modulated over a span of -0.5 kcal mol(-1) to +3.2 kcal mol(-1) by mutating the surrounding sequence context. Thus, these flanking residues play an active role in determining the transmembrane dimer stability. To assess the structural consequences of the thermodynamic effects of mutations, molecular models of mutant transmembrane domains were constructed, and a structure-based parameterization of the free energy change due to mutation was carried out. The changes in association free energy for glycophorin A mutants can be explained primarily by changes in packing interactions at the protein-protein interface. The energy cost of removing favorable van der Waals interactions was found to be 0.039 kcal mol(-1) per A2 of favorable occluded surface area. The value corresponds well with estimates for mutations in bacteriorhodopsin as well as for those mutations in the interiors of soluble proteins that create packing defects.     

Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers.

The aim of this study was to examine the differences between hydrophobicity and packing effects in specifying the three-dimensional structure and stability of proteins when mutating hydrophobes in the hydrophobic core. In DNA-binding proteins (leucine zippers), Leu residues are conserved at positions "d," and beta-branched amino acids, Ile and Val, often occur at positions "a" in the hydrophobic core. In order to discern what effect this selective distribution of hydrophobes has on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers, three Val or three Ile residues were simultaneously substituted for Leu at either positions "a" (9, 16, and 23) or "d" (12, 19, and 26) in both chains of a model coiled coil. The stability of the resulting coiled coils was monitored by CD in the presence of Gdn.HCl. The results of the mutations of Ile to Val at either positions "a" or "d" in the reduced or oxidized coiled coils showed a significant hydrophobic effect with the additional methylene group in Ile stabilizing the coiled coil (delta delta G values range from 0.45 to 0.88 kcal/mol/mutation). The results of mutations of Leu to Ile or Val at positions "a" in the reduced or oxidized coiled coils showed a significant packing effect in stabilizing the coiled coil (delta delta G values range from 0.59 to 1.03 kcal/mol/mutation). Our results also indicate the subtle control hydrophobic packing can have not only on protein stability but on the conformation adopted by the amphipathic alpha-helices. These structural findings correlate with the observation that in DNA-binding proteins, the conserved Leu residues at positions "d" are generally less tolerant of amino acid substitutions than the hydrophobic residues at positions "a."     

The role of packing interactions in stabilizing folded proteins.

In order to investigate the role of nonpolar side chains in determining protein stability, we have carried out a molecular dynamics simulation study of the thermodynamics of interconverting isoleucine and valine side chains in the core of ribonuclease T1. The free energy change in the unfolded state, which we take to be fully solvated, was small and agrees qualitatively with experimental studies of alkane solvation. In the two Ile----Val mutations studied, the protein was able to relax around the smaller side chains, while in the case of the two Val----Ile mutations, the ability of the core to accommodate the extra methylene group depended on where the mutation took place. We argue that the experimentally observed decrease in stability for mutating isoleucine into valine results from a loss of favorable packing interactions of the side chain in the folded form of the protein. This supports the view that packing interactions in the folded state are an important contributor to the overall stability of the folded protein and that the core of the native protein is packed efficiently and almost completely.     

Co-operative interactions during protein folding.

The theory for measuring co-operativity between interactions in proteins by protein engineering experiments is developed by introducing a procedure for analysing increasing orders of synergy in a protein with increasing numbers of residues. The (pairwise) interaction energy (delta 2Gint) between two side-chains may be measured experimentally by a double-mutant cycle consisting of the wild-type protein, the two single mutants and the double mutant. This procedure may be extended to three residues to give a value for delta 3Gint for a triple-mutant cube, and to higher orders using multi-dimensional mutant space. We now show that delta 3Gint is the excess energy of adding all three chains compared with the sum of all the pairwise values of delta 2Gint for each of the constituent double-mutant cycles and the sum of all the single addition energies. This physical interpretation extends to higher orders of mutation. delta nGint (i.e. the interaction energy for n residues), thus, reveals the layers of synergy in interactions as a protein is built up. This procedure is applied to measuring changes in synergy during the refolding of barnase for the triad of salt-linked residues Asp8, Asp12 and Arg110, which are mutated to alanine residues. The value of delta 3Gint in the folded structure is 0.77(+/- 0.06) kcal mol-1 (i.e. the triad is 0.77 kcal mol-1 more stable than expected from the sum of the individual pairwise interactions and single contributions). The value of delta 3Gint is still significant in the transition state for unfolding (0.60(+/- 0.07) kcal mol-1) and in the folding intermediate (0.60(+/- 0.13 kcal mol-1)). These results show that synergistic interactions exist in barnase, in its transition state for unfolding and in a refolding intermediate. A direct measurement of the change of co-operativity between the folded state and the transition state for unfolding shows a decrease of 0.17(+/- 0.04) kcal mol-1, suggesting that the initial stages of protein unfolding may be accompanied by some loosening of structure in parts that still interact. The similar extent of co-operativity in the transition state for unfolding and the intermediate in refolding suggests that the intermediate is homogeneous, at least in the region of the salt-linked triad, as heterogeneity would lower the co-operativity.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

Thermodynamic investigations of proteins. III. Thermodynamic description of lysozyme.

Standard functions of enthalpy, entropy and the Gibbs energy of native and denatured lysozyme in the range of 0-100 degrees C and pH 1.5-7.0 are represented in three-dimensional projections. The denaturational Gibbs energy change reaches 16 kcal mol-1 at conditions of maximal protein stability (0 degrees C, pH 4.5-7.0) and equals 14.5 kcal mol-1 at 25 degrees C and neutral pH. This result was found to be in agreement with the data reported from guanidine hydrochloride denaturation studies. Partial thermodynamic functions of the conformational and ionizational changes of the protein are obtained from entropy and Gibbs-energy changes in denaturation. The conformational partial entropy and Gibbs-energy change are found to be independent of pH. The pH-dependent partial ionizational entropy and Gibbs-energy changes are induced by normalization of the ionization behaviour of buried groups and cause a decrease of protein stability.     

Contribution of the hydrophobic effect to globular protein stability.

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.     

Effect of active site residues in barnase on activity and stability.

We have mutated residues in the active site of the ribonuclease, barnase, in order to determine their effects on both enzyme activity and protein stability. Mutation of several of the positively charged residues that interact with the negatively charged RNA substrate (Lys27----Ala, Arg59----Ala and His102----Ala) causes large decreases in activity. This is accompanied, however, by an increase in stability. There is presumably electrostatic strain in the active site where positively charged side-chains are clustered. Mutation of several residues that make hydrogen bonds (Ser57----Ala, Asn58----Asp and Tyr103----Phe) causes smaller decreases in activity, but increases or has no effect on stability. Deletion of hydrogen bonding groups elsewhere in proteins has been found previously to decrease stability by 0.5 to 1.5 kcal mol-1. Conversely, we find that two mutations (Asp54----Asn and Gln104----Ala) decrease stability and increase activity. Another mutation (Glu73----Ala) decreases both activity and stability. It is clear that many residues in the active site do not contribute to stability and that for some, but not all, of the residues there is a compromise between activity and stability. This suggests that certain types of local instability may be necessary for substrate binding and catalysis by barnase. This has implications for the understanding of enzyme activity and the design of enzymes.     

Novel folding and stability defects cause a deficiency of human glutathione transferase omega 1

The polymorphic deletion of Glu-155 from human glutathione transferase omega1 (GSTO1-1) occurs in most populations. Although the recombinant ΔGlu-155 enzyme expressed in Escherichia coli is active, the deletion causes a deficiency of the active enzyme in vivo. The crystal structure and the folding/unfolding kinetics of the ΔGlu-155 variant were determined in order to investigate the cause of the rapid loss of the enzyme in human cells. The crystal structure revealed altered packing around the Glu-155 deletion, an increase in the predicted solvent-accessible area and a corresponding reduction in the buried surface area. This increase in solvent accessibility was consistent with an elevated Stern-Volmer constant. The unfolding of both the wild type and ΔGlu-155 enzyme in urea is best described by a three-state model, and there is evidence for the more pronounced population of an intermediate state by the ΔGlu-155 enzymes. Studies using intrinsic fluorescence revealed a free energy change around 14.4 kcal/mol for the wild type compared with around 8.6 kcal/mol for the ΔGlu-155 variant, which indicates a decrease in stability associated with the Glu-155 deletion. Urea induced unfolding of the wild type GSTO1-1 was reversible through an initial fast phase followed by a second slow phase. In contrast, the ΔGlu-155 variant lacks the slow phase, indicating a refolding defect. It is possible that in some conditions in vivo, the increased solvent-accessible area and the low stability of the ΔGlu-155 variant may promote its unfolding, whereas the refolding defect limits its refolding, resulting in GSTO1-1 deficiency.     

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.     

Investigation of the structural stability of hUBF HMG box 5 by native-state hydrogen exchange

HMG box 5 of human upstream binding factor (hUBF) consists of three alpha-helices arranged in an L-shape with a hydrophobic core embraced by these helices and stabilized by extensive hydrophobic interactions between nonpolar residues around the core. The GdmCl-induced equilibrium unfolding transition of HMG box 5 of hUBF was monitored by both circular dichroism (CD) and fluorescence spectra. A cooperative two-state unfolding process was observed. The unfolding free energy, DeltaGU(D2O), and the cooperativity of the unfolding reaction, m, are 4.6 +/- 0.16 kcal x mol-1 and 1.62 +/- 0.06 kcal x mol-1 x M-1, respectively. Native-state hydrogen exchange (NHX) experiments under EX2 conditions were performed. NHX results clearly show that the hydrophobic core among the three helices is a slow-exchange core. The three helices would not contribute equally to the stability of the native protein. Helix 3 appears to contribute the least to the stability. The NHX data have also allowed the local, subglobal, and global unfolding structures of hUBF HMG box 5 to be dissected, and common global and subglobal unfolding units were successfully detected.     

Energetics of complementary side-chain packing in a protein hydrophobic core

The energetics of complementary packing of nonpolar side chains in the hydrophobic core of a protein were analyzed by protein engineering experiments. We have made the mutations Ile----Val, Ile----Ala, and Leu----Ala in a region of the small bacterial ribonuclease barnase where the major alpha-helix packs onto the central beta-sheet. The destabilization resulting from the creation of cavities was determined by measuring the decrease in free energy of folding from reversible denaturation induced by urea, guanidinium chloride, or heat. The different methods give consistent and reproducible results. The loss in free energy of folding for the mutant proteins is 1.0-1.6 kcal/mol per methylene group removed. This exceeds by severalfold the values obtained from model experiments of the partitioning of relevant side chains between aqueous and nonpolar solvents. Much of this discrepancy arises because two surfaces are buried when a protein folds--both the amino acid side chain in question and the portions of the protein into which it packs. These experiments directly demonstrate that the interior packing of a protein is crucial in stabilizing its three-dimensional structure: the conversion of leucine or isoleucine to alanine in the hydrophobic core loses half the net free energy of folding of barnase with a concomitant decrease in yield of the expressed recombinant protein.     

Structure-activity relationships in engineered proteins: characterization of disruptive deletions in the alpha-ammonium group binding site of tyrosyl-tRNA synthetase

Residues Asp-78 and Gln-173 of the tyrosyl-tRNA synthetase of Bacillus stearothermophilus form part of the binding site for tyrosine by making hydrogen bonds with the alpha-ammonium group. Asp-38 is close enough to the group to make an important electrostatic contribution. Unlike other residues in the active site that have been studied by site-directed mutagenesis, Asp-38, Asp-78, and Gln-173 are part of hydrogen-bonded networks. Each of these residues has been mutated to an alanine, and the resultant mutants have been studied by kinetics to construct the difference energy diagrams for the formation of tyrosyl adenylate. In each example, the binding of tyrosine is weakened by about 2.5 kcal mol-1. But, unlike previous mutants, the dissociation of the second substrate, in this case ATP, is also seriously affected, being weakened by some 2 kcal mol-1 for TyrTS(Ala-78) and TyrTS(Ala-173). The energy of the transition state for the formation of tyrosyl adenylate is raised by 7.8 kcal mol-1 for the former and 4.5 kcal mol-1 for the latter mutant. Addition of these mutants to linear free energy plots constructed for the nondisruptive mutants in the accompanying study [Fersht, A. R., Leatherbarrow, R. J., & Wells, T. N. C. (1987) Biochemistry (preceding paper in this issue)] reveals large deviations of the data for TyrTS(Ala-38) and TyrTS(Ala-78) from the regression line. These thus belong to a different class of mutations from previous nondisruptive examples. This observation combined with the structural evidence and difference energy diagrams strongly suggests that the mutations Asp----Ala-38 and Asp----Ala-78 are disruptive in nature.