Experimental and theoretical studies of the hydroxylation of a family of benzocycloarene compounds [benzocyclobutene, benzocyclopentene (indan), benzocyclohexene (tetralin), and benzocycloheptene] by wild type and Y96F mutant P450cam were performed in order to understand the factors affecting product distribution, catalytic rate and cofactor utilization. The products of all reactions except that of benzocycloheptene were regiospecifically hydroxylated in the 1-position. Reaction energetics predominated over active site steric constraints in this case so that quantum mechanical calculations (B3LYP/6-31G*) comparing the energetics of all possible radical intermediates successfully predicted hydroxylation at the 1- and 3-positions of benzocycloheptene, and at the 1-position for the other three compounds. However, the fact that the ratio of 1-alcohol to 3-alcohol changes significantly between wild type and Y96F mutant P450cam indicates that active site geometry and composition also play a significant role in determining BCA7 product regiospecificity. The indan and tetralin reaction products were stereoselective for the R enantiomer (88 and 94%, respectively). Steric constraints of the active site were confirmed by molecular dynamics calculations (locally enhanced sampling dynamics) to control enantiomer distribution for tetralin hydroxylation. NADH coupling, binding affinity, and product turnover rates were dramatically higher for Y96F P450cam, showing that the removal of the active site hydroxyl group on tyrosine makes the enzyme better suited for oxidation of these hydrophobic compounds. NADH coupling, binding affinity and product turnover rate for each enzyme generally increased with arene ring size. For both enzymes, NADH coupling and product turnover rates were correlated with the extent of high-spin shift upon substrate binding as determined by the shift in Soret absorption bands at 417 and 391 nm.