Development of In Silico Protocols to Predict Structural Insights into the Metabolic Activation Pathways of...

To establish in silico model to predict the structural insight into the metabolic bioactivation pathway of xenobiotics, we considered two specific and one non-specific mammary procarcinogen [e.g., dibenzo[a,l]pyrene (DBP), 7,12-dimethylbenz[a]anthracene (DMBA), and benzo[a]pyrene (BP)]. The CYP1A1, 1B1, 2C9, 1A2 and 2B6 reported in wet-lab studies to actively metabolize DBP also showed strong binding energies (kcal/mol) of -11.50, -10.67, -10.37, -9.76 and -9.72, respectively, with inhibition constants ranging between 0.01 and 0.08 µM. The CYP3A4 depicted minimum binding energy (-9.51 kcal/mol) which is in agreement with the wet-lab reports. Further, relatively better affinity of CYP1A1 and CYP1B1 with the dibenzo[a,l]pyrene-11,12-diol (DBPD) might be indicative of their involvement in carcinogenicity of parent compound. Like DBP, BP (-10.13 kcal/mol, Ki: 0.04 µM) and BP-diols (BPD) (-9.01 kcal/mol, Ki: 0.25 µM) observed plausible binding with CYP1A1 supporting to the reported data that emphasize the major contribution of CYP1A1 in the activation of similar procarcinogens and mutagens. Likewise, in silico results further highlighted the CYP1A1 as key player in bioactivation of DMBA to its carcinogenic metabolites. In case of PhIP metabolism, strong binding interaction predicted with CYP1A1 (-9.63 kcal/mol) rather than CYP1A2 (-8.84 kcal/mol). Dissimilarity in the binding affinity of PhIP might be due to its basic scaffold. Further, molecular dynamics (MD) simulation of 10 ns has been revealed that docked complexes of CYP1A1 with DBP, DMBA and BP are comparatively more stable than the complex of PhIP. Moreover, the current findings might be valuable as reference model in prediction and elucidation of the approximate metabolic pathway of xenobiotics.