Structure prediction of subtilisin BPN' mutants using molecular dynamics methods

In this paper we describe the achievements and pitfalls encounteredin doing structure predictions of protein mutants using moleculardynamics simulation techniques in which properties of atomsare slowly changed as a function of time. Basically the methodconsists of a thermodynamic integration (slow growth) calculationused for free energy determination, but aimed at structure prediction;this allows for a fast determination of the mutant structure.We compared the calculated structure of the mutants Met222Ala,Met222Phe and Met222Gln of subtilisin BPN' with the respectiveX-ray structures and found good agreement between predictedand X-ray structure. The conformation of the residue subjectto the mutation is relatively easy to predict and is mainlydetermined by packing criteria. When the side chain has polargroups its exact orientation may pose problems; long-range Coulombinteractions may generate a polarization feedback involvingsystem relaxation times beyond the simulation time. Changesinduced in the environment are harder to predict using thismethod. In particular, rearrangement of the hydration structurewas difficult to predict correctly, probably because of thelong relaxation times. In all conversions made the changes observedin the environment were found to be history-dependent and inparticular the hydrogen bonding patterns provided evidence formetastable substates. In all cases the structure predicted wascompared with available kinetic data and the reduced activitycould be explained in terms of changes in the configurationof the active site.