Analysis of Active Site Motions from a 175 picosecond Molecular Dynamics Simulation of Camphor-bound Cytochrome P450cam

Abstract
The structure and internal motions of the active site residues of camphor-bound cytochrome P450cam have been evaluated on the basis of a 175 psec molecular dynamics simulation. The active site residues generally show very small deviations away from their starting crystal positions. These residues also generally show much smaller fluctuations than for the enzyme as a whole. Phe 87 is dynamically very unusual and is suggested to play a role in substrate movement into and/or out of the active site. The average distance between the heme iron and atoms C5, C6, and C3 of camphor is 5.3,6.0, and 7.0 Å, respectively. This trend is consistent with the experimentally observed stereospecificity of the hydroxylation reaction. On the basis of distance and angle criteria, both 5-exo and 5-endo hydrogen abstraction are predicted to occur during the hydroxylation reaction; although the 5-exo pathway is expected to be 3-fold more likely.