Cytochrome P450cam-monoterpene interactions

The use of the molecular mechanics AMBER force field (FF) to predict product profiles for the hydroxylation of the monoterpenes 1R-camphor, 1S-camphor, 1R-norcamphor, 1S-norcamphor and camphane by the enzyme cytochrome P450cam from the soil bacterium Pseudomonas putida was investigated. Predictions were carried out by applying multiple substrate (starting) orientations in the enzyme pocket in two procedures: a procedure based on molecular dynamics (MD) and a procedure based on short MD simulations followed by geometry optimisations. The latter (GO) procedure is faster and enabled the use of more monoterpene starting orientations. Monoterpene orientations were transformed into product profiles by applying both energetic and geometrical criteria appropriate for the (monoterpene) hydrogen abstraction reaction. Good predictions compared to experimental data were obtained for most compounds in both the MD and GO procedures. Prior to the product profile calculations, the FF was calibrated by reproducing the experimental data for the binding energy of 1R-camphor and 1S-camphor to P450cam and the energy of vaporisation of water. Focus of the calibration was on the value for the scaling factor for the electrostatic interactions.