ON ACHIEVING BETTER THAN 1-ANGSTROM ACCURACY IN A SIMULATION OF A LARGE PROTEIN - STREPTOMYCES-GRISEUS PROTEASE-A

Computational methods are frequently used to simulate the properties of proteins. In these studies accuracy is clearly important, and the improvement of accuracy of protein simulation methodology is one of the major challenges in the application of theoretical methods, such as molecular dynamics, to structural studies of biological molecules. Much effort is being devoted to such improvements. Here, we present an analysis of a 187-ps molecular dynamics simulation of the serine protease Streptomyces griseus protease A in its crystal environment. The reproduction of the experimental structure is considerably better than has been achieved in earlier simulations-the root mean square deviation of the simulated structure from the x-ray structure being less than 1 angstrom, a significant step toward the goal of simulating proteins to within experimental error. The use of a longer cutoff with truncation rather than a switching function, inclusion of all crystalline water and the counterions in the crystallization medium, and use of the consistent valence force field characterize the differences in this calculation.