Structure-based statistical thermodynamic analysis of T4 lysozyme mutants: structural mapping of cooperative interactions.

The recent development of a structural parameterization of the energetics of protein folding has permitted the incorporation of the functions that describe the enthalpy, entropy and heat capacity changes, i.e. the individual components of the Gibbs energy, into a statistical thermodynamic formalism that describes the distribution of conformational states under equilibrium conditions. The goal of this approach is to construct with the computer a large ensemble of conformational states, and then to derive the most probable population distribution, i.e. the distribution of states that best accounts for a wide array of experimental observables. This analysis has been applied to four different mutants of T4 lysozyme (S44A, S44G, V131A, V131G). It is shown that the structural parameterization predicts well the stability of the protein and the effects of the mutations. The entire set of folding constants per residue has been calculated for the four mutants. In all cases, the effect of the mutations propagates beyond the mutation site itself through sequence and three-dimensional space. This phenomenon occurs despite the fact that the mutations are at solvent-exposed locations and do not directly affect other interactions in the protein. These results suggest that single amino acid mutations at solvent-exposed locations, or other locations that cause a minimal perturbation, can be used to identify the extent of cooperative interactions. The magnitude and extent of these effects and the accuracy of the algorithm can be tested by means of NMR-detected hydrogen exchange.

A novel GDP-mannose mannosyl hydrolase shares homology with the MutT family of enzymes.

The product of the Escherichia coli orf1.9, or yefc, gene (GenBank accession number L11721) has been expressed under the control of a T7 promoter, purified to apparent homogeneity, and identified as a novel enzyme that hydrolyzes GDP-mannose or GDP-glucose to GDP and the respective hexose. The enzyme has little or no activity on other nucleotides, dinucleotides, nucleotide sugars, or sugar phosphates. It has a pH optimum between 9.0 and 9.5, a Km of 0.3 mM, and a Vmax of 1.6 mumol min-1 mg-1 for GDP-mannose, and it requires divalent cations for activity. This enzyme of 160 amino acids (M(r) = 18, 405) contains the consensus sequence GX(I/L/V)(E/Q)(X)2ET(X)6R(X)4E(X)2(I/L), characteristic of the MutT family of proteins and previously shown to form part of the nucleotide-binding site of MutT (Frick, D. N., Weber, D. J., Abeygunawardana, C., Gittis, A. G., Bessman, M. J., and Mildvan, A. S. (1995) Biochemistry 34, 5577-5586). A comparison of the enzymatic reactions catalyzed by the GDP-mannose mannosyl hydrolase and the other enzymes of the MutT family suggests that the consensus signature sequence designates a novel nucleoside diphosphate binding site and catalytic motif.

Investigation of domain motions in bacteriophage T4 lysozyme.

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge foci were found in the simulation. One locus comprises residues 8-14 near the C-terminal of the A helix; the other site, residues 77-83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19-24) acted as a mobile 'flap' and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.