Predicting the product specificity and coupling of cytochrome P450cam

We present an analysis of several molecular dynamics trajectories of substrate-bound cytochrome P450cam. Trajectories were calculated for the native substrate, camphor, as well as for the alternative substrates, norcamphor and thiocamphor. The system modeled consisted of the crystallographically resolved amino acids, the heme group with a single oxygen atom as the distal ligand, the bound substrate, and the crystallographic waters. These trajectories of the presumptive ferryl oxygen intermediate were used to predict regiospecificity of hydroxylation and coupling between NADH consumption and product formation. Simple geometric criteria in combination with electronic considerations were used to calculate the probability of hydroxylation at specific sites on the substrate. We found that for all the cases examined, the predicted product ratios were in good agreement with the experimentally observed values. We also determined that these simple geometric criteria can be used to predict the degree of coupling between NADH consumption and product formation for a given substrate, which was in good agreement with the experimental values.