Abstract

To identify a possible pathway(s) for metyrapone egress from the active site of P450 3A4, a 5-ns conventional molecular dynamics simulation followed by steered molecular dynamics simulations was performed on the complex with metyrapone. The steered molecular dynamics simulations showed that metyrapone egress via channel 1, threading through the B-C loop, only required a relatively small rupture force and small displacement of residues, whereas egress via the third channel, between helix I and helices F′ and G′, required a relatively large force and perturbation of helices I, B′, and C. The conventional dynamics simulation indicated that channel 2, located between the β1 sheet, B-B′ loop, and F′-G′ region, is closed because of the movement of residues in the mouth of this channel. The findings suggest that channel 1 can be used for metyrapone egress, whereas both channel 2 and channel 3 have a low probability of serving as an exit channel for metyrapone. In addition, residues F108 and I120 appear to act as two gatekeepers to prevent the inhibitor from leaving the active site. These results are in agreement with previous site-directed mutagenesis experiments.