Exploring coumarin egress channels in human cytochrome P450 2A6 by random acceleration and steered molecular dynamics simulations

The kinetic analysis of coumarin oxidation by CYP2A6 suggested that substrate binding and release occurred in the multiple steps and such events proceeded rapidly. However, the crystal structure of the CYP2A6-coumarin complex reveals that no obvious channel is open enough to allow coumarin to pass through. Thus, an intriguing and important question arises: how coumarin enters and exits the active site, which is deeply buried at the center of CYP2A6 fold. In this study, geometric analysis of the potential openings was first performed on all the available crystal structures of CYP2A6. And then, random acceleration molecular dynamics simulations were used to explore the possible substrate egress channels in CYP2A6. Two channels were most frequently observed. Afterwards, steered molecular dynamics simulations were performed and potentials of mean force were constructed to compare the preference of the two channels serving as the substrate egress channel. The results showed that channel 2c, which is located between helices I and G and the helix B'-C region, was the most likely channel for coumarin egress. The opening of channel 2c was characterized by a rotation of Phe111 together with a bending of helix B'. Our findings will not only be helpful for understanding the unbinding mechanism of coumarin and for identifying structural determinants related to the biological function of CYP2A6, but also provide further insight into the channel selectivity of P450s.