The molecular basis for the inhibition of human cytochrome P450 1A2 by oroxylin and wogonin

In our previous kinetics studies the natural products oroxylin and wogonin were shown to have strong biological affinity for, and inhibitory effects against, human cytochrome P450 1A2, with IC₅₀ values of 579 and 248 nM, respectively; this might lead to the occurrence of drug–drug interactions when co-administered clinically. However, their inhibitory mechanisms against 1A2 remain elusive. In this study, molecular docking and molecular dynamics simulations were performed to better understand the molecular basis of their inhibitory mechanisms towards 1A2. Structural analysis revealed that oroxylin has a different binding pattern from wogonin and another very strongly binding inhibitor α-naphthoflavone (ANF, IC₅₀ = 49 nM). The O₇ atom of oroxylin forms hydrogen bonds with the OD1/OD2 atoms of Asp313, which is not observed in the 1A2–wogonin complex. Because of energetically unfavorable repulsions with the methoxy group at the 6 position of the oroxylin ring, significant conformational changes were observed for the sidechain of Thr118 in the MD simulated model. As a result, the larger and much more open binding-site architecture of the 1A2–oroxylin complex may account for its weaker inhibitory effect relative to the 1A2–ANF complex. Energy analysis indicated that oroxylin has a less negative predicted binding free energy of −19.8 kcal/mol than wogonin (−21.1 kcal/mol), which is consistent with our experimental assays. Additionally, our energy results suggest that van der Waals/hydrophobic and hydrogen-bonding interactions are important in the inhibitory mechanisms of oroxylin whereas the former is the underlying force responsible for strong inhibition by ANF and wogonin.