Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme.

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.

Docking molecules by families to increase the diversity of hits in database screens: computational strategy and experimental evaluation.

Molecular docking programs screen chemical databases for novel ligands that fit protein binding sites. When one compound fits the site well, close analogs typically do the same. Therefore, many of the compounds that are found in such screens resemble one another. This reduces the variety and novelty of the compounds suggested. In an attempt to increase the diversity of docking hit lists, the Available Chemicals Directory was grouped into families of related structures. All members of every family were docked and scored, but only the best scoring molecule of a high-ranking family was allowed in the hit list. The identity and scores of the other members of these families were recorded as annotations to the best family member, but they were not independently ranked. This family-based docking method was compared with molecule-by-molecule docking in screens against the structures of thymidylate synthase, dihydrofolate reductase (DHFR), and the cavity site of the mutant T4 lysozyme Leu99 --> Ala (L99A). In each case, the diversity of the hit list increased, and more families of known ligands were found. To investigate whether the newly identified hits were sensible, we tested representative examples experimentally for binding to L99A and DHFR. Of the six compounds tested against L99A, five bound to the internal cavity. Of the seven compounds tested against DHFR, six inhibited the enzyme with apparent K(i) values between 0.26 and 100 microM. The segregation of potential ligands into families of related molecules is a simple technique to increase the diversity of candidates suggested by database screens. The general approach should be applicable to most docking methods. Proteins 2001;42:279-293.

Use of a non-rigid region in T4 lysozyme to design an adaptable metal-binding site.

It is not easy to find candidate sites within a given protein where the geometry of the polypeptide chain matches that of metal-binding sites in known protein structures. By choosing a location in T4 lysozyme that is inherently flexible, it was possible to engineer a two-histidine site that binds different divalent cations. Crystallographic analysis shows that the geometry of binding of zinc is distorted tetrahedral while that of cobalt and nickel is octahedral. Insofar as spectroscopic data can be measured, they indicate that similar modes of coordination are retained in solution. The two substitutions, Thr21 --> His and Thr142 --> His, lie, respectively, on the surface of the N- and C-terminal domains on opposite sides of the active site cleft. The design takes advantage of hinge-bending motion which allows the binding site to adapt to the most favorable ligand geometry for the metal. Introduction of the two histidines increases the melting temperature of the protein by 2.0 degrees C at pH 7.4. Metal binding further increases the melting temperature, but only by a small amount (up to 1.5 degrees C). A third substitution, Gln141 --> His, which could act as a third ligand in principle, does not do so, demonstrating the difficulty in mimicking naturally occurring metal-binding sites.

Solid-state synthesis and mechanical unfolding of polymers of T4 lysozyme.

Recent advances in single molecule manipulation methods offer a novel approach to investigating the protein folding problem. These studies usually are done on molecules that are naturally organized as linear arrays of globular domains. To extend these techniques to study proteins that normally exist as monomers, we have developed a method of synthesizing polymers of protein molecules in the solid state. By introducing cysteines at locations where bacteriophage T4 lysozyme molecules contact each other in a crystal and taking advantage of the alignment provided by the lattice, we have obtained polymers of defined polarity up to 25 molecules long that retain enzymatic activity. These polymers then were manipulated mechanically by using a modified scanning force microscope to characterize the force-induced reversible unfolding of the individual lysozyme molecules. This approach should be general and adaptable to many other proteins with known crystal structures. For T4 lysozyme, the force required to unfold the monomers was 64 +/- 16 pN at the pulling speed used. Refolding occurred within 1 sec of relaxation with an efficiency close to 100%. Analysis of the force versus extension curves suggests that the mechanical unfolding transition follows a two-state model. The unfolding forces determined in 1 M guanidine hydrochloride indicate that in these conditions the activation barrier for unfolding is reduced by 2 kcal/mol.

The introduction of strain and its effects on the structure and stability of T4 lysozyme.

In order to try to better understand the role played by strain in the structure and stability of a protein a series of "small-to-large" mutations was made within the core of T4 lysozyme. Three different alanine residues, one involved in backbone contacts, one in side-chain contacts, and the third adjacent to a small cavity, were each replaced with subsets of the larger residues, Val, Leu, Ile, Met, Phe and Trp. As expected, the protein is progressively destabilized as the size of the introduced side-chain becomes larger. There does, however, seem to be a limit to the destabilization, suggesting that a protein of a given size may be capable of maintaining only a certain amount of strain. The changes in stability vary greatly from site to site. Substitution of larger residues for both Ala42 and Ala98 substantially destabilize the protein, even though the primary contacts in one case are predominantly with side-chain atoms and in the other with backbone. The results suggest that it is neither practical nor meaningful to try to separate the effects of introduced strain on side-chains from the effects on the backbone. Substitutions at Ala129 are much less destabilizing than at sites 42 or 98. This is most easily understood in terms of the pre-existing cavity, which provides partial space to accommodate the introduced side-chains. Crystal structures were obtained for a number of the mutants. These show that the changes in structure to accommodate the introduced side-chains usually consist of essentially rigid-body displacements of groups of linked atoms, achieved through relatively small changes in torsion angles. On rare occasions, a side-chain close to the site of substitution may change to a different rotamer. When such rotomer changes occur, they permit the structure to dissipate strain by a response that is plastic rather than elastic. In one case, a surface loop moves 1.2 A, not in direct response to a mutation, but in an interaction mediated via an intermolecular contact. It illustrates how the structure of a protein can be modified by crystal contacts.

Substitution with selenomethionine can enhance the stability of methionine-rich proteins.

The availability of a series of phage T4 lysozymes with up to 14 methionine residues incorporated within the protein has made it possible to systematically compare the effect on protein stability of selenomethionine relative to methionine. Wild-type lysozyme contains two fully buried methionine residues plus three more on the surface. The substitution of these methionine residues with selenomethionine slightly stabilizes the protein. As more and more methionine residues are substituted into the protein, there is a progressive loss of stability. This is, however, increasingly offset in the selenomethionine variants, ultimately resulting in a differential increase in melting temperature of about 7 degrees C. This increase, corresponding to about 0.25 kcal/mol per substitution, is in reasonable agreement with the difference in the solvent transfer free energy between the two amino acids.

Methionine and alanine substitutions show that the formation of wild-type-like structure in the carboxy-terminal domain of T4 lysozyme is a rate-limiting step in folding.

In an attempt to identify a systematic relation between the structure of a protein and its folding kinetics, the rate of folding was determined for 20 mutants of T4 lysozyme in which a bulky, buried, nonpolar wild-type residue (Leu, Ile, Phe, Val, or Met) was substituted with alanine. Methionine, which approximated the size of the original side chain but which is of different shape and flexibility, was also substituted at most of the same sites. Mutations that substantially destabilize the protein and are located in the carboxy-terminal domain generally slow the rate of folding. Destabilizing mutations in the amino-terminal domain, however, have little effect on the rate of folding. Mutations that have little effect on stability tend to have little effect on the rate, no matter where they are located. These results suggest that, at the rate-limiting step, elements of structure in the C-terminal domain are formed and have a structure similar to that of the fully folded protein. Consistent with this, two variants that somewhat increase the rate of folding (Phe104 --> Met and Val149 --> Met) are located within the carboxy-terminal domain and maintain or improve packing with very little perturbation of the wild-type structure.

Structural analysis of a non-contiguous second-site revertant in T4 lysozyme shows that increasing the rigidity of a protein can enhance its stability.

The mutation Glu108-->Val (E108V) in T4 lysozyme was previously isolated as a second-site revertant that specifically compensated for the loss of function associated with the destabilizing substitution Leu99-->Gly (L99G). Surprisingly, the two sites are 11 A apart, with Leu99 in the core and Glu108 on the surface of the protein. In order to better understand this result we have carried out a detailed thermodynamic, enzymatic and structural analysis of these mutant lysozymes as well as a related variant with the substitution Leu99-->Ala. It was found that E108V does increase the stability of L99G, but it also increases the stability of both the wild-type protein and L99A by essentially equal amounts. The effects of E108V on enzymatic activity are more complicated. The mutation slightly reduces the maximal rate of cell wall hydrolysis of wild-type, L99G and L99A. At the same time, L99G is an unstable protein and rapidly loses activity during the course of the assay, especially at temperatures above 20 degrees C. Thus, even though the double mutant L99G/E108V has a slightly lower maximal rate than L99G, over a period of 20-30 minutes it hydrolyzes more substrate. This decrease in the rate of thermal inactivation appears to be the basis of the action of E108V as a second-site revertant of L99G. Mutant L99A creates a cavity of volume 149 A(3). Instead of enlarging this cavity, mutant L99G results in a 4-5 A displacement of part of helix F (residues 108-113), creating a solvent-accessible declivity. In the double mutant, L99G/E108V, this helix returns to a position akin to wild-type, resulting in a cavity of volume 203 A(3). Whether the mutation Glu108-->Val is incorporated into either wild-type lysozyme, or L99A or L99G, it results in a decrease in crystallographic thermal factors, especially in the helices that include residues 99 and 108. This increase in rigidity, which appears to be due to a combination of increased hydrophobic stabilization plus a restriction of conformational fluctuation, provides a structural basis for the increase in thermostability.

Structural characterization of an engineered tandem repeat contrasts the importance of context and sequence in protein folding.

To test a different approach to understanding the relationship between the sequence of part of a protein and its conformation in the overall folded structure, the amino acid sequence corresponding to an alpha-helix of T4 lysozyme was duplicated in tandem. The presence of such a sequence repeat provides the protein with "choices" during folding. The mutant protein folds with almost wild-type stability, is active, and crystallizes in two different space groups, one isomorphous with wild type and the other with two molecules in the asymmetric unit. The fold of the mutant is essentially the same in all cases, showing that the inserted segment has a well-defined structure. More than half of the inserted residues are themselves helical and extend the helix present in the wild-type protein. Participation of additional duplicated residues in this helix would have required major disruption of the parent structure. The results clearly show that the residues within the duplicated sequence tend to maintain a helical conformation even though the packing interactions with the remainder of the protein are different from those of the original helix. It supports the hypothesis that the structures of individual alpha-helices are determined predominantly by the nature of the amino acids within the helix, rather than the structural environment provided by the rest of the protein.

Comparing the thermodynamic stabilities of a related thermophilic and mesophilic enzyme.

Several models have been proposed to explain the high temperatures required to denature enzymes from thermophilic organisms; some involve greater maximum thermodynamic stability for the thermophile, and others do not. To test these models, we reversibly melted two analogous protein domains in a two-state manner. E2cd is the isolated catalytic domain of cellulase E2 from the thermophile Thermomonospora fusca. CenAP30 is the analogous domain of the cellulase CenA from the mesophile Cellulomonas fimi. When reversibly denatured in a common buffer, the thermophilic enzyme E2cd had a temperature of melting (Tm) of 72.2 degrees C, a van't Hoff enthalpy of unfolding (DeltaHVH) of 190 kcal/mol, and an entropy of unfolding (DeltaSu) of 0.55 kcal/(mol*K); the mesophilic enzyme CenAP30 had a Tm of 56.4 degrees C, a DeltaHVH of 107 kcal/mol, and a DeltaSu of 0. 32 kcal/(mol*K). The higher DeltaHVH and DeltaSu values for E2cd suggest that its free energy of unfolding (DeltaGu) has a steeper dependence on temperature at the Tm than CenAP30. This result supports models that predict a greater maximum thermodynamic stability for thermophilic enzymes than for their mesophilic counterparts. This was further explored by urea denaturation. Under reducing conditions at 30 degrees C, E2cd had a concentration of melting (Cm) of 5.2 M and a DeltaGu of 11.2 kcal/mol; CenAP30 had a Cm of 2.6 M and a DeltaGu of 4.3 kcal/mol. Under nonreducing conditions, the Cm and DeltaGu of CenAP30 were increased to 4.5 M and 10.8 kcal/mol at 30 degrees C; the Cm for E2cd was increased to at least 7.4 M at 32 degrees C. We were unable to determine a DeltaGu value for E2cd under nonreducing conditions due to problems with reversibility. These data suggest that E2cd attains its greater thermal stability (DeltaTm = 15.8 degrees C) through a greater thermodynamic stability (DeltaDeltaGu = 6.9 kcal/mol) compared to its mesophilic analogue CenAP30.

Generation of ligand binding sites in T4 lysozyme by deficiency-creating substitutions.

Several variants of T4 lysozyme have been identified that sequester small organic ligands in cavities or clefts. To evaluate potential binding sites for non-polar molecules, we screened a number of hydrophobic large-to-small mutants for stabilization in the presence of benzene. In addition to Leu99-->Ala, binding was indicated for at least five other mutants. Variants Met102-->Ala and Leu133-->Gly, and a crevice mutant, Phe104-->Ala, were further characterized using X-ray crystallography and thermal denaturation. As predicted from the shape of the cavity in the benzene complex, mutant Leu133-->Gly also bound p-xylene. We attempted to enlarge the cavity of the Met102-->Ala mutant into a deep crevice through an additional substitution, but the double mutant failed to bind ligands because an adjacent helix rearranged into a non-helical structure, apparently due to the loss of packing interactions. In general, the protein structure contracted slightly to reduce the volume of the void created by truncating substitutions and expanded upon binding the non-polar ligand, with shifts similar to those resulting from the mutations.A polar molecule binding site was also created by truncating Arg95 to alanine. This creates a highly complementary buried polar environment that can be utilized as a specific "receptor" for a guanidinium ion. Our results suggest that creating a deficiency through truncating mutations of buried residues generates "binding potential" for ligands with characteristics similar to the deleted side-chain. Analysis of complex and apo crystal structures of binding and non-binding mutants suggests that ligand size and shape as well as protein flexibility and complementarity are all determinants of binding. Binding at non-polar sites is governed by hydrophobicity and steric interactions and is relatively permissive. Binding at a polar site is more restrictive and requires extensive complementarity between the ligand and the site.

Context-dependent protein stabilization by methionine-to-leucine substitution shown in T4 lysozyme.

The substitution of methionines with leucines within the interior of a protein is expected to increase stability both because of a more favorable solvent transfer term as well as the reduced entropic cost of holding a leucine side chain in a defined position. Together, these two terms are expected to contribute about 1.4 kcal/mol to protein stability for each Met --> Leu substitution when fully buried. At the same time, this expected beneficial effect may be offset by steric factors due to differences in the shape of leucine and methionine. To investigate the interplay between these factors, all methionines in T4 lysozyme except at the amino-terminus were individually replaced with leucine. Of these mutants, M106L and M120L have stabilities 0.5 kcal/mol higher than wild-type T4 lysozyme, while M6L is significantly destabilized (-2.8 kcal/mol). M102L, described previously, is also destabilized (-0.9 kcal/mol). Based on this limited sample it appears that methionine-to-leucine substitutions can increase protein stability but only in a situation where the methionine side chain is fully or partially buried, yet allows the introduction of the leucine without concomitant steric interference. The variants, together with methionine-to-lysine substitutions at the same sites, follow the general pattern that substitutions at rigid, internal sites tend to be most destabilizing, whereas replacements at more solvent-exposed sites are better tolerated.

The response of T4 lysozyme to large-to-small substitutions within the core and its relation to the hydrophobic effect.

To further examine the structural and thermodynamic basis of hydrophobic stabilization in proteins, all of the bulky non-polar residues that are buried or largely buried within the core of T4 lysozyme were substituted with alanine. In 25 cases, including eight reported previously, it was possible to determine the crystal structures of the variants. The structures of four variants with double substitutions were also determined. In the majority of cases the "large-to-small" substitutions lead to internal cavities. In other cases declivities or channels open to the surface were formed. In some cases the structural changes were minimal (mainchain shifts < or = 0.3 A); in other cases mainchain atoms moved up to 2 A. In the case of Ile 29 --> Ala the structure collapsed to such a degree that the volume of the putative cavity was zero. Crystallographic analysis suggests that the occupancy of the engineered cavities by solvent is usually low. The mutants Val 149 --> Ala (V149A) and Met 6 --> Ala (M6A), however, are exceptions and have, respectively, one and two well-ordered water molecules within the cavity. The Val 149 --> Ala substitution allows the solvent molecule to hydrogen bond to polar atoms that are occluded in the wild-type molecule. Similarly, the replacement of Met 6 with alanine allows the two solvent molecules to hydrogen bond to each other and to polar atoms on the protein. Except for Val 149 --> Ala the loss of stability of all the cavity mutants can be rationalized as a combination of two terms. The first is a constant for a given class of substitution (e.g., -2.1 kcal/mol for all Leu --> Ala substitutions) and can be considered as the difference between the free energy of transfer of leucine and alanine from solvent to the core of the protein. The second term can be considered as the energy cost of forming the cavity and is consistent with a numerical value of 22 cal mol(-1) A(-3). Physically, this term is due to the loss of van der Waal's interactions between the bulky sidechain that is removed and the atoms that form the wall of the cavity. The overall results are consistent with the prior rationalization of Leu --> Ala mutants in T4 lysozyme by Eriksson et al. (Eriksson et al., 1992, Science 255:178-183).

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.

A test of the "jigsaw puzzle" model for protein folding by multiple methionine substitutions within the core of T4 lysozyme.

To test whether the structure of a protein is determined in a manner akin to the assembly of a jigsaw puzzle, up to 10 adjacent residues within the core of T4 lysozyme were replaced by methionine. Such variants are active and fold cooperatively with progressively reduced stability. The structure of a seven-methionine variant has been shown, crystallographically, to be similar to wild type and to maintain a well ordered core. The interaction between the core residues is, therefore, not strictly comparable with the precise spatial complementarity of the pieces of a jigsaw puzzle. Rather, a certain amount of give and take in forming the core structure is permitted. A simplified hydrophobic core sequence, imposed without genetic selection or computer-based design, is sufficient to retain native properties in a globular protein.

Thermodynamic and structural compensation in "size-switch" core repacking variants of bacteriophage T4 lysozyme.

Previous analysis of randomly generated multiple mutations within the core of bacteriophage T4 lysozyme suggested that the "large-to-small" substitution Leu121 to Ala (L121A) and the spatially adjacent "small-to-large" substitution Ala129 to Met (A129M) might be mutually compensating. To test this hypothesis, the individual variants L121A and A129M were generated, as well as the double "size-switch" mutant L121A/A129M. To make the interchange symmetrical, the combination of L121A with A129L to give L121A/A129L was also constructed. The single mutations were all destabilizing. Somewhat surprisingly, the small-to-large substitutions, which increase hydrophobic stabilization but can also introduce strain, were less deleterious than the large-to-small replacements. Both Ala129 --> Leu and Ala129 --> Met offset the destabilization of L121A by about 50%. Also, in contrast to typical Leu --> Ala core substitutions, which destabilize by 2 to 5 kcal/mol, Leu121 --> Ala slightly stabilized A129L and A129M. Crystal structure analysis showed that a combination of side-chain and backbone adjustments partially accommodated changes in side-chain volume, but only to a limited degree. For example, the cavity that was created by the Leu121 to Ala replacement actually became larger in L121A/A129L. The results demonstrate that the destabilization associated with a change in volume of one core residue can be specifically compensated by an offsetting volume change in an adjacent residue. It appears, however, that complete compensation is unlikely because it is difficult to reconstitute an equivalent set of interactions. The relatively slow evolution of core relative to surface residues appears, therefore, to be due to two factors. First, a mutation in a single core residue that results in a substantial change in size will normally lead to a significant loss in stability. Such mutations will presumably be selected against. Second, if a change in bulk does occur in a buried residue, it cannot normally be fully compensated by a mutation of an adjacent residue. Thus, the most probable response will tend to be reversion to the parent protein.

Enhancement of protein stability by the combination of point mutations in T4 lysozyme is additive.

A number of mutations have been shown previously to stabilize T4 lysozyme. By combining up to seven such mutations in the same protein, the melting temperature was incrementally increased by up to 8.3 degrees C at pH 5.4 (delta delta G = 3.6 kcal/mol). This shows that it is possible to engineer a protein of enhanced thermostability by combining a series of rationally designed point mutations. It is also shown that this stabilization is achieved with only minor, localized changes in the structure of the protein. This is consistent with the observation that the change in stability of each of the multiple mutants is, in each case, additive, i.e. equal to the sum of the stability changes associated with the constituent single mutants. One of the seven substitutions, Asn116-->Asp, changes a residue that participates in substrate binding; not surprisingly, it causes a significant loss in activity. Ignoring this mutation, there is a gradual reduction in activity as successively more mutations are combined.

Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4 lysozyme.

To determine the constraints on interactions within the core of a folded protein, we have analyzed the binding of 91 different compounds to an internal cavity created in the interior of phage T4 lysozyme by site-directed mutagenesis [Eriksson et al. (1992a) Nature 355, 371-373]. The cavity is able to accommodate a variety of small, mainly nonpolar, ligands. Molecules which do not appear to bind include those that are very polar, those that are too large, and those that have appropriate volume and polarity but inappropriate shape. Calorimetric analysis of 16 of these ligands reveals that their free energies of binding are poorly correlated with their solvent-transfer free energies. In addition, their enthalpies of binding are much larger than expected on the basis of transfer of the ligands from an aqueous to a nonpolar liquid phase. The binding energetics were analyzed by dividing the reaction into three processes: desolvation, immobilization, and packing. This analysis indicates that all three processes contribute to binding specificity. For a subset of these ligands that are structurally related, however, packing interactions in the protein interior are well modeled by the interactions of the ligands with octanol.

Alanine scanning mutagenesis of the alpha-helix 115-123 of phage T4 lysozyme: effects on structure, stability and the binding of solvent.

A series of individual alanine mutations has been constructed in the helical region 115 to 123 in phage T4 lysozyme in order to evaluate the contribution to protein stability of the different side-chains within this region. Pairwise alanine mutations and a combination mutant with seven alanine substitutions were constructed to evaluate the additive effects upon structure and stability. Only three residues within this region (Ser117, Leu118 and Leu121) have a substantial influence upon stability (change in free energy of unfolding greater than 1.0 kcal/mol). Replacement of Ser117 with alanine results in an increase in protein stability of 1.27 kcal/mol, apparently due to the release of strain present in the wild-type protein. Replacement of the buried residues Leu118 and Leu121 is destabilizing. Substitution of the remaining six residues with alanine has relatively little effect on stability. This is consistent with prior studies showing that only 20 to 30% of the residues in amphipathic helices in T4 lysozyme are critical for stability. For some of the pairwise alanine mutants the effects on stability are additive. For most of these mutants, however, there is a slight (approximately 0.15 to 0.25 kcal/mol) non-additivity such that the double mutant is more stable than the sum of the constituent single mutants. This effect is consistently observed for residues with positions i, i +4; i.e. adjacent, but in consecutive turns of the helix, suggesting a weak but significant interaction between these amino acid residues. A more pronounced non-additivity (approximately 0.5 kcal/mol) is seen in the seven-alanine combination mutant. This non-additivity is due to a modest "collapse" or "repacking" that occurs for the combination mutants (especially the multiple alanine mutant) but is not possible for the single replacements. The truncation of some side-chains permits an increase in solvent accessibility of main-chain amide and carbonyl groups. This effect is most pronounced for the seven-alanine combination mutant, where two solvent molecules, not present in wild-type, hydrogen bond to main-chain carbonyl groups in the middle region of the helix. It has been suggested that the binding of such water molecules might represent the first step in solvent-mediated unfolding of an alpha-helix. The appearance of ordered solvent, however, appears to have very little effect on stability (approximately less than 0.2 kcal/mol).

A relationship between protein stability and protein function.

Enzymes are thought to use their ordered structures to facilitate catalysis. A corollary of this theory suggests that enzyme residues involved in function are not optimized for stability. We tested this hypothesis by mutating functionally important residues in the active site of T4 lysozyme. Six mutations at two catalytic residues, Glu-11 and Asp-20, abolished or reduced enzymatic activity but increased thermal stability by 0.7-1.7 kcal.mol-1. Nine mutations at two substrate-binding residues, Ser-117 and Asn-132, increased stability by 1.2-2.0 kcal.mol-1, again at the cost of reduced activity. X-ray crystal structures show that the substituted residues complement regions of the protein surface that are used for substrate recognition in the native enzyme. In two of these structures the enzyme undergoes a general conformational change, similar to that seen in an enzyme-product complex. These results support a relationship between stability and function for T4 lysozyme. Other evidence suggests that the relationship is general.