This study provides a comprehensive computational exploration of the inhibitory activity and metabolic pathways of 8-methoxypsoralen (8-MP), a furocoumarin derivative used for treating various skin disorders, on cytochrome P450 (P450). Employing quantum chemical DFT calculations, molecular docking, and molecular dynamics (MD) simulations analyses, the biotransformation mechanisms and the active site binding profile of 8-MP in CYP1B1 were investigated. Three plausible inactivation mechanisms were minutely scrutinized. Further analysis explored the formation of reactive metabolites in subsequent P450 metabolic processes, including covalent adduct formation through nucleophilic addition to the epoxide, 8-MP epoxide hydrolysis, and non-CYP-catalyzed epoxide ring opening. Special attention was paid to the catalytic effect of residue Phe268 on the mechanism-based inactivation (MBI) of P450 by 8-MP. Energetic profiles and facilitating conditions revealed a slight preference for the C4'=C5' epoxidation pathway, while recognizing a potential kinetic competition with the 8-OMe demethylation pathway due to comparable energy demands. The formation of covalent adducts via nucleophilic addition, particularly by phenylalanine, and the generation of potentially harmful reactive metabolites through autocatalyzed ring cleavage are likely to contribute significantly to P450 metabolism of 8-MP. Our findings highlight the key role of Phe268 in retaining 8-MP within the active site of CYP1B1, thereby facilitating initial oxygen addition transition states. This research offers crucial molecular-level insights that may guide the early stages of drug discovery and risk assessment related to the use of 8-MP.
Keywords: 8-methoxypsoralen; DFT; catalyzed biotransformation; cytochrome P450; molecular docking; molecular dynamics simulations.