Contributions of left-handed helical residues to the structure and stability of bacteriophage T4 lysozyme

Non-glycine residues in proteins are rarely observed to have "left-handed helical" conformations. For glycine, however, this conformation is common. To determine the contributions of left-handed helical residues to the stability of a protein, two such residues in phage T4 lysozyme, Asn55 and Lys124, were replaced with glycine. The mutant proteins fold normally and are fully active, showing that left-handed non-glycine residues, although rare, do not have an indispensable role in the folding of the protein or in its activity. The thermodynamic stability of the Lys124 to Gly variant is essentially identical with that of wild-type lysozyme. The Asn55 to Gly mutant protein is marginally less stable (0.5 kcal/mol). These results indicate that the conformational energy of a glycine and a non-glycine residue in the left-handed helical conformation are very similar. This is consistent with some theoretical energy distributions, but is inconsistent with others, which suggest that replacements of the sort described here might increase the stability of the protein by up to 5 kcal/mol. Crystallographic analysis of the mutant proteins shows that the backbone conformation of the Lys124 to Gly variant is essentially identical with that of the wild-type structure. In the case of the Asn55 to Gly replacement, however, the (phi, psi) values of residue 55 change by about 20 degrees. This suggests that the energy minimum for left-handed glycine residues is not the same as that for non-glycine residues. This is strongly indicated also by a survey of accurately determined protein crystal structures, which suggests that the energy minimum for left-handed glycine residues is near (phi = 90 degrees, psi = 0 degrees), whereas that for non-glycine residues is close to (phi = 60 degrees, psi = 30 degrees). This apparent energy minimum for glycine is not clearly predicted by any of the theoretical (phi, psi) energy contour maps.     

Thermodynamic stability and point mutations of bacteriophage T4 lysozyme

The thermodynamics of melting of bacteriophage T4 lysozyme and four of its mutants have been measured by van't Hoff methods. The effect of pH has been explored and utilized to obtain the dependence of the enthalpy on temperature as suggested by Privalov and co-workers. The enthalpy change is a steep linear function of temperature. ΔC_p is large and constant within experimental error. Changes in ΔH_u are as large as 30% for a single point mutation. Changes in enthalpy are largely compensated by changes in entropy. Changes in stability, as measured by the free energy of unfolding, are smaller than those of ΔH, but are very large in a relative sense, since ΔG is very much smaller than ΔH. Origins of the destabilization caused by mutations are discussed.     

Multiple methionine substitutions are tolerated in T4 lysozyme and have coupled effects on folding and stability

In order to further explore the tolerance of proteins to amino acid substitutions within the interior, a series of core residues was replaced by methionine within the C-terminal domain of T4 lysozyme. By replacing leucine, isoleucine, valine and phenylalanine residues a total of 10 methionines could be introduced, which corresponds to a third of the residues that are buried in this domain. As more methionines are incorporated the protein gradually loses stability. This is attributed in part to a reduction in hydrophobic stabilization, in part to the increased entropic cost of localizing the long, flexible methionine sidechains, and in part to steric clashes. The changes in structure of the mutants relative to the wildtype protein are modest but tend to increase in an additive fashion as more methionines are included. In the most extreme case, namely the 10-methionine mutant, much of the C-terminal domain remains quite similar to wildtype (root-mean-square backbone shifts of 0.56 Å), while the F and G helices undergo rotations of approximately 20° and center-of-mass shifts of approximately 1.4 Å. For up to six methionine substitutions the changes in stability are additive. Beyond this point, however, the multiple mutants are somewhat more stable than suggested from the sum of their constituents, especially for those including the replacement Val111→Met. This is interpreted in terms of the larger structural changes associated with this substitution. The substituted sidechains in the mutant structures have somewhat higher crystallographic thermal factors than their counterparts in WT*. Nevertheless, the interiors of the mutant proteins retain a well-defined structure with little suggestion of molten-globule characteristics. Lysozymes in which selenomethionine has been incorporated rather than methionine tend to have increased stability. At the same time they also fold faster. This provides further evidence that, at the rate-limiting step in folding, the structure of the C-terminal domain of T4 lysozyme is similar to that of the fully folded protein.     

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.     

Thermodynamic effects of mutations on the denaturation of T4 lysozyme.

We investigated the folding of substantially destabilized mutant forms of T4 lysozyme using differential scanning calorimetry and circular dichroism measurements. Three mutations in an alpha-helix in the protein's N-terminal region, the alanine insertion mutations S44[A] and K48[A], and the substitution A42K had previously been observed to result in unexpectedly low apparent enthalpy changes of melting, compared to a pseudo-wild-type reference protein. The pseudo-wild-type reference protein thermally unfolds in an essentially two-state manner. However, we found that the unfolding of the three mutant proteins has reduced cooperativity, which partially explains their lower apparent enthalpy changes. A three-state unfolding model including a discrete intermediate is necessary to describe the melting of the mutant proteins. The reduction in cooperativity must be considered for accurate calculation of the energy changes of folding. Unfolding in two stages reflects the underlying two-subdomain structure of the lysozyme protein family.     

Simulation analysis of the stability mutant R96H of T4 lysozyme

Free energy simulation methods are used to analyze the effects of the mutation Arg 96----His on the stability of T4 lysozyme. The calculated stability change and the lack of significant structural rearrangement in the folded state due to the mutation are in agreement with experimental studies [Kitamura, S., & Sturtevant, J. M. (1989) Biochemistry 28, 3788-3792; Weaver, L. H., et al. (1989) Biochemistry 28, 3793-3797]. By use of thermodynamic integration, the contributions of specific interactions to the free energy change are evaluated. It is shown that a number of contributions that stabilize the wild type or the mutant partially cancel in the overall free energy difference; some of these involve the unfolded state. Comparison of the results with conclusions based on structural and thermodynamic data leads to new insights into the origin of the stability difference between wild-type and mutant proteins. Of particular interest is the importance of the contributions of more distant residues, solvent water, and the covalent linkage of the mutated amino acid. Also, the analysis of the interactions of Arg/His 96 with the C-terminal end of a helix (residues 82-90) makes it clear that the nearby carbonyl groups (Tyr 88 and Asp 89) make the dominant contribution, that the amide groups do not contribute significantly, and that the helix-dipole model is inappropriate for this case.     

Theoretical studies on solvation contribution to the thermodynamic stability of mutants of lysozyme T4

Atomic solvation parameters (ASPs) are widely used to estimate the solvation contribution to the thermodynamic stability of proteins as well as the free energy of association for protein-ligand complexes. In view of discrepancies in the results of free energies of solvation of folding for various proteins obtained using different atomic solvation parameter sets, systematic studies have been carried out for the calculation of accessible surface area and the changes in free energy of solvation of folding (deltaG(s,f)) for mutants of lysozyme T4 where threonine 157 is replaced by amino acids: cysteine, aspartate, glutamate, phenylalanine, glycine, histidine, isoleucine, leucine, asparagine, arginine, serine and valine. The deviations of the calculated results from the experimental results are discussed to highlight the discrepancies in the atomic solvation parameter sets and possible reasons for them. The results are also discussed to throw light on the effect of chain free energy and hydrogen bonding on the stability of mutants. The octanol to water-based ASP sets 'Sch1' and 'EM' perform better than the vacuum to water-based ASP sets. The vacuum to water-based ASP sets 'Sch3' and 'WE' can be used to predict the stability of mutants if a proper method to calculate the hydrogen bond contribution to overall stability is in place.     

Role of medium- and long-range interactions to the stability of the mutants of T4 lysozyme

Inter-residue interactions play an important role to the folding and stability of protein molecules. In this work, we analyze the role of medium- and long-range interactions to the stability of T4 lysozyme mutants. We found that, in buried mutations, the increase in long-range contacts upon mutations destabilizes the protein, whereas, in surface mutations, the increase in long-range contacts increases the stability, indicating the importance of surrounding polar residues to the stability of surface mutations. Further, the increase in medium-range contacts decreases the stability of buried and surface mutations and a direct relationship is observed between the increase of medium-range contacts and increase in stability for partially buried/exposed mutations. Moreover, the relationship between amino acid properties and stability of T4 lysozyme mutants at positions Ile3, Phe53, and Leu99 showed that the effect of medium- and long-range contacts is less for buried mutations and the inter-residue contacts have significant correlation with the stability of partially buried mutations.     

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.     

Comparison of the predicted and observed secondary structure of T4 phage lysozyme.

Predictions of the secondary structure of T4 phage lysozyme, made by a number of investigators on the basis of the amino acid sequence, are compared with the structure of the protein determined experimentally by X-ray crystallography. Within the amino terminal half of the molecule the locations of helices predicted by a number of methods agree moderately well with the observed structure, however within the carboxyl half of the molecule the overall agreement is poor. For eleven different helix predictions, the coefficients giving the correlation between prediction and observation range from 0.14 to 0.42. The accuracy of the predictions for both beta-sheet regions and for turns are generally lower than for the helices, and in a number of instances the agreement between prediction and observation is no better than would be expected for a random selection of residues. The structural predictions for T4 phage lysozyme are much less successful than was the case for adenylate kinase (Schulz et al. (1974) Nature 250, 140-142). No one method of prediction is clearly superior to all others, and although empirical predictions based on larger numbers of known protein structure tend to be more accurate than those based on a limited sample, the improvement in accuracy is not dramatic, suggesting that the accuracy of current empirical predictive methods will not be substantially increased simply by the inclusion of more data from additional protein structure determinations.     

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.     

The tail lysozyme complex of bacteriophage T4

The tail baseplate of bacteriophage T4 contains a structurally essential, three-domain protein encoded by gene 5 in which the middle domain possesses lysozyme activity. The gene 5 product (gp5) undergoes post-translational cleavage, allowing the resultant N-terminal domain (gp5*) to assemble into the baseplate as a trimer. The lysozyme activity of the undissociated cleaved gp5 is inhibited until infection has been initiated, when the C-terminal portion of the molecule is detached and the rest of the molecule dissociates into monomers. The 3D structure of the undissociated cleaved gp5, complexed with gp27 (another component of the baseplate), shows that it is a cell-puncturing device that functions to penetrate the outer cell membrane and to locally dissolve the periplasmic cell wall.     

A systematic method for analysing the protein hydration structure of T4 lysozyme

A systematic method for the analysis of the hydration structure of proteins is demonstrated on the case study of lysozyme. The method utilises multiple structural data of the same protein deposited in the protein data bank. Clusters of high water occupancy are localised and characterised in terms of their interaction with protein. It is shown that they constitute a network of interconnected hydrogen bonds anchored to the protein molecule. The high occupancy of the clusters does not directly correlate with water–protein interaction energy as was originally hypothesised. The highly occupied clusters rather correspond to the nodes of the hydration network that have the maximum number of hydrogen bonds including both the protein atoms and the surrounding water clusters. Copyright © 2013 John Wiley & Sons, Ltd.     

Measured and calculated effects of mutations in bacteriophage T4 lysozyme on interactions in solution

Understanding the molecular determinants of protein interactions in solution has fundamental implications for understanding protein solution thermodynamics and, hence, processes as diverse as separations performance and cellular self-organization. Our earlier theoretical calculations indicate that the protein-protein interactions are dominated by a small number of configurations in which highly complementary surface regions are apposed, rather than by the overall colloidal interactions. To examine this paradigm more explicitly, we investigated the effects of protein structural modifications on protein-protein interactions. Experimental measurements are presented of B(22)(') values of a set of mutants of Ser44 in bacteriophage T4 lysozyme. Effects are seen with both charged and uncharged substitutions. The results with the charged substitutions follow the expected trends, whereas those with the uncharged substitutions may be explained by the impact of the mutations on the local protein geometry, which directly affects the complementarity of protein interactions. These effects are also captured well by molecular calculations that account for the mutations. The interaction energetics between protein pairs could provide information on the propensity for adventitious interactions, which can have important implications for separations and for normal and pathological self-assembly. Thus, protein structural data implicit in genomic information, coupled with appropriate calculational and experimental tools, can ultimately provide insights into protein interactions in vivo and in bioprocessing.