Abstract
Cytochrome P450 3A4 (CYP3A4) is a membrane bound monooxygenase. It metabolizes the largest proportion of all orally ingested drugs. Ligands can enter and exit the enzyme through flexible tunnels, which co-determine CYP3A4's ligand promiscuity. The flexibility can be represented by distinct conformational states of the enzyme. However, previous state definitions relied solely on crystal structures. We employed conventional molecular dynamics (cMD) simulations to sample the conformational space of CYP3A4. Five conformationally different crystal structures embedded in a membrane were simulated for 1 µs each. A Markov state model (MSM) coupled with spectral clustering (Robust Perron Cluster Analysis PCCA +) resulted in three distinct states: Two open conformations and an intermediate conformation. The tunnels inside CYP3A4 were calculated with CAVER3.0. Notably, we observed variations in bottleneck radii compared to those derived from crystallographic data. We want to point out the importance of simulations to characterize the dynamic behaviour. Moreover, we identified a mechanism, in which the membrane supports the opening of a tunnel. Therefore, CYP3A4 must be investigated in its membrane-bound state.