Aflatoxin B1 is a harmful hepatocarcinogen which is metabolized in our body by Cytochrome P450 enzymes, namely CYP1A2, CYP3A4, CYP3A5, and CYP3A7, into toxic (exo-8, 9-epoxide) and nontoxic (AFQ1, endo-epoxide) products. We have found from the literature that due to cooperativity, the rate of metabolic reactions increases in CYP1A2 and CYP3A4 involving more than one site of proteins to form two products at a given time, whereas the interaction of CYP3A5 and CYP3A7 is still unknown. Our work aims to study these four enzymes with AFB1 based on binding site pocket characterization and to find the probable resultant products at each binding site. We used computational approaches like homology modeling, molecular docking to form mono and double ligated systems, molecular dynamic simulations to analyze the potential energies (vdW & electrostatic), PCA, RMSF, and residue-wise interactions at the active as well as allosteric sites of these four enzymes. We found that CYP1A2, CYP3A4, and CYP3A5 were more hydrophobic at the first site and may induce epoxidation reaction to form toxic products, whereas the second site would be expected to be more polar and comprising charged interactions, thus enhancing non-toxic hydroxylated products. However, in CYP3A7, the first site favors hydroxylation, whereas the second site is involved in higher hydrophobic interactions. Thus, in the fetus where AFB1 is metabolized only by CYP3A7, a lower concentration of toxic metabolites will be expected, while in adults exhibiting CYP1A2, CYP3A4 and CYP3A5 may increase the concentration of the toxic metabolites due to the combined effect of these enzymes, consequently increasing liver toxicity. We believe that AFB1 binding characteristics will be helpful for medicinal chemists in the process of designing a new drug.