The understanding of mutagenic potency has been primarily approached using "quantitative structure-activity relationships" (QSAR). Often this method allows the prediction of mutagenic potency of the compound based on its structure. But it does not give the underlying reason why the mutagenic activities differ. We have taken a set of heterocyclic amine structures and used molecular dynamic calculations to dock these molecules into the active site of a computational model of the cytochrome P4501A2 enzyme. The calculated binding strength using Boltzman distribution constants was then compared to the QSAR value (HF/6-31G* optimized structures) and the Ames/Salmonella mutagenic potency. Further understanding will only come from knowing the complete set of mutagenic determinants. These include the nitrenium ion half-life, DNA adduct half-life, efficiency of repair of the adduct, and ultimately fixation of the mutation through cellular processes. For two isomers, PhIP and 3-Me-PhIP, we showed that for the 100-fold difference in the mutagenic potency a 5-fold difference can be accounted for by differences in the P450 oxidation. The other factor of 20 is not clearly understood but is downstream from the oxidation step. The application of QSAR (chemical characteristics) to biological principles related to mutagenesis is explored in this report.