A molecular dynamics simulation of bacteriophage T4 lysozyme

An analysis of a 400 ps molecular dynamics simulation of the164 amino acid enzyme T4 lysozyme is presented. The simulationwas carried out with all hydrogen atoms modeled explicitly,the inclusion of all 152 crystallographic waters and at a temperatureof 300 K. Temporal analysis of the trajectory versus energy,hydrogen bond stability, r.m.s. deviation from the startingcrystal structure and radius of gyration, demonstrates thatthe simulation was both stable and representative of the averageexperimental structure. Average structural properties were calculatedfrom the enzyme trajectory and compared with the crystal structure.The mean value of the C displacements of the average simulatedstructure from the X-ray structure was 1.1 ± 0.1 Å;differences of the backbone and angles between the averagesimulated structure and the crystal structure were also examined.Thermal-B factors were calculated from the simulation for heavyand backbone atoms and both were in good agreement with experimentalvalues. Relationships between protein secondary structure elementsand internal motions were studied by examining the positionalfluctuations of individual helix, sheet and turn structures.The structural integrity in the secondary structure units waspreserved throughout the simulation; however, the A helix didshow some unusually high atomic fluctuations. The largest backboneatom r.m.s. fluctuations were found in non-secondary structureregions; similar results were observed for r.m.s. fluctuationsof non-secondary structure and angles. In general, the calculatedvalues of r.m.s. fluctuations were quite small for the secondarystructure elements. In contrast, surface loops and turns exhibitedmuch larger values, being able to sample larger regions of conformationalspace. The C difference distance matrix and super-positioninganalyses comparing the X-ray structure with the average dynamicsstructure suggest that a ‘hinge-bending’ motionoccurs between the N- and C-terminal domains.