Considering the planar structure and nonpolar properity of benzo(a)pyrene(B(a)p) and the planar heme part of cytochrome P-450, stacking interaction is probable. MO calculation on B(a)P and heme part of cytochrome P-450 were carried out to dertermine probable stacking interaction models. In this case, orbital interaction is most important. Accordingly, the stacking positions have high eigen vector in frontier orbital and boning type between two molecules. In this way, five probate models were selected and examined by MN2 and MO method. The most probable .stacking interaction model which is the 4, 5, 6 positions of B(a)P overlap carbon atom and pyrrole ring of ring of heme group was determined.