Lessons from the lysozyme of phage T4

An overview is presented of some of the major insights that have come from studies of the structure, stability, and folding of T4 phage lysozyme. A major purpose of this review is to provide the reader with a complete tabulation of all of the variants that have been characterized, including melting temperatures, crystallographic data, Protein Data Bank access codes, and references to the original literature. The greatest increase in melting temperature (T_m) for any point mutant is 5.1°C for the mutant Ser 117 → Val. This is achieved in part not only by hydrophobic stabilization but also by eliminating an unusually short hydrogen bond of 2.48 Å that apparently has an unfavorable van der Waals contact. Increases in T_m of more than 3–4°C for point mutants are rare, whereas several different types of destabilizing substitutions decrease T_m by 20°C or thereabouts. The energetic cost of cavity creation and its relation to the hydrophobic effect, derived from early studies of “large‐to‐small” mutants in the core of T4 lysozyme, has recently been strongly supported by related studies of the intrinsic membrane protein bacteriorhodopsin. The L99A cavity in the C‐terminal domain of the protein, which readily binds benzene and many other ligands, has been the subject of extensive study. Crystallographic evidence, together with recent NMR analysis, suggest that these ligands are admitted by a conformational change involving Helix F and its neighbors. A total of 43 nonisomorphous crystal forms of different monomeric lysozyme mutants were obtained plus three more for synthetically‐engineered dimers. Among the 43 space groups, P2₁2₁2₁ and P2₁ were observed most frequently, consistent with the prediction of Wukovitz and Yeates.     

Hydrophobic core repacking and aromatic-aromatic interaction in the thermostable mutant of T4 lysozyme Ser 117-->Phe.

The T4 lysozyme mutant Ser 117-->Phe was isolated fortuitously and found to be more thermostable than wild-type by 1.1-1.4 kcal/mol. In the wild-type structure, the side chain of Ser 117 is in a sterically restricted region near the protein surface and forms a short hydrogen bond with Asn 132. The crystal structure of the S117F mutant shows that the introduced Phe side chain rotates by about 150 degrees about the C alpha-C beta bond relative to wild type and is buried in the hydrophobic core of the protein. Burial of Phe 117 is accommodated by rearrangements of the surrounding side chains of Leu 121, Leu 133, and Phe 153 and by main-chain shifts, which result in a minimal increase in packing density. The benzyl rings of Phe 117 and Phe 153 form a near-optimal edge-face interaction in the mutant structure. This aromatic-aromatic interaction, as well as increased hydrophobic stabilization and elimination of a close contact in the wild-type protein, apparently compensate for the loss of a hydrogen bond and the possible cost of structural rearrangements in the mutant. The structure illustrates the ability of a protein to accommodate a surprisingly large structural change in a manner that actually increases thermal stability. The mutant has activity about 10% that of wild-type, supportive of the prior hypothesis (Grütter, M.G. & Matthews, B.W., 1982, J. Mol. Biol. 154, 525-535) that the peptidoglycan substrate of T4 lysozyme makes extended contacts with the C-terminal domain in the vicinity of Ser 117.     

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.     

De novo design of the hydrophobic core of ubiquitin.

We have previously reported the development and evaluation of a computational program to assist in the design of hydrophobic cores of proteins. In an effort to investigate the role of core packing in protein structure, we have used this program, referred to as Repacking of Cores (ROC), to design several variants of the protein ubiquitin. Nine ubiquitin variants containing from three to eight hydrophobic core mutations were constructed, purified, and characterized in terms of their stability and their ability to adopt a uniquely folded native-like conformation. In general, designed ubiquitin variants are more stable than control variants in which the hydrophobic core was chosen randomly. However, in contrast to previous results with 434 cro, all designs are destabilized relative to the wild-type (WT) protein. This raises the possibility that beta-sheet structures have more stringent packing requirements than alpha-helical proteins. A more striking observation is that all variants, including random controls, adopt fairly well-defined conformations, regardless of their stability. This result supports conclusions from the cro studies that non-core residues contribute significantly to the conformational uniqueness of these proteins while core packing largely affects protein stability and has less impact on the nature or uniqueness of the fold. Concurrent with the above work, we used stability data on the nine ubiquitin variants to evaluate and improve the predictive ability of our core packing algorithm. Additional versions of the program were generated that differ in potential function parameters and sampling of side chain conformers. Reasonable correlations between experimental and predicted stabilities suggest the program will be useful in future studies to design variants with stabilities closer to that of the native protein. Taken together, the present study provides further clarification of the role of specific packing interactions in protein structure and stability, and demonstrates the benefit of using systematic computational methods to predict core packing arrangements for the design of proteins.     

Repacking the Core of T4 lysozyme by automated design

Automated protein redesign, as implemented in the program ORBIT, was used to redesign the core of phage T4 lysozyme. A total of 26 buried or partially buried sites in the C-terminal domain were allowed to vary both their sequence and side-chain conformation while the backbone and non-selected side-chains remained fixed. A variant with seven substitutions ("Core-7") was identified as having the most favorable energy. The redesign experiment was repeated with a penalty for the presence of methionine residues. In this case the redesigned protein ("Core-10") had ten amino acid changes. The two designed proteins, as well as the constituent single mutants, and several single-site revertants were over-expressed in Escherichia coli, purified, and subjected to crystallographic and thermal analyses. The thermodynamic and structural data show that some repacking was achieved although neither redesigned protein was more stable than the wild-type protein. The use of the methionine penalty was shown to be effective. Several of the side-chain rotamers in the predicted structure of Core-10 differ from those observed. Rather than changing to new rotamers predicted by the design process, side-chains tend to maintain conformations similar to those seen in the native molecule. In contrast, parts of the backbone change by up to 2.8A relative to both the designed structure and wild-type.Water molecules that are present within the lysozyme molecule were removed during the design process. In the redesigned protein the resultant cavities were, to some degree, re-occupied by side-chain atoms. In the observed structure, however, water molecules were still bound at or near their original sites. This suggests that it may be preferable to leave such water molecules in place during the design procedure. The results emphasize the specificity of the packing that occurs within the core of a typical protein. While point substitutions within the core are tolerated they almost always result in a loss of stability. Likewise, combinations of substitutions may also be tolerated but usually destabilize the protein. Experience with T4 lysozyme suggests that a general core repacking methodology with retention or enhancement of stability may be difficult to achieve without provision for shifts in the backbone.     

Subdomain interactions as a determinant in the folding and stability of T4 lysozyme.

The folding of large, multidomain proteins involves the hierarchical assembly of individual domains. It remains unclear whether the stability and folding of small, single-domain proteins occurs through a comparable assembly of small, autonomous folding units. We have investigated the relationship between two subdomains of the protein T4 lysozyme. Thermodynamically, T4 lysozyme behaves as a cooperative unit and the unfolding transition fits a two-state model. The structure of the protein, however, resembles a dumbbell with two potential subdomains: an N-terminal subdomain (residues 13-75), and a C-terminal subdomain (residues 76-164 and 1-12). To investigate the effect of uncoupling these two subdomains within the context of the native protein, we created two circular permutations, both at the subdomain interface (residues 13 and 75). Both variants adopt an active wild-type T4 lysozyme fold. The protein starting with residue 13 is 3 kcal/mol less stable than wild type, whereas the protein beginning at residue 75 is 9 kcal/mol less stable, suggesting that the placement of the termini has a major effect on protein stability while minimally affecting the fold. When isolated as protein fragments, the C-terminal subdomain folds into a marginally stable helical structure, whereas the N-terminal subdomain is predominantly unfolded. ANS fluorescence studies indicate that, at low pH, the C-terminal subdomain adopts a loosely packed acid state. An acid state intermediate is also seen for all of the full-length variants. We propose that this acid state is comprised of an unfolded N-terminal subdomain and a loosely folded C-terminal subdomain.     

Alteration of T4 lysozyme structure by second-site reversion of deleterious mutations.

Mutations that suppress the defects introduced into T4 lysozyme by single amino acid substitutions were isolated and characterized. Among 53 primary sites surveyed, 8 yielded second-site revertants; a total of 18 different mutants were obtained. Most of the restorative mutations exerted global effects, generally increasing lysozyme function in a number of primary mutant contexts. Six of them were more specific, suppressing only certain specific deleterious primary substitutions, or diminishing the function of lysozymes bearing otherwise nondeleterious primary substitutions. Some variants of proteins bearing primary substitutions at the positions of Asp 20 and Ala 98 are inferred to have significantly altered structures.     

Alanine-scanning mutagenesis of the beta-sheet region of phage T4 lysozyme suggests that tertiary context has a dominant effect on beta-sheet formation

In general, alpha-helical conformations in proteins depend in large part on the amino acid residues within the helix and their proximal interactions. For example, an alanine residue has a high propensity to adopt an alpha-helical conformation, whereas that of a glycine residue is low. The sequence preferences for beta-sheet formation are less obvious. To identify the factors that influence beta-sheet conformation, a series of scanning polyalanine mutations were made within the strands and associated turns of the beta-sheet region in T4 lysozyme. For each construct the stability of the folded protein was reduced substantially, consistent with removal of native packing interactions. However, the crystal structures showed that each of the mutants retained the beta-sheet conformation. These results suggest that the structure of the beta-sheet region of T4 lysozyme is maintained to a substantial extent by tertiary interactions with the surrounding parts of the protein. Such tertiary interactions may be important in determining the structures of beta-sheets in general.     

Activity losses among T4 lysozyme charge variants after adsorption to colloidal silica.

Enzyme structure and function depend to some extent on enzyme net charge and charge location. Altering the charge of even a single residue may affect the interaction between enzyme and substrate such that all catalytic activity is lost. In this study we investigated the effect of net charge and charge location on the enzymatic activity of synthetic mutants of bacteriophage T4 lysozyme in the presence of colloidal silica. Enzymatic activity decreased upon adsorption, and these changes were variant-specific. Results were interpreted with reference to differences in adsorbed enzyme structure and orientation, and electrostatic effects. By exploring the effects of enzyme charge on adsorption, it may be possible to gain a better understanding of how enzyme structure influences adsorption and function at an interface.     

Multiple alanine replacements within alpha-helix 126-134 of T4 lysozyme have independent, additive effects on both structure and stability.

In a systematic attempt to identify residues important in the folding and stability of T4 lysozyme, five amino acids within alpha-helix 126-134 were substituted by alanine, either singly or in selected combinations. Together with three alanines already present in the wild-type structure this provided a set of mutant proteins with up to eight alanines in sequence. All the variants behaved normally, suggesting that the majority of residues in the alpha-helix are nonessential for the folding of T4 lysozyme. Of the five individual alanine substitutions it is inferred that four result in slightly increased protein stability and one, the replacement of a buried leucine with alanine, substantially decreased stability. The results support the idea that alanine is a residue of high helix propensity. The change in protein stability observed for each of the multiple mutants is approximately equal to the sum of the energies associated with each of the constituent substitutions. All of the variants could be crystallized isomorphously with wild-type lysozyme, and, with one trivial exception, their structures were determined at high resolution. Substitution of the largely solvent-exposed residues Asp 127, Glu 128, and Val 131 with alanine caused essentially no change in structure except at the immediate site of replacement. Substitutions of the partially buried Asn 132 and the buried Leu 133 with alanine were associated with modest (< or = 0.4 A) structural adjustments. The structural changes seen in the multiple mutants were essentially a combination of those seen in the constituent single replacements. The different replacements therefore act essentially independently not only so far as changes in energy are concerned but also in their effect on structure. The destabilizing replacement Leu 133-->Ala made alpha-helix 126-134 somewhat less regular. Incorporation of additional alanine replacements tended to make the helix more uniform. For the penta-alanine variant a distinct change occurred in a crystal-packing contact, and the "hinge-bending angle" between the amino- and carboxy-terminal domains changed by 3.6 degrees. This tends to confirm that such hinge-bending in T4 lysozyme is a low-energy conformational change.     

Contributions of all 20 amino acids at site 96 to the stability and structure of T4 lysozyme

To try to resolve the loss of stability in the temperature‐sensitive mutant of T4 lysozyme, Arg 96 → His, all of the remaining 18 naturally occurring amino acids were substituted at site 96. Also, in response to suggestions that the charged residues Lys85 and Asp89, which are 5–8 Å away, may have important effects, each of these amino acids was replaced with alanine. Crystal structures were determined for many of the variants. With the exception of the tryptophan and valine mutants R96W and R96V, the crystallographic analysis shows that the substituted side chain following the path of Arg96 in wildtype (WT). The melting temperatures of the variants decrease by up to ～16°C with WT being most stable. There are two site 96 replacements, with lysine or glutamine, that leave the stability close to that of WT. The only element that the side chains of these residues have in common with the WT arginine is the set of three carbon atoms at the C^α, C^β, and C^γ positions. Although each side chain is long and flexible with a polar group at the distal position, the details of the hydrogen bonding to the rest of the protein differ in each case. Also, the glutamine replacement lacks a positive charge. This shows that there is some adaptability in achieving full stabilization at this site. At the other extreme, to be maximally destabilizing a mutation at site 96 must not only eliminate favorable interactions but also introduce an unfavorable element such as steric strain or a hydrogen‐bonding group that remains unsatisfied. Overall, the study highlights the essential need for atomic resolution site‐specific structural information to understand and to predict the stability of mutant proteins. It can be very misleading to simply assume that conservative amino acid substitutions cause small changes in stability, whereas large stability changes are associated with nonconservative replacements.     

Neutron crystallographic studies of T4 lysozyme at cryogenic temperature

Bacteriophage T4 lysozyme (T4L) has been used as a paradigm for seminal biophysical studies on protein structure, dynamics, and stability. Approximately 700 mutants of this protein and their respective complexes have been characterized by X‐ray crystallography; however, despite the high resolution diffraction limits attained in several studies, no hydrogen atoms were reported being visualized in the electron density maps. To address this, a 2.2 Å‐resolution neutron data set was collected at 80 K from a crystal of perdeuterated T4L pseudo‐wild type. We describe a near complete atomic structure of T4L, which includes the positions of 1737 hydrogen atoms determined by neutron crystallography. The cryogenic neutron model reveals explicit detail of the hydrogen bonding interactions in the protein, in addition to the protonation states of several important residues.     

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.     

Structural and thermodynamic consequences of burying a charged residue within the hydrophobic core of T4 lysozyme.

To determine the energetic and structural consequences of placing a charged group within the core of a protein, two "buried charge" mutants, Met 102----Lys (M102K) and Leu 133----Asp (L133D) were constructed in phage T4 lysozyme. Both proteins fold at neutral pH, although they are substantially less stable than wild type. The activity of M102K is about 35% that of wild type, while that of L133D is about 4%. M102K could be crystallized, and its structure was determined at high resolution. The crystal structure (at pH 6.8) of the mutant is very similar to that of wild type except for the alpha-helix that includes residues 108-113. In wild-type lysozyme, one side of this helix is exposed to solvent and the other contacts Met 102. In the M102K structure this alpha-helix becomes much more mobile, possibly allowing partial access of Lys 102 to solvent. The stability of M102K, determined by monitoring the unfolding of the protein with CD, is pH-dependent, consistent with the charged form of the substituted amino acid being more destabilizing than the uncharged form. The pKa of Lys 102 was estimated to be 6.5 both by differential titration and also by NMR analysis of isotopically labeled protein with 13C incorporated at the C epsilon position of all lysines. As the pH is lowered below pH 6.5, the overall three-dimensional structure of M102K at room temperature appears to be maintained to pH 3 or so, although there is evidence for some structural adjustment possibly allowing solvent accessibility to the protonated form of Lys 102.     

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.     

Structure and orientation of T4 lysozyme bound to the small heat shock protein alpha-crystallin

We have determined the structural changes that accompany the formation of a stable complex between a destabilized mutant of T4 lysozyme (T4L) and the small heat shock protein α-crystallin. Using pairs of fluorescence or spin label probes to fingerprint the T4L tertiary fold, we demonstrate that binding disrupts tertiary packing in the two domains as well as across the active-site cleft. Furthermore, increased distances between i and i + 4 residues of helices support a model in which the bound structure is not native-like but significantly unfolded. In the confines of the oligomer, T4L has a preferential orientation with residues in the more hydrophobic C-terminal domain sequestered in a buried environment, while residues in the N-terminal domain are exposed to the aqueous solvent. Furthermore, electron paramagnetic resonance spectral line shapes of sites in the N-terminal domain are narrower than in the folded, unbound T4L reflecting an unstructured backbone and an asymmetric pattern of contacts between T4L and α-crystallin. The net orientation is not affected by the location of the destabilizing mutation consistent with the notion that binding is not triggered by recognition of localized unfolding. Together, the structural and thermodynamic data indicate that the stably bound conformation of T4L is unfolded and support a model in which the two modes of substrate binding originate from two discrete binding sites on the chaperone.     

Improving the performance of Rosetta using multiple sequence alignment information and global measures of hydrophobic core formation

This study explores the use of multiple sequence alignment (MSA) information and global measures of hydrophobic core formation for improving the Rosetta ab initio protein structure prediction method. The most effective use of the MSA information is achieved by carrying out independent folding simulations for a subset of the homologous sequences in the MSA and then identifying the free energy minima common to all folded sequences via simultaneous clustering of the independent folding runs. Global measures of hydrophobic core formation, using ellipsoidal rather than spherical representations of the hydrophobic core, are found to be useful in removing non-native conformations before cluster analysis. Through this combination of MSA information and global measures of protein core formation, we significantly increase the performance of Rosetta on a challenging test set. Proteins 2001;43:1-11.     

Guanidinium derivatives bind preferentially and trigger long-distance conformational changes in an engineered T4 lysozyme

The binding of guanidinium ion has been shown to promote a large-scale translation of a tandemly duplicated helix in an engineered mutant of T4 lysozyme. The guanidinium ion acts as a surrogate for the guanidino group of an arginine side chain. Here we determine whether methyl- and ethylguanidinium provide better mimics. The results show that addition of the hydrophobic moieties to the ligand enhances the binding affinity concomitant with reduction in ligand solubility. Crystallographic analysis confirms that binding of the alternative ligands to the engineered site still drives the large-scale conformational change. Thermal analysis and NMR data show, in comparison to guanidinium, an increase in protein stability and in ligand affinity. This is presumably due to the successive increase in hydrophobicity in going from guanidinium to ethylguanidinium. A fluorescence-based optical method was developed to sense the ligand-triggered helix translation in solution. The results are a first step in the de novo design of a molecular switch that is not related to the normal function of the protein.     

Control of bacteriophage T4 tail lysozyme activity during the infection process

Bacteriophage T4 has an efficient mechanism for injecting the host Escherichiacoli cell with genomic DNA. Its gene product 5 (gp5) has a needle-like structure attached to the end of a tube through which the DNA passes on its way out of the head and into the host. The gp5 needle punctures the outer cell membrane and then digests the peptidoglycan cell wall in the periplasmic space. gp5 is normally post-translationally cleaved between residues 351 and 352. The function of this process in controlling the lysozyme activity of gp5 has now been investigated. When gp5 is over-expressed in E.coli, two mutants (S351H and S351A) showed a reduction of cleavage products and five other mutants (S351L, S351K, S351Y, S351Q, and S351T) showed no cleavage. Furthermore, in a complementation assay at 20 degrees C, the mutants that had no cleavage of gp5 produced a reduced number of plaques compared to wild-type T4. The crystal structure of the non-cleavage phenotype mutant of gp5, S351L, complexed with gene product 27, showed that the 18 residues in the vicinity of the potential cleavage site (disordered in the wild-type structure) had visible electron density. The polypeptide around the potential cleavage site is exposed, thus allowing access for an E.coli protease. The lysozyme activity is inhibited in the wild-type structure by a loop from the adjacent gp5 monomer that binds into the substrate-binding site. The same inhibition is apparent in the mutant structure, showing that the lysozyme is inhibited before gp5 is cleaved and, presumably, the lysozyme is activated only after gp5 has penetrated the outer membrane.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.